Search

We present the [OIII]52{\mu}m map of the dwarf galaxy IC10, obtained with the Field-Imaging Far-Infrared Line Spectrometer (FIFI-LS) on board the Stratospheric Observatory for Infrared Astronomy (SOFIA). We combine the [OIII]52{\mu}m map with Herschel and Spitzer observations, to estimate the electron density distribution of the brightest HII regions of IC10. We find that the line ratio [OIII]88{\mu}m/[OIII]52{\mu}m gives electron density (n_e) values (n_e_OIII) that cover a broad range, while the n_e values obtained using the line ratio [SIII]33{\mu}m/[SIII]18{\mu}m (n_e_SIII) are all similar within the uncertainties. n_e_OIII is similar to n_e_SIII for the M1, M2 and A1 regions, and it is higher than n_e_SIII for the two regions, A2 and M1b, which are the brightest in the 24{\mu}m continuum emission. These results suggest that for these regions the two ions, O++ and S++, trace two different ionised gas components, and that the properties of the ionised gas component traced by the O++ ion are more sensitive to the local physical conditions. In fact, while the gas layer traced by [SIII] does not keep track of the characteristics of the radiation field, the n_e_OIII, correlates with the star formation rate (SFR), the dust temperature and the 24{\mu}m. Therefore, n_e_OIII is an indicator of the evolutionary stage of the HII region and the radiation field, with higher n_e_OIII, found in younger SF regions and in more energetic environments., Comment: 19 pages, 13 figures, 5 tables. Paper accepted for publication in ApJ

Low-metallicity dwarf galaxies often show no or little CO emission, despite the intense star formation observed in local samples. Both simulations and resolved observations indicate that molecular gas in low-metallicity galaxies may reside in small dense clumps, surrounded by a substantial amount of more diffuse gas, not traced by CO. Constraining the relative importance of CO-bright versus CO-dark H2 star-forming reservoirs is crucial to understand how star formation proceeds at low metallicity. We put to the test classically used single component radiative transfer models and compare their results to those obtained assuming an increasingly complex structure of the interstellar gas, mimicking an inhomogeneous distribution of clouds with various physical properties. We compute representative models of the interstellar medium as combinations of several gas components, each with a specific set of physical parameters. We introduce physically-motivated models assuming power-law distributions for the density, ionization parameter, and the depth of molecular clouds. We confirm the presence of a predominantly CO-dark molecular reservoir in low-metallicity galaxies. The predicted total H2 mass is best traced by [C II]158um and, to a lesser extent, by [CI] 609um, rather than by CO(1-0). We examine the CO-to-H2 conversion factor vs. metallicity relation and find that its dispersion increases significantly when different geometries of the gas are considered. We define a clumpiness parameter that anti-correlates with [CII]/CO and explains the dispersion of the CO-to-H2 conversion factor vs. metallicity relation. We find that low-metallicity galaxies with high clumpiness may have CO-to-H2 conversion factor as low as the Galactic value. We identify the clumpiness of molecular gas as a key parameter to understand variations of geometry-sensitive quantities, such as CO-to-H2 conversion factor., Comment: Accepted for publication in A&A, 24 pages, 14 figures

(abridged) Quantifying the ISM porosity to ionizing photons in nearby galaxies may improve our understanding of the mechanisms leading to Lyman Continuum photons leakage from galaxies. Primitive galaxies with low metal and dust content have been shown to host a more patchy and porous ISM than their high-metallicity counterparts. To what extent this peculiar structure contributes to the leakage of ionizing photons remains to be quantitatively studied. To address these questions we build a refined grid of models including density-bounded regions and a possible contribution of an X-ray source. Using MULTIGRIS, a new Bayesian code based on Monte Carlo sampling, we combine the models as sectors under various assumptions to extract the probability density distributions of the parameters and infer the corresponding escape fractions from H II regions (fesc,HII). We apply this new code to a sample of 39 well-know local starbursting dwarf galaxies from the Dwarf Galaxy Survey. We confirm previous results hinting at an increased porosity to ionizing photons of the ISM in low-metallicity galaxies and provide, for the first time, quantitative predictions for fesc,HII. The predicted fesc,HII for low-metallicity objects span a large range of values, up to 60%, while the values derived for more metal-rich galaxies are globally lower. We also examine the influence of other parameters on the escape fractions, and find that the specific star-formation rate correlates best with fesc,HII . Finally, we provide observational line ratios which could be used as tracers of the photons escaping from density-bounded regions. Although this multi-sector modelling remains too simple to fully capture the ISM complexity, it can be used to preselect galaxy samples with potential leakage of ionizing photons based on current and up-coming spectral data in unresolved surveys of local and high-redshift galaxies., Comment: 39 pages, 24 figures, accepted for publication in A&A

While beta Pic is known to host silicates in ring-like structures, whether the properties of these silicate dust vary with stellocentric distance remains an open question. We re-analyze the beta Pictoris debris disk spectrum from the Spitzer Infrared Spectrograph (IRS) and a new IRTF/SpeX spectrum to investigate trends in Fe/Mg ratio, shape, and crystallinity in grains as a function of wavelength, a proxy for stellocentric distance. By analyzing a re-calibrated and re-extracted spectrum, we identify a new 18 micron forsterite emission feature and recover a 23 micron forsterite emission feature with a substantially larger line-to-continuum ratio than previously reported. We find that these prominent spectral features are primarily produced by small submicron-sized grains, which are continuously generated and replenished from planetesimal collisions in the disk and can elucidate their parent bodies' composition. We discover three trends about these small grains: as stellocentric distance increases, (1) small silicate grains become more crystalline (less amorphous), (2) they become more irregular in shape, and (3) for crystalline silicate grains, the Fe/Mg ratio decreases. Applying these trends to beta Pic's planetary architecture, we find that the dust population exterior to the orbits of beta Pic b and c differs substantially in crystallinity and shape. We also find a tentative 3-5 micron dust excess due to spatially unresolved hot dust emission close to the star. From our findings, we infer that the surfaces of large planetesimals are more Fe-rich and collisionally-processed closer to the star but more Fe-poor and primordial farther from the star., Comment: 19 pages, 12 figures, Accepted for Publication in ApJ

We present the Infrared Database of Extragalactic Observables from Spitzer (IDEOS), a homogeneous, publicly available, database of 77 fitted mid-infrared observables in the 5.4-36um range, comprising measurements for 3335 galaxies observed in the low-resolution staring mode of the Infrared Spectrometer onboard the Spitzer Space Telescope. Among the included observables are PAH fluxes and their equivalent widths, the strength of the 9.8um silicate feature, emission line fluxes, solid-state features, rest frame continuum fluxes, synthetic photometry, and a mid-infrared spectral classification. The IDEOS spectra were selected from the Cornell Atlas of Spitzer-IRS Sources. To our surprise we have detected at a >95% confidence level crystalline silicates in the spectra of 786 IDEOS galaxies. The detections range from single band detections to detections of all fitted crystalline bands (16, 19, 23, 28 and 33um). We find the strength of the crystalline silicate bands to correlate with the amorphous silicate strength, and the change from an emission to an absorption feature to occur at higher obscuration as the wavelength of the crystalline silicate band is longer. These observed characteristics are consistent with an origin for the amorphous and crystalline silicate features in a centrally heated dust geometry, either an edge-on disk or a cocoon. We find the 23 and 33um crystalline silicate bands to be well-suited to classify the obscuration level of galactic nuclei, even in the presence of strong circumnuclear star formation. Based on our detection statistics, we conclude that crystalline silicates are a common component of the interstellar medium of galactic nuclei., Comment: 46 pages, 34 figures, accepted for publication in ApJS. Associated spectra and data table can be downloaded from the IDEOS portal at http://ideos.astro.cornell.edu

The metallicity and gas density dependence of interstellar depletions, the dust-to-gas (D/G), and dust-to-metal (D/M) ratios have important implications for how accurately we can trace the chemical enrichment of the universe; either by using FIR dust emission as a tracer of the ISM; or by using spectroscopy of damped Lyman-alpha systems (DLAs) to measure chemical abundances over a wide range of redshifts. We collect and compare large samples of depletion measurements in the Milky Way (MW), LMC (Z=0.5 Zsun), and SMC (Z=0.2 Zsun). The relation between the depletions of different elements do not strongly vary between the three galaxies, implying that abundance ratios should trace depletions accurately down to 20% solar metallicity. From the depletions, we derive D/G and D/M. The D/G increases with density, consistent with the more efficient accretion of gas-phase metals onto dust grains in the denser ISM. For log N(H) > 21 cm^-2, the depletion of metallicity tracers (S, Zn) exceeds -0.5 dex, even at 20% solar metallicity. The gas fraction of metals increases from the MW to the LMC (factor 3) and SMC (factor 6), compensating the reduction in total heavy element abundances and resulting in those three galaxies having the same neutral gas-phase metallicities. The D/G derived from depletions are a factor of 2 (LMC) and 5 (SMC) higher than the D/G derived from FIR, 21 cm, and CO emission, likely due to the combined uncertainties on the dust FIR opacity and on the depletion of carbon and oxygen., Comment: 26 pages, 12 figures

In review of the White Papers from the Voyage 2050 process and after the public presentation of a number of these papers in October 2019 in Madrid, we as White Paper lead authors have identified a coherent science theme that transcends the divisions around which the Topical Teams are structured. This note aims to highlight this synergistic science theme and to make the Topical Teams and the Voyage 2050 Senior Committee aware of the wide importance of these topics and the broad support that they have across the worldwide science community., Comment: 5 pages, published in Frontiers in Astronomy and Space Science

While the CO(1-0) transition is often used to deduce the total molecular hydrogen in galaxies, it is challenging to detect in low metallicity galaxies, in spite of the star formation taking place. In contrast, the [CII] 158 micron line is relatively bright, highlighting a potentially important reservoir of H2 that is not traced by CO(1-0), but residing in the C+ - emitting regions. We explore a method to quantify the total H2 mass (MH2) in galaxies and learn what parameters control the CO-dark gas reservoir. We present Cloudy grids of density, radiation field and metallicity in terms of observed quantities, such as [OI], [CI], CO(1-0), [CII], total infrared luminosity and the total MH2 and provide recipes based on these models to derive total MH2 mass estimates from observations. The models are applied to the Herschel Dwarf Galaxy Survey, extracting the total MH2 for each galaxy which is compared to the H2 determined from the observed CO(1-0) line. While the H2 traced by CO(1-0) can be negligible, the [CII] 158 micron line can trace the total H2. 70% to 100% of the total H2 mass is not traced by CO(1-0) in the dwarf galaxies, but is well-traced by [CII] 158 micron line. The CO-dark gas mass fraction correlates with the observed L[CII]/LCO(1-0) ratio. A conversion factor for [CII] luminosity to total H2 and a new CO-to-total-MH2 conversion factor, as a function of metallicity, is presented. A recipe is provided to quantify the total mass of H2 in galaxies, taking into account the CO and [CII] observations. Accounting for this CO-dark H2 gas, we find that the star forming dwarf galaxies now fall on the Schmidt-Kennicutt relation. Their star-forming efficiency is rather normal, since the reservoir from which they form stars is now more massive when introducing the [CII] measures of the total H2, compared to the little amount of H2 in the CO-emitting region., Comment: Accepted for publication in A&A

Investigating star formation requires precise knowledge of the properties of the dense molecular gas. The low metallicity and wide range of star formation activity of the Large and Small Magellanic Clouds make them prime laboratories to study how local physical conditions impact the dense gas reservoirs. The aim of the Dense Gas Survey for the Magellanic Clouds (DeGaS-MC) project is to expand our knowledge of the relation between dense gas properties and star formation activity by targeting the LMC and SMC observed in the HCO+(2-1) and HCN(2-1) transitions. We carried out a pointing survey toward 30 LMC and SMC molecular clouds using the SEPIA180 instrument installed on the APEX telescope and a follow-up mapping campaign in 13 star-forming regions. This first paper provides line characteristic catalogs and integrated line-intensity maps of the sources. HCO+(2-1) is detected in 20 and HCN(2-1) in 8 of the 29 pointings observed. The dense gas velocity pattern follows the line-of-sight velocity field derived from the stellar population. The HCN emission is less extended than the HCO+ emission. The HCO+(2-1)/HCN(2-1) brightness temperature ratios range from 1 to 7, which is consistent with the large ratios commonly observed in low-metallicity environments. A larger number of young stellar objects are found at high HCO+ intensities and lower HCO+/HCN flux ratios, and thus toward denser lines of sight. The dense gas luminosities correlate with the star formation rate traced by the total infrared luminosity over the two orders of magnitude covered by our observations, although substantial region-to-region variations are observed., Comment: Accepted for publication in A&A; 25 pages, 13 figures, 5 tables

(abridged) The ambiguous origin of [CII] 158um in the interstellar medium complicates its use for diagnostics concerning the star-formation rate and physical conditions in photodissociation regions (PDRs). We observed the giant HII region N11 in the Large Magellanic Cloud with SOFIA/GREAT in order to investigate the origin of [CII] to obtain the total H2 gas content, the fraction of CO-dark H2 gas, and the influence of environmental effects such as stellar feedback. We present an innovative spectral decomposition method that allows statistical trends to be derived. The [CII] line is resolved in velocity and compared to HI and CO, using a Bayesian approach to decompose the profiles. A simple model accounting for collisions in the neutral atomic and molecular gas was used in order to derive the H2 column density traced by C+. The profile of [CII] most closely resembles that of CO, but the integrated [CII] line width lies between that of CO and that of HI. Using various methods, we find that [CII] mostly originates from the neutral gas. We show that [CII] mostly traces the CO-dark H2 gas but there is evidence of a weak contribution from neutral atomic gas preferentially in the faintest components. Most of the molecular gas is CO-dark. The fraction of CO-dark H2 gas decreases with increasing CO column density, with a slope that seems to depend on the impinging radiation field from nearby massive stars. Finally we extend previous measurements of the photoelectric-effect heating efficiency, which we find is constant across regions probed with Herschel, with [CII] and [OI] being the main coolants in faint and diffuse, and bright and compact regions, respectively, and with PAH emission tracing the CO-dark H2 gas heating where [CII] and [OI] emit. Our study highlights the importance of velocity-resolved PDR diagnostics and higher spatial resolution for HI observations., Comment: Accepted for publication in A&A

We present Herschel Space Observatory Photodetector Array Camera and Spectrometer (PACS) and Spectral and Photometric Imaging Receiver Fourier Transform Spectrometer (SPIRE FTS) spectroscopy of a sample of twenty massive Young Stellar Objects (YSOs) in the Large and Small Magellanic Clouds (LMC and SMC). We analyse the brightest far infrared (far-IR) emission lines, that diagnose the conditions of the heated gas in the YSO envelope and pinpoint their physical origin.We compare the properties of massive Magellanic and Galactic YSOs.We find that [OI] and [CII] emission, that originates from the photodissociation region associated with the YSOs, is enhanced with respect to the dust continuum in the Magellanic sample. Furthermore the photoelectric heating efficiency is systematically higher for Magellanic YSOs, consistent with reduced grain charge in low metallicity environments. The observed CO emission is likely due to multiple shock components. The gas temperatures, derived from the analysis of CO rotational diagrams, are similar to Galactic estimates. This suggests a common origin to the observed CO excitation, from low-luminosity to massive YSOs, both in the Galaxy and the Magellanic Clouds. Bright far-IR line emission provides a mechanism to cool the YSO environment. We find that, even though [OI], CO and [CII] are the main line coolants, there is an indication that CO becomes less important at low metallicity, especially for the SMC sources. This is consistent with a reduction in CO abundance in environments where the dust is warmer due to reduced ultraviolet-shielding. Weak H$_2$O and OH emission is detected, consistent with a modest role in the energy balance of wider massive YSO environments., Comment: 38 pages; accepted for publication by MNRAS; full integrated version; journal version will include appendices on-line only

(abridged) The far-UV wavelength range (912-2000A) provides access to atomic and molecular transitions of many species the interstellar medium (ISM), circumgalactic medium (CGM), and intergalactic medium, within phases spanning a wide range of ionization, density, temperature, and molecular gas fraction. Far-UV space telescopes have enabled detailed studies of the ISM in the Milky Way thanks to absorption features appearing in the UV spectra of hot stars and yielding fundamental insights into the composition and physical characteristics of all phases of the ISM along with the processes that influence them. However, we have yet to design a spectrometer able to observe the full UV domain at resolving power R>10^5 with a signal-to-noise ratio SNR>500. Such a resolution is necessary to resolve lines from both the cold molecular hydrogen and the warm metal ions with a turbulent velocity of about 1 km s-1, and to differentiate distinct velocity components. Future UV spectroscopic studies of the Milky Way ISM must revolutionize our understanding of the ISM as a dynamical, unstable, and magnetized medium, and rise to the challenge brought forward by current theories. Another interesting prospect is to transpose the same level of details that has been reached for the Milky Way to the ISM in external galaxies, in particular in metal-poor galaxies, where the ISM chemical composition, physical conditions, and topology change dramatically, with significant consequences on the star-formation properties. Finally, we need to be able to perform statistical analyses of background quasar lines of sight intersecting the CGM of galaxies at various redshifts and to comprehend the role of gas exchanges and flows for galaxy evolution., Comment: This ESA white paper draws from and expands on the following white papers previously submitted for the Astro 2020 Decadal Survey, Gry et al. (2019), Lebouteiller et al. (2019), and Yan et al. (2019). The original ESA white paper was reformatted as a 2-column paper for the present document. (arXiv admin note: substantial text overlap with arXiv:1903.07295, arXiv:1903.06250, arXiv:1903.11065)

The sensitive infrared telescopes, Spitzer and Herschel, have been used to target low-metallicity star-forming galaxies, allowing us to investigate the properties of their interstellar medium (ISM) in unprecedented detail. Interpretation of the observations in physical terms relies on careful modeling of those properties. We have employed a multiphase approach to model the ISM phases (HII region and photodissociation region) with the spectral synthesis code Cloudy. Our goal is to characterize the physical conditions (gas densities, radiation fields, etc.) in the ISM of the galaxies from the Herschel Dwarf Galaxy Survey. We are particularly interested in correlations between those physical conditions and metallicity or star-formation rate. Other key issues we have addressed are the contribution of different ISM phases to the total line emission, especially of the [CII]157um line, and the characterization of the porosity of the ISM. We find that the lower-metallicity galaxies of our sample tend to have higher ionization parameters and galaxies with higher specific star-formation rates have higher gas densities. The [CII] emission arises mainly from PDRs and the contribution from the ionized gas phases is small, typically less than 30% of the observed emission. We also find correlation - though with scatter - between metallicity and both the PDR covering factor and the fraction of [CII] from the ionized gas. Overall, the low metal abundances appear to be driving most of the changes in the ISM structure and conditions of these galaxies, and not the high specific star-formation rates. These results demonstrate in a quantitative way the increase of ISM porosity at low metallicity. Such porosity may be typical of galaxies in the young Universe., Comment: Accepted in A&A. 20 pages, 12 figures, 5 tables plus appendices

The wavelength range 912-2000A (hereafter far-UV) provides access to absorption lines of the interstellar medium (ISM), circumgalactic medium (CGM), and intergalactic medium (IGM) in phases spanning a wide range of ionization, density, temperature, and molecular gas fraction. Far-UV space telescopes have enabled detailed studies of the ISM in the Milky Way, paving the way to understand in particular depletion of elements by condensation onto dust grains, molecular gas formation pathways, ionization processes, and the dynamics of shocks. An interesting prospect is to transpose the same level of details to the ISM in external galaxies, in particular in metal-poor galaxies, where the ISM chemical composition, physical conditions, and topology (arrangement of various phases) change dramatically, with significant consequences on the star-formation properties and on the overall galaxy evolution. To circumvent current systematic biases in column density determinations and to examine the ISM enrichment as a function of the environment, we advocate for a versatile far-UV space mission able to observe individual O/B stars or stellar clusters in galaxies up to few 100s Mpc away with a spectral resolution power R~10^{4-5}. Such requirements would also make it possible to probe multiple quasar lines of sight intersecting the CGM of galaxies at various redshifts, making it possible to reconcile the various physical scales and phases and to comprehend the role of gas exchanges and flows for galaxy evolution., Comment: Astro2020 Science White Paper

Our objective is to investigate the physical properties of the ionised gas of the low-metallicity dwarf galaxy, IC 10, at various spatial scales: from individual HII regions to the entire galaxy scale and examine whether diagnostics for integrated measurements introduce bias in the results. We modeled the ionised gas combining the mid- and far-infrared fine-structure cooling lines observed with Spitzer/IRS and Herschel/PACS, with the photoionisation code Cloudy. The free parameters of the models are the age of the stellar cluster, the density and the ionisation parameter of the ionised gas as well as the depth of the cloud. The latter is used to investigate the leakage of the ionising photons from the analysed regions of IC 10. We investigate HII regions in the main star-forming body, on scales of ~25 pc, three in the main star-forming region in the center of the galaxy and two on the first arc. We then consider larger sizes on the scale of ~200 pc. We find that most clumps have nearly identical properties, density ~10$^{2.} $ - 10$^{2.6}$ cm$^{-3}$, ionisation parameter between 10$^{-2.2}$ and 10$^{-1.6}$ and age of the stellar cluster ~5.5 Myr. All of them are matter-bounded regions, allowing ionising photons to leak. The relatively uniform physical properties of the clumps suggest a common origin for their star formation activity, which could be related to the feedback from stellar winds or supernovae of a previous generation of stars. The properties derived for ~200 pc size "zones" have similar properties as the HII regions they encompass, but with the larger regions tending to be more radiation-bounded. Finally, we investigate the fraction of [CII] 157.7 {\mu}m, [SiII] 34.8 {\mu}m, and [FeII] 25.9 {\mu}m, emission arising from the ionised gas phase and we find that most of the emission originates from the neutral gas, not from the ionised gas., Comment: accepted for publication in A&A

We present the Herschel-SPIRE photometric atlas for a complete flux limited sample of 43 local Ultraluminous Infrared Galaxies (ULIRGs), selected at 60${\mu}$m by IRAS, as part of the HERschel ULIRG Survey (HERUS). Photometry observations were obtained using the SPIRE instrument at 250, 350 and 500${\mu}$m. We describe these observations, present the results, and combine the new observations with data from IRAS to examine the far-IR spectral energy distributions (SEDs) of these sources. We fit the observed SEDs of HERUS objects with a simple parameterised modified black body model where temperature and emissivity $\beta$ are free parameters. We compare the fitted values to those of non-ULIRG local galaxies, and find, in agreement with earlier results, that HERUS ULIRGs have warmer dust (median temperature T = 37.9+/-4.7 K compared to 21.3+/-3.4 K) but a similar $\beta$ distribution (median $\beta$ = 1.7 compared to 1.8) to the Herschel reference sample (HRS, Cortese et al., 2014) galaxies. Dust masses are found to be in the range of 10^7.5 to 10^9 Msun significantly higher than that of Herschel Reference Sample (HRS) sources. We compare our results for local ULIRGs with higher redshift samples selected at 250 and 850${\mu}$m. These latter sources generally have cooler dust and/or redder 100-to-250 ${\mu}$m colours than our 60${\mu}$m-selected ULIRGs. We show that this difference may in part be the result of the sources being selected at different wavelengths rather than being a simple indication of rapid evolution in the properties of the population., Comment: Accepted for publication in MNRAS

The SPICA mid and far-infrared telescope will address fundamental issues in our understanding of star formation and ISM physics in galaxies. A particular hallmark of SPICA is the outstanding sensitivity enabled by the cold telescope, optimized detectors, and wide instantaneous bandwidth throughout the mid- and far-infrared. The spectroscopic, imaging and polarimetric observations that SPICA will be able to collect will help in clarifying the complex physical mechanisms which underlie the baryon cycle of galaxies. In particular: (i) The access to a large suite of atomic and ionic fine-structure lines for large samples of galaxies will shed light on the origin of the observed spread in star formation rates within and between galaxies. (ii) Observations of HD rotational lines (out to $\sim$10 Mpc) and fine structure lines such as [CII] 158 $\mu$m (out to $\sim$100 Mpc) will clarify the main reservoirs of interstellar matter in galaxies, including phases where CO does not emit. (iii) Far-infrared spectroscopy of dust and ice features will address uncertainties in the mass and composition of dust in galaxies, and the contributions of supernovae to the interstellar dust budget will be quantified by photometry and monitoring of supernova remnants in nearby galaxies. (iv) Observations of far-infrared cooling lines such as [OI] 63 $\mu$m from star-forming molecular clouds in our Galaxy will evaluate the importance of shocks to dissipate turbulent energy. The paper concludes with requirements for the telescope and instruments, and recommendations for the observing strategy., Comment: Accepted for publication in PASA; 19 pages, 8 figures, 1 table

The abundance of metals in galaxies is a key parameter which permits to distinguish between different galaxy formation and evolution models. Most of the metallicity determinations are based on optical line ratios. However, the optical spectral range is subject to dust extinction and, for high-z objects (z > 3), some of the lines used in optical metallicity diagnostics are shifted to wavelengths not accessible to ground based observatories. For this reason, we explore metallicity diagnostics using far-infrared (IR) line ratios which can provide a suitable alternative in such situations. To investigate these far-IR line ratios, we modeled the emission of a starburst with the photoionization code CLOUDY. The most sensitive far-IR ratios to measure metallicities are the [OIII]52$\mu$m and 88$\mu$m to [NIII]57$\mu$m ratios. We show that this ratio produces robust metallicities in the presence of an AGN and is insensitive to changes in the age of the ionizing stellar. Another metallicity sensitive ratio is the [OIII]88$\mu$m/[NII]122$\mu$m ratio, although it depends on the ionization parameter. We propose various mid- and far-IR line ratios to break this dependency. Finally, we apply these far-IR diagnostics to a sample of 19 local ultraluminous IR galaxies (ULIRGs) observed with Herschel and Spitzer. We find that the gas-phase metallicity in these local ULIRGs is in the range 0.7 < Z_gas/Z_sun < 1.5, which corresponds to 8.5 < 12 + log (O/H) < 8.9. The inferred metallicities agree well with previous estimates for local ULIRGs and this confirms that they lie below the local mass-metallicity relation., Comment: Accepted for publication in MNRAS; 25 pages, 19 figures

Properties of the cold interstellar medium of low-metallicity galaxies are not well-known due to the faintness and extremely small scale on which emission is expected. We present deep ALMA band 6 (230GHz) observations of the nearby, low-metallicity (12 + log(O/H) = 7.25) blue compact dwarf galaxy SBS0335-052 at an unprecedented resolution of 0.2 arcsec (52 pc). The 12CO J=2-1 line is not detected and we report a 3-sigma upper limit of LCO(2-1) = 3.6x10^4 K km/s pc^2. Assuming that molecular gas is converted into stars with a given depletion time, ranging from 0.02 to 2 Gyr, we find lower limits on the CO-to-H2 conversion factor alpha_CO in the range 10^2-10^4 Msun pc^-2 (K km/s)^-1. The continuum emission is detected and resolved over the two main super star clusters. Re-analysis of the IR-radio spectral energy distribution suggests that the mm-fluxes are not only free-free emission but are most likely also associated with a cold dust component coincident with the position of the brightest cluster. With standard dust properties, we estimate its mass to be as large as 10^5 Msun. Both line and continuum results suggest the presence of a large cold gas reservoir unseen in CO even with ALMA., Comment: corrected typos

(abridged) The dominant thermal mechanisms in the neutral interstellar medium, which acts as a star-forming gas reservoir, are uncertain in extremely metal-poor galaxies. Our objective is to identify the heating mechanisms in one such galaxy, IZw18, and assess the diagnostic value of fine-structure cooling lines. We also seek to constrain the mass of H$_2$, which, despite being an important catalyst and tracer of star formation, remains elusive in this object. Building on a previous photoionization model within a multi-sector topology, we provide additional constraints from the [CII] and [OI] lines and the dust mass recently measured with Herschel. The heating of the HI region appears to be mainly due to photoionization by radiation from a bright X-ray binary source, while photoelectric effect (PE) is negligible. The [CII] and [OI] lines imply an average X-ray luminosity of $4\times10^{40}$ erg s$^{-1}$, while the [NeV] upper limits bring strong constraints to the soft X-ray flux arising from the binary. A negligible amount of H$_2$ is predicted, but $\lesssim10^7$ M$_\odot$ of H$_2$ may be hidden in sufficiently dense clouds of order $\lesssim10$ pc in size. Regardless of the presence of significant amounts of H$_2$, [CII] and [OI] do not trace the so-called CO-dark gas, but the almost purely atomic medium. Although the [CII]+[OI]/TIR ratio is close to values found in more metal-rich sources, it cannot be safely used as a PE heating efficiency proxy. This ratio seems to be kept stable due to a correlation between the X-ray luminosity and the star-formation rate. We propose that X-ray heating is an important process in extremely metal-poor sources. The weak PE heating due to the low dust-to-gas ratio could be compensated for by the larger occurrence and power of X-ray binaries in low-metallicity galaxies. We speculate that X-ray heating may quench star formation., Comment: Accepted for publication in A&A. Language editing completed in version 2

We report on the energetics of molecular outflows in 14 local Ultraluminous Infrared Galaxies (ULIRGs) that show unambiguous outflow signatures (P-Cygni profiles or high-velocity absorption wings) in the far-infrared lines of OH measured with the Herschel/PACS spectrometer. Detection of both ground-state (at 119 and 79 um) and one or more radiatively-excited (at 65 and 84 um) lines allows us to model the nuclear gas (<~300 pc) as well as the more extended components using spherically symmetric radiative transfer models. The highest molecular outflow velocities are found in buried sources, in which slower but massive expansion of the nuclear gas is also observed. With the exception of a few outliers, the outflows have momentum fluxes of (2-5)xL_IR/c and mechanical luminosities of (0.1-0.3)% of L_IR. The moderate momentum boosts in these sources (<~3) suggest that the outflows are mostly momentum-driven by the combined effects of AGN and nuclear starbursts, as a result of radiation pressure, winds, and supernovae remnants. In some sources (~20%), however, powerful (10^{10.5-11} Lsun) AGN feedback and (partially) energy-conserving phases are required, with momentum boosts in the range 3-20. These outflows appear to be stochastic strong-AGN feedback events that occur throughout the merging process. In a few sources, the outflow activity in the innermost regions has subsided in the last ~1 Myr. While OH traces the molecular outflows at sub-kpc scales, comparison of the masses traced by OH with those previously inferred from tracers of more extended outflowing gas suggests that most mass is loaded (with loading factors of Mdot/SFR=1-10) from the central galactic cores (a few x 100 pc). Outflow depletion timescales are <10^8 yr, shorter than the gas consumption timescales by factors of 1.1-15, and are anti-correlated with the AGN luminosity.

Aims. The magnetic star HD148937 is the only Galactic Of?p star surrounded by a nebula. The structure of this nebula is particularly complex and is composed, from the center out outwards, of a close bipolar ejecta nebula (NGC6164/5), an ellipsoidal wind-blown shell, and a spherically symmetric Stromgren sphere. The exact formation process of this nebula and its precise relation to the star's evolution remain unknown. Methods. We analyzed infrared Spitzer IRS and far-infrared Herschel/PACS observations of the NGC6164/5 nebula. The Herschel imaging allowed us to constrain the global morphology of the nebula. We also combined the infrared spectra with optical spectra of the central star to constrain its evolutionary status. We used these data to derive the abundances in the ejected material. To relate this information to the evolutionary status of the star, we also determined the fundamental parameters of HD148937 using the CMFGEN atmosphere code. Results. The H$\alpha$ image displays a bipolar or "8"-shaped ionized nebula, whilst the infrared images show dust to be more concentrated around the central object. We determine nebular abundance ratios of N/O = 1.06 close to the star, and N/O = 1.54 in the bright lobe constituting NGC6164. Interestingly, the parts of the nebula located further from HD148937 appear more enriched in stellar material than the part located closer to the star. Evolutionary tracks suggest that these ejecta have occured $\sim$1.2-1.3 and $\sim$0.6 Myrs ago, respectively. In addition, we derive abundances of argon for the nebula compatible with the solar values and we find a depletion of neon and sulfur. The combined analyses of the known kinematics and of the new abundances of the nebula suggest either a helical morphology for the nebula, possibly linked to the magnetic geometry, or the occurrence of a binary merger., Comment: 15 pages, 3 pages of appendix, 12 figures, A&A in press

More complete knowledge of galaxy evolution requires understanding the process of star formation and interaction between the interstellar radiation field and the interstellar medium in galactic environments traversing a wide range of physical parameter space. Here we focus on the impact of massive star formation on the surrounding low metallicity ISM in 30 Doradus in the Large Magellanic Cloud. A low metal abundance, as is the case of some galaxies of the early universe, results in less ultra-violet shielding for the formation of the molecular gas necessary for star formation to proceed. The half-solar metallicity gas in this region is strongly irradiated by the super star cluster R136, making it an ideal laboratory to study the structure of the ISM in an extreme environment. Our spatially resolved study investigates the gas heating and cooling mechanisms, particularly in the photo-dissociation regions where the chemistry and thermal balance are regulated by far-ultraviolet photons (6 eV< h\nu <13.6 eV). We present Herschel observations of far-infrared fine-structure lines obtained with PACS and SPIRE/FTS. We have combined atomic fine-structure lines from Herschel and Spitzer observations with ground-based CO data to provide diagnostics on the properties and the structure of the gas by modeling it with the Meudon PDR code. We derive the spatial distribution of the radiation field, the pressure, the size, and the filling factor of the photodissociated gas and molecular clouds. We find a range of pressure of ~ 10^5 - 1.7x10^6 cm^{-3} K and a range of incident radiation field G_UV ~ 10^2 - 2.5x10^4 through PDR modeling. Assuming a plane-parallel geometry and a uniform medium, we find a total extinction of 1-3 mag , which correspond to a PDR cloud size of 0.2 to 3pc, with small CO depth scale of 0.06 to 0.5pc. We also determine the three dimensional structure of the gas. (Abridged), Comment: 20 pages, 23 figures, accepted in A&A

We combine Spitzer and Herschel data of the star-forming region N11 in the Large Magellanic Cloud to produce detailed maps of the dust properties in the complex and study their variations with the ISM conditions. We also compare APEX/LABOCA 870um observations with our model predictions in order to decompose the 870um emission into dust and non-dust (free-free emission and CO(3-2) line) contributions. We find that in N11, the 870um can be fully accounted for by these 3 components. The dust surface density map of N11 is combined with HI and CO observations to study local variations in the gas-to-dust mass ratios. Our analysis leads to values lower than those expected from the LMC low-metallicity as well as to a decrease of the gas-to-dust mass ratio with the dust surface density. We explore potential hypotheses that could explain the low observed gas-to-dust mass ratios (variations in the XCO factor, presence of CO-dark gas or of optically thick HI or variations in the dust abundance in the dense regions). We finally decompose the local SEDs using a Principal Component Analysis (i.e. with no a priori assumption on the dust composition in the complex). Our results lead to a promising decomposition of the local SEDs in various dust components (hot, warm, cold) coherent with that expected for the region. Further analysis on a larger sample of galaxies will follow in order to understand how unique this decomposition is or how it evolves from one environment to another., Comment: 24 pages, 16 figures, accepted for publication in MNRAS

The surroundings of HII regions can have a profound influence on their development, morphology, and evolution. This paper explores the effect of the environment on H II regions in the MonR2 molecular cloud. We aim to investigate the density structure of envelopes surrounding HII regions and to determine their collapse and ionisation expansion ages. The Mon R2 molecular cloud is an ideal target since it hosts an H II region association. Column density and temperature images derived from Herschel data were used together to model the structure of HII bubbles and their surrounding envelopes. The resulting observational constraints were used to follow the development of the Mon R2 ionised regions with analytical calculations and numerical simulations. The four hot bubbles associated with H II regions are surrounded by dense, cold, and neutral gas envelopes. The radial density profiles are reminiscent of those of low-mass protostellar envelopes. The inner parts of envelopes of all four HII regions could be free-falling because they display shallow density profiles. As for their outer parts, the two compact HII regions show a density profile, which is typical of the equilibrium structure of an isothermal sphere. In contrast, the central UCHii region shows a steeper outer profile, that could be interpreted as material being forced to collapse. The size of the heated bubbles, the spectral type of the irradiating stars, and the mean initial neutral gas density are used to estimate the ionisation expansion time, texp, 0.1Myr,for the dense UCHII and compact HII regions and 0.35 Myr for the extended HII region. The envelope transition radii between the shallow and steeper density profiles are used to estimate the time elapsed since the formation of the first proto stellar embryo, Tinf : 1Myr, for the ultra-compact, 1.5 / 3Myr for the compact, and greater than 6Myr for the extended HII regions., Comment: 25pages A&A sous presse

In this work, we aim at providing a consistent analysis of the dust properties from metal-poor to metal-rich environments by linking them to fundamental galactic parameters. We consider two samples of galaxies: the Dwarf Galaxy Survey (DGS) and KINGFISH, totalling 109 galaxies, spanning almost 2 dex in metallicity. We collect infrared (IR) to submillimetre (submm) data for both samples and present the complete data set for the DGS sample. We model the observed spectral energy distributions (SED) with a physically-motivated dust model to access the dust properties. Using a different SED model (modified blackbody), dust composition (amorphous carbon), or wavelength coverage at submm wavelengths results in differences in the dust mass estimate of a factor two to three, showing that this parameter is subject to non-negligible systematic modelling uncertainties. For eight galaxies in our sample, we find a rather small excess at 500 microns (< 1.5 sigma). We find that the dust SED of low-metallicity galaxies is broader and peaks at shorter wavelengths compared to more metal-rich systems, a sign of a clumpier medium in dwarf galaxies. The PAH mass fraction and the dust temperature distribution are found to be driven mostly by the specific star-formation rate, SSFR, with secondary effects from metallicity. The correlations between metallicity and dust mass or total-IR luminosity are direct consequences of the stellar mass-metallicity relation. The dust-to-stellar mass ratios of metal-rich sources follow the well-studied trend of decreasing ratio for decreasing SSFR. The relation is more complex for highly star-forming low-metallicity galaxies and depends on the chemical evolutionary stage of the source (i.e., gas-to-dust mass ratio). Dust growth processes in the ISM play a key role in the dust mass build-up with respect to the stellar content at high SSFR and low metallicity. (abridged), Comment: 44 pages (20 pages main body plus 5 Appendices), 11 figures, 9 tables, accepted for publication in A&A

The Infrared Spectrograph (IRS) on board the Spitzer Space Telescope observed about 15,000 objects during the cryogenic mission lifetime. Observations provided low-resolution (R~60-127) spectra over ~5-38um and high-resolution (R~600) spectra over ~10-37um. The Cornell Atlas of Spitzer/IRS Sources (CASSIS) was created to provide publishable quality spectra to the community. Low-resolution spectra have been available in CASSIS since 2011, and we present here the addition of the high-resolution spectra. The high-resolution observations represent approximately one third of all staring observations performed with the IRS instrument. While low-resolution observations are adapted to faint objects and/or broad spectral features (e.g., dust continuum, molecular bands), high-resolution observations allow more accurate measurements of narrow features (e.g., ionic emission lines) as well as a better sampling of the spectral profile of various features. Given the narrow aperture of the two high-resolution modules, cosmic ray hits and spurious features usually plague the spectra. Our pipeline is designed to minimize these effects through various improvements. A super sampled point-spread function was created in order to enable the optimal extraction in addition to the full aperture extraction. The pipeline selects the best extraction method based on the spatial extent of the object. For unresolved sources, the optimal extraction provides a significant improvement in signal-to-noise ratio over a full aperture extraction. We have developed several techniques for optimal extraction, including a differential method that eliminates low-level rogue pixels (even when no dedicated background observation was performed). The updated CASSIS repository now includes all the spectra ever taken by the IRS, with the exception of mapping observations., Comment: Accepted on May 11th 2015 in ApJS

The far-infrared (FIR) lines are key tracers of the physical conditions of the interstellar medium (ISM) and are becoming workhorse diagnostics for galaxies throughout the universe. Our goal is to explain the differences and trends observed in the FIR line emission of dwarf galaxies compared to more metal-rich galaxies. We present Herschel PACS spectroscopic observations of the CII157um, OI63 and 145um, OIII88um, NII122 and 205um, and NIII57um fine-structure cooling lines in a sample of 48 low-metallicity star-forming galaxies of the guaranteed time key program Dwarf Galaxy Survey. We correlate PACS line ratios and line-to-LTIR ratios with LTIR, LTIR/LB, metallicity, and FIR color, and interpret the observed trends in terms of ISM conditions and phase filling factors with Cloudy radiative transfer models. We find that the FIR lines together account for up to 3 percent of LTIR and that star-forming regions dominate the overall emission in dwarf galaxies. Compared to metal-rich galaxies, the ratios of OIII/NII122 and NIII/NII122 are high, indicative of hard radiation fields. In the photodissociation region (PDR), the CII/OI63 ratio is slightly higher than in metal-rich galaxies, with a small increase with metallicity, and the OI145/OI63 ratio is generally lower than 0.1, demonstrating that optical depth effects should be small on the scales probed. The OIII/OI63 ratio can be used as an indicator of the ionized gas/PDR filling factor, and is found ~4 times higher in the dwarfs than in metal-rich galaxies. The high CII/LTIR, OI/LTIR, and OIII/LTIR ratios, which decrease with increasing LTIR and LTIR/LB, are interpreted as a combination of moderate FUV fields and low PDR covering factor. Harboring compact phases of low filling factor and a large volume filling factor of diffuse gas, the ISM of low-metallicity dwarf galaxies has a more porous structure than that in metal-rich galaxies., Comment: Accepted for publication in A&A

Context: The north-west photo-dissociation region (PDR) in the reflection nebula NGC 7023 displays a complex structure. Filament-like condensations at the edge of the cloud can be traced via the emission of the main cooling lines, offering a great opportunity to study the link between the morphology and energetics of these regions. Aims: We study the spatial variation of the far-infrared fine-structure lines of [C II] (158 um) and [O I] (63 and 145 um). These lines trace the local gas conditions across the PDR. Methods: We used observations from the Herschel/PACS instrument to map the spatial distribution of these fine-structure lines. The observed region covers a square area of about 110" x 110" with an angular resolution that varies from 4" to 11". We compared this emission with ground-based and Spitzer observations of H2 lines, Herschel/SPIRE observations of CO lines, and Spitzer/IRAC 3.6 um images that trace the emission of polycyclic aromatic hydrocarbons. Results: The [C II] (158 um) and [O I] (63 and 145 um) lines arise from the warm cloud surface where the PDR is located and the gas is warm, cooling the region. We find that although the relative contribution to the cooling budget over the observed region is dominated by [O I]63 um (>30%), H2 contributes significantly in the PDR (35%), as does [C II]158 um outside the PDR (30%). Other species contribute little to the cooling ([O I]145 um 9%, and CO 4%). The [O I] maps resolve these condensations into two structures and show that the peak of [O I] is slightly displaced from the molecular H2 emission. The size of these structures is about 8" (0.015 pc) and in surface cover about 9% of the PDR emission. Finally, we did not detect emission from [N II]122 um, suggesting that the cavity is mostly filled with non-ionised gas., Comment: 10 pages, 7 figures, 1 Table

The rate at which interstellar gas is converted into stars, and its dependence on environment, is one of the pillars on which our understanding of the visible Universe is build. We present a comparison of the surface density of young stars (Sigma_*) and dust surface density (Sigma_d) across NGC346 (N66) in 115 independent pixels of 6x6 pc^2. We find a correlation between Sigma_* and Sigma_d with a considerable scatter. A power law fit to the data yields a steep relation with an exponent of 2.6+-0.2. We convert Sigma_d to gas surface density (Sigma_g) and Sigma_* to star formation rate (SFR) surface densities (Sigma_SFR), using simple assumptions for the gas-to-dust mass ratio and the duration of star formation. The derived total SFR (4+-1 10^-3 M_sun/yr) is consistent with SFR estimated from the Ha emission integrated over the Ha nebula. On small scales the Sigma_SFR derived using Ha systematically underestimates the count-based Sigma_SFR, by up to a factor of 10. This is due to ionizing photons escaping the area, where the stars are counted. We find that individual 36 pc^2 pixels fall systematically above integrated disc-galaxies in the Schmidt-Kennicutt diagram by on average a factor of ~7. The NGC346 average SFR over a larger area (90 pc radius) lies closer to the relation but remains high by a factor of ~3. The fraction of the total mass (gas plus young stars) locked in young stars is systematically high (~10 per cent) within the central 15 pc and systematically lower outside (2 per cent), which we interpret as variations in star formation efficiency. The inner 15 pc is dominated by young stars belonging to a centrally condensed cluster, while the outer parts are dominated by a dispersed population. Therefore, the observed trend could reflect a change of star formation efficiency between clustered and non-clustered star-formation., Comment: Accepted for publication in MNRAS

We examine the relation between polycyclic aromatic hydrocarbon (PAH) emission at 8 microns and far-infrared emission from hot dust grains at 24 microns and from large dust grains at 160 and 250 microns in the nearby spiral galaxies NGC 2403 and M83 using data from the Spitzer Space Telescope and Herschel Space Observatory. We find that the PAH emission in NGC 2403 is better correlated with emission at 250 microns from dust heated by the diffuse interstellar radiation field (ISRF) and that the 8/250 micron surface brightness ratio is well-correlated with the stellar surface brightness as measured at 3.6 microns. This implies that the PAHs in NGC 2403 are intermixed with cold large dust grains in the diffuse interstellar medium (ISM) and that the PAHs are excited by the diffuse ISRF. In M83, the PAH emission appears more strongly correlated with 160 micron emission originating from large dust grains heated by star forming regions. However, the PAH emission in M83 is low where the 24 micron emission peaks within star forming regions, and enhancements in the 8/160 micron surface brightness ratios appear offset relative to the dust and the star forming regions within the spiral arms. This suggests that the PAHs observed in the 8 micron band are not excited locally within star forming regions but either by light escaping non-axisymmetrically from star forming regions or locally by young, non-photoionising stars that have migrated downstream from the spiral density waves. The results from just these two galaxies show that PAHs may be excited by different stellar populations in different spiral galaxies., Comment: 21 pages, 19 figures, accepted for publication in MNRAS

We present new Herschel photometric and spectroscopic observations of Supernova 1987A, carried out in 2012. Our dedicated photometric measurements provide new 70 micron data and improved imaging quality at 100 and 160 micron compared to previous observations in 2010. Our Herschel spectra show only weak CO line emission, and provide an upper limit for the 63 micron [O I] line flux, eliminating the possibility that line contaminations distort the previously estimated dust mass. The far-infrared spectral energy distribution (SED) is well fitted by thermal emission from cold dust. The newly measured 70 micron flux constrains the dust temperature, limiting it to nearly a single temperature. The far-infrared emission can be fitted by 0.5+-0.1 Msun of amorphous carbon, about a factor of two larger than the current nucleosynthetic mass prediction for carbon. The observation of SiO molecules at early and late phases suggests that silicates may also have formed and we could fit the SED with a combination of 0.3 Msun of amorphous carbon and 0.5 Msun of silicates, totalling 0.8 Msun of dust. Our analysis thus supports the presence of a large dust reservoir in the ejecta of SN 1987A. The inferred dust mass suggests that supernovae can be an important source of dust in the interstellar medium, from local to high-redshift galaxies., Comment: ApJ accepted, 8 pages

We present Herschel PACS and SPIRE spectroscopy of the most important far-infrared cooling lines in the nearby edge-on spiral galaxy, NGC 891: [CII] 158 $\mu$m, [NII] 122, 205 $\mu$m, [OI] 63, 145 $\mu$m, and [OIII] 88 $\mu$m. We find that the photoelectric heating efficiency of the gas, traced via the ([CII]+[OII]63)/$F_{\mathrm{TIR}}$ ratio, varies from a mean of 3.5$\times$10$^{-3}$ in the centre up to 8$\times$10$^{-3}$ at increasing radial and vertical distances in the disc. A decrease in ([CII]+[OII]63)/$F_{\mathrm{TIR}}$ but constant ([CII]+[OI]63)/$F_{\mathrm{PAH}}$ with increasing FIR colour suggests that polycyclic aromatic hydrocarbons (PAHs) may become important for gas heating in the central regions. We compare the observed flux of the FIR cooling lines and total IR emission with the predicted flux from a PDR model to determine the gas density, surface temperature and the strength of the incident far-ultraviolet (FUV) radiation field, $G_{0}$. Resolving details on physical scales of ~0.6 kpc, a pixel-by-pixel analysis reveals that the majority of the PDRs in NGC 891's disc have hydrogen densities of 1 < log ($n$/cm$^{-3}$) < 3.5 experiencing an incident FUV radiation field with strengths of 1.7 < log $G_0$ < 3. Although these values we derive for most of the disc are consistent with the gas properties found in PDRs in the spiral arms and inter-arm regions of M51, observed radial trends in $n$ and $G_0$ are shown to be sensitive to varying optical thickness in the lines, demonstrating the importance of accurately accounting for optical depth effects when interpreting observations of high inclination systems. With an empirical relationship between the MIPS 24 $\mu$m and [NII] 205 $\mu$m emission, we estimate an enhancement of the FUV radiation field strength in the far north-eastern side of the disc., Comment: Accepted for publication in A&A. 25 pages, including 17 figures and 3 tables, abstract abridged for arXiv

We search for variations in the disk of Centaurus A of the emission from atomic fine structure lines using Herschel PACS and SPIRE spectroscopy. In particular we observe the [C II](158 $\mu$m), [N II](122 and 205 $\mu$m), [O I](63 and 145 $\mu$m) and [O III](88 $\mu$m) lines, which all play an important role in cooling the gas in photo-ionized and photodissociation regions. We determine that the ([C II]+[O I]$_{63}$)/$F_{TIR}$ line ratio, a proxy for the heating efficiency of the gas, shows no significant radial trend across the observed region, in contrast to observations of other nearby galaxies. We determine that 10 - 20% of the observed [C II] emission originates in ionized gas. Comparison between our observations and a PDR model shows that the strength of the far-ultraviolet radiation field, $G_0$, varies between $10^{1.75}$ and $10^{2.75}$ and the hydrogen nucleus density varies between $10^{2.75}$ and $10^{3.75}$ cm$^{-3}$, with no significant radial trend in either property. In the context of the emission line properties of the grand-design spiral galaxy M51 and the elliptical galaxy NGC 4125, the gas in Cen A appears more characteristic of that in typical disk galaxies rather than elliptical galaxies., Comment: Accepted for publication in the Astrophysical Journal. 22 pages, 10 figures, 5 tables

We aim to quantify the molecular gas reservoir in a subset of 6 low-metallicity galaxies from the Herschel Dwarf Galaxy Survey with newly acquired CO data, and link this reservoir to the observed star formation activity. We present CO(1-0), CO(2-1), and CO(3-2) observations obtained at the ATNF Mopra 22-m, APEX, and IRAM 30-m telescopes, as well as [CII] 157um and [OI] 63um observations obtained with the Herschel/PACS spectrometer in the 6 galaxies: Haro11, Mrk1089, Mrk930, NGC4861, NGC625, and UM311. We derive molecular gas mass from several methods including the use of the CO-to-H2 conversion factor Xco (both Galactic and metallicity-scaled values) and of dust measurements. The molecular and atomic gas reservoirs are compared to the star formation activity. We also constrain the physical conditions of the molecular clouds using the non-LTE code RADEX and the spectral synthesis code Cloudy. We detect CO in 5 of the 6 galaxies, including first detections in Haro11 (Z~0.4 Zsun), Mrk930 (0.2 Zsun), and UM311 (0.5 Zsun), but CO remains undetected in NGC4861 (0.2 Zsun). The CO luminosities are low while [CII] is bright in these galaxies, resulting in [CII]/CO(1-0)>10000. Our dwarf galaxies are in relatively good agreement with the Schmidt-Kennicutt relation for total gas. They show short molecular depletion time scales, even when considering metallicity-scaled Xco factors. Those galaxies are dominated by their HI gas, except Haro11 which has high star formation efficiency and is dominated by ionized and molecular gas. We determine the mass of each ISM phase in Haro11 using Cloudy and estimate an equivalent Xco factor which is 10 times higher than the Galactic value. Overall, our results confirm the emerging picture that CO suffers from significant selective photodissociation in low-metallicity dwarf galaxies., Comment: 19 pages, 14 figures. Accepted for publication in A&A

This paper analyses the behaviour of the gas-to-dust mass ratio (G/D) of local Universe galaxies over a large metallicity range. We combine three samples: the Dwarf Galaxy Survey, the KINGFISH survey and a subsample from Galametz et al. (2011) totalling 126 galaxies, covering a 2 dex metallicity range, with 30% of the sample with 12+log(O/H) < 8.0. The dust masses are homogeneously determined with a semi-empirical dust model, including submm constraints. The atomic and molecular gas masses are compiled from the literature. Two XCO are used to estimate molecular gas masses: the Galactic XCO, and a XCO depending on the metallicity (as Z^{-2}). Correlations with morphological types, stellar masses, star formation rates and specific star formation rates are discussed. The trend between G/D and metallicity is empirically modelled using power-laws (slope of -1 and free) and a broken power-law. We compare the evolution of the G/D with predictions from chemical evolution models. We find that out of the five tested galactic parameters, metallicity is the galactic property driving the observed G/D. The G/D versus metallicity relation cannot be represented by a power-law with a slope of -1 over the whole metallicity range. The observed trend is steeper for metallicities lower than ~ 8.0. A large scatter is observed in the G/D for a given metallicity, with a dispersion of 0.37 dex in metallicity bins of ~0.1 dex. The broken power-law reproduces best the observed G/D and provides estimates of the G/D that are accurate to a factor of 1.6. The good agreement of the G/D and its scatter with the three tested chemical evolution models shows that the scatter is intrinsic to galactic properties, reflecting the different star formation histories, dust destruction efficiencies, dust grain size distributions and chemical compositions across the sample. (abriged), Comment: 23 pages, 12 figures, accepted in Astronomy & Astrophysics

The mid-infrared (MIR) spectra observed with the \textit{Spitzer} Infrared Spectrograph (IRS) provide a valuable dataset for untangling the physical processes and conditions within galaxies. This paper presents the first attempt to blindly learn fundamental spectral components of MIR galaxy spectra, using non-negative matrix factorisation (NMF). NMF is a recently developed multivariate technique shown to be successful in blind source separation problems. Unlike the more popular multivariate analysis technique, principal component analysis, NMF imposes the condition that weights and spectral components are non-negative. This more closely resembles the physical process of emission in the mid-infrared, resulting in physically intuitive components. By applying NMF to galaxy spectra in the Cornell Atlas of Spitzer/IRS sources (CASSIS), we find similar components amongst different NMF sets. These similar components include two for AGN emission and one for star formation. [... ABBREVIATED...] We show an NMF set with seven components can reconstruct the general spectral shape of a wide variety of objects, though struggle to fit the varying strength of emission lines. We also show that the seven components can be used to separate out different types of objects. We model this separation with Gaussian Mixtures modelling and use the result to provide a classification tool. We also show the NMF components can be used to separate out the emission from AGN and star formation regions and define a new star formation/AGN diagnostic which is consistent with all mid-infrared diagnostics already in use but has the advantage that it can be applied to mid-infrared spectra with low signal to noise or with limited spectral range. The 7 NMF components and code for classification are made public on arxiv and are available at: \url{https://github.com/pdh21/NMF_software/}, Comment: 23 pages, 19 figures, accepted for publication in MNRAS

We present Herschel photometry and spectroscopy, carried out as part of the Herschel ULIRG survey (HERUS), and a model for the infrared to submillimetre emission of the ultraluminous infrared galaxy IRAS 08572+3915. This source shows one of the deepest known silicate absorption features and no polycyclic aromatic hydrocarbon (PAH) emission. The model suggests that this object is powered by an active galactic nucleus (AGN) with a fairly smooth torus viewed almost edge-on and a very young starburst. According to our model the AGN contributes about 90% of the total luminosity of 1.1 x 10^13 Lo, which is about a factor of five higher than previous estimates. The large correction of the luminosity is due to the anisotropy of the emission of the best fit torus. Similar corrections may be necessary for other local and high-z analogs. This correction implies that IRAS 08572+3915 at a redshift of 0.05835 may be the nearest hyperluminous infrared galaxy and probably the most luminous infrared galaxy in the local (z < 0.2) Universe. IRAS 08572+3915 shows a low ratio of [CII] to IR luminosity (log L_[CII]/L_{IR} < -3.8) and a [OI]63um to [CII]158um line ratio of about 1 that supports the model presented in this letter., Comment: MNRAS Letters accepted

Data from the Herschel Space Observatory have revealed an unusual elliptical galaxy, NGC 4125, which has strong and extended submillimeter emission from cold dust but only very strict upper limits to its CO and HI emission. Depending on the dust emissivity, the total dust mass is 2-5x10^6 Msun. While the neutral gas-to-dust mass ratio is extremely low (< 12-30), including the ionized gas traced by [CII] emission raises this limit to < 39-100. The dust emission follows a similar r^{1/4} profile to the stellar light and the dust to stellar mass ratio is towards the high end of what is found in nearby elliptical galaxies. We suggest that NGC 4125 is currently in an unusual phase where evolved stars produced in a merger-triggered burst of star formation are pumping large amounts of gas and dust into the interstellar medium. In this scenario, the low neutral gas-to-dust mass ratio is explained by the gas being heated to temperatures >= 10^4 K faster than the dust is evaporated. If galaxies like NGC 4125, where the far-infrared emission does not trace neutral gas in the usual manner, are common at higher redshift, this could have significant implications for our understanding of high redshift galaxies and galaxy evolution., Comment: Accepted to ApJ Letters

We present new photometric data from our Herschel Key Programme, the Dwarf Galaxy Survey (DGS), dedicated to the observation of the gas and dust in 48 low-metallicity environments. They were observed with PACS and SPIRE onboard Herschel at 70,100,160,250,350, and 500 microns. We focus on a systematic comparison of the derived FIR properties (FIR luminosity, dust mass, dust temperature and emissivity index) with more metal-rich galaxies and investigate the detection of a potential submm excess. The data reduction method is adapted for each galaxy to derive the most reliable photometry from the final maps. PACS flux densities are compared with the MIPS 70 and 160 microns bands. We use colour-colour diagrams and modified blackbody fitting procedures to determine the dust properties of the DGS galaxies. We also include galaxies from the Herschel KINGFISH sample, containing more metal-rich environments, totalling 109 galaxies. The location of the DGS galaxies on Herschel colour-colour diagrams highlights the differences in global environments of low-metallicity galaxies. The dust in DGS galaxies is generally warmer than in KINGFISH galaxies (T_DGS~32 K, T_KINGFISH~23 K). The emissivity index, beta, is ~1.7 in the DGS, but metallicity does not make a strong effect on beta. The dust-to-stellar mass ratio is lower in low-metallicity galaxies: M_dust/M_star~0.02% for the DGS vs 0.1% for KINGFISH. Per unit dust mass, dwarf galaxies emit ~6 times more in the FIR than higher metallicity galaxies. Out of the 22 DGS galaxies detected at 500 micron, 41% present an excess in the submm not explained by our dust SED model. The excess mainly appears in lower metallicity galaxies (12+log(O/H) < 8.3), and the strongest excesses are detected in the most metal-poor galaxies. We stress the need for observations longwards of the Herschel wavelengths to detect any submm excess appearing beyond 500 micron., Comment: 32 pages, 20 figures, accepted in A&A

We present Herschel PACS and SPIRE spectroscopy of the most important far-infrared cooling lines in M51, [CII](158 \mu m), [NII](122 & 205 \mu m), [OI](63 and 145 \mu m) and [OIII](88 \mu m). We compare the observed flux of these lines with the predicted flux from a photon dominated region model to determine characteristics of the cold gas such as density, temperature and the far-ultraviolet radiation field, G_0, resolving details on physical scales of roughly 600 pc. We find an average [CII]/F_TIR of 4 x 10^{-3}, in agreement with previous studies of other galaxies. A pixel-by-pixel analysis of four distinct regions of M51 shows a radially decreasing trend in both the far-ultraviolet (FUV) radiation field, G_0 and the hydrogen density, n, peaking in the nucleus of the galaxy, then falling off out to the arm and interarm regions. We see for the first time that the FUV flux and gas density are similar in the differing environments of the arm and interarm regions, suggesting that the inherent physical properties of the molecular clouds in both regions are essentially the same., Comment: 24 pages, 7 figures, 9 tables; accepted for publication in the Astrophysical Journal

We report on our observations of the 79 and 119um doublet transitions of OH for 24 local (z<0.262) ULIRGs observed with Herschel-PACS as part of the Herschel ULIRG Survey (HERUS). Some OH119 profiles display a clear P-Cygni shape and therefore imply outflowing OH gas, other profiles are predominantly in absorption or are completely in emission. We find that the relative strength of the OH emission component decreases as the silicate absorption increases. This locates the OH outflows inside the obscured nuclei. The maximum outflow velocities for our sources range from less than 100 to 2000 km/s, with 15/24 (10/24) sources showing OH absorption at velocities exceeding 700 km/s (1000 km/s). Three sources show maximum OH outflow velocities exceeding that of Mrk231. Since outflow velocities above 500-700 km/s are thought to require an active galactic nucleus (AGN) to drive them, about 2/3 of our ULIRG sample may host AGN-driven molecular outflows. This finding is supported by the correlation we find between the maximum OH outflow velocity and the IR-derived bolometric AGN luminosity. No such correlation is found with the IR-derived star formation rate. The highest outflow velocities are found among sources which are still deeply embedded. We speculate that the molecular outflows in these sources may be in an early phase of disrupting the nuclear dust veil before these sources evolve into less obscured AGN. Four of our sources show high-velocity wings in their [C II] fine-structure line profiles implying neutral gas outflow masses of at least 2-4.5 x 10^8 Msun., Comment: 23 pages, 12 figures, 4 tables

The Dwarf Galaxy Survey (DGS) program is studying low-metallicity galaxies using 230h of far-infrared (FIR) and submillimetre (submm) photometric and spectroscopic observations of the Herschel Space Observatory and draws to this a rich database of a wide range of wavelengths tracing the dust, gas and stars. This sample of 50 galaxies includes the largest metallicity range achievable in the local Universe including the lowest metallicity (Z) galaxies, 1/50 Zsun, and spans 4 orders of magnitude in star formation rates. The survey is designed to get a handle on the physics of the interstellar medium (ISM) of low metallicity dwarf galaxies, especially on their dust and gas properties and the ISM heating and cooling processes. The DGS produces PACS and SPIRE maps of low-metallicity galaxies observed at 70, 100, 160, 250, 350, and 500 mic with the highest sensitivity achievable to date in the FIR and submm. The FIR fine-structure lines, [CII] 158 mic, [OI] 63 mic, [OI] 145 mic, [OIII] 88 mic, [NIII] 57 mic and [NII] 122 and 205 mic have also been observed with the aim of studying the gas cooling in the neutral and ionized phases. The SPIRE FTS observations include many CO lines (J=4-3 to J=13-12), [NII] 205 mic and [CI] lines at 370 and 609 mic. This paper describes the sample selection and global properties of the galaxies, the observing strategy as well as the vast ancillary database available to complement the Herschel observations. The scientific potential of the full DGS survey is described with some example results included., Comment: Article accepted for PASP (april 25 2013)

We investigate star formation and dust heating in the compact FIR bright sources detected in the Herschel maps of M83. We use the source extraction code GETSOURCES to detect and extract sources in the FIR, as well as their photometry in the MIR and H{\alpha}. By performing infrared SED fitting and applying an H{\alpha} based star formation rate (SFR) calibration, we derive the dust masses and temperatures, SFRs, gas masses and star formation efficiencies (SFEs). The detected sources lie exclusively on the spiral arms and represent giant molecular associations (GMAs), with gas masses and sizes of 10^6-10^8 solar masses and 200-300 pc, respectively. The inferred parameters show little to no radial dependence and there is only a weak correlation between the SFRs and gas masses, which suggests that more massive clouds are less efficient at forming stars. Dust heating is mainly due to local star formation. However, although the sources are not optically thick, the total intrinsic young stellar population luminosity can almost completely account for the dust luminosity. This suggests that other radiation sources contribute to the dust heating as well and approximately compensate for the unabsorbed fraction of UV light., Comment: 26 pages, 20 figures, accepted to MNRAS

[abridged] Aims. The aim of the present paper is to provide new and more detailed relations at the kpc scale between the gas surface density and the face-on optical depth directly calibrated on galaxies, in order to compute the attenuation not only for semi-analytic models but also observationally as new and upcoming radio observatories are able to trace gas ever farther in the Universe. Methods. We have selected a sample of 4 nearby resolved galaxies and a sample of 27 unresolved galaxies from the Herschel Reference Survey and the Very Nearby Galaxies Survey, for which we have a large set of multi-wavelength data from the FUV to the FIR including metallicity gradients for resolved galaxies, along with radio HI and CO observations. For each pixel in resolved galaxies and for each galaxy in the unresolved sample, we compute the face-on optical depth from the attenuation determined with the CIGALE SED fitting code and an assumed geometry. We determine the gas surface density from HI and CO observations with a metallicity-dependent XCO factor. Results. We provide new, simple to use, relations to determine the face-on optical depth from the gas surface density, taking the metallicity into account, which proves to be crucial for a proper estimate. The method used to determine the gas surface density or the face-on optical depth has little impact on the relations except for galaxies that have an inclination over 50d. Finally, we provide detailed instructions on how to compute the attenuation practically from the gas surface density taking into account possible information on the metallicity. Conclusions. Examination of the influence of these new relations on simulated FUV and IR luminosity functions shows a clear impact compared to older oft-used relations, which in turn could affect the conclusions drawn from studies based on large scale cosmological simulations., Comment: 24 pages, 21 figures, accepted for publication in A&A

[Abridged] We present an integrated photometric spectral energy distribution (SED) of the Magellanic-type galaxy NGC 4449 from the far-ultraviolet (UV) to the submillimetre, including new observations acquired by the Herschel Space Observatory. We include integrated UV photometry from the Swift Ultraviolet and Optical Telescope using a measurement technique which is appropriate for extended sources with coincidence loss. In this paper, we examine the available multiwavelength data to infer a range of ages, metallicities and star formation rates for the underlying stellar populations, as well as the composition and the total mass of dust in NGC 4449. We present an iterative scheme, which allows us to build an in-depth and multicomponent representation of NGC 4449 `bottom-up', taking advantage of the broad capabilities of the photoionization and radiative transfer code MOCASSIN (MOnte CArlo SimulationS of Ionized Nebulae). We fit the observed SED, the global ionization structure and the emission line intensities, and infer a recent SFR of 0.4 Msolar/yr and a total stellar mass of approximately 1e9 Msolar emitting with a bolometric luminosity of 5.7e9 Lsolar. Our fits yield a total dust mass of 2.9e6 Msolar including 2 per cent attributed to polycyclic aromatic hydrocarbons. We deduce a dust to gas mass ratio of 1/190 within the modelled region. While we do not consider possible additional contributions from even colder dust, we note that including the extended HI envelope and the molecular gas is likely to bring the ratio down to as low as ~ 1/800., Comment: 22 pages, 11 figures, 8 tables, published in MNRAS, 431, 2493 (2013)

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., Comment: Accepted for publication in MNRAS; 23 pages, 11 figures, 2 tables

(Abridged) The Magellanic Clouds provide a nearby laboratory for metal-poor dwarf galaxies. The low dust abundance enhances the penetration of UV photons into the interstellar medium (ISM), resulting in a relatively larger filling factor of the ionized gas. Furthermore, there is likely a hidden molecular gas reservoir probed by the [CII]157um line. We present Herschel/PACS maps in several tracers, [CII], [OI]63um,145um, [NII]122um, [NIII]57um, and [OIII]88um in the HII region N11B in the Large Magellanic Cloud. Halpha and [OIII]5007A images were used as complementary data to investigate the effect of dust extinction. Observations were interpreted with photoionization models to infer the gas conditions and estimate the ionized gas contribution to the [CII] emission. Photodissociation regions (PDRs) are probed through polycyclic aromatic hydrocarbons (PAHs). We first study the distribution and properties of the ionized gas. We then constrain the origin of [CII]157um by comparing to tracers of the low-excitation ionized gas and of PDRs. [OIII] is dominated by extended emission from the high-excitation diffuse ionized gas; it is the brightest far-infrared line, ~4 times brighter than [CII]. The extent of the [OIII] emission suggests that the medium is rather fragmented, allowing far-UV photons to permeate into the ISM to scales of >30pc. Furthermore, by comparing [CII] with [NII], we find that 95% of [CII] arises in PDRs, except toward the stellar cluster for which as much as 15% could arise in the ionized gas. We find a remarkable correlation between [CII]+[OI] and PAH emission, with [CII] dominating the cooling in diffuse PDRs and [OI] dominating in the densest PDRs. The combination of [CII] and [OI] provides a proxy for the total gas cooling in PDRs. Our results suggest that PAH emission describes better the PDR gas heating as compared to the total infrared emission., Comment: Accepted for publication in Astronomy and Astrophysics. Fixed inverted line ratio in Sect. 5.2

We report the detection of oxygen-rich circumstellar envelopes in stars of the nearby (700 kpc) starburst galaxy IC 10. The star formation history and the chemical environment of this galaxy makes it an ideal target to observe dust production by high-mass stars in a low-metallicity environment. The goal of this study is to identify oxygen-rich stars in IC 10 and to constrain their nature between asymptotic giant branch stars (AGBs), red supergiants (RSGs), and other infrared bright sources. We examine the mass-loss rate of the stars and compare to results obtained for the Magellanic Clouds. Our objectives are to (1) assess whether RSGs can be significant dust producers in IC 10, and (2), solve the discrepancy between the star formation history of IC 10 and the relatively low number of RSGs detected in the optical. We search for silicate dust in emission by using the spectral map observed with the Infrared Spectrograph on board the Spitzer Space Telescope. The optical (UBVRI) and infrared (JHK, Spitzer/IRAC and Spitzer/MIPS) photometry is used to assert the membership of the stars to IC 10 and disentangle between AGBs and RSGs. Radiative models are used to infer mass-loss rates and stellar luminosities. The luminosity and colors of at least 9 silicate emission sources are consistent with stars within IC 10. Furthermore, the photometry of 2 of these sources is consistent with RSGs. We derive dust mass-loss rates similar to the values found in the Magellanic Clouds. Accounting for the sample completeness, RSGs are not important contributors to the dust mass budget in IC 10., Comment: Accepted for publication in A&A

(abridged) Our goal is to describe the multi-phase ISM of the IR bright low-metallicity galaxy Haro 11, dissecting the photoionised and photodissociated gas components. We present observations of the mid- and far-IR fine-structure cooling lines obtained with the Spitzer/IRS and Herschel/PACS spectrometers. We use the spectral synthesis code Cloudy to methodically model the ionised and neutral gas from which these lines originate. We find that the mid- and far-IR lines account for ~1% of the total IR luminosity L_TIR. Haro 11 is undergoing a phase of intense star formation, as traced by the brightest line [OIII] 88um, with L_[OIII]/L_TIR ~0.3%, and high ratios of [NeIII]/[NeII] and [SIV]/[SIII]. Due to their different origins, the observed lines require a multi-phase modeling comprising: a compact HII region, dense fragmented photodissociation regions (PDRs), a diffuse extended low-ionisation/neutral gas which has a volume filling factor of at least 90%, and porous warm dust in proximity to the stellar source. For a more realistic picture of the ISM of Haro 11 we would need to model the clumpy source and gas structures. We combine these 4 model components to explain the emission of 17 spectral lines, investigate the global energy balance of the galaxy through its spectral energy distribution, and establish a phase mass inventory. While the ionic emission lines of Haro 11 essentially originate from the dense HII region component, a diffuse low-ionisation gas is needed to explain the [NeII], [NII], and [CII] line intensities. The [OIII] 88um line intensity is not fully reproduced by our model, hinting towards the possible presence of yet another low-density high-ionisation medium. The [OI] emission is consistent with a dense PDR of low covering factor, and we find no evidence for an X-ray dominated component. The PDR component accounts for only 10% of the [CII] emission., Comment: Accepted for publication in A&A, 23 pages, 13 figures