We present the first spatially resolved study of molecular gas in the vicinity of a Gamma Ray Burst, using CO(2-1) emission line observations with the Atacama Large Millimetre Array (ALMA) at ~50 pc scales. The host galaxy of GRB 980425 contains a ring of high column density HI gas which is likely to have formed due to a collision between the GRB host and its companion galaxy, within which the GRB is located. We detect eleven molecular gas clumps in the galaxy, seven of which are within the gas ring. The clump closest to the GRB position is at a projected separation of ~280 pc. Although it is plausible that the GRB progenitor was ejected from clusters formed in this clump, we argue that the in situ formation of the GRB progenitor is the most likely scenario. We measure the molecular gas masses of the clumps and find them to be sufficient for forming massive star clusters. The molecular gas depletion times of the clumps show a variation of ~2 dex, comparable with the large variation in depletion times found in starburst galaxies in the nearby Universe. This demonstrates the presence of starburst modes of star formation on local scales in the galaxy, even while the galaxy as a whole cannot be categorised as a starburst based on its global properties. Our findings suggest that the progenitor of GRB 9802425 was originated in a young massive star cluster formed in the starburst mode of star formation.

Determining the efficiency with which gas is converted into stars in galaxies requires an accurate determination of the total reservoir of molecular gas mass. However, despite being the most abundant molecule in the Universe, H$_2$ is challenging to detect through direct observations and indirect methods have to be used to estimate the total molecular gas reservoir. These are often based on scaling relations from tracers such as CO or dust, and are generally calibrated in the Milky Way. Yet, evidence that these scaling relations are environmentally dependent is growing. In particular, the commonly used CO-to-H$_2$ conversion factor (X$_{\rm CO}$) is expected to be higher in metal-poor and/or strongly UV-irradiated environments. We use new SOFIA/FIFI-LS observations of far-infrared fine structure lines from the ionised and neutral gas and the Meudon photodissociation region model to constrain the physical properties and the structure of the gas in the massive star-forming region of 30 Doradus in the Large Magellanic Cloud, and determine the spatially resolved distribution of the total reservoir of molecular gas in the proximity of the young massive cluster R136. We compare this value with the molecular gas mass inferred from ground-based CO observations and dust-based estimates to quantify the impact of this extreme environment on commonly used tracers of the molecular gas. We find that the strong radiation field combined with the half-solar metallicity of the surrounding gas are responsible for a large reservoir of "CO-dark" molecular gas, leaving a large fraction of the total H$_2$ gas (> 75%) undetected when adopting a standard X$_{\rm CO}$ factor in this massive star-forming region.

We present the highest resolution (1") $^{12}$CO(1-0) observations of molecular gas in the dwarf starburst galaxy NGC 625, to date, obtained with ALMA. Molecular gas is distributed in discrete clouds within an area of $0.4$ kpc$^2$ and does not have well-ordered large-scale motions. We measure a total molecular mass in NGC 625 of $5.3\times 10^6 M_\odot$, assuming a Milky Way CO-to-H$_2$ conversion factor. We use the CPROPS package to identify molecular clouds and measure their properties. The 19 resolved CO clouds have a median radius of 20 pc, a median line width 2.5 km s$^{-1}$, and a median surface density of 169 $M_\odot$pc$^{-2}$. Larson scaling relations suggest that molecular clouds in NGC 625 are mostly in virial equilibrium. Comparison of our high-resolution CO observations with a star formation rate map, inferred from ancillary optical observations observations, suggests that about 40% of the molecular clouds coincide with the brightest HII regions. These bright HII regions have a range of molecular gas depletion timescales, all within a factor of $\sim 3$ of the global depletion time in NGC 625 of 106-134 Myr. The highest surface density molecular clouds towards the southwest of the galaxy, in a region we call the Butterfly, do not show strong star formation activity and suggest a depletion time scale longer than 5 Gyr.

We present new observations, carried out with IRAM NOEMA, of the atomic neutral carbon transitions [CI](1-0) at 492 GHz and [CI](2-1) at 809 GHz of GN20, a well-studied star-bursting galaxy at $z=4.05$. The high luminosity line ratio [CI](2-1)/[CI](1-0) implies an excitation temperature of $48^{+14}_{-9}$ K, which is significantly higher than the apparent dust temperature of $T_{\rm d}=33\pm2$ K ($\beta=1.9$) derived under the common assumption of an optically thin far-infrared dust emission, but fully consistent with $T_{\rm d}=52\pm5$ K of a general opacity model where the optical depth ($\tau$) reaches unity at a wavelength of $\lambda_0=170\pm23$ $\mu$m. Moreover, the general opacity solution returns a factor of $\sim 2\times$ lower dust mass and, hence, a lower molecular gas mass for a fixed gas-to-dust ratio, than with the optically thin dust model. The derived properties of GN20 thus provide an appealing solution to the puzzling discovery of starbursts appearing colder than main-sequence galaxies above $z>2.5$, in addition to a lower dust-to-stellar mass ratio that approaches the physical value predicted for starburst galaxies.

(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.

Cloud-scale surveys of molecular gas reveal the link between molecular clouds properties and star formation (SF) across a range of galactic environments. Cloud populations in galaxy disks are considered to be representative of the `normal' SF. At high resolution, however, clouds with exceptional gas properties and SF activity may also be observed in normal disk environments. In this paper, we study the brightest cloud traced in CO emission in the disk of NGC628. The cloud is spatially coincident with an extremely bright HII region. We characterize its molecular gas properties and investigate how feedback and large-scale processes influence the properties of the molecular gas. High resolution CO ALMA observations are used to characterize its mass and dynamical state, which are compared to other clouds in NGC628. A LVG analysis is used to constrain the beam-diluted density and temperature of the molecular gas. We analyze the MUSE spectrum using Starburst99 to characterize the young stellar population associated with the HII region. The cloud is massive ($1-2\times10^7$M$_{\odot}$), with a beam-diluted density of $n_{\rm H_2}=5\times10^4$ cm$^{-3}$. It has a low virial parameter, suggesting that its CO emission may be overluminous due to heating by the HII region. A young ($2-4$ Myr), massive $3\times10^{5}$ M$_{\odot}$ stellar population is associated. We argue that the cloud is currently being destroyed by feedback from young massive stars. Due to the cloud's large mass, this phase of the cloud's evolution is long enough for the impact of feedback on the excitation of the gas to be observed. Its high mass may be related to its location at a spiral co-rotation radius, where gas experiences reduced galactic shear compared to other regions of the disk, and receives a sustained inflow of gas that can promote the cloud's mass growth.

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.

We present EMPIRE, an IRAM 30-m large program that mapped $\lambda = 3{-}4$ mm dense gas tracers at $\sim 1{-}2\,$kpc resolution across the whole star-forming disk of nine nearby, massive, spiral galaxies. We describe the EMPIRE observing and reduction strategies and show new whole-galaxy maps of HCN(1-0), HCO$^+$(1-0), HNC(1-0) and CO(1-0). We explore how the HCN-to-CO and IR-to-HCN ratios, observational proxies for the dense gas fraction and dense gas star formation efficiency, depend on host galaxy and local environment. We find that the fraction of dense gas correlates with stellar surface density, gas surface density, molecular-to-atomic gas ratio, and dynamical equilibrium pressure. In EMPIRE, the star formation rate per unit dense gas anti-correlates with these same environmental parameters. Thus, although dense gas appears abundant the central regions of many spiral galaxies, this gas appears relatively inefficient at forming stars. These results qualitatively agree with previous work on nearby galaxies and the Milky Way's Central Molecular Zone. To first order, EMPIRE demonstrates that the conditions in a galaxy disk set the gas density distribution and that the dense gas traced by HCN shows an environment-dependent relation to star formation. However, our results also show significant ($\pm 0.2$ dex) galaxy-to-galaxy variations. We suggest that gas structure below the scale of our observations and dynamical effects likely also play an important role.

We report the Atacama Large Millimeter/submillimeter Array (ALMA) observations of the metal rich host galaxy of superluminous supernova (SLSN) PTF10tpz, a barred spiral galaxy at z=0.03994. We find the CO(1-0) emission to be confined within the bar of the galaxy. The distribution and kinematics of molecular gas in the host galaxy resemble gas flows along two lanes running from the tips of the bar towards the galaxy center. These gas lanes end in a gaseous structure in the inner region of the galaxy, likely associated with an inner Lindblad resonance. The interaction between the large-scale gas flows in the bar and the gas in the inner region plausibly leads to the formation of massive molecular clouds and consequently massive clusters. This in turn can result in formation of massive stars, and thus the likely progenitor of the SLSN in a young, massive cluster. This picture is consistent with SLSN PTF10tpz being located near the inner structure. We find the molecular gas in the vicinity of the SLSN to have high surface densities, comparable with those in interacting galaxies or starburst regions in nearby galaxies. This lends support to high densities being favorable conditions for formation of SLSNe progenitors.

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.

We use ALMA and IRAM 30-m telescope data to investigate the relationship between the spectroscopically-traced dense gas fraction and the cloud-scale (120 pc) molecular gas surface density in five nearby, star-forming galaxies. We estimate the dense gas mass fraction at 650 pc and 2800 pc scales using the ratio of HCN (1-0) to CO (1-0) emission. We then use high resolution (120 pc) CO (2-1) maps to calculate the mass-weighted average molecular gas surface density within 650 pc or 2770 pc beam where the dense gas fraction is estimated. On average, the dense gas fraction correlates with the mass-weighted average molecular gas surface density. Thus, parts of a galaxy with higher mean cloud-scale gas surface density also appear to have a larger fraction of dense gas. The normalization and slope of the correlation do vary from galaxy to galaxy and with the size of the regions studied. This correlation is consistent with a scenario where the large-scale environment sets the gas volume density distribution, and this distribution manifests in both the cloud-scale surface density and the dense gas mass fraction.

We present the first results of an ALMA survey of the lower fine structure line of atomic carbon [C I]$(^3P_1\,-\,^{3}P_0)$ in far infrared-selected galaxies on the main sequence at $z\sim1.2$ in the COSMOS field. We compare our sample with a comprehensive compilation of data available in the literature for local and high-redshift starbursting systems and quasars. We show that the [C I]($^3P_1$$\rightarrow$$^3P_0$) luminosity correlates on global scales with the infrared luminosity $L_{\rm IR}$ similarly to low-$J$ CO transitions. We report a systematic variation of $L'_{\rm [C\,I]^3P_1\,-\, ^3P_0}$/$L_{\rm IR}$ as a function of the galaxy type, with the ratio being larger for main-sequence galaxies than for starbursts and sub-millimeter galaxies at fixed $L_{\rm IR}$. The $L'_{\rm [C\,I]^3P_1\,-\, ^3P_0}$/$L'_{\rm CO(2-1)}$ and $M_{\rm{[C I]}}$/$M_{\rm dust}$ mass ratios are similar for main-sequence galaxies and for local and high-redshift starbursts within a 0.2 dex intrinsic scatter, suggesting that [C I] is a good tracer of molecular gas mass as CO and dust. We derive a fraction of $f_{\rm{[C\,I]}} = M_{\rm{[C\,I]}} / M_{\rm{C}}\sim3-13$% of the total carbon mass in the atomic neutral phase. Moreover, we estimate the neutral atomic carbon abundance, the fundamental ingredient to calibrate [C I] as a gas tracer, by comparing $L'_{\rm [C\,I]^3P_1\,-\, ^3P_0}$ and available gas masses from CO lines and dust emission. We find lower [C I] abundances in main-sequence galaxies than in starbursting systems and sub-millimeter galaxies, as a consequence of the canonical $\alpha_{\rm CO}$ and gas-to-dust conversion factors. This argues against the application to different galaxy populations of a universal standard [C I] abundance derived from highly biased samples.

This review describes where we are today in light of the dust and gas properties and their relation to star formation, in low metallicity galaxies of the local universe following recent surveys from sensitive infrared space telescopes, mainly Spitzer and Herschel space observatories as well as ground-based observations of the molecular gas reservoir. Models to interpret the ISM properties are gaining sophistication in order to account for the wide range of valuable observational diagnostics that we have today to trace the different gas phases, the broad range of photometry we have, from mid-infrared to submillimetre dust emission and the various galactic size scales that we can sample today. This review summarizes the rich multi-phase observations we can exploit today, and the multi-phase modeling approach to interpret the observations.

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.

We present ALMA Band 9 observations of the [C II]158um emission for a sample of 10 main-sequence galaxies at redshift z ~ 2, with typical stellar masses (log M*/Msun ~ 10.0 - 10.9) and star formation rates (~ 35 - 115 Msun/yr). Given the strong and well understood evolution of the interstellar medium from the present to z = 2, we investigate the behaviour of the [C II] emission and empirically identify its primary driver. We detect [C II] from six galaxies (four secure, two tentative) and estimate ensemble averages including non detections. The [C II]-to-infrared luminosity ratio (L[C II]/LIR) of our sample is similar to that of local main-sequence galaxies (~ 2 x 10^-3), and ~ 10 times higher than that of starbursts. The [C II] emission has an average spatial extent of 4 - 7 kpc, consistent with the optical size. Complementing our sample with literature data, we find that the [C II] luminosity correlates with galaxies' molecular gas mass, with a mean absolute deviation of 0.2 dex and without evident systematics: the [C II]-to-H2 conversion factor (alpha_[C II] ~ 30 Msun/Lsun) is largely independent of galaxies' depletion time, metallicity, and redshift. [C II] seems therefore a convenient tracer to estimate galaxies' molecular gas content regardless of their starburst or main-sequence nature, and extending to metal-poor galaxies at low and high redshifts. The dearth of [C II] emission reported for z > 6 - 7 galaxies might suggest either a high star formation efficiency or a small fraction of UV light from star formation reprocessed by dust.

We use new ALMA observations to investigate the connection between dense gas fraction, star formation rate, and local environment across the inner region of four local galaxies showing a wide range of molecular gas depletion times. We map HCN (1-0), HCO$^+$ (1-0), CS (2-1), $^{13}$CO (1-0), and C$^{18}$O (1-0) across the inner few kpc of each target. We combine these data with short spacing information from the IRAM large program EMPIRE, archival CO maps, tracers of stellar structure and recent star formation, and recent HCN surveys by Bigiel et al. and Usero et al. We test the degree to which changes in the dense gas fraction drive changes in the SFR. $I_{HCN}/I_{CO}$ (tracing the dense gas fraction) correlates strongly with $I_{CO}$ (tracing molecular gas surface density), stellar surface density, and dynamical equilibrium pressure, $P_{DE}$. Therefore, $I_{HCN}/I_{CO}$ becomes very low and HCN becomes very faint at large galactocentric radii, where ratios as low as $I_{HCN}/I_{CO} \sim 0.01$ become common. The apparent ability of dense gas to form stars, $\Sigma_{SFR}/\Sigma_{dense}$ (where $\Sigma_{dense}$ is traced by the HCN intensity and the star formation rate is traced by a combination of H$\alpha$ and 24$\mu$m emission), also depends on environment. $\Sigma_{SFR}/\Sigma_{dense}$ decreases in regions of high gas surface density, high stellar surface density, and high $P_{DE}$. Statistically, these correlations between environment and both $\Sigma_{SFR}/\Sigma_{dense}$ and $I_{HCN}/I_{CO}$ are stronger than that between apparent dense gas fraction ($I_{HCN}/I_{CO}$) and the apparent molecular gas star formation efficiency $\Sigma_{SFR}/\Sigma_{mol}$. We show that these results are not specific to HCN.

Carbon monoxide (CO) provides crucial information about the molecular gas properties of galaxies. While $^{12}$CO has been targeted extensively, isotopologues such as $^{13}$CO have the advantage of being less optically thick and observations have recently become accessible across full galaxy discs. We present a comprehensive new dataset of $^{13}$CO(1-0) observations with the IRAM 30-m telescope of the full discs of 9 nearby spiral galaxies from the EMPIRE survey at a spatial resolution of $\sim$1.5kpc. $^{13}$CO(1-0) is mapped out to $0.7-1r_{25}$ and detected at high signal-to-noise throughout our maps. We analyse the $^{12}$CO(1-0)-to-$^{13}$CO(1-0) ratio ($\Re$) as a function of galactocentric radius and other parameters such as the $^{12}$CO(2-1)-to-$^{12}$CO(1-0) intensity ratio, the 70-to-160$\mu$m flux density ratio, the star-formation rate surface density, the star-formation efficiency, and the CO-to-H$_2$ conversion factor. We find that $\Re$ varies by a factor of 2 at most within and amongst galaxies, with a median value of 11 and larger variations in the galaxy centres than in the discs. We argue that optical depth effects, most likely due to changes in the mixture of diffuse/dense gas, are favored explanations for the observed $\Re$ variations, while abundance changes may also be at play. We calculate a spatially-resolved $^{13}$CO(1-0)-to-H$_2$ conversion factor and find an average value of $1.0\times10^{21}$ cm$^{-2}$ (K.km/s)$^{-1}$ over our sample with a standard deviation of a factor of 2. We find that $^{13}$CO(1-0) does not appear to be a good predictor of the bulk molecular gas mass in normal galaxy discs due to the presence of a large diffuse phase, but it may be a better tracer of the mass than $^{12}$CO(1-0) in the galaxy centres where the fraction of dense gas is larger.

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.

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.

(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.

We use the IRAM Large Program EMPIRE and new high-resolution ALMA data to measure 13CO(1-0)/C18O(1-0) intensity ratios across nine nearby spiral galaxies. These isotopologues of CO are typically optically thin across most of the area in galaxy disks, and this ratio allows us to gauge their relative abundance due to chemistry or stellar nucleosynthesis effects. Resolved 13CO/C18O gradients across normal galaxies have been rare due to the faintness of these lines. We find a mean 13CO/C18O ratio of 6.0$\pm$0.9 for the central regions of our galaxies. This agrees well with results in the Milky Way, but differs from results for starburst galaxies (3.4$\pm$0.9) and ultraluminous infrared galaxies (1.1$\pm$0.4). In our sample, the 13CO/C18O ratio consistently increases with increasing galactocentric radius and decreases with increasing star formation rate surface density. These trends qualitatively agree with expectations for carbon and oxygen isotopic abundance variations due to stellar nucleosynthesis, with a possible effect of fractionation.

High critical density molecular lines like HCN(1-0) or HCO+(1-0) represent our best tool to study currently star-forming, dense molecular gas at extragalactic distances. The optical depth of these lines is a key ingredient to estimate the effective density required to excite emission. However, constraints on this quantity are even scarcer in the literature than measurements of the high density tracers themselves. Here, we combine new observations of HCN, HCO+ and HNC(1-0) and their optically thin isotopologues H13CN, H13CO+ and HN13C(1-0) to measure isotopologue line ratios. We use IRAM 30-m observations from the large program EMPIRE and new ALMA observations, which together target 6 nearby star-forming galaxies. Using spectral stacking techniques, we calculate or place strong upper limits on the HCN/H13CN, HCO+/H13CO+ and HNC/HN13C line ratios in the inner parts of these galaxies. Under simple assumptions, we use these to estimate the optical depths of HCN(1-0) and HCO+(1-0) to be \tau ~2-11 in the active, inner regions of our targets. The critical densities are consequently lowered to values between 5-20$\times 10^5$, 1-3$\times 10^5$ and 9$\times 10^4$ cm-3 for HCN, HCO+ and HNC, respectively. We study the impact of having different beam-filling factors, $\eta$, on these estimates and find that the effective critical densities decrease by a factor of $\frac{\eta_{12}}{\eta_{13}}\,\tau_{12}$. A comparison to existing work in NGC 5194 and NGC 253 shows HCN/H13CN and HCO+/H13CO+ ratios in agreement with our measurements within the uncertainties. The same is true for studies in other environments such as the Galactic Centre or nuclear regions of AGN-dominated nearby galaxies.

We explore the use of mm-wave emission line ratios to trace molecular gas density when observations integrate over a wide range of volume densities within a single telescope beam. For observations targeting external galaxies, this case is unavoidable. Using a framework similar to that of Krumholz and Thompson (2007), we model emission for a set of common extragalactic lines from lognormal and power law density distributions. We consider the median density of gas producing emission and the ability to predict density variations from observed line ratios. We emphasize line ratio variations, because these do not require knowing the absolute abundance of our tracers. Patterns of line ratio variations have the prospect to illuminate the high-end shape of the density distribution, and to capture changes in the dense gas fraction and median volume density. Our results with and without a high density power law tail differ appreciably; we highlight better knowledge of the PDF shape as an important area. We also show the implications of sub-beam density distributions for isotopologue studies targeting dense gas tracers. Differential excitation often implies a significant correction to the naive case. We provide tabulated versions of many of our results, which can be used to interpret changes in mm-wave line ratios in terms of changes in the underlying density distributions.

We make an inventory of the interstellar medium material in three low-metallicity dwarf spheroidal galaxies of the Local Group (NGC147, NGC185 and NGC205). Ancillary HI, CO, Spitzer IRS spectra, H{\alpha} and X-ray observations are combined to trace the atomic, cold and warm molecular, ionised and hot gas phases. We present new Nobeyama CO(1-0) observations and Herschel SPIRE FTS [CI] observations of NGC205 to revise its molecular gas content. We derive total gas masses of M_gas = 1.9-5.5x10^5 Msun for NGC185 and M_gas = 8.6-25.0x10^5 Msun for NGC205. Non-detections combine to an upper limit on the gas mass of M_gas =< 0.3-2.2x10^5 Msun for NGC147. The observed gas reservoirs are significantly lower compared to the expected gas masses based on a simple closed-box model that accounts for the gas mass returned by planetary nebulae and supernovae. The gas-to-dust mass ratios GDR~37-107 and GDR~48-139 are also considerably lower compared to the expected GDR~370 and GDR~520 for the low metal abundances in NGC 185 (0.36 Zsun) and NGC205 (0.25 Zsun), respectively. To simultaneously account for the gas deficiency and low gas-to-dust ratios, we require an efficient removal of a large gas fraction and a longer dust survival time (~1.6 Gyr). We believe that efficient galactic winds (combined with heating of gas to sufficiently high temperatures in order for it to escape from the galaxy) and/or environmental interactions with neighbouring galaxies are responsible for the gas removal from NGC147, NGC185 and NGC205.

The [CII] 158 um fine structure line is one of the dominant cooling lines in the interstellar medium (ISM) and is an important tracer of star formation. Recent velocity-resolved studies with Herschel/HIFI and SOFIA/GREAT showed that the [CII] line can constrain the properties of the ISM phases in star-forming regions. The [CII] line as a tracer of star formation is particularly important in low-metallicity environments where CO emission is weak because of the presence of large amounts of CO-dark gas. The nearby irregular dwarf galaxy NGC 4214 offers an excellent opportunity to study an actively star-forming ISM at low metallicity. We analyzed the spectrally resolved [CII] line profiles in three distinct regions at different evolutionary stages of NGC 4214 with respect to ancillary HI and CO data in order to study the origin of the [CII] line. We used SOFIA/GREAT [CII] 158 um observations, HI data from THINGS, and CO(2-1) data from HERACLES to decompose the spectrally resolved [CII] line profiles into components associated with neutral atomic and molecular gas. We use this decomposition to infer gas masses traced by [CII] under different ISM conditions. Averaged over all regions, we associate about 46% of the [CII] emission with the HI emission. However, we can assign only around 9% of the total [CII] emission to the cold neutral medium (CNM). We found that about 79% of the total molecular hydrogen mass is not traced by CO emission. On average, the fraction of CO-dark gas dominates the molecular gas mass budget. The fraction seems to depend on the evolutionary stage of the regions: it is highest in the region covering a super star cluster in NGC 4214, while it is lower in a more compact, more metal-rich region.

We present the Herschel SPIRE Fourier Transform Spectroscopy (FTS) atlas for a complete flux limited sample of local Ultra-Luminous Infra-Red Galaxies as part of the HERschel ULIRG Survey (HERUS). The data reduction is described in detail and was optimized for faint FTS sources with particular care being taken with the subtraction of the background which dominates the continuum shape of the spectra. Special treatment in the data reduction has been given to any observation suffering from artefacts in the data caused by anomalous instrumental effects to improve the final spectra. Complete spectra are shown covering $200 - 671\mu$m with photometry in the SPIRE bands at 250$\mu$m, 350$\mu$m and 500$\mu$m. The spectra include near complete CO ladders for over half of our sample, as well as fine structure lines from [CI] 370 $\mu$m, [CI] 609 $\mu$m, and [NII] 205 $\mu$m. We also detect H$_{2}$O lines in several objects. We construct CO Spectral Line Energy Distributions (SLEDs) for the sample, and compare their slopes with the far-infrared colours and luminosities. We show that the CO SLEDs of ULIRGs can be broadly grouped into three classes based on their excitation. We find that the mid-J (5$<$J$<$8) lines are better correlated with the total far-infrared luminosity, suggesting that the warm gas component is closely linked to recent star-formation. The higher J transitions do not linearly correlate with the far-infrared luminosity, consistent with them originating in hotter, denser gas unconnected to the current star-formation. {\bf We conclude that in most cases more than one temperature components are required to model the CO SLEDs.}

One of the key goals of the Bluedisk survey is to characterize the impact of gas accretion in disc galaxies in the context of galaxy evolution. It contains 50 disc galaxies in the stellar mass range 10^10-10^11 Msun, of which half are bluer and more HI-rich galaxies than their HI-normal (control) counterparts. In this paper, we investigate how ongoing disc growth affects the molecular gas distribution and the star-formation efficiency in these galaxies. We present 12CO observations from the IRAM 30-m telescope in 26 galaxies of the Bluedisk survey. We compare the amount and spatial distribution of the molecular gas to key quantities such as atomic gas, stellar mass and surface density, star-formation rate and metallicity. We analyse the star-formation rate per unit gas (SFR/HI and SFR/H2) and relate all those parameters to general galaxy properties (HI-rich/control disc, morphology, etc.). We find that the HI-rich galaxies have similar H2 masses as the control galaxies. In their centres, HI-rich galaxies have lower H2/HI ratios and marginally shorter molecular gas depletion times. However, the main differences between the two samples occur in the outer parts of the discs, with the HI-rich galaxies having slightly smaller CO discs (relative to the optical radius R25) and steeper CO and metallicity gradients than the control galaxies. The ongoing accretion of HI at large radii has thus not led to an appreciable growth of the CO discs in our sample. Based on depletion times, we estimate that this gas will contribute to star formation on time-scales of at least 5 Gyr.

We present Herschel SPIRE Fourier Transform Spectrometer (FTS) observations of N159W, an active star-forming region in the Large Magellanic Cloud (LMC). In our observations, a number of far-infrared cooling lines including CO(4-3) to CO(12-11), [CI] 609 and 370 micron, and [NII] 205 micron are clearly detected. With an aim of investigating the physical conditions and excitation processes of molecular gas, we first construct CO spectral line energy distributions (SLEDs) on 10 pc scales by combining the FTS CO transitions with ground-based low-J CO data and analyze the observed CO SLEDs using non-LTE radiative transfer models. We find that the CO-traced molecular gas in N159W is warm (kinetic temperature of 153-754 K) and moderately dense (H2 number density of (1.1-4.5)e3 cm-3). To assess the impact of the energetic processes in the interstellar medium on the physical conditions of the CO-emitting gas, we then compare the observed CO line intensities with the models of photodissociation regions (PDRs) and shocks. We first constrain the properties of PDRs by modelling Herschel observations of [OI] 145, [CII] 158, and [CI] 370 micron fine-structure lines and find that the constrained PDR components emit very weak CO emission. X-rays and cosmic-rays are also found to provide a negligible contribution to the CO emission, essentially ruling out ionizing sources (ultraviolet photons, X-rays, and cosmic-rays) as the dominant heating source for CO in N159W. On the other hand, mechanical heating by low-velocity C-type shocks with ~10 km/s appears sufficient enough to reproduce the observed warm CO.

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.

Dwarf spheroidal galaxies are among the most numerous galaxy population in the Universe, but their main formation and evolution channels are still not well understood. The three dwarf spheroidal satellites (NGC147, NGC185, and NGC205) of the Andromeda galaxy are characterised by very different interstellar medium (ISM) properties, which might suggest them being at different galaxy evolutionary stages. While the dust content of NGC205 has been studied in detail by De Looze et al. (2012), we present new Herschel dust continuum observations of NGC147 and NGC185. The non-detection of NGC147 in Herschel SPIRE maps puts a strong constraint on its dust mass (< 128 Msun). For NGC185, we derive a total dust mass M_d = 5.1 x 10^3 Msun, which is a factor of ~2-3 higher than that derived from ISO and Spitzer observations and confirms the need for longer wavelength observations to trace more massive cold dust reservoirs. We, furthermore, estimate the dust production by asymptotic giant branch (AGB) stars and supernovae (SNe). For NGC147, the upper limit on the dust mass is consistent with expectations of the material injected by the evolved stellar population. In NGC185 and NGC205, the observed dust content is one order of magnitude higher compared to the estimated dust production by AGBs and SNe. Efficient grain growth, and potentially longer dust survival times (3-6 Gyr) are required to account for their current dust content. Our study confirms the importance of grain growth in the gas phase to account for the current dust reservoir in galaxies.

We present a method to simultaneously model the dust far-infrared spectral energy distribution (SED) and the total infrared $-$ carbon monoxide (CO) integrated intensity $(S_{\rm IR}-I_{\rm CO})$ relationship. The modelling employs a hierarchical Bayesian (HB) technique to estimate the dust surface density, temperature ($T_{\rm eff}$), and spectral index at each pixel from the observed far-infrared (FIR) maps. Additionally, given the corresponding CO map, the method simultaneously estimates the slope and intercept between the FIR and CO intensities, which are global properties of the observed source. The model accounts for correlated and uncorrelated uncertainties, such as those present in Herschel observations. Using synthetic datasets, we demonstrate the accuracy of the HB method, and contrast the results with common non-hierarchical fitting methods. As an initial application, we model the dust and gas on 100 pc scales in the Magellanic Clouds from Herschel FIR and NANTEN CO observations. The slopes of the $\log S_{\rm IR}-\log I_{\rm CO}$ relationship are similar in both galaxies, falling in the range 1.1$-$1.7. However, in the SMC the intercept is nearly 3 times higher, which can be explained by its lower metallicity than the LMC, resulting in a larger $S_{\rm IR}$ per unit $I_{\rm CO}$. The HB modelling evidences an increase in $T_{\rm eff}$ in regions with the highest $I_{\rm CO}$ in the LMC. This may be due to enhanced dust heating in the densest molecular regions from young stars. Such simultaneous dust and gas modelling may reveal variations in the properties of the ISM and its association with other galactic characteristics, such as star formation rates and/or metallicities.

We present the first results from the EMPIRE survey, an IRAM large program that is mapping tracers of high density molecular gas across the disks of nine nearby star-forming galaxies. Here, we present new maps of the 3-mm transitions of HCN, HCO+, and HNC across the whole disk of our pilot target, M51. As expected, dense gas correlates with tracers of recent star formation, filling the "luminosity gap" between Galactic cores and whole galaxies. In detail, we show that both the fraction of gas that is dense, f_dense traced by HCN/CO, and the rate at which dense gas forms stars, SFE_dense traced by IR/HCN, depend on environment in the galaxy. The sense of the dependence is that high surface density, high molecular gas fraction regions of the galaxy show high dense gas fractions and low dense gas star formation efficiencies. This agrees with recent results for individual pointings by Usero et al. 2015 but using unbiased whole-galaxy maps. It also agrees qualitatively with the behavior observed contrasting our own Solar Neighborhood with the central regions of the Milky Way. The sense of the trends can be explained if the dense gas fraction tracks interstellar pressure but star formation occurs only in regions of high density contrast.

Comparison of the ISM properties of a wide range of metal-poor galaxies with normal metal-rich galaxies reveals striking differences. We find that the combination of the low dust abundance and the active star formation results in a very porous ISM filled with hard photons, heating the dust in dwarf galaxies to overall higher temperatures than their metal-rich counterparts. This results in photodissociation of molecular clouds to greater depths, leaving relatively large PDR envelopes and difficult-to-detect CO cores. From detailed modeling of the low-metallicity ISM, we find significant fractions of CO-dark H2 - a reservoir of molecular gas not traced by CO, but present in the [CII] and [CI]-emitting envelopes. Self-consistent analyses of the neutral and ionized gas diagnostics along with the dust SED is the necessary way forward in uncovering the multiphase structure of galaxies

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)

A correlation between the 24 {\mu}m continuum and the [NII] 205 {\mu}m line emission may arise if both quantities trace the star formation activity on spatially-resolved scales within a galaxy, yet has so far only been observed in the nearby edge-on spiral galaxy NGC 891. We therefore assess whether the [NII] 205 - 24 {\mu}m emission correlation has some physical origin or is merely an artefact of line-of-sight projection effects in an edge-on disc. We search for the presence of a correlation in Herschel and Spitzer observations of two nearby face-on galaxies, M51 and M83, and the interacting Antennae galaxies NGC 4038 and 4039. We show that not only is this empirical relationship also observed in face-on galaxies, but also that the correlation appears to be governed by the star formation rate (SFR). Both the nuclear starburst in M83 and the merger-induced star formation in NGC 4038/9 exhibit less [NII] emission per unit SFR surface density than the normal star-forming discs. These regions of intense star formation exhibit stronger ionization parameters, as traced by the 70/160 {\mu}m far-infrared colour, that suggest the presence of higher ionization lines that may become more important for gas cooling, thereby reducing the observed [NII] 205 {\mu}m line emission in regions with higher star formation rates. Finally, we present a general relation between the [NII] 205 {\mu}m line flux density and SFR density for normal star-forming galaxies, yet note that future studies should extend this analysis by including observations with wider spatial coverage for a larger sample of galaxies.

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.

We present far-infrared spectral line observations of five galaxies from the LITTLE THINGS sample: DDO 69, DDO 70, DDO 75, DDO 155, and WLM. While most studies of dwarfs focus on bright systems or starbursts due to observational constraints, our data extend the observed parameter space into the regime of low surface brightness dwarf galaxies with low metallicities and moderate star formation rates. Our targets were observed with Herschel at the [CII] 158um, [OI] 63um, [OIII] 88um, and NII 122um emission lines using the PACS Spectrometer. These high-resolution maps allow us for the first time to study the far-infrared properties of these systems on the scales of larger star-forming complexes. The spatial resolution in our maps, in combination with star formation tracers, allows us to identify separate PDRs in some of the regions we observed. Our systems have widespread [CII] emission that is bright relative to continuum, averaging near 0.5% of the total infrared budget - higher than in solar-metallicity galaxies of other types. [NII] is weak, suggesting that the [CII] emission in our galaxies comes mostly from PDRs instead of the diffuse ionized ISM. These systems exhibit efficient cooling at low dust temperatures, as shown by ([OI]+[CII])/TIR in relation to 60um/100um, and low [OI]/[CII] ratios which indicate that [CII] is the dominant coolant of the ISM. We observe [OIII]/[CII] ratios in our galaxies that are lower than those published for other dwarfs, but similar to levels noted in spirals.

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)

We present a model for the interstellar medium of NGC4214 with the objective to probe the physical conditions in the two main star-forming regions and their connection with the star formation activity of the galaxy. We used the spectral synthesis code Cloudy to model an HII region and the associated photodissociation region (PDR) to reproduce the emission of mid- and far-infrared fine-structure cooling lines from the Spitzer and Herschel space telescopes for these two regions. Input parameters of the model, such as elemental abundances and star formation history, are guided by earlier studies of the galaxy, and we investigated the effect of the mode in which star formation takes place (bursty or continuous) on the line emission. Furthermore, we tested the effect of adding pressure support with magnetic fields and turbulence on the line predictions. We find that this model can satisfactorily predict (within a factor of ~2) all observed lines that originate from the ionized medium ([SIV] 10.5um, [NeIII] 15.6um, [SIII] 18.7um, [SIII] 33.5um, and [OIII] 88um), with the exception of [NeII] 12.8um and [NII] 122um, which may arise from a lower ionization medium. In the PDR, the [OI] 63um, [OI] 145um, and [CII] 157um lines are matched within a factor of ~5 and work better when weak pressure support is added to the thermal pressure or when the PDR clouds are placed farther away from the HII regions and have covering factors lower than unity. Our models of the HII region agree with different evolutionary stages found in previous studies, with a more evolved, diffuse central region, and a younger, more compact southern region. However, the local PDR conditions are averaged out on the 175 pc scales that we probe and do not reflect differences observed in the star formation properties of the two regions.

As part of the Bluedisk survey we analyse the radial gas-phase metallicity profiles of 50 late-type galaxies We compare the metallicity profiles of a sample of HI-rich galaxies against a control sample of HI-'normal' galaxies. We find the metallicity gradient of a galaxy to be strongly correlated with its HI mass fraction (M(HI) / Mstar). We note that some galaxies exhibit a steeper metallicity profile in the outer disc than in the inner disc. These galaxies are found in both the HI-rich and control samples. This contradicts a previous indication that these outer drops are exclusive to HI-rich galaxies. These effects are not driven by bars, although we do find some indication that barred galaxies have flatter metallicity profiles. By applying a simple analytical model we are able to account for the variety of metallicity profiles that the two samples present. The success of this model implies that the metallicity in these isolated galaxies may be in a local equilibrium, regulated by star formation. This insight could provide an explanation of the observed local mass-metallicity relation.

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.

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.

Since the launch of the Herschel Space Observatory, our understanding about the photo-dissociation regions (PDR) has taken a step forward. In the bandwidth of the Fourier Transform Spectrometer (FTS) of the Spectral and Photometric Imaging REceiver (SPIRE) on board Herschel, ten CO rotational transitions, including J=4-3 to J=13-12, and three fine structure lines, including [CI] 609, [CI] 370, and [NII] 250 micron, are covered. In this paper, we present our findings from the FTS observations at the nuclear region of M83, based on the spatially resolved physical parameters derived from the CO spectral line energy distribution (SLED) map and the comparisons with the dust properties and star-formation tracers. We discuss (1) the potential of using [NII] 250 and [CI] 370 micron as star-formation tracers; (2) the reliability of tracing molecular gas with CO; (3) the excitation mechanisms of warm CO; (4) the possibility of studying stellar feedback by tracing the thermal pressure of molecular gas in the nuclear region of M83.

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.

Context: Dust reprocesses about half of the stellar radiation in galaxies. The thermal re-emission by dust of absorbed energy is considered driven merely by young stars and, consequently, often applied to trace the star formation rate in galaxies. Recent studies have argued that the old stellar population might anticipate a non-negligible fraction of the radiative dust heating. Aims: In this work, we aim to analyze the contribution of young (< 100 Myr) and old (~ 10 Gyr) stellar populations to radiative dust heating processes in the nearby grand-design spiral galaxy M51 using radiative transfer modeling. High-resolution 3D radiative transfer (RT) models are required to describe the complex morphologies of asymmetric spiral arms and clumpy star-forming regions and model the propagation of light through a dusty medium. Methods: In this paper, we present a new technique developed to model the radiative transfer effects in nearby face-on galaxies. We construct a high-resolution 3D radiative transfer model with the Monte-Carlo code SKIRT accounting for the absorption, scattering and non-local thermal equilibrium (NLTE) emission of dust in M51. The 3D distribution of stars is derived from the 2D morphology observed in the IRAC 3.6 {\mu}m, GALEX FUV, H{\alpha} and MIPS 24 {\mu}m wavebands, assuming an exponential vertical distribution with an appropriate scale height. The dust geometry is constrained through the far-ultraviolet (FUV) attenuation, which is derived from the observed total-infrared-to-far-ultraviolet luminosity ratio. The stellar luminosity, star formation rate and dust mass have been scaled to reproduce the observed stellar spectral energy distribution (SED), FUV attenuation and infrared SED. (abridged)

We analyze the applicability of far-infrared fine-structure lines [CII] 158 micron, [OI] 63 micron and [OIII] 88 micron to reliably trace the star formation rate (SFR) in a sample of low-metallicity dwarf galaxies from the Herschel Dwarf Galaxy Survey and compare with a broad sample of galaxies of various types and metallicities in the literature. We study the trends and scatter in the relation between the SFR (as traced by GALEX FUV and MIPS 24 micron) and far-infrared line emission, on spatially resolved and global galaxy scales, in dwarf galaxies. We assemble far-infrared line measurements from the literature and infer whether the far-infrared lines can probe the SFR (as traced by the total-infrared luminosity) in a variety of galaxy populations. In metal-poor dwarfs, the [OI] and [OIII] lines show the strongest correlation with the SFR with an uncertainty on the SFR estimates better than a factor of 2, while the link between [CII] emission and the SFR is more dispersed (uncertainty factor of 2.6). The increased scatter in the SFR-L([CII]) relation towards low metal abundances, warm dust temperatures, large filling factors of diffuse, highly ionized gas suggests that other cooling lines start to dominate depending on the density and ionization state of the gas. For the literature sample, we evaluate the correlations for a number of different galaxy populations. The [CII] and [OI] lines are considered to be reliable SFR tracers in starburst galaxies, recovering the star formation activity within an uncertainty of factor 2. [Abridged]

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.

We present Herschel/PACS 100 and 160 micron integrated photometry for the 323 galaxies in the Herschel Reference Survey (HRS), a K-band-, volume-limited sample of galaxies in the local Universe. Once combined with the Herschel/SPIRE observations already available, these data make the HRS the largest representative sample of nearby galaxies with homogeneous coverage across the 100-500 micron wavelength range. In this paper, we take advantage of this unique dataset to investigate the properties and shape of the far-infrared/sub-millimeter spectral energy distribution in nearby galaxies. We show that, in the stellar mass range covered by the HRS (8<log(M*/Msun)<12), the far-infrared/sub-millimeter colours are inconsistent with a single modified black-body having the same dust emissivity index beta for all galaxies. In particular, either beta decreases, or multiple temperature components are needed, when moving from metal-rich/gas-poor to metal-poor/gas-rich galaxies. We thus investigate how the dust temperature and mass obtained from a single modified black-body depend on the assumptions made on beta. We show that, while the correlations between dust temperature, galaxy structure and star formation rate are strongly model dependent, the dust mass scaling relations are much more reliable, and variations of beta only change the strength of the observed trends.

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.

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)

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.

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.

We present Herschel observations of six fine-structure lines in 25 Ultraluminous Infrared Galaxies at z<0.27. The lines, [O III]52, [N III]57, [O I]63, [N II]122, [O I]145, and [C II]158, are mostly single gaussians with widths <600 km s-1 and luminosities of 10^7 - 10^9 Solar. There are deficits in the [O I]63/L_IR, [N II]/L_IR, [O I]145/L_IR, and [C II]/L_IR ratios compared to lower luminosity systems. The majority of the line deficits are consistent with dustier H II regions, but part of the [C II] deficit may arise from an additional mechanism, plausibly charged dust grains. This is consistent with some of the [C II] originating from PDRs or the ISM. We derive relations between far-IR line luminosities and both IR luminosity and star formation rate. We find that [N II] and both [O I] lines are good tracers of IR luminosity and star formation rate. In contrast, [C II] is a poor tracer of IR luminosity and star formation rate, and does not improve as a tracer of either quantity if the [C II] deficit is accounted for. The continuum luminosity densities also correlate with IR luminosity and star formation rate. We derive ranges for the gas density and ultraviolet radiation intensity of 10^1 < n < 10^2.5 and 10^2.2 < G_0 < 10^3.6, respectively. These ranges depend on optical type, the importance of star formation, and merger stage. We do not find relationships between far-IR line properties and several other parameters; AGN activity, merger stage, mid-IR excitation, and SMBH mass. We conclude that these far-IR lines arise from gas heated by starlight, and that they are not strongly influenced by AGN activity.

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.

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.

We present ALMA observations of 30 Doradus -- the highest resolution view of molecular gas in an extragalactic star formation region to date (~0.4pc x 0.6pc). The 30Dor-10 cloud north of R136 was mapped in 12CO 2-1, 13CO 2-1, C18O 2-1, 1.3mm continuum, the H30alpha recombination line, and two H2CO 3-2 transitions. Most 12CO emission is associated with small filaments and clumps (<1pc, ~1000 Msun at the current resolution). Some clumps are associated with protostars, including "pillars of creation" photoablated by intense radiation from R136. Emission from molecular clouds is often analyzed by decomposition into approximately beam-sized clumps. Such clumps in 30 Doradus follow similar trends in size, linewidth, and surface density to Milky Way clumps. The 30 Doradus clumps have somewhat larger linewidths for a given size than predicted by Larson's scaling relation, consistent with pressure confinement. They extend to higher surface density at a given size and linewidth compared to clouds studied at 10pc resolution. These trends are also true of clumps in Galactic infrared-dark clouds; higher resolution observations of both environments are required. Consistency of clump masses calculated from dust continuum, CO, and the virial theorem reveals that the CO abundance in 30 Doradus clumps is not significantly different from the LMC mean, but the dust abundance may be reduced by ~2. There are no strong trends in clump properties with distance from R136; dense clumps are not strongly affected by the external radiation field, but there is a modest trend towards lower dense clump filling fraction deeper in the cloud.

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.

[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.

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.

We present the sub-millimeter spectra from 450 GHz to 1550 GHz of eleven nearby active galaxies observed with the SPIRE Fourier Transform Spectrometer (SPIRE/FTS) onboard Herschel. We detect CO transitions from J_up = 4 to 12, as well as the two [CI] fine structure lines at 492 and 809 GHz and the [NII] 461 GHz line. We used radiative transfer models to analyze the observed CO spectral line energy distributions (SLEDs). The FTS CO data were complemented with ground-based observations of the low-J CO lines. We found that the warm molecular gas traced by the mid-J CO transitions has similar physical conditions (n_H2 ~ 10^3.2 - 10^3.9 cm^-3 and T_kin ~ 300 - 800 K) in most of our galaxies. Furthermore, we found that this warm gas is likely producing the mid-IR rotational H2 emission. We could not determine the specific heating mechanism of the warm gas, however it is possibly related to the star-formation activity in these galaxies. Our modeling of the [CI] emission suggests that it is produced in cold (T_kin < 30 K) and dense (n_H2 > 10^3 cm^-3) molecular gas. Transitions of other molecules are often detected in our SPIRE/FTS spectra. The HF J=1-0 transition at 1232 GHz is detected in absorption in UGC05101 and in emission in NGC7130. In the latter, near-infrared pumping, chemical pumping, or collisional excitation with electrons are plausible excitation mechanisms likely related to the AGN of this galaxy. In some galaxies few H2O emission lines are present. Additionally, three OH+ lines at 909, 971, and 1033 GHz are identified in NGC7130.

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.

We study for the first time the complete sub-millimeter spectra (450 GHz to 1550 GHz) of a sample of nearby active galaxies observed with the SPIRE Fourier Transform Spectrometer (SPIRE/FTS) onboard Herschel. The CO ladder (from Jup = 4 to 12) is the most prominent spectral feature in this range. These CO lines probe warm molecular gas that can be heated by ultraviolet photons, shocks, or X-rays originated in the active galactic nucleus or in young star-forming regions. In these proceedings we investigate the physical origin of the CO emission using the averaged CO spectral line energy distribution (SLED) of six Seyfert galaxies. We use a radiative transfer model assuming an isothermal homogeneous medium to estimate the molecular gas conditions. We also compare this CO SLED with the predictions of photon and X-ray dominated region (PDR and XDR) models.

(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.

(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.

The low metallicity interstellar medium of dwarf galaxies gives a different picture in the far infrared(FIR)/submillimetre(submm)wavelengths than the more metal-rich galaxies. Excess emission is often found in the submm beginning at or beyond 500 mu. Even without taking this excess emission into account as a possible dust component, higher dust-to-gas mass ratios (DGR) are often observed compared to that expected from their metallicity for moderately metal-poor galaxies. The SEDs of the lowest metallicity galaxies, however, give very low dust masses and excessively low values of DGR, inconsistent with the amount of metals expected to be captured into dust if we presume the usual linear relationship holding for all metallicities, including the more metal-rich galaxies. This transition seems to appear near metalllicities of 12 + log(O/H) ~ 8.0 - 8.2. These results rely on accurately quantifying the total molecular gas reservoir, which is uncertain in low metallicity galaxies due to the difficulty in detecting CO(1-0) emission. Dwarf galaxies show an exceptionally high [CII] 158 mu/CO (1-0) ratio which may be indicative of a significant reservoir of 'CO-free' molecular gas residing in the photodissociated envelope, and not traced by the small CO cores.

We present Herschel dust continuum, James Clerk Maxwell Telescope CO(3-2) observations and a search for [CII] 158 micron and [OI] 63 micron spectral line emission for the brightest early-type dwarf satellite of Andromeda, NGC 205. While direct gas measurements (Mgas ~ 1.5e+6 Msun, HI + CO(1-0)) have proven to be inconsistent with theoretical predictions of the current gas reservoir in NGC 205 (> 1e+7 Msun), we revise the missing interstellar medium mass problem based on new gas mass estimates (CO(3-2), [CII], [OI]) and indirect measurements of the interstellar medium content through dust continuum emission. Based on Herschel observations, covering a wide wavelength range from 70 to 500 micron, we are able to probe the entire dust content in NGC 205 (Mdust ~ 1.1-1.8e+4 Msun at Tdust ~ 18-22 K) and rule out the presence of a massive cold dust component (Mdust ~ 5e+5 Msun, Tdust ~ 12 K), which was suggested based on millimeter observations from the inner 18.4 arcsec. Assuming a reasonable gas-to-dust ratio of ~ 400, the dust mass in NGC 205 translates into a gas mass Mgas ~ 4-7e+6 Msun. The non-detection of [OI] and the low L_[CII]-to-L_CO(1-0) line intensity ratio (~ 1850) imply that the molecular gas phase is well traced by CO molecules in NGC 205. We estimate an atomic gas mass of 1.5e+4 Msun associated with the [CII] emitting PDR regions in NGC 205. From the partial CO(3-2) map of the northern region in NGC 205, we derive a molecular gas mass of M_H2 ~ 1.3e+5 Msun. [abridged]

Molecular globules and pillars are spectacular features, found only in the interface region between a molecular cloud and an HII-region. Impacting Far-ultraviolet (FUV) radiation creates photon dominated regions (PDRs) on their surfaces that can be traced by typical cooling lines. With the GREAT receiver onboard SOFIA we mapped and spectrally resolved the [CII] 158 micron atomic fine-structure line and the highly excited 12CO J=11-10 molecular line from three objects in Cygnus X (a pillar, a globule, and a strong IRAS source). We focus here on the globule and compare our data with existing Spitzer data and recent Herschel Open-Time PACS data. Extended [CII] emission and more compact CO-emission was found in the globule. We ascribe this emission mainly to an internal PDR, created by a possibly embedded star-cluster with at least one early B-star. However, external PDR emission caused by the excitation by the Cyg OB2 association cannot be fully excluded. The velocity-resolved [CII] emission traces the emission of PDR surfaces, possible rotation of the globule, and high-velocity outflowing gas. The globule shows a velocity shift of ~2 km/s with respect to the expanding HII-region, which can be understood as the residual turbulence of the molecular cloud from which the globule arose. This scenario is compatible with recent numerical simulations that emphazise the effect of turbulence. It is remarkable that an isolated globule shows these strong dynamical features traced by the [CII]-line, but it demands more observational studies to verify if there is indeed an embedded cluster of B-stars.

We present photometry of the nearby galaxy NGC 5128 (Centaurus A) observed with the PACS and SPIRE instruments on board the Herschel Space Observatory, at 70, 160, 250, 350 and 500 {\mu}m, as well as new CO J = 3-2 observations taken with the HARP-B instrument on the JCMT. Using a single component modified blackbody, we model the dust spectral energy distribution within the disk of the galaxy using all five Herschel wavebands, and find dust temperatures of ~30 K towards the centre of the disk and a smoothly decreasing trend to ~20 K with increasing radius. We find a total dust mass of (1.59 \pm 0.05) \times 10^7 M\odot, and a total gas mass of (2.7 \pm 0.2) \times 10^9 M\odot. The average gas-to-dust mass ratio is 103 \pm 8 but we find an interesting increase in this ratio to approximately 275 toward the centre of Cen A. We discuss several possible physical processes that may be causing this effect, including dust sputtering, jet entrainment and systematic variables such as the XCO factor. Dust sputtering by X-rays originating in the AGN or the removal of dust by the jets are our most favoured explanations.

The Herschel Dwarf Galaxy Survey investigates the interplay of star formation activity and the the metal-poor gas and dust of dwarf galaxies using FIR and submillimetre imaging spectroscopic and photometric observations in the 50 to 550mu window of the Herschel Space Observatory. The dust SEDs are well constrained with the new Herschel and MIR Spitzer data. A submillimetre excess is often found in low metallicity galaxies, which,if tracing very cold dust, would highlight large dust masses not easily reconciled in some cases, given the low metallicities and expected gas-to-dust mass ratios. The galaxies are also mapped in the FIR fine-structure lines (63 and 145mu OI, 158mu CII, 122 and 205mu NII, 88mu OIII) probing the low density ionised gas, the HII regions and photodissociation regions. While still early in the Herschel mission we can already see, along with earlier studies, that line ratios in the metal-poor ISM differ remarkably from those in the metal-rich starburst environments. In dwarf galaxies, L[CII]/L(CO) (>10^4) is at least an order of magnitude greater than in the most metal-rich starburst galaxies. The enhanced [CII] arises from the larger photodissociation region where H2, not traced by the CO, can exist. The 88mu [OIII] line usually dominates the FIR line emission over galaxy-wide scales in dwarf galaxies, not the 158mu [CII] line which is the dominant FIR cooling line in metal-rich galaxies. All of the FIR lines together can contribute 1% to 2% of the L(TIR). The Herschel Dwarf Galaxy survey will provide statistical information on the nature of the dust and gas in low metallicity galaxies, elucidating the origin of the submm excess in dwarf galaxies, and help determine a ([CII] +CO) to H2 conversion factor, thus providing observational constraints on chemical evolution models of galaxies.

Galactic cirrus emission at far-infrared wavelengths affects many extragalactic observations. Separating this emission from that associated with extragalactic objects is both important and difficult. In this paper we discuss a particular case, the M81 group, and the identification of diffuse structures prominent in the infrared, but also detected at optical wavelengths. The origin of these structures has previously been controversial, ranging from them being the result of a past interaction between M81 and M82 or due to more local Galactic emission. We show that over of order a few arcminute scales the far-infrared (Herschel 250 &\mu&m) emission correlates spatially very well with a particular narrow velocity (2-3 km/s) component of the Galactic HI. We find no evidence that any of the far-infrared emission associated with these features actually originates in the M81 group. Thus we infer that the associated diffuse optical emission must be due to galactic light back scattered off dust in our galaxy. Ultra-violet observations pick out young stellar associations around M81, but no detectable far-infrared emission. We consider in detail one of the Galactic cirrus features, finding that the far-infrared HI relation breaks down below arc minute scales and that at smaller scales there can be quite large dust temperature variations.

We study a 24\,$\mu$m selected sample of 330 galaxies observed with the Infrared Spectrograph for the 5\,mJy Unbiased Spitzer Extragalactic Survey. We estimate accurate total infrared luminosities by combining mid-IR spectroscopy and mid-to-far infrared photometry, and by utilizing new empirical spectral templates from {\em Spitzer} data. The infrared luminosities of this sample range mostly from 10$^9$L$_\odot$ to $10^{13.5}$L$_\odot$, with 83% in the range 10$^{10}$L$_\odot$$<$L$_{\rm IR}$$<10^{12}$L$_\odot$. The redshifts range from 0.008 to 4.27, with a median of 0.144. The equivalent widths of the 6.2\,$\mu$m aromatic feature have a bimodal distribution. We use the 6.2\,$\mu$m PAH EW to classify our objects as SB-dominated (44%), SB-AGN composite (22%), and AGN-dominated (34%). The high EW objects (SB-dominated) tend to have steeper mid-IR to far-IR spectral slopes and lower L$_{\rm IR}$ and redshifts. The low EW objects (AGN-dominated) tend to have less steep spectral slopes and higher L$_{\rm IR}$ and redshifts. This dichotomy leads to a gross correlation between EW and slope, which does not hold within either group. AGN dominated sources tend to have lower log(L$_{\rm PAH 7.7\mu m}$/L$_{\rm PAH 11.3\mu m}$) ratios than star-forming galaxies, possibly due to preferential destruction of the smaller aromatics by the AGN. The log(L$_{\rm PAH 7.7\mu m}$/L$_{\rm PAH 11.3\mu m}$) ratios for star-forming galaxies are lower in our sample than the ratios measured from the nuclear spectra of nearby normal galaxies, most probably indicating a difference in the ionization state or grain size distribution between the nuclear regions and the entire galaxy. Finally, we provide a calibration relating the monochromatic 5.8, 8, 14 and 24um continuum or Aromatic Feature luminosity to L$_{\rm IR}$ for different types of objects.

By combining Herschel-SPIRE data with archival Spitzer, HI, and CO maps, we investigate the spatial distribution of gas and dust in the two famous grand-design spirals M99 and M100 in the Virgo cluster. Thanks to the unique resolution and sensitivity of the Herschel-SPIRE photometer, we are for the first time able to measure the distribution and extent of cool, submillimetre (submm)-emitting dust inside and beyond the optical radius. We compare this with the radial variation in both the gas mass and the metallicity. Although we adopt a model-independent, phenomenological approach, our analysis provides important insights. We find the dust extending to at least the optical radius of the galaxy and showing breaks in its radial profiles at similar positions as the stellar distribution. The colour indices f350/f500 and f250/f350 decrease radially consistent with the temperature decreasing with radius. We also find evidence of an increasing gas to dust ratio with radius in the outer regions of both galaxies.

We present Herschel-SPIRE observations of the perturbed galaxy NGC4438 in the Virgo cluster. These images reveal the presence of extra-planar dust up to ~4-5 kpc away from the galaxy's disk. The dust closely follows the distribution of the stripped atomic and molecular hydrogen, supporting the idea that gas and dust are perturbed in a similar fashion by the cluster environment. Interestingly, the extra-planar dust lacks a warm temperature component when compared to the material still present in the disk, explaining why it was missed by previous far-infrared investigations. Our study provides evidence for dust stripping in clusters of galaxies and illustrates the potential of Herschel data for our understanding of environmental effects on galaxy evolution.

We present Herschel SPIRE and PACS photometeric observations of the low metallicity (Z ~ 0.35 solar) nearby dwarf galaxy, NGC 1705, in six wavelength bands as part of the Dwarf Galaxy Survey guaranteed time Herschel Key Program. We confirm the presence of two dominant circumnuclear IR-bright regions surrounding the central super star cluster that had been previously noted at mid-IR wavelengths and in the sub-mm by LABOCA. On constructing a global spectral energy distribution using the SPIRE and PACS photometry, in conjunction with archival IR measurements, we note the presence of an excess at sub-mm wavelengths. This excess suggests the presence of a significant cold dust component within NGC 1705 and was modeled as an additional cold component in the SED. Although alternative explanations for the sub-mm excess beyond 350 microns, such as changes to the dust emissivity cannot be ruled out, the most likely explanation for the observed submillimetre excess is that of an additional cold dust component.

The standard method of mapping the interstellar medium in a galaxy, by observing the molecular gas in the CO 1-0 line and the atomic gas in the 21-cm line, is largely limited with current telescopes to galaxies in the nearby universe. In this letter, we use SPIRE observations of the galaxies M99 and M100 to explore the alternative approach of mapping the interstellar medium using the continuum emission from the dust. We have compared the methods by measuring the relationship between the star-formation rate and the surface density of gas in the galaxies. We find the two methods give relationships with a similar dispersion, confirming that observing the continuum emission from the dust is a promising method of mapping the interstellar medium in galaxies.

With appropriate spatial resolution, images of spiral galaxies in thermal infrared (~10 micron and beyond) often reveal a bright central component, distinct from the stellar bulge, superimposed on a disk with prominent spiral arms. ISO and Spitzer studies have shown that much of the scatter in the mid-infrared colors of spiral galaxies is related to changes in the relative importance of these two components, rather than to other modifications, such as the morphological type or star formation rate, that affect the properties of the galaxy as a whole. With the Herschel imaging capability from 70 to 500 micron, we revisit this two-component approach at longer wavelengths, to see if it still provides a working description of the brightness distribution of galaxies, and to determine its implications on the interpretation of global far-infrared properties of galaxies.

We present Herschel/PACS spectroscopic maps of the dwarf galaxy NC4214 observed in 6 far infrared fine-structure lines: [C II] 158mu, [O III] 88mu, [O I] 63mu, [O I] 146mu, [N II] 122mu, and [N II] 205mu. The maps are sampled to the full telescope spatial resolution and reveal unprecedented detail on ~ 150 pc size scales. We detect [C II] emission over the whole mapped area, [O III] being the most luminous FIR line. The ratio of [O III]/[C II] peaks at about 2 toward the sites of massive star formation, higher than ratios seen in dusty starburst galaxies. The [C II]/CO ratios are 20 000 to 70 000 toward the 2 massive clusters, which are at least an order of magnitude larger than spiral or dusty starbursts, and cannot be reconciled with single-slab PDR models. Toward the 2 massive star-forming regions, we find that L[CII] is 0.5 to 0.8% of the LTIR . All of the lines together contribute up to 2% of LTIR . These extreme findings are a consequence of the lower metallicity and young, massive-star formation commonly found in dwarf galaxies. These conditions promote large-scale photodissociation into the molecular reservoir, which is evident in the FIR line ratios. This illustrates the necessity to move to multiphase models applicable to star-forming clusters or galaxies as a whole.

We use Herschel Space Observatory data to place observational constraints on the peak and Rayleigh-Jeans slope of dust emission observed at 70-500 microns in the nearby spiral galaxy M81. We find that the ratios of wave bands between 160 and 500 microns are primarily dependent on radius but that the ratio of 70 to 160 micron emission shows no clear dependence on surface brightness or radius. These results along with analyses of the spectral energy distributions imply that the 160-500 micron emission traces 15-30 K dust heated by evolved stars in the bulge and disc whereas the 70 micron emission includes dust heated by the active galactic nucleus and young stars in star forming regions.

We present the observations of the starburst galaxy M82 taken with the Herschel SPIRE Fourier Transform Spectrometer. The spectrum (194-671 {\mu}m) shows a prominent CO rotational ladder from J = 4-3 to 13-12 emitted by the central region of M82. The fundamental properties of the gas are well constrained by the high J lines observed for the first time. Radiative transfer modeling of these high-S/N 12CO and 13CO lines strongly indicates a very warm molecular gas component at ~500 K and pressure of ~3x10^6 K cm^-3, in good agreement with the H_2 rotational lines measurements from Spitzer and ISO. We suggest that this warm gas is heated by dissipation of turbulence in the interstellar medium (ISM) rather than X-rays or UV flux from the straburst. This paper illustrates the promise of the SPIRE FTS for the study of the ISM of nearby galaxies.

We present the first Herschel PACS and SPIRE images of the low-metallicity galaxy NGC6822 observed from 70 to 500 mu and clearly resolve the HII regions with PACS and SPIRE. We find that the ratio 250/500 is dependent on the 24 mu surface brightness in NGC6822, which would locally link the heating processes of the coldest phases of dust in the ISM to the star formation activity. We model the SEDs of some regions HII regions and less active regions across the galaxy and find that the SEDs of HII regions show warmer ranges of dust temperatures. We derive very high dust masses when graphite is used in our model to describe carbon dust. Using amorphous carbon, instead, requires less dust mass to account for submm emission due to its lower emissivity properties. This indicates that SED models including Herschel constraints may require different dust properties than commonly used.

We present infrared colours (in the 25-500 mic spectral range) and UV to radio continuum spectral energy distributions of a sample of 51 nearby galaxies observed with SPIRE on Herschel. The observed sample includes all morphological classes, from quiescent ellipticals to active starbursts. Active galaxies have warmer colour temperatures than normal spirals. In ellipticals hosting a radio galaxy, the far-infrared (FIR) emission is dominated bynthe synchrotron nuclear emission. The colour temperature of the cold dust is higher in quiescent E-S0a than in star-forming systems probably because of the different nature of their dust heating sources (evolved stellar populations, X-ray, fast electrons) and dust grain properties. In contrast to the colour temperature of the warm dust, the f350/f500 index sensitive to the cold dust decreases with star formation and increases with metallicity, suggesting an overabundance of cold dust or an emissivity parameter beta<2 in low metallicity, active systems.

We present Herschel-SPIRE observations at 250-500um of the giant elliptical galaxy M86 and examine the distribution of the resolved cold dust emission and its relation with other galactic tracers. The SPIRE images reveal three dust components: emission from the central region; a dust lane extending north-south; and a bright emission feature 10kpc to the south-east. We estimate that approximately 10^6 solar masses of dust is spatially coincident with atomic and ionized hydrogen, originating from stripped material from the nearby spiral NGC4438 due to recent tidal interactions with M86. The gas-to-dust ratio of the cold gas component ranges from ~20-80. We discuss the different heating mechanisms for the dust features.

M82 is a unique representative of a whole class of galaxies, starbursts with superwinds, in the Very Nearby Galaxy Survey with Herschel. In addition, its interaction with the M81 group has stripped a significant portion of its interstellar medium from its disk. SPIRE maps now afford better characterization of the far-infrared emission from cool dust outside the disk, and sketch a far more complete picture of its mass distribution and energetics than previously possible. They show emission coincident in projection with the starburst wind and in a large halo, much more extended than the PAH band emission seen with Spitzer. Some complex substructures coincide with the brightest PAH filaments, and others with tidal streams seen in atomic hydrogen. We subtract the far-infrared emission of the starburst and underlying disk from the maps, and derive spatially-resolved far-infrared colors for the wind and halo. We interpret the results in terms of dust mass, dust temperature, and global physical conditions. In particular, we examine variations in the dust physical properties as a function of distance from the center and the wind polar axis, and conclude that more than two thirds of the extraplanar dust has been removed by tidal interaction, and not entrained by the starburst wind.

The Herschel Reference Survey is a guaranteed time Herschel key project and will be a benchmark study of dust in the nearby universe. The survey will complement a number of other Herschel key projects including large cosmological surveys that trace dust in the distant universe. We will use Herschel to produce images of a statistically-complete sample of 323 galaxies at 250, 350 and 500 micron. The sample is volume-limited, containing sources with distances between 15 and 25 Mpc and flux limits in the K-band to minimize the selection effects associated with dust and with young high-mass stars and to introduce a selection in stellar mass. The sample spans the whole range of morphological types (ellipticals to late-type spirals) and environments (from the field to the centre of the Virgo Cluster) and as such will be useful for other purposes than our own. We plan to use the survey to investigate (i) the dust content of galaxies as a function of Hubble type, stellar mass and environment, (ii) the connection between the dust content and composition and the other phases of the interstellar medium and (iii) the origin and evolution of dust in galaxies. In this paper, we describe the goals of the survey, the details of the sample and some of the auxiliary observing programs that we have started to collect complementary data. We also use the available multi-frequency data to carry out an analysis of the statistical properties of the sample.

Coupled, multi-field models of inflation can provide several attractive features unavailable in the case of a single inflaton field. These models have a rich dynamical structure resulting from the interaction of the fields and their associated fluctuations. We present a formalism to study the nonequilibrium dynamics of coupled scalar fields. This formalism solves the problem of renormalizing interacting models in a transparent way using dimensional regularization. The evolution is generated by a renormalized effective Lagrangian which incorporates the dynamics of the mean fields and their associated fluctuations at one-loop order. We apply our method to two problems of physical interest: (i) a simple two-field model which exemplifies applications to reheating in inflation, and (ii) a supersymmetric hybrid inflation model. This second case is interesting because inflation terminates via a smooth phase transition which gives rise to a spinodal instability in one of the fields. We study the evolution of the zero mode of the fields and the energy density transfer to the fluctuations from the mean fields. We conclude that back reaction effects can be significant over a wide parameter range. In particular for the supersymmetric hybrid model we find that particle production can be suppressed due to these effects.

We investigate the out of equilibrium dynamics of global chiral supersymmetry at finite energy density. We concentrate on two specific models. The first is the massive Wess-Zumino model which we study in a selfconsistent one-loop approximation. We find that for energy densities above a certain threshold, the fields are driven dynamically to a point in field space at which the fermionic component of the superfield is massless. The state, however is found to be unstable, indicating a breakdown of the one-loop approximation. To investigate further, we consider an O(N) massive chiral model which is solved exactly in the large $N$ limit. For sufficiently high energy densities, we find that for late times the fields reach a nonperturbative minimum of the effective potential degenerate with the perturbative minimum. This minimum is a true attractor for O(N) invariant states at high energy densities, and this provides a mechanism for determining which of the otherwise degenerate vacua is chosen by the dynamics. The final state for large energy density is a cloud of massless particles (both bosons and fermions) around this new nonperturbative supersymmetric minimum. By introducing boson masses which softly break the supersymmetry, we demonstrate a see-saw mechanism for generating small fermion masses. We discuss some of the cosmological implications of our results.

The accelerated expansion of the Universe measured by high redshift Type Ia Supernovae observations is explained using the non-equilibrium dynamics of naturally soft boson fields. Spinodal instabilities inevitably present in such systems yield a very attractive mechanism for arriving at the required equation of state at late times, while satisfying all the known constraints on models of quintessence.

We study the out of equilibrium dynamics of inflationary phase transitions and compute the resulting spectrum of metric perturbations relevant to observation. We show that simple single field models of inflation may produce an adiabatic perturbation spectrum with a blue spectral tilt and that the precise spectrum depends on initial conditions at the outset of inflation.

We analyse the dynamics of spinodal decomposition in inflationary cosmology using the closed time path formalism of out of equilibrium quantum field theory combined with the non-perturbative Hartree approximation. In addition to a general analysis, we compute the detailed evolution of two inflationary models of particular importance: lambda Phi^4 new inflation and natural inflation. We compute the metric fluctuations resulting from inflationary phase transitions in the slow roll approximation, showing that there exists a regime for which quantum fluctuations of the inflaton field result in a significant deviation in the predictions of the spectrum of primordial density perturbations from standard results. We provide case examples for which a blue tilt to the power spectrum (i.e. n_s > 1) results from the evolution of a single inflaton field, and demonstrate that field fluctuations may result in a scalar amplitude of fluctuations significantly below standard predictions, resulting in a slight alleviation of the inflationary fine tuning problem. We show explicitly that the metric perturbation spectrum resulting from inflation depends upon the state at the outset of the inflationary phase.

Out-of-equilibrium, non-perturbative, quantum effects significantly modify the standard picture of inflation in a wide class of models including new, natural, and hybrid inflation. We find that the quantum evolution of a single real inflaton field may be modeled by a classical theory of two homogeneous scalar fields. We briefly discuss the important observational consequences that are expected to result.

For simple inflationary models, we provide a consistent and complete scheme by which the macro-physical details of early universe inflation may be determined explicitly from the underlying micro-physical theory. We examine inflationary dynamics within the context of a fully renormalized, non-perturbative, and non-equilibrium quantum field theory using the closed time path formalism. We study the non-perturbative dynamics using the self-consistent Hartree approximation and the large N limit. In addition, we include gravitation dynamically by means of the semi-classical approximation to Einstein gravity, allowing the cosmological geometry to be determined self-consistently by the evolution of the particle fields. We develop a simple and numerically implementable regularization and renormalization scheme for which the renormalization counterterms contain no explicit dependence on the initial state. We concern ourselves first with the dynamics of new inflation models, showing how the dynamics become quantum fluctuation driven. We then describe how these quantum fluctuations may be reassembled into an effective field which behaves classically and acts as the source for the gravitational background. This allows us to determine the spectrum of primordial density perturbations in a consistent manner. We also examine explosive particle production after inflation through the processes of spinodal decomposition and parametric amplification, showing that the non-perturbative backreaction on the inflaton field prevents catastrophic particle production. We find a strong dependence of particle production on the expansion rate and on the symmetry properties of the inflaton.

The energy and time scales during the inflationary stage of the universe calls for an out of equilibrium quantum field treatment. Moreover, the high energy densities involved make necessary the use of non-perturbative approaches as large N and Hartree methods. We start these lectures by introducing such non-perturbative out of equilibrium methods in cosmological universes. We discuss the renormalization procedure and the choice of initial conditions. We then study the nonlinear dynamics of quantum fields in matter and radiation dominated FRW and de Sitter universes. For a variety of initial conditions, we compute the evolution of the inflaton,its quantum fluctuations and the equation of state. We investigate the explosive particle production due to spinodal unstabilities and parametric amplification in FRW and de Sitter universes with and without symmetry breaking.We find that the particle production is sensitive to the expansion of the universe.For symmetry breaking scenarios, we determine generic late time fields behavior for FRW and deSitter cosmologies.We find that quantum fluctuations damp in FRW as the square of the scale factor while the order parameter approaches a minimum of the potential at the same rate.We con- sider an O(N) inflaton model coupled self-consistently to gravity in the semi- classical approximation for a `new inflation' scenario. We find that spinodal instabilities drive the growth of non-perturbatively large quantum fluctuations which shut off inflation. These fluctuations assemble with the inflaton zero mode yielding a new effective field that actually rolls down behaving classi- cally. We compute the amplitude and index for scalar density and tensor perturbations. In all models of this type the spinodal instabilities produce a `red' spectrum of primordial scalar density perturbations.

We consider an O(N) model coupled self-consistently to gravity in the semiclassical approximation, where the field is subject to `new inflation' type initial conditions. We study the dynamics self-consistently and non-perturbatively with non-equilibrium field theory methods in the large N limit. We find that spinodal instabilities drive the growth of non-perturbatively large quantum fluctuations which shut off the inflationary growth of the scale factor. We find that a very specific combination of these large fluctuations plus the inflaton zero mode assemble into a new effective field. This new field behaves classically and it is the object which actually rolls down. We show how this reinterpretation saves the standard picture of how metric perturbations are generated during inflation and that the spinodal growth of fluctuations dominates the time dependence of the Bardeen variable for superhorizon modes during inflation. We compute the amplitude and index for the spectrum of scalar density and tensor perturbations and argue that in all models of this type the spinodal instabilities are responsible for a `red' spectrum of primordial scalar density perturbations. A criterion for the validity of these models is provided and contact with the reconstruction program is established validating some of the results within a non-perturbative framework. The decoherence aspects and the quantum to classical transition through inflation are studied in detail by following the full evolution of the density matrix and relating the classicality of cosmological perturbations to that of long-wavelength matter fluctuations.

The physics of the inflationary universe requires the study of the out of equilibrium evolution of quantum fields in curved spacetime. We present the evolution for both the geometry and the matter (described by the quantum inflaton field) by means of the non-perturbative large N limit combined with semi-classical gravitational dynamics including the back-reaction of quantum fluctuations self-consistently for a new inflation scenario. We provide a criterion for the validity of the classical approximation and a full analysis of the case in which spinodal quantum fluctuations drive the evolution of the scale factor. Under carefully determined conditions, we show that the full field equations may be well approximated by those of a single composite field which obeys the classical equation of motion in all cases. The de Sitter stage is found to be followed by a matter dominated phase. We compute the spectrum of scalar density perturbations and argue that the spinodal instabilities are responsible for a `red' spectrum with more power at longer wavelengths. A criterion for the validity of these models is provided and contact with the reconstruction program is established.

We study the non-linear dynamics of quantum fields in matter and radiation dominated universes, using the non-equilibrium field theory approach combined with the non-perturbative Hartree and the large N approximations. We examine the phenomenon of explosive particle production due to spinodal instabilities and parametric amplification in expanding universes with and without symmetry breaking. For a variety of initial conditions, we compute the evolution of the inflaton, its quantum fluctuations, and the equation of state. We find explosive growth of quantum fluctuations, although particle production is somewhat sensitive to the expansion of the universe. In the large N limit for symmetry breaking scenarios, we determine generic late time solutions for any flat Friedman-Robertson-Walker cosmology. We also present a complete and numerically implementable renormalization scheme for the equation of motion and the energy momentum tensor in flat FRW cosmologies. In this scheme the renormalization constants are independent of time and of the initial conditions.

The nonlinear time evolution of the quantum fields is studied in the O(N) model for large N in a radiation dominated FRW universe, with a view towards the phenomenon of explosive particle production due to either spinodal instabilities or parametric amplification, i.e. preheating. Quantum backreaction effects due to the produced particles are included consistently within the large N approximation. We find that preheating persists when the expansion is included, although the amount of particle production is reduced compared to the values found in Minkowski space. We also see that the behavior of the evolving zero mode is very different from that in Minkowski space, though the late time behavior in all cases is determined by a sum rule that implies the existence of Goldstone bosons in the final state.