We present the first nonparametric morphological analysis of a set of spiral galaxies from UV to submm wavelengths. Our study is based on high-quality multi-wavelength imaging for nine well-resolved spiral galaxies from the DustPedia database, combined with nonparametric morphology indicators calculated in a consistent way using the {\tt{StatMorph}} package. We measure the half-light radius, the concentration index, the asymmetry index, the smoothness index, the Gini coefficient and the $M_{20}$ indicator in various wavebands from UV to submm wavelengths, as well as in stellar mass, dust mass and star formation rate maps. We find that the interstellar dust in galaxies is distributed in a more extended, less centrally concentrated, more asymmetric, and more clumpy way than the stars. This is particularly evident when comparing morphological indicators based on the stellar mass and dust mass maps. This should serve as a warning sign against treating the dust in galaxies as a simple smooth component. We argue that the nonparametric galaxy morphology of galaxies from UV to submm wavelengths is an interesting test for cosmological hydrodynamics simulations.

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.

We present the detection of multiple CO line transitions with ALMA in a few tens of infrared-selected galaxies on and above the main sequence at z=1.1-1.7. We reliably detected the emission of CO(5-4), CO(2-1), and CO(7-6)+[CI](2-1) in 50, 33, and 13 galaxies, respectively, and we complemented this information with available CO(4-3) and [CI](1-0) fluxes for part of the sample, and modeling of the optical-to-mm SEDs. We retrieve a quasi-linear relation between LIR and CO(5-4) or CO(7-6) for main-sequence galaxies and starbursts, corroborating the hypothesis that these transitions can be used as SFR tracers. We find the CO excitation to steadily increase as a function of the star formation efficiency, the mean intensity of the radiation field warming the dust, the surface density of SFR, and, less distinctly, with the distance from the main sequence. This adds to the tentative evidence for higher excitation of the CO+[CI] SLED of starbursts relative to that for main-sequence objects, where the dust opacities play a minor role in shaping the high-J CO transitions in our sample. However, the distinction between the average SLED of upper main-sequence and starburst galaxies is blurred, driven by a wide variety of intrinsic shapes. LVG radiative transfer modeling demonstrates the existence of a highly excited component that elevates the CO SLED of high-redshift main-sequence and starbursting galaxies above the typical values observed in the disk of the Milky Way. This excited component is dense and it encloses ~50% of the total molecular gas mass in main-sequence objects. We interpret the observed trends involving the CO excitation as mainly driven by a combination of large SFRs and compact sizes, naturally connected with enhanced dense molecular gas fractions and higher dust and gas temperatures, due to increasing UV radiation fields, cosmic ray rates, and dust/gas coupling. [Abridged]

Dust emission, an important diagnostic of star formation and ISM mass throughout the Universe, can be powered by sources unrelated to ongoing star formation. In the framework of the DustPedia project we have set out to disentangle the radiation of the ongoing star formation from that of the older stellar populations. This is done through detailed, 3D radiative transfer simulations of face-on spiral galaxies. In this particular study, we focus on NGC 1068, which contains an active galactic nucleus (AGN). The effect of diffuse dust heating by AGN (beyond the torus) was so far only investigated for quasars. This additional dust heating source further contaminates the broadband fluxes on which classic galaxy modelling tools rely to derive physical properties. We aim to fit a realistic model to the observations of NGC 1068 and quantify the contribution of the several dust heating sources. Our model is able to reproduce the global spectral energy distribution of the galaxy. It matches the resolved optical and infrared images fairly well, but deviates in the UV and the submm. We find a strong wavelength dependency of AGN contamination to the broadband fluxes. It peaks in the MIR, drops in the FIR, but rises again at submm wavelengths. We quantify the contribution of the dust heating sources in each 3D dust cell and find a median value of 83% for the star formation component. The AGN contribution is measurable at the percentage level in the disc, but quickly increases in the inner few 100 pc, peaking above 90%. This is the first time the phenomenon of an AGN heating the diffuse dust beyond its torus is quantified in a nearby star-forming galaxy. NGC 1068 only contains a weak AGN, meaning this effect can be stronger in galaxies with a more luminous AGN. This could significantly impact the derived star formation rates and ISM masses for such systems.

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.

Interstellar dust absorbs stellar light very efficiently and thus shapes the energetic output of galaxies. Studying the impact of different stellar populations on the dust heating remains hard because it requires decoupling the relative geometry of stars and dust, and involves complex processes as scattering and non-local dust heating. We aim to constrain the relative distribution of dust and stellar populations in the spiral galaxy M81 and create a realistic model of the radiation field that describes the observations. Investigating the dust-starlight interaction on local scales, we want to quantify the contribution of young and old stellar populations to the dust heating. We aim to standardise the setup and model selection of such inverse radiative transfer simulations so this can be used for comparable modelling of other nearby galaxies. We present a semi-automated radiative transfer modelling pipeline that implements the necessary steps such as the geometric model construction and the normalisation of the components through an optimisation routine. We use the Monte Carlo radiative transfer code SKIRT to calculate a self-consistent, panchromatic model of the interstellar radiation field. By looking at different stellar populations independently, we can quantify to what extent different stellar age populations contribute to the dust heating. Our method takes into account the effects of non-local heating. We obtain a realistic 3D radiative transfer model of the face-on galaxy M81. We find that only 50.2\% of the dust heating can be attributed to young stellar populations. We confirm a tight correlation between the specific star formation rate and the heating fraction by young stellar populations, both in sky projection and in 3D, also found for radiative transfer models of M31 and M51. We conclude that... (abridged)

Context: Dust in late-type galaxies in the local Universe is responsible for absorbing approximately one third of the energy emitted by stars. It is often assumed that dust heating is mainly attributable to the absorption of UV and optical photons emitted by the youngest (<= 100 Myr) stars. Consequently, thermal re-emission by dust at FIR wavelengths is often linked to the star-formation activity of a galaxy. However, several studies argue that the contribution to dust heating by much older stars might be more significant. Advances in radiation transfer (RT) simulations finally allow us to actually quantify the heating mechanisms of diffuse dust by the stellar radiation field. Aims: As one of the main goals in the DustPedia project, we have constructed detailed 3D stellar and dust RT models for nearby galaxies. We analyse the contribution of the different stellar populations to the dust heating in four face-on barred galaxies: NGC1365, M83, M95, and M100. We aim to quantify the fraction directly related to young stars, both globally and on local scales, and to assess the influence of the bar on the heating fraction. Results: We derive global attenuation laws for each galaxy and confirm that galaxies of high sSFR have shallower attenuation curves and weaker UV bumps. On average, 36.5% of the bolometric luminosity is absorbed by dust. We report a clear effect of the bar structure on the radial profiles of the dust-heating fraction by the young stars, and the dust temperature. We find that the young stars are the main contributors to the dust heating, donating, on average ~59% of their luminosity to this purpose throughout the galaxy. This dust-heating fraction drops to ~53% in the bar region and ~38% in the bulge region where the old stars are the dominant contributors to the dust heating. We also find a strong link between the heating fraction by the young stars and the sSFR.

We compare the spectral energy distributions (SEDs) and inferred physical properties for simulated and observed galaxies at low redshift. We exploit UV-submillimetre mock fluxes of ~7000 z=0 galaxies from the EAGLE suite of cosmological simulations, derived using the radiative transfer code SKIRT. We compare these to ~800 observed galaxies in the UV-submillimetre range, from the DustPedia sample of nearby galaxies. To derive global properties, we apply the SED fitting code CIGALE consistently to both data sets, using the same set of ~80 million models. The results of this comparison reveal overall agreement between the simulations and observations, both in the SEDs and in the derived physical properties, with a number of discrepancies. The optical and far-infrared regimes, and the scaling relations based upon the global emission, diffuse dust and stellar mass, show high levels of agreement. However, the mid-infrared fluxes of the EAGLE galaxies are overestimated while the far-UV domain is not attenuated enough, compared to the observations. We attribute these discrepancies to a combination of galaxy population differences between the samples, and limitations in the subgrid treatment of star-forming regions in the EAGLE-SKIRT post-processing recipe. Our findings show the importance of detailed radiative transfer calculations and consistent comparison, and provide suggestions for improved numerical models.

We report ALMA observations of the neutral atomic carbon transitions [CI] and multiple CO lines in a sample of $\sim30$ main sequence galaxies at $z\sim1$, including novel information on [CI](2-1) and CO(7-6) for 7 of such normal objects. We complement our observations with a collection of $>200$ galaxies with coverage of similar transitions, spanning the $z=0-4$ redshift interval and a variety of ambient conditions from local to high-redshift starbursts. We find systematic variations in the [CI]/IR and [CI]/high-$J$ ($J=7$) CO luminosity ratios among the various samples. We interpret these differences as increased dense molecular gas fractions and star formation efficiencies in the strongest high-redshift starbursts with respect to normal main sequence galaxies. We further report constant $L'_{\rm [CI]2-1}$/$L'_{\rm [CI]1-0}$ ratios across the galaxy populations and redshifts, suggesting that gas temperatures $T_{\rm exc}$ traced by [CI] do not strongly vary. We find only a mild correlation with $T_{\rm dust}$ and that, generally, $T_{\rm exc} \lesssim T_{\rm dust}$. We fit the line ratios with classical PDR models, retrieving consistently larger densities and intensities of the UV radiation fields in submm galaxies than in main sequence and local objects. However, these simple models fall short in representing the complexity of a multiphase interstellar medium and should be treated with caution. Finally, we compare our observations with the Santa Cruz semi-analytical model of galaxy evolution, recently extended to simulate submm emission. While we confirm the success in reproducing the CO lines, we find systematically larger [CI] luminosities at fixed IR luminosity than predicted theoretically. This highlights the necessity of improving our understanding of the mechanisms regulating the [CI] emission on galactic scales. We release our data compilation to the community.

High-resolution ($<$0.1pc) ALMA observations of the 30Dor-10 molecular cloud 15pc north of R136 are presented. The $^{12}$CO 2-1 emission morphology contains clumps near the locations of known mid-infrared massive protostars, as well as a series of parsec-long filaments oriented almost directly towards R136. There is elevated kinetic energy (linewidths at a given size scale) in 30Dor-10 compared to other LMC and Galactic star formation regions, consistent with large scale energy injection to the region. Analysis of the cloud substructures is performed by segmenting emission into disjoint approximately round "cores" using clumpfind, by considering the hierarchical structures defined by isointensity contours using dendrograms, and by segmenting into disjoint long thin "filaments" using Filfinder. Identified filaments have widths $\sim$0.1pc. The inferred balance between gravity and kinematic motions depends on the segmentation method: Entire objects identified with clumpfind are consistent with free-fall collapse or virial equilibrium with moderate external pressure, whereas many dendrogram-identified parts of hierarchical structures have higher mass surface densities $\Sigma_{LTE}$ than if gravitational and kinetic energies were in balance. Filaments have line masses that vary widely compared to the critical line mass calculated assuming thermal and nonthermal support. Velocity gradients in the region do not show any strong evidence for accretion of mass along filaments. The upper end of the "core" mass distribution is consistent with a power-law with the same slope as the stellar initial mass function.

We have conducted ALMA CO isotopes and 1.3 mm continuum observations toward filamentary molecular clouds of the N159W-South region in the Large Magellanic Cloud with an angular resolution of $\sim$0"25 ($\sim$0.07 pc). Although the previous lower-resolution ($\sim$1") ALMA observations revealed that there is a high-mass protostellar object at an intersection of two line-shaped filaments in $^{13}$CO with the length scale of $\sim$10 pc, the spatially resolved observations, in particular, toward the highest column density part traced by the 1.3 mm continuum emission, the N159W-South clump, show complicated hub-filamentary structures. We also discovered that there are multiple protostellar sources with bipolar outflows along the massive filament. The redshifted/blueshifted components of the $^{13}$CO emission around the massive filaments/protostars have complementary distributions, which is considered to be a possible piece of evidence for a cloud-cloud collision. We propose a new scenario in which the supersonically colliding gas flow triggers the formation of both the massive filament and protostars. This is a modification of the earlier scenario of cloud-cloud collision, by Fukui et al., that postulated the two filamentary clouds occur prior to the high-mass star formation. A recent theoretical study of the shock compression in colliding molecular flows by Inoue et al. demonstrates that the formation of filaments with hub structure is a usual outcome of the collision, lending support for the present scenario. The theory argues that the filaments are formed as dense parts in a shock compressed sheet-like layer, which resembles $"$an umbrella with pokes.$"$

We present the ALMA observations of CO isotopes and 1.3 mm continuum emission toward the N159E-Papillon Nebula in the Large Magellanic Cloud (LMC). The spatial resolution is 0"25-0"28 (0.06-0.07 pc), which is a factor of 3 higher than the previous ALMA observations in this region. The high resolution allowed us to resolve highly filamentary CO distributions with typical widths of $\sim$0.1 pc (full width half maximum) and line masses of a few 100 $M_{\odot}$ pc$^{-1}$. The filaments (more than ten in number) show an outstanding hub-filament structure emanating from the nebular center toward the north. We identified for the first time two massive protostellar outflows of $\sim$10$^4$ yr dynamical age along one of the most massive filaments. The observations also revealed several pillar-like CO features around the Nebula. The H II region and the pillars have a complementary spatial distribution and the column density of the pillars is an order of magnitude higher than that of the pillars in the Eagle nebula (M16) in the Galaxy, suggesting an early stage of pillar formation with an age younger than $\sim$10$^5$ yr. We suggest that a cloud-cloud collision triggered the formation of the filaments and protostar within the last $\sim$2 Myr. It is possible that the collision is more recent, as part of the kpc-scale H I flows come from the tidal interaction resulting from the close encounter between the LMC and SMC $\sim$200 Myr ago as suggested for R136 by Fukui et al.

The purpose of this work is the characterization of the main scaling relations between all the ISM components (dust, atomic/molecular/total gas), gas-phase metallicity, and other galaxy properties, such as Mstar and galaxy morphology, for late-type galaxies in the Local Universe. This study is performed by extracting late-type galaxies from the entire DustPedia sample and by exploiting the large and homogeneous dataset available thanks to the DustPedia project. The sample consists of 436 galaxies with morphological stage from T = 1 to 10, Mstar from 6 x 10^7 to 3 x 10^11 Msun, SFR from 6 x 10^(-4) to 60 Msun/yr, and 12 + log(O/H) from 8 to 9.5. The scaling relations involving the molecular gas are studied by assuming both a constant and a metallicity-dependent XCO. The analysis has been performed by means of the survival analysis technique. We confirm that the dust mass correlates very well with the total gas mass, and find -- for the first time -- that the dust mass correlates better with the atomic gas mass than the molecular one. We characterize important mass ratios such as gas fraction, molecular-to-atomic gas mass ratio, dust-to-total gas mass ratio (DGR), and dust-to-stellar mass ratio. Only the assumption of a metallicity-dependent XCO reproduces the expected decrease of the DGR with increasing morphological stage and decreasing gas-phase metallicity, with a slope of about 1. DGR, gas-phase metallicity, and the dust-to-stellar mass ratio are, for our galaxy sample, directly linked to the galaxy morphology. The molecular-to-atomic gas mass ratio and the DGR show a positive correlation for low molecular gas fractions, but for molecular gas rich galaxies this trend breaks down. This trend has never been found previously, to our knowledge. It provides new constraints for theoretical models of galaxy evolution and a reference for high-redshift studies.

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

Observations of evolution in the dust-to-metal ratio allow us to constrain the dominant dust processing mechanisms. In this work, we present a study of the dust-to-metal and dust-to-gas ratios in a sub-sample of ~500 DustPedia galaxies. Using literature and MUSE emission line fluxes, we derived gas-phase metallicities (oxygen abundances) for over 10000 individual regions and determine characteristic metallicities for each galaxy. We study how the relative dust, gas, and metal contents of galaxies evolve by using metallicity and gas fraction as proxies for evolutionary state. The global oxygen abundance and nitrogen-to-oxygen ratio are found to increase monotonically as galaxies evolve. Additionally, unevolved galaxies (gas fraction > 60%, metallicity 12 + log(O/H) < 8.2) have dust-to-metal ratios that are about a factor of 2.1 lower (a factor of six lower for galaxies with gas fraction > 80%) than the typical dust-to-metal ratio (Md/MZ ~ 0.214) for more evolved sources. However, for high gas fractions, the scatter is larger due to larger observational uncertainties as well as a potential dependence of the dust grain growth timescale and supernova dust yield on local conditions and star formation histories. We find chemical evolution models with a strong contribution from dust grain growth describe these observations reasonably well. The dust-to-metal ratio is also found to be lower for low stellar masses and high specific star formation rates (with the exception of some sources undergoing a starburst). Finally, the metallicity gradient correlates weakly with the HI-to-stellar mass ratio, the effective radius and the dust-to-stellar mass ratio, but not with stellar mass.

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.

The dust mass absorption coefficient, $\kappa_{d}$, is the conversion function used to infer physical dust masses from observations of dust emission. However, it is notoriously poorly constrained, and it is highly uncertain how it varies, either between or within galaxies. Here we present the results of a proof-of concept study, using the DustPedia data for two nearby face-on spiral galaxies M74 (NGC 628) and M83 (NGC 5236), to create the first ever maps of $\kappa_{d}$ in galaxies. We determine $\kappa_{d}$ using an empirical method that exploits the fact that the dust-to-metals ratio of the interstellar medium is constrained by direct measurements of the depletion of gas-phase metals. We apply this method pixel-by-pixel within M74 and M83, to create maps of $\kappa_{d}$. We also demonstrate a novel method of producing metallicity maps for galaxies with irregularly-sampled measurements, using the machine learning technique of Gaussian process regression. We find strong evidence for significant variation in $\kappa_{d}$. We find values of $\kappa_{d}$ at 500 $\mu$m spanning the range 0.11-0.25 ${\rm m^{2}\,kg^{-1}}$ in M74, and 0.15-0.80 ${\rm m^{2}\,kg^{-1}}$ in M83. Surprisingly, we find that $\kappa_{d}$ shows a distinct inverse correlation with the local density of the interstellar medium. This inverse correlation is the opposite of what is predicted by standard dust models. However, we find this relationship to be robust against a large range of changes to our method - only the adoption of unphysical or highly unusual assumptions would be able to suppress it.

The zodiacal dust complex, a population of dust and small particles that pervades the Solar System, provides important insight into the formation and dynamics of planets, comets, asteroids, and other bodies. Here we present a new set of data obtained using a novel technique: direct measurements of momentum transfer to a spacecraft from individual particle impacts. This technique is made possible by the extreme precision of the instruments flown on the LISA Pathfinder spacecraft, a technology demonstrator for a future space-based gravitational wave observatory that operated near the first Sun-Earth Lagrange point from early 2016 through Summer of 2017. Using a simple model of the impacts and knowledge of the control system, we show that it is possible to detect impacts and measure properties such as the transferred momentum (related to the particle's mass and velocity), direction of travel, and location of impact on the spacecraft. In this paper, we present the results of a systematic search for impacts during 4348 hours of Pathfinder data. We report a total of 54 candidates with momenta ranging from 0.2$\,\mu\textrm{Ns}$ to 230$\,\mu\textrm{Ns}$. We furthermore make a comparison of these candidates with models of micrometeoroid populations in the inner solar system including those resulting from Jupiter-family comets, Oort-cloud comets, Hailey-type comets, and Asteroids. We find that our measured population is consistent with a population dominated by Jupiter-family comets with some evidence for a smaller contribution from Hailey-type comets. This is in agreement with consensus models of the zodiacal dust complex in the momentum range sampled by LISA Pathfinder.

Aims: We compare the far-infrared to sub-millimetre dust emission properties measured in high Galactic latitude cirrus with those determined in a sample of 204 late-type DustPedia galaxies. The aim is to verify if it is appropriate to use Milky Way dust properties to derive dust masses in external galaxies. Methods: We used Herschel observations and atomic and molecular gas masses to estimate the disc-averaged dust emissivity at 250 micrometres, and from this, the absorption cross section per H atom and per dust mass. The emissivity requires one assumption, which is the CO-to-H_2 conversion factor, and the dust temperature is additionally required for the absorption cross section per H atom; yet another constraint on the dust-to-hydrogen ratio D/H, depending on metallicity, is required for the absorption cross section dust mass. Results: We find epsilon(250) = 0.82 +/- 0.07 MJy sr^-1 (1E20 H cm^-2)^-1 for galaxies with 4 < F(250)/F(500) < 5. This depends only weakly on the adopted CO-to-H_2 conversion factor. The value is almost the same as that for the Milky Way at the same colour ratio. Instead, for F(250)/F(500) > 6, epsilon(250) is lower than predicted by its dependence on the heating conditions. The reduction suggests a variation in dust emission properties for spirals of earlier type, higher metallicity, and with a higher fraction of molecular gas. When the standard emission properties of Galactic cirrus are used for these galaxies, their dust masses might be underestimated by up to a factor of two. Values for the absorption cross sections at the Milky Way metallicity are also close to those of the cirrus. Mild trends of the absorption cross sections with metallicity are found, although the results depend on the assumptions made.

The future of exoplanet detection lies in the mid-infrared (MIR). The MIR region contains the blackbody peak of both hot and habitable zone exoplanets, making the contrast between starlight and planet light less extreme. It is also the region where prominent chemical signatures indicative of life exist, such as ozone at 9.7 microns. At a wavelength of 4 microns the difference in emission between an Earth-like planet and a star like our own is 80 dB. However a jovian planet, at the same separation exhibits 60 dB of contrast, or only 20 dB if it is hot due to its formation energy or being close to its host star. A two dimensional nulling interferometer, made with chalcogenide glass, has been measured to produce a null of 20 dB, limited by scattered light. Measures to increase the null depth to the theoretical limit of 60 dB are discussed.

With an aim of probing the physical conditions and excitation mechanisms of warm molecular gas in individual star-forming regions, we performed Herschel SPIRE FTS observations of 30 Doradus in the LMC. In our FTS observations, important FIR cooling lines in the ISM, including CO J=4-3 to 13-12, [CI] 370 micron, and [NII] 205 micron, were clearly detected. In combination with ground-based CO data, we then constructed CO spectral line energy distributions (SLEDs) on 10 pc scales over a 60 pc x 60 pc area and found that the shape of the observed CO SLEDs considerably changes across 30 Doradus, e.g., the peak transition varies from J=6-5 to 10-9, while the slope characterized by the high-to-intermediate J ratio ranges from 0.4 to 1.8. To examine the source(s) of these variations in CO transitions, we analyzed the CO observations, along with [CII] 158 micron, [CI] 370 micron, [OI] 145 micron, H2 0-0 S(3), and FIR luminosity data, using state-of-the-art models of PDRs and shocks. Our detailed modeling showed that the observed CO emission likely originates from highly-compressed (thermal pressure ~ 1e7-1e9 K cm-3) clumps on 0.7-2 pc scales, which could be produced by either UV photons (UV radiation field ~ 1e3-1e5 Mathis fields) or low-velocity C-type shocks (pre-shock medium density ~ 1e4-1e6 cm-3 and shock velocity ~ 5-10 km s-1). Considering the stellar content in 30 Doradus, however, we tentatively excluded the stellar origin of CO excitation and concluded that low-velocity shocks driven by kpc scale processes (e.g., interaction between the Milky Way and the Magellanic Clouds) are likely the dominant source of heating for CO. The shocked CO-bright medium was then found to be warm (temperature ~ 100-500 K) and surrounded by a UV-regulated low pressure component (a few (1e4-1e5) K cm-3) that is bright in [CII] 158 micron, [CI] 370 micron, [OI] 145 micron, and FIR dust continuum emission.

We use a sub-set of the DustPedia galaxy sample (461 galaxies) to investigate the effect the environment has had on galaxies. We consider Virgo cluster and field samples and also assign a density contrast parameter to each galaxy, as defined by the local density of SDSS galaxies. We consider their chemical evolution (using M_{Dust}/M_{Baryon} and M_{Gas}/M_{Baryon}), their specific star formation rate (SFR/M_{Stars}), star formation efficiency (SFR/M_{Gas}), stars-to-dust mass ratio (M_{Stars}/M_{Dust}), gas-to-dust mass ratio (M_{Gas}/M_{Dust}) and the relationship between star formation rate per unit mass of dust and dust temperature (SFR/M_{Dust} and T_{Dust}). Late type galaxies (later than Sc) in all of the environments can be modelled using simple closed box chemical evolution and a simple star formation history (SFR(t) \propto t\exp{-t/\tau}). For earlier type galaxies the physical mechanisms that give rise to their properties are clearly much more varied and require a more complicated model (mergers, gas in or outflow). However, we find little or no difference in the properties of galaxies of the same morphological type within the cluster, field or with different density contrasts. It appears that it is morphology, how and whenever this is laid down, and consistent internal physical processes that primarily determine the derived properties of galaxies in the DustPedia sample and not processes related to differences in the local environment.

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.

Within the framework of the DustPedia project we investigate the properties of cosmic dust and its interaction with the stellar radiation (originating from different stellar populations) for 814 galaxies in the nearby Universe, all observed by the Herschel Space Observatory. We take advantage of the widely used galaxy SED fitting code CIGALE, properly adapted to include the state-of-the-art dust model THEMIS. Using the DustPedia photometry we determine the physical properties of the galaxies, such as, the dust and stellar mass, the star-formation rate, the bolometric luminosity as well as the unattenuated and the absorbed by dust stellar light, for both the old (> 200 Myr) and young (<= 200 Myr) stellar populations. We show how the mass of stars, dust, and atomic gas, as well as the star-formation rate and the dust temperature vary between galaxies of different morphologies and provide recipes to estimate these parameters given their Hubble stage (T). We find a mild correlation between the mass fraction of the small a-C(:H) grains with the specific star-formation rate. On average, young stars are very efficient in heating the dust, with absorption fractions reaching as high as ~77% of the total, unattenuated luminosity of this population. On the other hand, the maximum absorption fraction of old stars is ~24%. Dust heating in early-type galaxies is mainly due to old stars, up to a level of ~90%. Young stars progressively contribute more for `typical' spiral galaxies and they become the dominant source of dust heating for Sm type and irregular galaxies, donating up to ~60% of their luminosity to this purpose. Finally, we find a strong correlation of the dust heating fraction by young stars with morphology and the specific star-formation rate.

Galaxy formation depends on a complex interplay between gravitational collapse, gas accretion, merging, and feedback processes. Yet, after many decades of investigation, these concepts are poorly understood. This paper presents the argument that warm H$_2$ can be used as a tool to unlock some of these mysteries. Turbulence, shocks and outflows, driven by star formation, AGN activity or inflows, may prevent the rapid buildup of star formation in galaxies. Central to our understanding of how gas is converted into stars is the process by which gas can dissipate its mechanical energy through turbulence and shocks in order to cool. H$_2$ lines provide direct quantitative measurements of kinetic energy dissipation in molecular gas in galaxies throughout the Universe. Based on the detection of very powerful H$_2$ lines from z = 2 galaxies and proto-clusters at the detection limits of {\it Spitzer}, we are confident that future far-IR and UV H$_2$ observations will provide a wealth of new information and insight into galaxy evolution to high-z. Finally, at the very earliest epoch of star and galaxy formation, warm H$_2$ may also provide a unique glimpse of molecular gas collapse at 7 $<$ z $<$ 12 in massive dark matter (DM) halos on their way to forming the very first galaxies. Such measurements are beyond the reach of existing and planned observatories.

Most radiative transfer models assume that dust in spiral galaxies is distributed exponentially. In this paper our goal is to verify this assumption by analysing the two-dimensional large-scale distribution of dust in galaxies from the DustPedia sample. For this purpose, we make use of Herschel imaging in five bands, from 100 to 500{\mu}m, in which the cold dust constituent is primarily traced and makes up the bulk of the dust mass in spiral galaxies. For a subsample of 320 disc galaxies, we successfully perform a simultaneous fitting with a single S\'ersic model of the Herschel images in all five bands using the multiband modelling code GALFITM. We report that the S\'ersic index $n$, which characterises the shape of the S\'ersic profile, lies systematically below 1 in all Herschel bands and is almost constant with wavelength. The average value at 250{\mu}m is $0.67\pm0.37$ (187 galaxies are fitted with $n_{250}\leq0.75$, 87 galaxies have $0.75<n_{250}\leq1.25$, and 46 - with $n_{250}>1.25$). Most observed profiles exhibit a depletion in the inner region (at $r<0.3-0.4$ of the optical radius $r_{25}$ ) and are more or less exponential in the outer part. We also find breaks in the dust emission profiles at longer distances $(0.5-0.6)r_{25}$ which are associated with the breaks in the optical and near-infrared. We assume that the observed deficit of dust emission in the inner galaxy region is related to the depression in the radial profile of the HI surface density in the same region because the atomic gas reaches high enough surface densities there to be transformed into molecular gas. If a galaxy has a triggered star formation in the inner region (for example, because of a strong bar instability, which transfers the gas inwards to the centre, or a pseudobulge formation), no depletion or even an excess of dust emission in the centre is observed.

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 study the fraction of stellar radiation absorbed by dust, f_abs, in 814 galaxies of different morphological types. The targets constitute the vast majority (93%) of the DustPedia sample, including almost all large (optical diameter larger than 1'), nearby (v <= 3000 km/s) galaxies observed with the Herschel Space Observatory. For each object, we model the spectral energy distribution from the ultraviolet to the sub-millimetre using the dedicated, aperture-matched DustPedia photometry and the fitting code CIGALE. The value of f_abs is obtained from the total luminosity emitted by dust and from the bolometric luminosity, which are estimated by the fit. On average, 19% of the stellar radiation is absorbed by dust in DustPedia galaxies. The fraction rises to 25% if only late-type galaxies are considered. The dependence of f_abs on morphology, showing a peak for Sb-Sc galaxies, is weak; it reflects a stronger, yet broad, positive correlation with the bolometric luminosity, which is identified for late-type, disk-dominated, high-specific-star-formation rate, gas-rich objects. We find no variation of f_abs with inclination, at odds with radiative transfer models of edge-on galaxies. These results call for a self-consistent modelling of the evolution of the dust mass and geometry along the build-up of the stellar content. We also provide template spectral energy distributions in bins of morphology and luminosity and study the variation of f_abs with stellar mass and specific star formation rate. We confirm that the local Universe is missing the high f_abs}, luminous and actively star-forming objects necessary to explain the energy budget in observations of the extragalactic background light.

With the recent discovery of a dozen dusty star-forming galaxies and around 30 quasars at z>5 that are hyper-luminous in the infrared ($\mu$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $\mu$ is a lensing magnification factor), the possibility has opened up for SPICA, the proposed ESA M5 mid-/far-infrared mission, to extend its spectroscopic studies toward the epoch of reionization and beyond. In this paper, we examine the feasibility and scientific potential of such observations with SPICA's far-infrared spectrometer SAFARI, which will probe a spectral range (35-230 $\mu$m) that will be unexplored by ALMA and JWST. Our simulations show that SAFARI is capable of delivering good-quality spectra for hyper-luminous infrared galaxies (HyLIRGs) at z=5-10, allowing us to sample spectral features in the rest-frame mid-infrared and to investigate a host of key scientific issues, such as the relative importance of star formation versus AGN, the hardness of the radiation field, the level of chemical enrichment, and the properties of the molecular gas. From a broader perspective, SAFARI offers the potential to open up a new frontier in the study of the early Universe, providing access to uniquely powerful spectral features for probing first-generation objects, such as the key cooling lines of low-metallicity or metal-free forming galaxies (fine-structure and H2 lines) and emission features of solid compounds freshly synthesized by Population III supernovae. Ultimately, SAFARI's ability to explore the high-redshift Universe will be determined by the availability of sufficiently bright targets (whether intrinsically luminous or gravitationally lensed). With its launch expected around 2030, SPICA is ideally positioned to take full advantage of upcoming wide-field surveys such as LSST, SKA, Euclid, and WFIRST, which are likely to provide extraordinary targets for SAFARI.

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.

Measurements in the infrared wavelength domain allow us to assess directly the physical state and energy balance of cool matter in space, thus enabling the detailed study of the various processes that govern the formation and early evolution of stars and planetary systems in galaxies over cosmic time. Previous infrared missions, from IRAS to Herschel, have revealed a great deal about the obscured Universe, but sensitivity has been limited because up to now it has not been possible to fly a telescope that is both large and cold. SPICA is a mission concept aimed at taking the next step in mid- and far-infrared observational capability by combining a large and cold telescope with instruments employing state-of-the-art ultra-sensitive detectors. The mission concept foresees a 2.5-meter diameter telescope cooled to below 8 K. With cooling provided by mechanical coolers instead of depending on a limited cryogen supply, the mission lifetime can extend significantly beyond the required three years. SPICA offers instrumentation with spectral resolving powers ranging from R ~50 through 11000 in the 17-230 $\mu$m domain as well as R~28.000 spectroscopy between 12 and 18 $\mu$m. Additionally SPICA will provide efficient 30-37 $\mu$m broad band mapping, and polarimetric imaging in the 100-350 $\mu$m range. SPICA will provide unprecedented spectroscopic sensitivity of ~5 x $10^{-20}$ W/m$^2$ (5$\sigma$/1hr) - at least two orders of magnitude improvement over what has been attained to date. With this exceptional leap in performance, new domains in infrared astronomy will become accessible, allowing us, for example, to unravel definitively galaxy evolution and metal production over cosmic time, to study dust formation and evolution from very early epochs onwards, and to trace the formation history of planetary systems.

We present the molecular cloud properties of N55 in the Large Magellanic Cloud using $^{12}$CO(1-0) and $^{13}$CO(1-0) observations obtained with Atacama Large Millimeter Array. We have done a detailed study of molecular gas properties, to understand how the cloud properties of N55 differ from Galactic clouds. Most CO emission appears clumpy in N55, and molecular cores that have YSOs show larger linewidths and masses. The massive clumps are associated with high and intermediate mass YSOs. The clump masses are determined by local thermodynamic equilibrium and virial analysis of the $^{12}$CO and $^{13}$CO emissions. These mass estimates lead to the conclusion that, (a) the clumps are in self-gravitational virial equilibrium, and (b) the $^{12}$CO(1-0)-to-H$_2$ conversion factor, X$_{\rm CO}$, is 6.5$\times$10$^{20}$cm$^{-2}$(K km s$^{-1}$)$^{-1}$. This CO-to-H$_2$ conversion factor for N55 clumps is measured at a spatial scale of $\sim$0.67 pc, which is about two times higher than the X$_{\rm CO}$ value of Orion cloud at a similar spatial scale. The core mass function of N55 clearly show a turnover below 200M$_{\odot}$, separating the low-mass end from the high-mass end. The low-mass end of the $^{12}$CO mass spectrum is fitted with a power law of index 0.5$\pm$0.1, while for $^{13}$CO it is fitted with a power law index 0.6$\pm$0.2. In the high-mass end, the core mass spectrum is fitted with a power index of 2.0$\pm$0.3 for $^{12}$CO, and with 2.5$\pm$0.4 for $^{13}$CO. This power-law behavior of the core mass function in N55 is consistent with many Galactic clouds.

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.

IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Large Millimeter Array (ALMA) at longer wavelengths. SPICA, with its 2.5-m telescope actively-cooled to below 8K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei (AGN), sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. AGN and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. SPICA's large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z~6. Furthermore, SPICA spectroscopy will uncover the most luminous galaxies in the first few hundred million years of the Universe, through their characteristic dust and molecular hydrogen features.

The physical processes driving the chemical evolution of galaxies in the last $\sim 11\, \rm{Gyr}$ cannot be understood without directly probing the dust-obscured phase of star-forming galaxies and active galactic nuclei. This phase, hidden to optical tracers, represents the bulk of star formation and black hole accretion activity in galaxies at $1 < z < 3$. Spectroscopic observations with a cryogenic infrared (IR) observatory like SPICA will be sensitive enough to peer through the dust-obscured regions of galaxies and access the rest-frame mid- to far-IR range in galaxies at high-$z$. This wavelength range contains a unique suite of spectral lines and dust features that serve as proxies for the abundances of heavy elements and the dust composition, providing tracers with a feeble response to both extinction and temperature. In this work, we investigate how SPICA observations could be exploited to understand key aspects in the chemical evolution of galaxies: the assembly of nearby galaxies based on the spatial distribution of heavy element abundances, the global content of metals in galaxies reaching the knee of the luminosity function up to $z \sim 3$, and the dust composition of galaxies at high-$z$. Possible synergies with facilities available in the late 2020s are also discussed.

The DustPedia project is capitalising on the legacy of the Herschel Space Observatory, using cutting-edge modelling techniques to study dust in the 875 DustPedia galaxies - representing the vast majority of extended galaxies within 3000 km s$^{-1}$ that were observed by Herschel. This work requires a database of multiwavelength imagery and photometry that greatly exceeds the scope (in terms of wavelength coverage and number of galaxies) of any previous local-Universe survey. We constructed a database containing our own custom Herschel reductions, along with standardised archival observations from GALEX, SDSS, DSS, 2MASS, WISE, Spitzer, and Planck. Using these data, we performed consistent aperture-matched photometry, which we combined with external supplementary photometry from IRAS and Planck. We present our multiwavelength imagery and photometry across 42 UV-microwave bands for the 875 DustPedia galaxies. Our aperture-matched photometry, combined with the external supplementary photometry, represents a total of 21,857 photometric measurements. A typical DustPedia galaxy has multiwavelength photometry spanning 25 bands. We also present the Comprehensive & Adaptable Aperture Photometry Routine (CAAPR), the pipeline we developed to carry out our aperture-matched photometry. CAAPR is designed to produce consistent photometry for the enormous range of galaxy and observation types in our data. In particular, CAAPR is able to determine robust cross-compatible uncertainties, thanks to a novel method for reliably extrapolating the aperture noise for observations that cover a very limited amount of background. Our rich database of imagery and photometry is being made available to the community

The purpose of this work is the characterization of the radial distribution of dust, stars, gas, and star-formation rate (SFR) in a sub-sample of 18 face-on spiral galaxies extracted from the DustPedia sample. This study is performed by exploiting the multi-wavelength, from UV to sub-mm bands, DustPedia database, in addition to molecular (12CO) and atomic (HI) gas maps and metallicity abundance information available in the literature. We fitted the surface brightness profiles of the tracers of dust and stars, the mass surface density profiles of dust, stars, molecular gas, and total gas, and the SFR surface density profiles with an exponential curve and derived their scale-lengths. We also developed a method to solve for the CO-to-H2 conversion factor (alpha_CO) per galaxy by using dust and gas mass profiles. Although each galaxy has its own peculiar behaviour, we identified a common trend of the exponential scale-lengths vs. wavelength. On average, the scale-lengths normalized to the B-band 25 mag/arcsec^2 radius decrease from UV to 70 micron, from 0.4 to 0.2, and then increase back up to 0.3 at 500 microns. The main result is that, on average, the dust mass surface density scale-length is about 1.8 times the stellar one derived from IRAC data and the 3.6 micron surface brightness, and close to that in the UV. We found a mild dependence of the scale-lengths on the Hubble stage T: the scale-lengths of the Herschel bands and the 3.6 micron scale-length tend to increase from earlier to later types, the scale-length at 70 micron tends to be smaller than that at longer sub-mm wavelength with ratios between longer sub-mm wavelengths and 70 micron that decrease with increasing T. The scale-length ratio of SFR and stars shows a weak increasing trend towards later types.

We use the Spitzer SAGE survey of the Magellanic Clouds to evaluate the relationship between the 8-mic PAH emission, 24-mic hot dust emission, and HII region radiative transfer. We confirm that in the higher-metallicity Large Magellanic Cloud, PAH destruction is sensitive to optically thin conditions in the nebular Lyman continuum: objects identified as optically thin candidates based on nebular ionization structure show 6 times lower median 8-mic surface brightness (0.18 mJy arcsec^-2) than their optically thick counterparts (1.2 mJy arcsec^-2). The 24-mic surface brightness also shows a factor of 3 offset between the two classes of objects (0.13 vs 0.44 mJy arcsec^-2, respectively), which is driven by the association between the very small dust grains and higher density gas found at higher nebular optical depths. In contrast, PAH and dust formation in the low-metallicity Small Magellanic Cloud is strongly inhibited such that we find no variation in either 8-mic or 24-mic emission between our optically thick and thin samples. This is attributable to extremely low PAH and dust production together with high, corrosive UV photon fluxes in this low-metallicity environment. The dust mass surface densities and gas-to-dust ratios determined from dust maps using Herschel HERITAGE survey data support this interpretation.

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 report on electrostatic measurements made on board the European Space Agency mission LISA Pathfinder. Detailed measurements of the charge-induced electrostatic forces exerted on free-falling test masses (TMs) inside the capacitive gravitational reference sensor are the first made in a relevant environment for a space-based gravitational wave detector. Employing a combination of charge control and electric-field compensation, we show that the level of charge-induced acceleration noise on a single TM can be maintained at a level close to 1.0 fm/s^2/sqrt(Hz) across the 0.1-100 mHz frequency band that is crucial to an observatory such as LISA. Using dedicated measurements that detect these effects in the differential acceleration between the two test masses, we resolve the stochastic nature of the TM charge build up due to interplanetary cosmic rays and the TM charge-to-force coupling through stray electric fields in the sensor. All our measurements are in good agreement with predictions based on a relatively simple electrostatic model of the LISA Pathfinder instrument.

Understanding exoplanet formation and finding potentially habitable exoplanets is vital to an enhanced understanding of the universe. The use of nulling interferometry to strongly attenuate the central starlight provides the opportunity to see objects closer to the star than ever before. Given that exoplanets are usually warm, the 4 microns Mid-Infrared region is advantageous for such observations. The key performance parameters for a nulling interferometer are the extinction ratio it can attain and how well that is maintained across the operational bandwidth. Both parameters depend on the design and fabrication accuracy of the subcomponents and their wavelength dependence. Via detailed simulation it is shown in this paper that a planar chalcogenide photonic chip, consisting of three highly fabrication tolerant multimode interference couplers, can exceed an extinction ratio of 60 dB in double nulling operation and up to 40 dB for a single nulling operation across a wavelength window of 3.9 to 4.2 microns. This provides a beam combiner with sufficient performance, in theory, to image exoplanets.

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.

The radiation of stars heats dust grains in the diffuse interstellar medium and in star-forming regions in galaxies. Modelling this interaction provides information on dust in galaxies, a vital ingredient for their evolution. It is not straightforward to identify the stellar populations heating the dust, and to link attenuation to emission on a sub-galactic scale. Radiative transfer models are able to simulate this dust-starlight interaction in a realistic, three-dimensional setting. We investigate the dust heating mechanisms on a local and global galactic scale, using the Andromeda galaxy (M31) as our laboratory. We perform a series of panchromatic radiative transfer simulations of Andromeda with our code SKIRT. The high inclination angle of M31 complicates the 3D modelling and causes projection effects. However, the observed morphology and flux density are reproduced fairly well from UV to sub-millimeter wavelengths. Our model reveals a realistic attenuation curve, compatible with previous, observational estimates. We find that the dust in M31 is mainly (91 % of the absorbed luminosity) heated by the evolved stellar populations. The bright bulge produces a strong radiation field and induces non-local heating up to the main star-forming ring at 10 kpc. The relative contribution of unevolved stellar populations to the dust heating varies strongly with wavelength and with galactocentric distance.The dust heating fraction of unevolved stellar populations correlates strongly with NUV-r colour and specific star formation rate. These two related parameters are promising probes for the dust heating sources at a local scale.

The European Space Agency has invested heavily in two cornerstones missions; Herschel and Planck. The legacy data from these missions provides us with an unprecedented opportunity to study cosmic dust in galaxies so that we can answer fundamental questions about, for example: the origin of the chemical elements, physical processes in the interstellar medium (ISM), its effect on stellar radiation, its relation to star formation and how this relates to the cosmic far infrared background. In this paper we describe the DustPedia project, which is enabling us to develop tools and computer models that will help us relate observed cosmic dust emission to its physical properties (chemical composition, size distribution, temperature), to its origins (evolved stars, super novae, growth in the ISM) and the processes that destroy it (high energy collisions and shock heated gas). To carry out this research we will combine the Herschel/Planck data with that from other sources of data, providing observations at numerous wavelengths (< 41) across the spectral energy distribution, thus creating the DustPedia database. To maximise our spatial resolution and sensitivity to cosmic dust we limit our analysis to 4231 local galaxies (v < 3000 km/s) selected via their near infrared luminosity (stellar mass). To help us interpret the data we have developed a new physical model for dust (THEMIS), a new Bayesian method of fitting and interpreting spectral energy distributions (HerBIE) and a state-of-the-art Monte Carlo photon tracing radiative transfer model (SKIRT). In this the first of the DustPedia papers we describe the project objectives, data sets used and provide an insight into the new scientific methods we plan to implement.

We study the correlation between far-infared/submm dust emission and atomic gas column density in order to derive the properties of the high Galactic latitude, low density, Milky Way cirrus in the foreground of the Virgo cluster of galaxies. Dust emission maps from 60 to 850 um are obtained from SPIRE observations carried out within the Herschel Virgo Cluster Survey, complemented by IRAS-IRIS and Planck-HFI maps. Data from the Arecibo legacy Fast ALFA Survey is used to derive atomic gas column densities for two broad velocity components, low and intermediate velocity clouds. Dust emissivities are derived for each gas component and each far-infared/submm band. For the low velocity clouds, we measure an average emissivity 0.79 +/- 0.08 times 1E-20 MJy sr^-1 cm^2 at 250 um. After fitting a modified blackbody to the available bands, we estimated a dust absorption cross-section 0.49 +/- 0.13 times 1E-25 cm^2 H^-1 at 250 um (with dust temperature T = 20.4 +/- 1.5 K and spectral index beta = 1.53 +/- 0.17). The results are in excellent agreement with those obtained by Planck over a much larger coverage of the high Galactic latitude cirrus (50% of the sky vs 0.2% in our work). For dust associated with intermediate velocity gas, we confirm earlier Planck results and find a higher temperature and lower emissivity and cross-section. After subtracting the modelled components, we find regions at scales smaller than 20' where the residuals deviate significantly from the average, cosmic-infrared-background dominated, scatter. These large residuals are most likely due to local variations in the cirrus dust properties (and/or the dust/atomic-gas correlation) or to high-latitude molecular clouds with average N_H2 <~ 1E20 cm^-2. We find no conclusive evidence for intracluster dust emission in Virgo.

Photonic integrated circuits are established as the technique of choice for a number of astronomical processing functions due to their compactness, high level of integration, low losses, and stability. Temperature control, mechanical vibration and acoustic noise become controllable for such a device enabling much more complex processing than can realistically be considered with bulk optics. To date the benefits have mainly been at wavelengths around 1550 nm but in the important Mid-Infrared region, standard photonic chips absorb light strongly. Chalcogenide glasses are well known for their transparency to beyond 10000 nm, and the first results from coupler devices intended for use in an interferometric nuller for exoplanetary observation in the Mid-Infrared L band (3800-4200 nm) are presented here showing that suitable performance can be obtained both theoretically and experimentally for the first fabricated devices operating at 4000 nm.

We present a line survey of the ultra-luminous infrared galaxy Arp 220, taken with the newly installed SEPIA Band 5 instrument on APEX. We illustrate the capacity of SEPIA to detect the 183.3 GHz H2O 31,3-22,0 line against the atmospheric H2O absorption feature. We confirm the previous detection of the HCN(2-1) line, and detect new transitions of standard dense gas tracers such as HNC(2-1), HCO+(2-1), CS(4-3), C34S(4-3), HC3N(20-19). We also detect HCN(2-1) v2=1 and the 193.5 GHz methanol (4-3) group for the first time. The absence of time variations in the megamaser water line compared to previous observations seems to rule out an AGN nuclear origin for the line. It could, on the contrary, favor a thermal origin instead, but also possibly be a sign that the megamaser emission is associated with star-forming cores washed-out in the beam. We finally discuss how the new transitions of HCN, HNC, HCO+ refine our knowledge of the ISM physical conditions in Arp 220.

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

Several dedicated surveys focusing on early-type galaxies (ETGs) reveal that significant fractions of them are detectable in all interstellar medium phases studied to date. We select ETGs from the Herschel Reference Survey that have both far-infrared Herschel and either HI or CO detection (or both). We derive their star formation rates (SFR), stellar masses and dust masses via modelling their spectral energy distributions. We combine these with literature information on their atomic and molecular gas properties, in order to relate their star formation, total gas mass and dust mass on global scales. The ETGs deviate from the dust mass-SFR relation and the Schmidt-Kennicutt relation that SDSS star forming galaxies define: compared to SDSS galaxies, ETGs have more dust at the same SFR, or less SFR at the same dust mass. When placing them in the M*-SFR plane, ETGs show a much lower specific SFR as compared to normal star-forming galaxies. ETGs show a large scatter compared to the Schmidt-Kennicutt relation found locally within our Galaxy, extending to lower SFRs and gas mass surface densities. Using an ETG's SFR and the Schmidt-Kennicutt law to predict its gas mass leads to an underestimate. ETGs have similar observed-gas-to-modelled-dust mass ratios to star forming-galaxies of the same stellar mass, as well as they exhibit a similar scatter.

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 the ALMA Band 3 and Band 6 results of 12CO(2-1), 13$CO(2-1), H30alpha recombination line, free-free emission around 98 GHz, and the dust thermal emission around 230 GHz toward the N159 East Giant Molecular Cloud (N159E) in the Large Magellanic Cloud (LMC). LMC is the nearest active high-mass star forming face-on galaxy at a distance of 50 kpc and is the best target for studing high-mass star formation. ALMA observations show that N159E is the complex of filamentary clouds with the width and length of ~1 pc and 5 pc - 10 pc, respectively. The total molecular mass is 0.92 x 10^5 Msun from the 13CO(2-1) intensity. N159E harbors the well-known Papillon Nebula, a compact high-excitation HII region. We found that a YSO associated with the Papillon Nebula has the mass of 35 Msun and is located at the intersection of three filamentary clouds. It indicates that the formation of the high-mass YSO was induced by the collision of filamentary clouds. Fukui et al. 2015 reported a similar kinematic structure toward a YSO in the N159 West region which is another YSO that has the mass larger than 35 Msun in these two regions. This suggests that the collision of filamentary clouds is a primary mechanism of high-mass star formation. We found a small molecular hole around the YSO in Papillon Nebula with sub-pc scale. It is filled by free-free and H30alpha emission. Temperature of the molecular gas around the hole reaches ~ 80 K. It indicates that this YSO has just started the distruction of parental molecular cloud.

We present observations of the dense molecular gas tracers HCN, HNC, and HCO+ in the J=1-0 transition using ALMA. We supplement our datasets with previous observations of CO J=1-0, which traces the total molecular gas content. We separate the Antennae into 7 bright regions in which we detect emission from all three molecules, including the nuclei of NGC 4038 and NGC 4039, 5 super giant molecular complexes in the overlap region and 2 additional bright clouds. We find that the ratio of L(HCN)/L(CO), which traces the dense molecular gas fraction, is greater in the two nuclei (L(HCN)/L(CO) ~ 0.07 - 0.08) than in the overlap region (L(HCN)/L(CO) <0.05). We attribute this to an increase in pressure due to the stellar potential within the nuclei, similar to what has been seen previously in the Milky Way and nearby spiral galaxies. Furthermore, the ratio of L(HNC)/L(HCN) ~ 0.3-0.4 does not vary by more than a factor of 1.5 between regions. By comparing our measured ratios to PDR models including mechanical heating, we find that the ratio of L(HNC)/L(HCN) is consistent with mechanical heating contributing >= 5% - 10% of the PDR surface heating to the total heating budget. Finally, the ratio of L(HCN)/L(HCO+) varies from ~1 in the nucleus of NGC 4038 down to ~0.5 in the overlap region. The lower ratio in the overlap region may be due to an increase in the cosmic ray rate from the increased supernova rate within this region.

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)

We present $^{12}$CO(1-0) and $^{12}$CO(2-1) observations of a sample of 20 star-forming dwarfs selected from the Herschel Virgo Cluster Survey, with oxygen abundances ranging from 12 + log(O/H) ~ 8.1 to 8.8. CO emission is observed in ten galaxies and marginally detected in another one. CO fluxes correlate with the FIR 250 $\mu$m emission, and the dwarfs follow the same linear relation that holds for more massive spiral galaxies extended to a wider dynamical range. We compare different methods to estimate H2 molecular masses, namely a metallicity-dependent CO-to-H2 conversion factor and one dependent on H-band luminosity. The molecular-to-stellar mass ratio remains nearly constant at stellar masses <~ 10$^9$ M$_{\odot}$, contrary to the atomic hydrogen fraction, M$_{HI}$/M$_*$, which increases inversely with M$_*$. The flattening of the M$_{H_2}$/M$_*$ ratio at low stellar masses does not seem to be related to the effects of the cluster environment because it occurs for both HI-deficient and HI-normal dwarfs. The molecular-to-atomic ratio is more tightly correlated with stellar surface density than metallicity, confirming that the interstellar gas pressure plays a key role in determining the balance between the two gaseous components of the interstellar medium. Virgo dwarfs follow the same linear trend between molecular gas mass and star formation rate as more massive spirals, but gas depletion timescales, $\tau_{dep}$, are not constant and range between 100 Myr and 6 Gyr. The interaction with the Virgo cluster environment is removing the atomic gas and dust components of the dwarfs, but the molecular gas appears to be less affected at the current stage of evolution within the cluster. However, the correlation between HI deficiency and the molecular gas depletion time suggests that the lack of gas replenishment from the outer regions of the disc is lowering the star formation activity.

The dust in nearby galaxies absorbs a fraction of the UV-optical-near-infrared radiation produced by stars. This energy is consequently re-emitted in the infrared. We investigate the portion of the stellar radiation absorbed by spiral galaxies from the HRS by modelling their UV-to-submillimetre spectral energy distributions. Our models provide an attenuated and intrinsic SED from which we find that on average 32 % of all starlight is absorbed by dust. We define the UV heating fraction as the percentage of dust luminosity that comes from absorbed UV photons and find that this is 56 %, on average. This percentage varies with morphological type, with later types having significantly higher UV heating fractions. We find a strong correlation between the UV heating fraction and specific star formation rate and provide a power-law fit. Our models allow us to revisit the IRX-AFUV relations, and derive these quantities directly within a self-consistent framework. We calibrate this relation for different bins of NUV-r colour and provide simple relations to relate these parameters. We investigated the robustness of our method and we conclude that the derived parameters are reliable within the uncertainties which are inherent to the adopted SED model. This calls for a deeper investigation on how well extinction and attenuation can be determined through panchromatic SED modelling.

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.

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

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.

We have carried out 13CO(J=2-1) observations of the active star-forming region N159 West in the LMC with ALMA. We have found that the CO distribution at a sub-pc scale is highly elongated with a small width. These elongated clouds called "filaments" show straight or curved distributions with a typical width of 0.5-1.0pc and a length of 5-10pc. All the known infrared YSOs are located toward the filaments. We have found broad CO wings of two molecular outflows toward young high-mass stars in N159W-N and N159W-S, whose dynamical timescale is ~10^4 yrs. This is the first discovery of protostellar outflow in external galaxies. For N159W-S which is located toward an intersection of two filaments we set up a hypothesis that the two filaments collided with each other ~10^5 yrs ago and triggered formation of the high-mass star having ~37 Mo. The colliding clouds show significant enhancement in linewidth in the intersection, suggesting excitation of turbulence in the shocked interface layer between them as is consistent with the magneto-hydro-dynamical numerical simulations (Inoue & Fukui 2013). This turbulence increases the mass accretion rate to ~6x10^-4 Mo yr^-1, which is required to overcome the stellar feedback to form the high-mass star.

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.

The LISA Pathfinder mission will demonstrate the technology of drag-free test masses for use as inertial references in future space-based gravitational wave detectors. To accomplish this, the Pathfinder spacecraft will perform drag-free flight about a test mass while measuring the acceleration of this primary test mass relative to a second reference test mass. Because the reference test mass is contained within the same spacecraft, it is necessary to apply forces on it to maintain its position and attitude relative to the spacecraft. These forces are a potential source of acceleration noise in the LISA Pathfinder system that are not present in the full LISA configuration. While LISA Pathfinder has been designed to meet it's primary mission requirements in the presence of this noise, recent estimates suggest that the on-orbit performance may be limited by this `suspension noise'. The drift-mode or free-flight experiments provide an opportunity to mitigate this noise source and further characterize the underlying disturbances that are of interest to the designers of LISA-like instruments. This article provides a high-level overview of these experiments and the methods under development to analyze the resulting data.

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

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.

To assess the effects of the cluster environment on the different components of the interstellar medium, we analyse the FIR-submm properties of a sample of star-forming dwarf (SFD) galaxies detected by the Herschel Virgo Cluster Survey (HeViCS). We determine dust masses and dust temperatures by fitting a modified black body (MBB) function to the spectral energy distributions (SEDs). Stellar and gas masses, star formation rates (SFRs), and metallicities are obtained from the analysis of a set of ancillary data. Dust is detected in 49 out of 140 optically identified dwarfs covered by the HeViCS field; considering only dwarfs brighter than $m_B$ = 18 mag, this gives a detection rate of 43%. After evaluating different emissivity indices, we find that the FIR-submm SEDs are best-fit by $\beta$=1.5, with a median dust temperature $T_d$ = 22.4 K. Assuming $\beta$=1.5, 67% of the 23 galaxies detected in all five Herschel bands show emission at 500 $\mu$m in excess of the MBB model. The excess is inversely correlated with SFR and stellar masses. To study the variations in the global properties of our sample due to environmental effects, we compare the Virgo SFDs to other Herschel surveys, such as KINGFISH, the Dwarf Galaxy Survey (DGS), and the HeViCS bright galaxy catalogue (BGC). We explore the relations between stellar mass and HI fraction, specific SFR, dust fraction, gas-to-dust ratio over a wide range of stellar masses. Highly HI-deficient Virgo dwarf galaxies are mostly characterised by quenched star formation activity and lower dust fractions giving hints for dust stripping in cluster dwarfs. However, we find that the fraction of dust removed has to be less than that of the HI component. Since the Virgo SFDs are likely to be crossing the cluster for the first time, a longer timescale might be necessary to strip the more centrally concentrated dust distribution.

The spatial variations of the gas-to-dust ratio (GDR) provide constraints on the chemical evolution and lifecycle of dust in galaxies. We examine the relation between dust and gas at 10-50 pc resolution in the Large and Small Magellanic Clouds (LMC and SMC) based on Herschel far-infrared (FIR), H I 21 cm, CO, and Halpha observations. In the diffuse atomic ISM, we derive the gas-to-dust ratio as the slope of the dust-gas relation and find gas-to-dust ratios of 380+250-130 in the LMC, and 1200+1600-420 in the SMC, not including helium. The atomic-to-molecular transition is located at dust surface densities of 0.05 Mo pc-2 in the LMC and 0.03 Mo pc-2 in the SMC, corresponding to AV ~ 0.4 and 0.2, respectively. We investigate the range of CO-to-H2 conversion factor to best account for all the molecular gas in the beam of the observations, and find upper limits on XCO to be 6x1020 cm-2 K-1 km-1 s in the LMC (Z=0.5Zo) at 15 pc resolution, and 4x 1021 cm-2 K-1 km-1 s in the SMC (Z=0.2Zo) at 45 pc resolution. In the LMC, the slope of the dust-gas relation in the dense ISM is lower than in the diffuse ISM by a factor ~2, even after accounting for the effects of CO-dark H2 in the translucent envelopes of molecular clouds. Coagulation of dust grains and the subsequent dust emissivity increase in molecular clouds, and/or accretion of gas-phase metals onto dust grains, and the subsequent dust abundance (dust-to-gas ratio) increase in molecular clouds could explain the observations. In the SMC, variations in the dust-gas slope caused by coagulation or accretion are degenerate with the effects of CO-dark H2. Within the expected 5--20 times Galactic XCO range, the dust-gas slope can be either constant or decrease by a factor of several across ISM phases. Further modeling and observations are required to break the degeneracy between dust grain coagulation, accretion, and CO-dark H2.

We present the detection and analysis of molecular hydrogen emission toward ten interstellar regions in the Large Magellanic Cloud. We examined low-resolution infrared spectral maps of twelve regions obtained with the Spitzer infrared spectrograph (IRS). The pure rotational 0--0 transitions of H$_2$ at 28.2 and 17.1${\,\rm \mu m}$ are detected in the IRS spectra for ten regions. The higher level transitions are mostly upper limit measurements except for three regions, where a 3$\sigma$ detection threshold is achieved for lines at 12.2 and 8.6${\,\rm \mu m}$. The excitation diagrams of the detected H$_2$ transitions are used to determine the warm H$_2$ gas column density and temperature. The single-temperature fits through the lower transition lines give temperatures in the range $86-137\,{\rm K}$. The bulk of the excited H$_2$ gas is found at these temperatures and contributes $\sim$5-17% to the total gas mass. We find a tight correlation of the H$_2$ surface brightness with polycyclic aromatic hydrocarbon and total infrared emission, which is a clear indication of photo-electric heating in photodissociation regions. We find the excitation of H$_2$ by this process is equally efficient in both atomic and molecular dominated regions. We also present the correlation of the warm H$_2$ physical conditions with dust properties. The warm H$_2$ mass fraction and excitation temperature show positive correlations with the average starlight intensity, again supporting H$_2$ excitation in photodissociation regions.

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.

With the newly available SPIRE images at 250 and 500 micron from Herschel Space Observatory, we study quantitative correlations over a sub-kpc scale among three distinct emission components in the interstellar medium of the nearby spiral galaxy M 81 (NGC 3031): (a) $I_{8}$ or $I_{24}$, the surface brightness of the mid-infrared emission observed in the Spitzer IRAC 8 or MIPS 24 micron band, with $I_8$ and $I_{24}$ being dominated by the emissions from Polycyclic Aromatic Hydrocarbons (PAHs) and very small grains (VSGs) of dust, respectively; (b) $I_{500}$, that of the cold dust continuum emission in the Herschel SPIRE 500 micron band, dominated by the emission from large dust grains heated by evolved stars, and (c) $I_{{\rm H}\alpha}$, a nominal surface brightness of the H$\alpha$ line emission, from gas ionized by newly formed massive stars. The results from our correlation study, free from any assumption on or modeling of dust emissivity law or dust temperatures, present solid evidence for significant heating of PAHs and VSGs by evolved stars. In the case of M 81, about 67% (48%) of the 8 micron (24 micron) emission derives its heating from evolved stars, with the remainder attributed to radiation heating associated with ionizing stars.

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)

NGC1569 has some of the most vigorous star formation among nearby galaxies. It hosts two super star clusters (SSCs) and has a higher star formation rate (SFR) per unit area than other starburst dwarf galaxies. Extended emission beyond the galaxy's optical body is observed in warm and hot ionised and atomic hydrogen gas; a cavity surrounds the SSCs. We aim to understand the impact of the massive star formation on the surrounding interstellar medium in NGC1569 through a study of its stellar and dust properties. We use Herschel and ancillary multiwavelength observations, from the ultraviolet to the submillimeter regime, to construct its spectral energy distribution, which we model with magphys on ~300pc scales at the SPIRE250 {\mu}m resolution. The multiwavelength morphology shows low levels of dust emission in the cavity, and a concentration of several dust knots in its periphery. The extended emission seen in the ionised and neutral hydrogen observations is also present in the far-infrared emission. The dust mass is higher in the periphery of the cavity, driven by ongoing star formation and dust emission knots. The SFR is highest in the central region, while the specific SFR is more sensitive to the ongoing star formation. The region encompassing the cavity and SSCs contains only 12 per cent of the dust mass of the central starburst, in accord with other tracers of the interstellar medium. The gas-to-dust mass ratio is lower in the cavity and fluctuates to higher values in its periphery.

In this paper we present a catalogue of Giant Molecular Clouds (GMCs) in the Andromeda (M31) galaxy extracted from the Hershel Exploitation of Local Galaxy Andromeda (HELGA) dataset. GMCs are identified from the Herschel maps using a hierarchical source extraction algorithm. We present the results of this new catalogue and characterise the spatial distribution and spectral energy properties of its clouds based on the radial dust/gas properties found by Smith et al (2012). 236 GMCs in the mass range 10^4-10^7 M_sol are identified, their cumulative mass distribution is found to be proportional to M^-1.45 in agreement with earlier studies. The GMCs appear to follow the same cloud mass to L_CO correlation observed in the Milky Way. However, comparison between this catalogue and interferometry studies also shows that the GMCs are substructured below the Herschel resolution limit suggesting that we are observing associations of GMCs. Following Gordon et al. (2006), we study the spatial structure of M31 by splitting the observed structure into a set of spiral arms and offset rings. We fit radii of 10.5 and 15.5 kpc to the two most prominent rings. We then fit a logarithmic spiral with a pitch angle of 8.9 deg to the GMCs not associated with either ring. Lastly, we comment upon the effects of deprojection on our results and investigate the effect different models for M31's inclination will have upon the projection of an unperturbed spiral arm system.

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

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.

The Magellanic Bridge is the nearest low-metallicity, tidally stripped environment, offering a unique high-resolution view of physical conditions in merging and forming galaxies. In this paper we present analysis of candidate massive young stellar objects (YSOs), i.e., {\it in situ, current} massive star formation (MSF) in the Bridge using {\it Spitzer} mid-IR and complementary optical and near-IR photometry. While we definitely find YSOs in the Bridge, the most massive are $\sim10 M_\odot$, $\ll45 M_\odot$ found in the Large Magellanic Cloud (LMC). The intensity of MSF in the Bridge also appears decreasing, as the most massive YSOs are less massive than those formed in the past. To investigate environmental effects on MSF, we have compared properties of massive YSOs in the Bridge to those in the LMC. First, YSOs in the Bridge are apparently less embedded than in the LMC: 81% of Bridge YSOs show optical counterparts, compared to only 56% of LMC sources with the same range of mass, circumstellar dust mass, and line-of-sight extinction. Circumstellar envelopes are evidently more porous or clumpy in the Bridge's low-metallicity environment. Second, we have used whole samples of YSOs in the LMC and the Bridge to estimate the probability of finding YSOs at a given \hi\ column density, N(HI). We found that the LMC has $\sim3\times$ higher probability than the Bridge for N(HI) $>10\times10^{20}$ cm$^{-2}$, but the trend reverses at lower N(HI). Investigating whether this lower efficiency relative to HI is due to less efficient molecular cloud formation, or less efficient cloud collapse, or both, will require sensitive molecular gas observations.

We investigate the connection between dust and gas in the nearby edge-on spiral galaxy NGC 891. High resolution Herschel PACS and SPIRE 70, 100, 160, 250, 350, and 500 $\mu$m images are combined with JCMT SCUBA 850 $\mu$m observations to trace the far-infrared/submillimetre spectral energy distribution (SED). Maps of the HI 21 cm line and CO(J=3-2) emission trace the atomic and molecular hydrogen gas, respectively. We fit one-component modified blackbody models to the integrated SED, finding a global dust mass of 8.5$\times$10$^{7}$ M$_{\odot}$ and an average temperature of 23$\pm$2 K. We also fit the pixel-by-pixel SEDs to produce maps of the dust mass and temperature. The dust mass distribution correlates with the total stellar population as traced by the 3.6 $\mu$m emission. The derived dust temperature, which ranges from approximately 17 to 24 K, is found to correlate with the 24 $\mu$m emission. Allowing the dust emissivity index to vary, we find an average value of $\beta$ = 1.9$\pm$0.3. We confirm an inverse relation between the dust emissivity spectral index and dust temperature, but do not observe any variation of this relationship with vertical height from the mid-plane of the disk. A comparison of the dust properties with the gaseous components of the ISM reveals strong spatial correlations between the surface mass densities of dust and the molecular hydrogen and total gas surface densities. Observed asymmetries in the dust temperature, and the H$_{2}$-to-dust and total gas-to-dust ratios hint that an enhancement in the star formation rate may be the result of larger quantities of molecular gas available to fuel star formation in the NE compared to the SW. Whilst the asymmetry likely arises from dust obscuration due to the geometry of the line-of-sight projection of the spiral arms, we cannot exclude an enhancement in the star formation rate in the NE side of the disk.

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 gather infrared (IR) photometric data from 8 to 500 microns (Spitzer, WISE, IRAS and Herschel) for all of the HRS galaxies. Draine & Li (2007) models are fit to the data from which the stellar contribution has been carefully removed. We find that our photometric coverage is sufficient to constrain all of the models parameters and that a strong constraint on the 20-60 microns range is mandatory to estimate the relative contribution of the photo-dissociation regions to the IR SED. The SED models tend to systematically under-estimate the observed 500 microns flux densities, especially for low mass systems. We provide the output parameters for all of the galaxies: the minimum intensity of the interstellar radiation field (ISRF), the fraction of PAH, the relative contribution of PDR and evolved stellar population to the dust heating, the $M_{dust}$ and the $L_{IR}$. For a subsample of gas-rich galaxies, we analyze the relations between these parameters and the integrated properties of galaxies, such as $M_*$, SFR, metallicity, H$\alpha$ and H-band surface brightness, and the FUV attenuation. A good correlation between the fraction of PAH and the metallicity is found implying a weakening of the PAH emission in galaxies with low metallicities. The intensity of the IRSF and the H-band and H$\alpha$ surface brightnesses are correlated, suggesting that the diffuse dust component is heated by both the young stars in star forming regions and the diffuse evolved population. We use these results to provide a new set of IR templates calibrated with Herschel observations on nearby galaxies and a mean SED template to provide the z=0 reference for cosmological studies. For the same purpose, we put our sample on the SFR-$M_*$ diagram. The templates are compared to the most popular IR SED libraries, enlightening a large discrepancy between all of them in the 20-100 microns range.

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 SPIRE-FTS observations of the Antennae (NGC 4038/39), a well studied, nearby ($22$ Mpc) ongoing merger between two gas rich spiral galaxies. We detect 5 CO transitions ($J=4-3$ to $J=8-7$), both [CI] transitions and the [NII]$205\mu m$ transition across the entire system, which we supplement with ground based observations of the CO $J=1-0$, $J=2-1$ and $J=3-2$ transitions, and Herschel PACS observations of [CII] and [OI]$63\mu m$. Using the CO and [CI] transitions, we perform both a LTE analysis of [CI], and a non-LTE radiative transfer analysis of CO and [CI] using the radiative transfer code RADEX along with a Bayesian likelihood analysis. We find that there are two components to the molecular gas: a cold ($T_{kin}\sim 10-30$ K) and a warm ($T_{kin} \gtrsim 100$ K) component. By comparing the warm gas mass to previously observed values, we determine a CO abundance in the warm gas of $x_{CO} \sim 5\times 10^{-5}$. If the CO abundance is the same in the warm and cold gas phases, this abundance corresponds to a CO $J=1-0$ luminosity-to-mass conversion factor of $\alpha_{CO} \sim 7 \ M_{\odot}{pc^{-2} \ (K \ km \ s^{-1})^{-1}}$ in the cold component, similar to the value for normal spiral galaxies. We estimate the cooling from H$_2$, [CII], CO and [OI]$63\mu m$ to be $\sim 0.01 L_{\odot}/M_{\odot}$. We compare PDR models to the ratio of the flux of various CO transitions, along with the ratio of the CO flux to the far-infrared flux in NGC 4038, NGC 4039 and the overlap region. We find that the densities recovered from our non-LTE analysis are consistent with a background far-ultraviolet field of strength $G_0\sim 1000$. Finally, we find that a combination of turbulent heating, due to the ongoing merger, and supernova and stellar winds are sufficient to heat the molecular gas.

Herschel FIR observations are used to construct Virgo cluster galaxy luminosity functions and to show that the cluster lacks the very bright and the numerous faint sources detected in field galaxy surveys. The far-infrared SEDs are fitted to obtain dust masses and temperatures and the dust mass function. The cluster is over dense in dust by about a factor of 100 compared to the field. The same emissivity (beta) temperature relation applies for different galaxies as that found for different regions of M31. We use optical and HI data to show that Virgo is over dense in stars and atomic gas by about a factor of 100 and 20 respectively. Metallicity values are used to measure the mass of metals in the gas phase. The mean metallicity is about 0.7 solar and 50% of the metals are in the dust. For the cluster as a whole the mass density of stars in galaxies is 8 times that of the gas and the gas mass density is 130 times that of the metals. We use our data to consider the chemical evolution of the individual galaxies, inferring that the measured variations in effective yield are due to galaxies having different ages, being affected to varying degrees by gas loss. Four galaxy scaling relations are considered: mass-metallicity, mass-velocity, mass-star formation rate and mass-radius - we suggest that initial galaxy mass is the prime driver of a galaxy's ultimate destiny. Finally, we use X-ray observations and galaxy dynamics to assess the dark and baryonic matter content compared to the cosmological model.

The old, red stars which constitute the bulges of galaxies, and the massive black holes at their centres, are the relics of a period in cosmic history when galaxies formed stars at remarkable rates and active galactic nuclei (AGN) shone brightly from accretion onto black holes. It is widely suspected, but unproven, that the tight correlation in mass of the black hole and stellar components results from the AGN quenching the surrounding star formation as it approaches its peak luminosity. X-rays trace emission from AGN unambiguously, while powerful star-forming galaxies are usually dust-obscured and are brightest at infrared to submillimetre wavelengths. Here we report observations in the submillimetre and X-ray which show that rapid star formation was common in the host galaxies of AGN when the Universe was 2-6 Gyrs old, but that the most vigorous star formation is not observed around black holes above an X-ray luminosity of 10^44 erg/s. This suppression of star formation in the host galaxies of powerful AGN is a key prediction of models in which the AGN drives a powerful outflow, expelling the interstellar medium of its host galaxy and transforming the galaxy's properties in a brief period of cosmic time.

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

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 use dust scaling relations to investigate the hypothesis that Virgo cluster transition-type dwarfs are infalling star-forming field galaxies, which is argued based on their optical features (e.g. disks, spiral arms, bars) and kinematic properties similar to late-type galaxies. After their infall, environmental effects gradually transform them into early-type galaxies through the removal of their interstellar medium and quenching of all star formation activity. In this paper, we aim to verify whether this hypothesis holds using far-infrared diagnostics based on Herschel observations of the Virgo cluster taken as part of the Herschel Virgo Cluster Survey (HeViCS). We select transition-type objects in the nearest cluster, Virgo, based on spectral diagnostics indicative for their residual or ongoing star formation. We detect dust Md ~ 10^{5-6} Msun in 36% of the transition-type dwarfs located on the high end of the stellar mass distribution. This suggests that the dust reservoirs present in non-detections fall just below the Herschel detection limit (< 1.1x10^5 Msun). Dust scaling relations support the hypothesis of a transformation between infalling late-type galaxies to quiescent low-mass spheroids governed by environmental effects, with dust-to-stellar mass fractions for transition-type dwarfs in between values characteristic for late-type objects and the lower dust fractions observed in early-type galaxies. Several transition-type dwarfs demonstrate blue central cores, hinting at the radially outside-in removal of gas and quenching of star formation activity. The fact that dust is also confined to the inner regions suggests that metals are stripped in the outer regions along with the gas. In the scenario of most dust being stripped from the galaxy along with the gas, we argue that... (abridged)

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 exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8<z<2.2 range. Sources with mass in the 10< log(M/Msun) <11 range are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the SMBH grows and is obscured by dust, is followed by an AGN dominated phase, then evolving toward a local elliptical. On the other hand, moderately starforming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.

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.

LISA Pathfinder is a mission planned by the European Space Agency to test the key technologies that will allow the detection of gravitational waves in space. The instrument on-board, the LISA Technology package, will undergo an exhaustive campaign of calibrations and noise characterisation campaigns in order to fully describe the noise model. Data analysis plays an important role in the mission and for that reason the data analysis team has been developing a toolbox which contains all the functionalities required during operations. In this contribution we give an overview of recent activities, focusing on the improvements in the modelling of the instrument and in the data analysis campaigns performed both with real and simulated data.

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

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

We study the emission by dust and stars in the Large and Small Magellanic Clouds, a pair of low-metallicity nearby galaxies, as traced by their spatially resolved spectral energy distributions (SEDs). This project combines Herschel Space Observatory PACS and SPIRE far-infrared photometry with other data at infrared and optical wavelengths. We build maps of dust and stellar luminosity and mass of both Magellanic Clouds, and analyze the spatial distribution of dust/stellar luminosity and mass ratios. These ratios vary considerably throughout the galaxies, generally between the range $0.01\leq L_{\rm dust}/L_\ast\leq 0.6$ and $10^{-4}\leq M_{\rm dust}/M_\ast\leq 4\times10^{-3}$. We observe that the dust/stellar ratios depend on the interstellar medium (ISM) environment, such as the distance from currently or previously star-forming regions, and on the intensity of the interstellar radiation field (ISRF). In addition, we construct star formation rate (SFR) maps, and find that the SFR is correlated with the dust/stellar luminosity and dust temperature in both galaxies, demonstrating the relation between star formation, dust emission and heating, though these correlations exhibit substantial scatter.

We present Herschel observations of the Fornax cluster at 100, 160, 250, 350 and 500u with a spatial resolution of 7 - 36 arc sec (10" = 1 kpc at d_Fornax=17.9 Mpc). We define a sample of 11 bright galaxies, selected at 500u, directly comparable with our past work on Virgo. We find good agreement with previous observations made by IRAS and Planck. The FIR luminosity density is higher (factor of three) in Fornax compared to Virgo. The 100u (42.5-122.5u) luminosity is two orders of magnitude larger in Fornax than in the local field as measured by IRAS. Using stellar (L_{0.4-2.5}) and FIR (L_{100-500}) luminosities we estimate a mean optical depth of tau=0.4+/-0.1 - the same value as Virgo. For 10 of the 11 galaxies (NGC1399 excepted) we fit a modified blackbody curve (beta=2.0) to the SEDs to derive dust masses and temperatures of 10^{6.54-8.35} M_0 and T=14.6-24.2K respectively, comparable to Virgo. The derived stars-to-gas(atomic) and gas(atomic)-to-dust ratios vary from 1.1-67.6 and 9.8-436.5 respectively, again consistent with Virgo. Fornax is a mass overdensity in stars and dust of about 120 compared to the local field (30 for Virgo). Fornax and Virgo are both a factor of 6 lower over densities in gas(atomic) than in stars and dust indicating loss of gas, but not dust and stars, in the cluster environment. As the brightest FIR source in either Fornax and Virgo, NGC1365 is detected by Planck. The Planck data fit the PACS/SPIRE SED out to 1382u with no evidence of other sources of emission ('spinning dust', free-free, synchrotron). At the opposite end of the scale NGC1399 is detected only at 500$\mu$m with the emission probably arising from the nuclear radio source rather than inter-stellar dust.

The Herschel Virgo Cluster Survey (HeViCS) is the deepest, confusion-limited survey of the Virgo Cluster at far-infrared (FIR) wavelengths. The entire survey at full depth covers $\sim$55 sq. deg. in 5 bands (100-500 \micron), encompassing the areas around the central dominant elliptical galaxies (M87, M86 & M49) and extends as far as the NW cloud, the W cloud and the Southern extension. The survey extends beyond this region with lower sensitivity so that the total area covered is 84 sq. deg. In this paper we describe the data, the data acquisition techniques and present the detection rates of the optically selected Virgo Cluster Catalogue (VCC). We detect 254 (34%) of 750 VCC galaxies found within the survey boundary in at least one band and 171 galaxies are detected in all five bands. For the remainder of the galaxies we have measured strict upper limits for their FIR emission. The population of detected galaxies contains early- as well as late-types although the latter dominate the detection statistics. We have modelled 168 galaxies, showing no evidence of a strong synchrotron component in their FIR spectra, using a single-temperature modified blackbody spectrum with a fixed emissivity index ($\beta = 2$). A study of the $\chi^2$ distribution indicates that this model is not appropriate in all cases, and this is supported by the FIR colours which indicate a spread in $\beta$=1--2. Statistical comparison of the dust mass and temperature distributions from 140 galaxies with $\chi^2_{dof=3} < 7.8$ (95% confidence level) shows that late-types have typically colder, more massive dust reservoirs; the early-type dust masses have a mean of ${\rm log}(<M> / M_{\sun}) = 6.3 \pm 0.3 $, while for late-types ${\rm log}(<M> / M_{\sun}) =7.1 \pm 0.1$... (abridged)

We combine new dust continuum observations of the edge-on spiral galaxy NGC 4565 in all Herschel/SPIRE (250, 350, 500 micron) wavebands, obtained as part of the Herschel Reference Survey, and a large set of ancillary data (Spitzer, SDSS, GALEX) to analyze its dust energy balance. We fit a radiative transfer model for the stars and dust to the optical maps with the fitting algorithm FitSKIRT. To account for the observed UV and mid-infrared emission, this initial model was supplemented with both obscured and unobscured star-forming regions. Even though these star-forming complexes provide an additional heating source for the dust, the far-infrared/submillimeter emission long wards of 100 micron is underestimated by a factor of 3-4. This inconsistency in the dust energy budget of NGC 4565 suggests that a sizable fraction (two-thirds) of the total dust reservoir (Mdust ~ 2.9e+8 Msun) consists of a clumpy distribution with no associated young stellar sources. The distribution of those dense dust clouds would be in such a way that they remain unresolved in current far-infrared/submillimeter observations and hardly comtribute to the attenuation at optical wavelengths. More than two-thirds of the dust heating in NGC 4565 is powered by the old stellar population, with localized embedded sources supplying the remaining dust heating in NGC 4565. The results from this detailed dust energy balance study in NGC 4565 is consistent with that of similar analyses of other edge-on spirals.

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

We have obtained Herschel images at five wavelengths from 100 to 500 micron of a ~5.5x2.5 degree area centred on the local galaxy M31 (Andromeda), our nearest neighbour spiral galaxy, as part of the Herschel guaranteed time project "HELGA". The main goals of HELGA are to study the characteristics of the extended dust emission, focusing on larger scales than studied in previous observations of Andromeda at an increased spatial resolution, and the obscured star formation. In this paper we present data reduction and Herschel maps, and provide a description of the far-infrared morphology, comparing it with features seen at other wavelengths. We use high--resolution maps of the atomic hydrogen, fully covering our fields, to identify dust emission features that can be associated to M31 with confidence, distinguishing them from emission coming from the foreground Galactic cirrus. Thanks to the very large extension of our maps we detect, for the first time at far-infrared wavelengths, three arc-like structures extending out to ~21, ~26 and ~31 kpc respectively, in the south-western part of M31. The presence of these features, hosting ~2.2e6 Msol of dust, is safely confirmed by their detection in HI maps. Overall, we estimate a total dust mass of ~5.8e7 Msol, about 78% of which is contained in the two main ring-like structures at 10 and 15 kpc, at an average temperature of 16.5 K. We find that the gas-to-dust ratio declines exponentially as a function of the galacto-centric distance, in agreement with the known metallicity gradient, with values ranging from 66 in the nucleus to ~275 in the outermost region. [Abridged]

The first complete submillimetre spectrum (190-670um) of the Seyfert 2 galaxy NGC1068 has been observed with the SPIRE Fourier Transform Spectrometer onboard the {\it Herschel} Space Observatory. The sequence of CO lines (Jup=4-13), lines from water, the fundamental rotational transition of HF, two o-H_2O+ lines and one line each from CH+ and OH+ have been detected, together with the two [CI] lines and the [NII]205um line. The observations in both single pointing mode with sparse image sampling and in mapping mode with full image sampling allow us to disentangle two molecular emission components, one due to the compact circum-nuclear disk (CND) and one from the extended region encompassing the star forming ring (SF-ring). Radiative transfer models show that the two CO components are characterized by density of n(H_2)=10^4.5 and 10^2.9 cm^-3 and temperature of T=100K and 127K, respectively. The comparison of the CO line intensities with photodissociation region (PDR) and X-ray dominated region (XDR) models, together with other observational constraints, such as the observed CO surface brightness and the radiation field, indicate that the best explanation for the CO excitation of the CND is an XDR with density of n(H_2) 10^4 cm^-3 and X-ray flux of 9 erg s^-1 cm^-2, consistent with illumination by the active galactic nucleus, while the CO lines in the SF-ring are better modeled by a PDR. The detected water transitions, together with those observed with the \her \sim PACS Spectrometer, can be modeled by an LVG model with low temperature (T_kin \sim 40K) and high density (n(H_2) in the range 10^6.7-10^7.9 cm^-3).

We present an analysis of the dust and gas in Andromeda, using Herschel images sampling the entire far-infrared peak. We fit a modified-blackbody model to ~4000 quasi-independent pixels with spatial resolution of ~140pc and find that a variable dust-emissivity index (beta) is required to fit the data. We find no significant long-wavelength excess above this model suggesting there is no cold dust component. We show that the gas-to-dust ratio varies radially, increasing from ~20 in the center to ~70 in the star-forming ring at 10kpc, consistent with the metallicity gradient. In the 10kpc ring the average beta is ~1.9, in good agreement with values determined for the Milky Way (MW). However, in contrast to the MW, we find significant radial variations in beta, which increases from 1.9 at 10kpc to ~2.5 at a radius of 3.1kpc and then decreases to 1.7 in the center. The dust temperature is fairly constant in the 10kpc ring (ranging from 17-20K), but increases strongly in the bulge to ~30K. Within 3.1kpc we find the dust temperature is highly correlated with the 3.6 micron flux, suggesting the general stellar population in the bulge is the dominant source of dust heating there. At larger radii, there is a weak correlation between the star formation rate and dust temperature. We find no evidence for 'dark gas' in M31 in contrast to recent results for the MW. Finally, we obtained an estimate of the CO X-factor by minimising the dispersion in the gas-to-dust ratio, obtaining a value of (1.9+/-0.4)x10^20 cm^-2 [K kms^-1]^-1.

We present GALEX far-ultraviolet (FUV) and near-ultraviolet (NUV) as well as SDSS g, r, i photometry and structural parameters for the Herschel Reference Survey, a magnitude-, volume-limited sample of nearby galaxies in different environments. We use this unique dataset to investigate the ultraviolet (UV) structural scaling relations of nearby galaxies and to determine how the properties of the UV disk vary with atomic hydrogen content and environment. We find a clear change of slope in the stellar mass vs. effective surface brightness relation when moving from the optical to the UV, with more massive galaxies having brighter optical but fainter UV surface brightnesses than smaller systems. A similar change of slope is also seen in the radius vs. surface brightness relation. By comparing our observations with the predictions of a simple multi-zone chemical model of galaxy evolution, we show that these findings are a natural consequence of a much more efficient inside-out growth of the stellar disk in massive galaxies. We confirm that isophotal radii are always a better proxy for the size of the stellar/star-forming disk than effective quantities and we show that the extent of the UV disk (normalized to the optical size) is strongly correlated to the integrated HI gas fraction. This relation still holds even when cluster spirals are considered, with HI-deficient systems having less extended star-forming disks than HI-normal galaxies. Interestingly, the star formation in the inner part of HI-deficient galaxies is significantly less affected by the removal of the atomic hydrogen, as expected in a simple ram-pressure stripping scenario. These results suggest that it is the amount of HI that regulates the growth of the star-forming disk in the outskirts of galaxies.

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 new Herschel-SPIRE imaging spectroscopy (194-671 microns) of the bright starburst galaxy M82. Covering the CO ladder from J=4-3 to J=13-12, spectra were obtained at multiple positions for a fully sampled ~ 3 x 3 arcminute map, including a longer exposure at the central position. We present measurements of 12CO, 13CO, [CI], [NII], HCN, and HCO+ in emission, along with OH+, H2O+ and HF in absorption and H2O in both emission and absorption, with discussion. We use a radiative transfer code and Bayesian likelihood analysis to model the temperature, density, column density, and filling factor of multiple components of molecular gas traced by 12CO and 13CO, adding further evidence to the high-J lines tracing a much warmer (~ 500 K), less massive component than the low-J lines. The addition of 13CO (and [CI]) is new and indicates that [CI] may be tracing different gas than 12CO. No temperature/density gradients can be inferred from the map, indicating that the single-pointing spectrum is descriptive of the bulk properties of the galaxy. At such a high temperature, cooling is dominated by molecular hydrogen. Photon-dominated region (PDR) models require higher densities than those indicated by our Bayesian likelihood analysis in order to explain the high-J CO line ratios, though cosmic-ray enhanced PDR models can do a better job reproducing the emission at lower densities. Shocks and turbulent heating are likely required to explain the bright high-J emission.

Using the far-infrared emission, as observed by the Herschel Virgo Cluster Survey (HeViCS), and the integrated HI and CO brightness, we infer the dust and total gas mass for a magnitude limited sample of 35 metal rich spiral galaxies in Virgo. The CO flux correlates tightly and linearly with far-infrared fluxes observed by Herschel. Molecules in these galaxies are more closely related to cold dust rather than to dust heated by star formation or to optical/NIR brightness. We show that dust mass establishes a stronger correlation with the total gas mass than with the atomic or molecular component alone. The dust-to-gas ratio increases as the HI deficiency increases, but in highly HI deficient galaxies it stays constant. Dust is in fact less affected than atomic gas by weak cluster interactions, which remove most of the HI gas from outer and high latitudes regions. Highly disturbed galaxies, in a dense cluster environment, can instead loose a considerable fraction of gas and dust from the inner regions of the disk keeping constant the dust-to-gas ratio. There is evidence that the molecular phase is also quenched. This quencing becomes evident by considering the molecular gas mass per unit stellar mass. Its amplitude, if confirmed by future studies, highlights that molecules are missing in Virgo HI deficient spirals, but to a somewhat lesser extent than dust.

The Herschel Reference Survey (HRS) is a guaranteed time Herschel key project aimed at studying the physical properties of the interstellar medium in galaxies of the nearby universe. This volume limited, K-band selected sample is composed of galaxies spanning the whole range of morphological types (from ellipticals to late-type spirals) and environments (from the field to the centre of the Virgo Cluster). We present flux density measurements of the whole sample of 323 galaxies of the HRS in the three bands of the Spectral and Photometric Imaging Receiver (SPIRE), at 250, 350 and 500 microns. Aperture photometry is performed on extended galaxies and point spread function (PSF) fitting on timeline data for unresolved objects; we carefully estimate errors and upper limits. The flux densities are found to be in good agreement with those of the HeViCS and KINGFISH key projects in all SPIRE bands, and of the Planck consortium at 350 and 550 microns, for the galaxies in common. This submillimetre catalogue of nearby galaxies is a benchmark for the study of the dust properties in the local universe, giving the zero redshift reference for any cosmological survey.

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]

The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 \mu m), and Herschel-PACS (at 100 and 160 \mu m), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the re-processed optical and ultra-violet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multi-wavelength understanding of galaxy formation and evolution. The survey will detect of order 100,000 galaxies at 5\sigma in some of the best studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; rapidly identify unusual objects; and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include: the total infrared emission of galaxies; the evolution of the luminosity function; the clustering properties of dusty galaxies; and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques. This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.

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.

We combine new Herschel/SPIRE sub-millimeter observations with existing multiwavelength data to investigate the dust scaling relations of the Herschel Reference Survey, a magnitude-, volume-limited sample of ~300 nearby galaxies in different environments. We show that the dust-to-stellar mass ratio anti-correlates with stellar mass, stellar mass surface density and NUV-r colour across the whole range of parameters covered by our sample. Moreover, the dust-to-stellar mass ratio decreases significantly when moving from late- to early-type galaxies. These scaling relations are similar to those observed for the HI gas-fraction, supporting the idea that the cold dust is tightly coupled to the cold atomic gas component in the interstellar medium. We also find a weak increase of the dust-to-HI mass ratio with stellar mass and colour but no trend is seen with stellar mass surface density. By comparing galaxies in different environments we show that, although these scaling relations are followed by both cluster and field galaxies, HI-deficient systems have, at fixed stellar mass, stellar mass surface density and morphological type systematically lower dust-to-stellar mass and higher dust-to-HI mass ratios than HI-normal/field galaxies. This provides clear evidence that dust is removed from the star-forming disk of cluster galaxies but the effect of the environment is less strong than what is observed in the case of the HI disk. Such effects naturally arise if the dust disk is less extended than the HI and follows more closely the distribution of the molecular gas phase, i.e., if the dust-to-atomic gas ratio monotonically decreases with distance from the galactic center.

We provide an overview of ancillary 24, 70, and 160 micron data from the Multiband Imaging Photometer for Spitzer (MIPS) that are intended to complement the 70-500 micron Herschel Space Observatory photometry data for nearby galaxies obtained by the Herschel-SPIRE Local Galaxies Guaranteed Time Programs and the Herschel Virgo Cluster Survey. The MIPS data can be used to extend the photometry to wave bands that are not observed in these Herschel surveys and to check the photometry in cases where Herschel performs observations at the same wavelengths. Additionally, we measured globally-integrated 24-160 micron flux densities for the galaxies in the sample that can be used for the construction of spectral energy distributions. Using MIPS photometry published by other references, we have confirmed that we are obtaining accurate photometry for these galaxies.

We present Herschel observations of 62 Early-Type Galaxies (ETGs), including 39 galaxies morphologically classified as S0+S0a and 23 galaxies classified as ellipticals using SPIRE at 250, 350 and 500 microns (and PACS 100 and 160 microns for 19 sources) as part of the volume-limited Herschel Reference Survey. We detect dust emission in 24% of the ellipticals and 62% of the S0s. The mean temperature of the dust is 23.9+/-0.8 K, warmer than that found for late-type galaxies in the Virgo Cluster. Including the non-detections, the mean dust mass is log(Mdust) = 5.9+/-0.1 and 5.2+/-0.1 Msun for the S0s and elliptical galaxies respectively. The mean dust-to-stellar mass is log(Mdust/Mstar) = -4.4+/-0.1 (S0s) and -5.8+/-0.1 (ellipticals). Virtually all the galaxies lie close to the red sequence yet the large number of detections of cool dust, the gas-to-dust ratios and the ratios of far-infrared to radio emission all suggest that many ETGs contain a cool interstellar medium similar to that in late-type galaxies. The mean dust-to-stellar mass ratio for S0s is approximatly a factor of ten less than for early-type spirals and the sizes of the dust sources in the S0s are also much smaller. We show that the difference cannot be explained by either the different bulge-to-disk ratios or environmental effects such as ram-pressure stripping. The wide range in the dust-to-stellar mass ratio for ETGs and the lack of a correlation between dust mass and optical luminosity suggest that much of the dust in the ETGs detected by Herschel has been acquired as the result of gravitational interactions; these interactions are unlikely to have had a major effect on the stellar masses of the ETGs. The Herschel observations tentatively suggest that in the most massive ETGs, the mass of the interstellar medium is unconnected to the evolution of the stellar populations.

It has often been suggested that an alternative to the standard CO/21-cm method for estimating the mass of the interstellar medium (ISM) in a galaxy might be to estimate the mass of the ISM from the continuum dust emission. In this paper, we investigate the potential of this technique using Herschel observations of ten galaxies in the Herschel Reference Survey and in the Herschel Virgo Cluster Survey. We show that the emission detected by Herschel is mostly from dust that has a temperature and emissivity index similar to that of dust in the local ISM in our galaxy, with the temperature generally increasing towards the centre of each galaxy. We calibrate the dust method using the CO and 21-cm observations to provide an independent estimate of the mass of hydrogen in each galaxy, solving the problem of the uncertain `X factor' for the molecular gas by minimizing the dispersion in the ratio of the masses estimated using the two methods. With the calibration for the dust method and the estimate of the X-factor produced in this way, the dispersion in the ratio of the two gas masses is 30%, which gives an upper limit on the fundamental accuracy of the dust method. The calibration we obtain for the dust method is very similar to an independent Herschel measurement for M31 and to the calibration for the Milky Way from Planck measurements.

We study the far infrared (60-500 $\mu$m) colours of late-type galaxies in the $Herschel$ Reference Survey, a K-band selected, volume limited sample of nearby galaxies. The far infrared colours are correlated with each other, with tighter correlations for the indices that are closer in wavelength. We also compare the different colour indices to various tracers of the physical properties of the target galaxies, such as the surface brightness of the ionising and non-ionising stellar radiation, the dust attenuation and the metallicity. The emission properties of the cold dust dominating the far infrared spectral domain are regulated by the properties of the interstellar radiation field. Consistent with that observed in nearby, resolved galaxies, our analysis shows that the ionising and the non-ionising stellar radiation, including that emitted by the most evolved, cold stars, both contribute to the heating of the cold dust component. This work also shows that metallicity is another key parameter characterising the cold dust emission of normal, late-type galaxies. A single modified black body with a grain emissivity index $\beta$=1.5 better fits the observed SPIRE flux density ratios $S250/S350$ vs. $S350/S500$ than $\beta$=2, although values of $\beta$ $\simeq$ 2 are possible in metal rich, high surface brightness galaxies. Values of $\beta$ $\lesssim$ 1.5 better represent metal poor, low surface brightness objects. This observational evidence provides strong constraints for dust emission models of normal, late type galaxies.

We examine the dust and gas properties of the nearby, barred galaxy M83, which is part of the Very Nearby Galaxy Survey. Using images from the PACS and SPIRE instruments of Herschel, we examine the dust temperature and dust mass surface density distribution. We find that the nuclear, bar and spiral arm regions exhibit higher dust temperatures and masses compared to interarm regions. However, the distribution of dust temperature and mass are not spatially coincident. Assuming a trailing spiral structure, the dust temperature peaks in the spiral arms lie ahead of the dust surface density peaks. The dust mass surface density correlates well with the distribution of molecular gas as traced by CO (J=3-2) images (JCMT) and the star formation rate as traced by H?2 with a correction for obscured star formation using 24 micron emission. Using HI images from THINGS to trace the atomic gas component, we make total gas mass surface density maps and calculate the gas-to-dust ratio. We find a mean gas-to-dust ratio of 84 \pm 4 with higher values in the inner region assuming a constant CO-to-H2 conversion factor. We also examine the gas-to-dust ratio using CO-to-H2 conversion factor that varies with metallicity.

We use the photoionisation and dust radiative transfer code MOCASSIN to create a model of the dwarf irregular galaxy NGC 4449. The best-matching model reproduces the global optical emission line fluxes and the observed spectral energy distribution (SED) spanning wavelengths from the UV to sub-mm, and requires the bolometric luminosity of 6.25e9 Lsolar for the underlying stellar component, M_d/M_g of 1/680 and M_d of 2.2e6 Msolar.

We investigate the potential of submm-mm and submm-mm-radio photometric redshifts using a sample of mm-selected sources as seen at 250, 350 and 500 {\mu}m by the SPIRE instrument on Herschel. From a sample of 63 previously identified mm-sources with reliable radio identifications in the GOODS-N and Lockman Hole North fields 46 (73 per cent) are found to have detections in at least one SPIRE band. We explore the observed submm/mm colour evolution with redshift, finding that the colours of mm-sources are adequately described by a modified blackbody with constant optical depth {\tau} = ({\nu}/{\nu}0)^{\beta} where {\beta} = +1.8 and {\nu}0 = c/100 {\mu}m. We find a tight correlation between dust temperature and IR luminosity. Using a single model of the dust temperature and IR luminosity relation we derive photometric redshift estimates for the 46 SPIRE detected mm-sources. Testing against the 22 sources with known spectroscopic, or good quality optical/near-IR photometric, redshifts we find submm/mm photometric redshifts offer a redshift accuracy of |z|/(1+z) = 0.16 (< |z| >= 0.51). Including constraints from the radio-far IR correlation the accuracy is improved to |z|/(1 + z) = 0.15 (< |z| >= 0.45). We estimate the redshift distribution of mm-selected sources finding a significant excess at z > 3 when compared to ~ 850 {\mu}m selected samples.

We describe the Herschel Virgo Cluster Survey (HeViCS) and the first data that cover the complete survey area (four 4 x 4 deg2 regions). We use these data to measure and compare the global far infrared properties of 78 optically bright galaxies that are selected at 500 \mum and detected in all five far-infrared bands. We show that our measurements and calibration are broadly consistent with previous data obtained by IRAS, ISO, Spitzer and Planck. We use SPIRE and PACS photometry data to produce 100, 160, 250, 350 and 500 \mum cluster luminosity distributions. These luminosity distributions are not power laws, but peaked, with small numbers of both faint and bright galaxies. We measure a cluster 100-500 micron far-infrared luminosity density of 1.6(7.0) \pm 0.2 x 10^9 Lsun/Mpc3. This compares to a cluster 0.4-2.5 \mum optical luminosity density of 5.0(20.0) x 10^9 Lsun/Mpc3, some 3.2(2.9) times larger than the far-infrared. A typical photon originates from an optical depth of 0.4\pm0.1. Most of our sample galaxies are well fitted by a single modified blackbody (beta=2), leading to a mean dust mass of log Mdust = 7.31 Msun and temperature of 20.0 K. We also derive both stellar and atomic hydrogen masses from which we calculate mean values for the stars:gas(atomic) and gas(atomic): dust mass ratios of 15.1 and 58.2 respectively. Using our derived dust, atomic gas and stellar masses we estimate cluster mass densities of 8.6(27.8) x 10^6, 4.6(13.9) x 10^8, 7.8(29.7) x 10^9 Msun/Mpc3, respectively for dust, atomic gas and stars. These values are higher than those derived for field galaxies by factors of 39(126), 6(18) and 34(129) respectively. In the above luminosity/mass densities are given using the whole sample with values in brackets using just those galaxies that lie between 17 and 23 Mpc. We provide a data table of flux densities in all the Herschel bands for all 78 bright Virgo cluster galaxies.

Aims: In this paper, we perform detailed modelling of the Spitzer and Herschel observations of the LMC, in order to: (i) systematically study the uncertainties and biases affecting dust mass estimates; and to (ii) explore the peculiar ISM properties of the LMC. Methods: To achieve these goals, we have modelled the spatially resolved SEDs with two alternate grain compositions, to study the impact of different submillimetre opacities on the dust mass. We have rigorously propagated the observational errors (noise and calibration) through the entire fitting process, in order to derive consistent parameter uncertainties. Results: First, we show that using the integrated SED leads to underestimating the dust mass by ~50 % compared to the value obtained with sufficient spatial resolution, for the region we studied. This might be the case, in general, for unresolved galaxies. Second, we show that Milky Way type grains produce higher gas-to-dust mass ratios than what seems possible according to the element abundances in the LMC. A spatial analysis shows that this dilemma is the result of an exceptional property: the grains of the LMC have on average a larger intrinsic submm opacity (emissivity index beta~1.7 and opacity kappa_abs(160 microns)=1.6 m2/kg) than those of the Galaxy. By studying the spatial distribution of the gas-to-dust mass ratio, we are able to constrain the fraction of unseen gas mass between ~10, and ~100 % and show that it is not sufficient to explain the gas-to-dust mass ratio obtained with Milky Way type grains. Finally, we confirm the detection of a 500 microns extended emission excess with an average relative amplitude of ~15 %, varying up to 40 %. This excess anticorrelates well with the dust mass surface density. Although we do not know the origin of this excess, we show that it is unlikely the result of very cold dust, or CMB fluctuations.

We present Herschel SPIRE-FTS observations of Arp~220, a nearby ULIRG. The FTS continuously covers 190 -- 670 microns, providing a good measurement of the continuum and detection of several molecular and atomic species. We detect luminous CO (J = 4-3 to 13-12) and water ladders with comparable total luminosity; very high-J HCN absorption; OH+, H2O+, and HF in absorption; and CI and NII. Modeling of the continuum yields warm dust, with T = 66 K, and an unusually large optical depth of ~5 at 100 microns. Non-LTE modeling of the CO shows two temperature components: cold molecular gas at T ~ 50 K and warm molecular gas at T ~1350 K. The mass of the warm gas is 10% of the cold gas, but dominates the luminosity of the CO ladder. The temperature of the warm gas is in excellent agreement with H2 rotational lines. At 1350 K, H2 dominates the cooling (~20 L_sun/M_sun) in the ISM compared to CO (~0.4 L_sun/M_sun). We found that only a non-ionizing source such as the mechanical energy from supernovae and stellar winds can excite the warm gas and satisfy the energy budget of ~20 L_sun/M_sun. We detect a massive molecular outflow in Arp 220 from the analysis of strong P-Cygni line profiles observed in OH+, H2O+, and H2O. The outflow has a mass > 10^{7} M_sun and is bound to the nuclei with velocity < 250 km/s. The large column densities observed for these molecular ions strongly favor the existence of an X-ray luminous AGN (10^{44} ergs/s) in Arp 220.

Star formation is arguably the most important physical process in the cosmos. It is a fundamental driver of galaxy evolution and the ultimate source of most of the energy emitted by galaxies. A correct interpretation of star formation rate (SFR) measures is therefore essential to our understanding of galaxy formation and evolution. Unfortunately, however, no single SFR estimator is universally available or even applicable in all circumstances: the numerous galaxies found in deep surveys are often too faint (or too distant) to yield significant detections with most standard SFR measures, and until now there have been no global, multi-band observations of nearby galaxies that span all the conditions under which star-formation is taking place. To address this need in a systematic way, we have undertaken a multi-band survey of all types of star-forming galaxies in the local Universe. This project, the Star Formation Reference Survey (SFRS), is based on a statistically valid sample of 369 nearby galaxies that span all existing combinations of dust temperature, SFR, and specific SFR. Furthermore, because the SFRS is blind with respect to AGN fraction and environment it serves as a means to assess the influence of these factors on SFR. Our panchromatic global flux measurements (including GALEX FUV+NUV, SDSS ugriz, 2MASS JHKs, Spitzer 3-8{\mu}m, and others) furnish uniform SFR measures and the context in which their reliability can be assessed. This paper describes the SFRS survey strategy, defines the sample, and presents the multi-band photometry collected to date.

We use Spitzer Space Telescope and Herschel Space Observatory far-infrared data along with ground-based optical and near-infrared data to understand how dust heating in the nearby face-on spiral galaxies M81, M83, and NGC 2403 is affected by the starlight from all stars and by the radiation from star forming regions. We find that 70/160 micron surface brightness ratios tend to be more strongly influenced by star forming regions. However, the 250/350 micron and 350/500 micron surface brightness ratios are more strongly affected by the light from the total stellar populations, suggesting that the dust emission at >250 microns originates predominantly from a component that is colder than the dust seen at <160 microns and that is relatively unaffected by star formation activity. We conclude by discussing the implications of this for modelling the spectral energy distributions of both nearby and more distant galaxies and for using far-infrared dust emission to trace star formation.

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

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

Using Herschel data from the Open Time Key Project the Herschel Virgo Cluster Survey (HeViCS), we investigated the relationship between the metallicity gradients expressed by metal abundances in the gas phase as traced by the chemical composition of HII regions, and in the solid phase, as traced by the dust-to-gas mass ratio. We derived the radial gradient of the dust-to-gas mass ratio for all galaxies observed by HeViCS whose metallicity gradients are available in the literature. They are all late type Sbc galaxies, namely NGC4254, NGC4303, NGC4321, and NGC4501. We examined different dependencies on metallicity of the CO-to-H$_2$ conversion factor (\xco), used to transform the $^{12}$CO observations into the amount of molecular hydrogen. We found that in these galaxies the dust-to-gas mass ratio radial profile is extremely sensitive to choice of the \xco\ value, since the molecular gas is the dominant component in the inner parts. We found that for three galaxies of our sample, namely NGC4254, NGC4321, and NGC4501, the slopes of the oxygen and of the dust-to-gas radial gradients agree up to $\sim$0.6-0.7R$_{25}$ using \xco\ values in the range 1/3-1/2 Galactic \xco. For NGC4303 a lower value of \xco$\sim0.1\times$ 10$^{20}$ is necessary. We suggest that such low \xco\ values might be due to a metallicity dependence of \xco (from close to linear for NGC4254, NGC4321, and NGC4501 to superlinear for NGC4303), especially in the radial regions R$_G<$0.6-0.7R$_{25}$ where the molecular gas dominates. On the other hand, the outer regions, where the atomic gas component is dominant, are less affected by the choice of \xco, and thus we cannot put constraints on its value.

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.

We present a galaxy (SAGE1CJ053634.78-722658.5) at a redshift of 0.14 of which the IR is entirely dominated by emission associated with the AGN. We present the 5-37 um Spitzer/IRS spectrum and broad wavelength SED of SAGE1CJ053634, an IR point-source detected by Spitzer/SAGE (Meixner et al 2006). The source was observed in the SAGE-Spec program (Kemper et al., 2010) and was included to determine the nature of sources with deviant IR colours. The spectrum shows a redshifted (z=0.14+-0.005) silicate emission feature with an exceptionally high feature-to-continuum ratio and weak polycyclic aromatic hydrocarbon (PAH) bands. We compare the source with models of emission from dusty tori around AGNs from Nenkova et al. (2008). We present a diagnostic diagram that will help to identify similar sources based on Spitzer/MIPS and Herschel/PACS photometry. The SED of SAGE1CJ053634 is peculiar because it lacks far-IR emission and a clear stellar counterpart. We find that the SED and the IR spectrum can be understood as emission originating from the inner ~10 pc around an accreting black hole. There is no need to invoke emission from the host galaxy, either from the stars or from the interstellar medium, although a possible early-type host galaxy cannot be excluded based on the SED analysis. The hot dust around the accretion disk gives rise to a continuum, which peaks at 4 um, whereas the strong silicate features may arise from optically thin emission of dusty clouds within ~10 pc around the black hole. The weak PAH emission does not appear to be linked to star formation, as star formation templates strongly over-predict the measured far-IR flux levels. The SED of SAGE1CJ053634 is rare in the local universe but may be more common in the more distant universe. The conspicuous absence of host-galaxy IR emission places limits on the far-IR emission arising from the dusty torus alone.

We are studying the effects of submm observations on the total dust mass and thus dust-to-gas mass ratio measurements. We gather a wide sample of galaxies that have been observed at submm wavelengths to model their Spectral Energy Distributions using submm observations and then without submm observational constraints in order to quantify the error on the dust mass when submm data are not available. Our model does not make strong assumptions on the dust temperature distribution to precisely avoid submm biaises in the study. Our sample includes 52 galaxies observed at submm wavelengths. Out of these, 9 galaxies show an excess in submm which is not accounted for in our fiducial model, most of these galaxies being low- metallicity dwarfs. We chose to add an independant very cold dust component (T=10K) to account for this excess. We find that metal-rich galaxies modelled with submm data often show lower dust masses than when modelled without submm data. Indeed, these galaxies usually have dust SEDs that peaks at longer wavelengths and require constraints above 160 um to correctly position the peak and sample the submillimeter slope of their SEDs and thus correctly cover the dust temperature distribution. On the other hand, some metal-poor dwarf galaxies modelled with submm data show higher dust masses than when modelled without submm data. Using submm constraints for the dust mass estimates, we find a tightened correlation of the dust-to-gas mass ratio with the metallicity of the galaxies. We also often find that when there is a submm excess present, it occurs preferentially in low-metallicity galaxies. We analyse the conditions for the presence of this excess and find a relation between the 160/850 um ratio and the submm excess.

The GALEX Ultraviolet Virgo Cluster Survey (GUViCS) is a complete blind survey of the Virgo cluster covering about 40 sq. deg. in the far UV (FUV, lambda_eff=1539A, Delta-lambda=442A) and about 120 sq. deg. in the near UV (NUV, lambda_eff=2316A, Delta-lambda=1060A). The goal of the survey is to study the ultraviolet (UV) properties of galaxies in a rich cluster environment, spanning a wide luminosity range from giants to dwarfs, and regardless of prior knowledge of their star formation activity. The UV data will be combined with those in other bands (optical: NGVS; far-infrared - submm: HeViCS; HI: ALFALFA) and with our multizone chemo-spectrophotometric models of galaxy evolution to make a complete and exhaustive study of the effects of the environment on the evolution of galaxies in high density regions. We present here the scientific objectives of the survey, describing the observing strategy and briefly discussing different data reduction techniques. Using UV data already in-hand for the central 12 sq. deg. we determine the FUV and NUV luminosity functions of the Virgo cluster core for all cluster members and separately for early- and late-type galaxies and compare it to the one obtained in the field and other nearby clusters (Coma, A1367). This analysis shows that the FUV and NUV luminosity functions of the core of the Virgo clusters are flatter (alpha about -1.1) than those determined in Coma and A1367. We discuss the possible origin of this difference

This paper presents a quantitative assessment of the performance of the upcoming LISA Pathfinder geodesic explorer mission. The findings are based on the results of extensive ground testing and simulation campaigns using flight hardware and flight control and operations algorithms. The results show that, for the central experiment of measuring the stray differential acceleration between the LISA test masses, LISA Pathfinder will be able to verify the overall acceleration noise to within a factor two of the LISA requirement at 1 mHz and within a factor 6 at 0.1 mHz. We also discuss the key elements of the physical model of disturbances, coming from LISA Pathfinder and ground measurement, that will guarantee the LISA performance.

The integrated Spectral Energy Distributions of the Large and Small Magellanic Cloud appear significantly flatter than expected from dust models based on their FIR and radio emission. The origin of this millimetre excess is still unexplained, and is here investigated using the Planck data. The background CMB contribution is subtracted using an ILC method performed locally around the galaxies. The foreground emission from the Milky Way is subtracted. After subtraction, the emission of both galaxies correlates closely with the gas emission of the LMC and SMC. The millimetre excess in the LMC can be explained by CMB fluctuations, but a significant excess is still present in the SMC SED. The Planck and IRIS data at 100 micron are combined to produce thermal dust temperature and optical depth maps of the two galaxies. The LMC temperature map shows the presence of a warm inner arm already found with the Spitzer data, but also shows the existence of a previously unidentified cold outer arm. Several cold regions are found along this arm, some of which are associated with known molecular clouds. The average emissivity spectral index is found to be consistent with beta=1.5 and beta=1.2 below 500 microns for the LMC and SMC respectively, significantly flatter than the values observed in the Milky Way. Furthermore, there is evidence in the SMC for a further flattening of the SED in the sub-mm. The spatial distribution of the millimetre dust excess in the SMC follows the gas and thermal dust distribution. Different models are explored in order to fit the dust emission in the SMC. It is concluded that the millimetre excess is unlikely to be caused by very cold dust emission and that it could be due to a combination of spinning dust. emission and thermal dust emission by more amorphous dust grains than those present in our Galaxy.

The all-sky coverage of the Planck Early Release Compact Source Catalogue (ERCSC) provides an unsurpassed survey of galaxies at submillimetre (submm) wavelengths, representing a major improvement in the numbers of galaxies detected, as well as the range of far-IR/submm wavelengths over which they have been observed. We here present the first results on the properties of nearby galaxies using these data. We match the ERCSC catalogue to IRAS-detected galaxies in the Imperial IRAS Faint Source Redshift Catalogue (IIFSCz), so that we can measure the spectral energy distributions (SEDs) of these objects from 60 to 850 microns. This produces a list of 1717 galaxies with reliable associations between Planck and IRAS, from which we select a subset of 468 for SED studies, namely those with strong detections in the three highest frequency Planck bands and no evidence of cirrus contamination. The SEDs are fitted using parametric dust models to determine the range of dust temperatures and emissivities. We find evidence for colder dust than has previously been found in external galaxies, with T<20K. Such cold temperatures are found using both the standard single temperature dust model with variable emissivity beta, or a two dust temperature model with beta fixed at 2. We also compare our results to studies of distant submm galaxies (SMGs) which have been claimed to contain cooler dust than their local counterparts. We find that including our sample of 468 galaxies significantly reduces the distinction between the two populations. Fits to SEDs of selected objects using more sophisticated templates derived from radiative transfer models confirm the presence of the colder dust found through parameteric fitting. We thus conclude that cold (T<20K) dust is a significant and largely unexplored component of many nearby galaxies.

The extragalactic background light at far-infrared wavelengths originates from optically-faint, dusty, star-forming galaxies in the universe with star-formation rates at the level of a few hundred solar masses per year. Due to the relatively poor spatial resolution of far-infrared telescopes, the faint sub-millimetre galaxies are challenging to study individually. Instead, their average properties can be studied using statistics such as the angular power spectrum of the background intensity variations. A previous attempt at measuring this power spectrum resulted in the suggestion that the clustering amplitude is below the level computed with a simple ansatz based on a halo model. Here we report a clear detection of the excess clustering over the linear prediction at arcminute angular scales in the power spectrum of brightness fluctuations at 250, 350, and 500 microns. From this excess, we find that sub-millimetre galaxies are located in dark matter halos with a minimum mass of log[M_min/M_sun ]= 11.5^+0.7_-0.2 at 350 microns. This minimum dark matter halo mass corresponds to the most efficient mass scale for star formation in the universe, and is lower than that predicted by semi-analytical models for galaxy formation.

We investigate the emission properties of polycyclic aromatic hydrocarbons (PAHs) in various metallicity environments with the Infrared Spectrograph on board Spitzer. Local giant HII regions are used as references as they enable access to the distinct interstellar medium components that contribute to the mid-infrared spectrum of star-forming galaxies: photodissociation regions (PDRs), photoionized gas, stellar clusters, and embedded regions. Three objects are considered, NGC3603 in the Milky Way, 30Doradus in the Large Magellanic Cloud, and N66 in the Small Magellanic Cloud. From the variations of the PAH/14um ratio, we find that PAHs are destroyed in the ionized gas for a radiation field such that [NeIII]/[NeII]>3. From the variations of the PAH/Hu-alpha ratio, we find that the PAH emission sources in the giant HII regions follow the same photodestruction law regardless of metallicity. We then compare these results with observations of starburst galaxies, HII galaxies, and blue compact dwarf galaxies (BCDs). While the integrated mid-infrared spectra of BCDs are reminiscent of a warm dusty ionized gas, we observe a significant contribution to the PAH emission in starburst galaxies that is not arising from PDRs.

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.

The objective of this paper is to construct a catalog providing the dust properties and the star formation efficiency (SFE) of the molecular clouds in the Large Magellanic Cloud (LMC). We use the infrared (IR) data obtained with the Spitzer telescope as part of the ``Surveying the Agents of a Galaxy's Evolution'' (SAGE) Legacy survey as well as the IRAS data. We also work with extinction (Av) maps of the LMC. A total of 272 molecular clouds have been detected in the LMC in a previous molecular survey, accounting for 230 giant molecular clouds and 42 smaller clouds. We perform correlations between the IR emission/extinction, and atomic and molecular gas tracers. We compare the atomic gas that surrounds the molecular cloud with the molecular gas in the cloud. Using a dust emission model, we derive the physical properties of dust in and outside the clouds: equilibrium temperature, emissivity and extinction. We also determine the luminosity of the interstellar radiation field intercepted by the cloud, and the total IR luminosity. Statistically, we do not find any significant difference in the dust properties between the atomic and the molecular phases. In particular we do not find evidence for a systematic decrease of the dust temperature in the molecular phase, with respect to the surrounding, presumably atomic gas. This is probably because giant molecular clouds are the sites of star formation, which heat the dust, while the smallest clouds are unresolved. The ratio between the IR luminosity and the cloud mass (LDust/Mgas) does not seem to correlate with Mgas. The highest value of the ratio we derived is 18.1 Lsol/Msol in the 30 Doradus region, which is known to be the most prominent star formation region of the LMC, while the most likely value is 0.5 and is representative of quiescent clouds. We provide a prescription to associate the various stages of star formation with its LDust/Mgas.

The origin of the far-infrared emission from the nearby radio galaxy M87 remains a matter of debate. Some studies find evidence of a far-infrared excess due to thermal dust emission, whereas others propose that the far-infrared emission can be explained by synchrotron emission without the need for an additional dust emission component. We observed M87 with PACS and SPIRE as part of the Herschel Virgo Cluster Survey (HeViCS). We compare the new Herschel data with a synchrotron model based on infrared, submm and radio data to investigate the origin of the far-infrared emission. We find that both the integrated SED and the Herschel surface brightness maps are adequately explained by synchrotron emission. At odds with previous claims, we find no evidence of a diffuse dust component in M87.

We present the classification of 197 point sources observed with the Infrared Spectrograph in the SAGE-Spec Legacy program on the Spitzer Space Telescope. We introduce a decision-tree method of object classification based on infrared spectral features, continuum and spectral energy distribution shape, bolometric luminosity, cluster membership, and variability information, which is used to classify the SAGE-Spec sample of point sources. The decision tree has a broad application to mid-infrared spectroscopic surveys, where supporting photometry and variability information are available. We use these classifications to make deductions about the stellar populations of the Large Magellanic Cloud and the success of photometric classification methods. We find 90 asymptotic giant branch (AGB) stars, 29 young stellar objects, 23 post-AGB objects, 19 red supergiants, eight stellar photospheres, seven background galaxies, seven planetary nebulae, two HII regions and 12 other objects, seven of which remain unclassified.

Dusty, star forming galaxies contribute to a bright, currently unresolved cosmic far-infrared background. Deep Herschel-SPIRE images designed to detect and characterize the galaxies that comprise this background are highly confused, such that the bulk lies below the classical confusion limit. We analyze three fields from the HerMES programme in all three SPIRE bands (250, 350, and 500 microns); parameterized galaxy number count models are derived to a depth of ~2 mJy/beam, approximately 4 times the depth of previous analyses at these wavelengths, using a P(D) (probability of deflection) approach for comparison to theoretical number count models. Our fits account for 64, 60, and 43 per cent of the far-infrared background in the three bands. The number counts are consistent with those based on individually detected SPIRE sources, but generally inconsistent with most galaxy number counts models, which generically overpredict the number of bright galaxies and are not as steep as the P(D)-derived number counts. Clear evidence is found for a break in the slope of the differential number counts at low flux densities. Systematic effects in the P(D) analysis are explored. We find that the effects of clustering have a small impact on the data, and the largest identified systematic error arises from uncertainties in the SPIRE beam.

We present first results of a study of the submillimetre (rest frame far-infrared) properties of z~3 Lyman Break Galaxies (LBGs) and their lower-redshift counterparts BX/BM galaxies, based on Herschel-SPIRE observations of the Northern field of the Great Observatories Origins Deep Survey (GOODS-N). We use stacking analysis to determine the properties of LBGs well below the current limit of the survey. Although LBGs are not detected individually, stacking the infrared luminous LBGs (those detected with Spitzer at 24 microns yields a statistically significant submm detection with mean flux <S_{250}>= 5.9+/-1.4 mJy confirming the power of SPIRE in detecting UV-selected high-redshift galaxies at submillimetre wavelengths. In comparison, the Spitzer 24 microns detected BX/BM galaxies appear fainter with a stacked value of <S_{250}> = 2.7 +/-0.8 mJy. By fitting the Spectral Energy Distributions (SEDs) we derive median infrared luminosities, L_{IR}, of 2.8x10^{12} Lsun and 1.5x10^{11} Lsun for z~3 LBGs and BX/BMs, respectively. We find that $L_{IR} estimates derived from present measurements are in good agreement with those based on UV data for z~2 BX/BM galaxies, unlike the case for z~3 infrared luminous LBGs where the UV underestimates the true $L_{IR}. Although sample selection effects may influence this result we suggest that differences in physical properties (such as morphologies, dust distribution and extent of star-forming regions) between z ~3 LBGs and z~2 BX/BMs may also play a significant role.

The Spectral and Photometric Imaging Receiver (SPIRE) on Herschel has been carrying out deep extragalactic surveys, one of whose aims is to establish spectral energy distributions (SED)s of individual galaxies spanning the infrared/submillimeter (IR/SMM) wavelength region. We report observations of the (IR/SMM) emission from the Lockman North field (LN) and Great Observatories Origins Deep Survey field North (GOODS-N). Because galaxy images in the wavelength range covered by Herschel generally represent a blend with contributions from neighboring galaxies, we present sets of galaxies in each field especially free of blending at 250, 350, and 500 microns. We identify the cumulative emission of these galaxies and the fraction of the far infrared cosmic background radiation they contribute. Our surveys reveal a number of highly luminous galaxies at redshift z ~< 3 and a novel relationship between infrared and visible emission that shows a dependence on luminosity and redshift.

We present the cross-identification and source photometry techniques used to process Herschel SPIRE imaging taken as part of the Herschel Multi-Tiered Extragalactic Survey (HerMES). Cross-identifications are performed in map-space so as to minimise source blending effects. We make use of a combination of linear inversion and model selection techniques to produce reliable cross-identification catalogues based on Spitzer MIPS 24 micron source positions. Testing on simulations and real Herschel observations show that this approach gives robust results for even the faintest sources S250~10 mJy. We apply our new technique to HerMES SPIRE observations taken as part of the science demostration phase of Herschel. For our real SPIRE observations we show that, for bright unconfused sources, our flux density estimates are in good agreement with those produced via more traditional point source detection methods (SussExtractor; Savage & Oliver et al. 2006) by Smith et al. 2010. When compared to the measured number density of sources in the SPIRE bands, we show that our method allows the recovery of a larger fraction of faint sources than these traditional methods. However this completeness is heavily dependant on the relative depth of the existing 24 micron catalogues and SPIRE imaging. Using our deepest multi-wavelength dataset in GOODS-N, we estimate that the use of shallow 24 micron in our other fields introduces an incompleteness at faint levels of between 20-40 per cent at 250 micron.

We study the dust properties of galaxies in the redshift range 0.1<z<2.8 observed by the Herschel Space Observatory in the field of the Great Observatories Origins Deep Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity (L_IR) and dust temperature (T_dust) of galaxies are derived from the spectral energy distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR luminosities and dust temperatures of local galaxies at z<0.1 using AKARI and IRAS data in the field of the Sloan Digital Sky Survey. We compare the L_IR-T_dust relation between the two samples and find that: the median T_dust of Herschel-selected galaxies at z>0.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs.

The reliability of infrared (IR) and ultraviolet (UV) emissions to measure star formation rates in galaxies is investigated for a large sample of galaxies observed with the SPIRE and PACS instruments on Herschel as part of the HerMES project. We build flux-limited 250 micron samples of sources at redshift z<1, cross-matched with the Spitzer/MIPS and GALEX catalogues. About 60 % of the Herschel sources are detected in UV. The total IR luminosities, L_IR, of the sources are estimated using a SED-fitting code that fits to fluxes between 24 and 500 micron. Dust attenuation is discussed on the basis of commonly-used diagnostics: the L_IR/L_UV ratio and the slope, beta, of the UV continuum. A mean dust attenuation A_UV of ~ 3 mag is measured in the samples. L_IR/L_UV is found to correlate with L_IR. Galaxies with L_IR > 10 ^{11} L_sun and 0.5< z<1 exhibit a mean dust attenuation A_UV about 0.7 mag lower than that found for their local counterparts, although with a large dispersion. Our galaxy samples span a large range of beta and L_IR/L_UV values which, for the most part, are distributed between the ranges defined by the relations found locally for starburst and normal star-forming galaxies. As a consequence the recipe commonly applied to local starbursts is found to overestimate the dust attenuation correction in our galaxy sample by a factor ~2-3 .

The HERschel Inventory of The Agents of Galaxy Evolution (HERITAGE) of the Magellanic Clouds will use dust emission to investigate the life cycle of matter in both the Large and Small Magellanic Clouds (LMC and SMC). Using the Herschel Space Observatory's PACS and SPIRE photometry cameras, we imaged a 2x8 square degree strip through the LMC, at a position angle of ~22.5 degrees as part of the science demonstration phase of the Herschel mission. We present the data in all 5 Herschel bands: PACS 100 and 160 {\mu}m and SPIRE 250, 350 and 500 {\mu}m. We present two dust models that both adequately fit the spectral energy distribution for the entire strip and both reveal that the SPIRE 500 {\mu}m emission is in excess of the models by 6 to 17%. The SPIRE emission follows the distribution of the dust mass, which is derived from the model. The PAH-to-dust mass (f_PAH) image of the strip reveals a possible enhancement in the LMC bar in agreement with previous work. We compare the gas mass distribution derived from the HI 21 cm and CO J=1-0 line emission maps to the dust mass map from the models and derive gas-to-dust mass ratios (GDRs). The dust model, which uses the standard graphite and silicate optical properties for Galactic dust, has a very low GDR = 65(+15,-18) making it an unrealistic dust model for the LMC. Our second dust model, which uses amorphous carbon instead of graphite, has a flatter emissivity index in the submillimeter and results in a GDR = 287(+25,-42) that is more consistent with a GDR inferred from extinction.

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.

The Spectral and Photometric Imaging Receiver (SPIRE), is the Herschel Space Observatory`s submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 microns, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 microns (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4' x 8', observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6'. The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the Herschel cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5-2.

The Photodetector Array Camera and Spectrometer (PACS) is one of the three science instruments on ESA's far infrared and submillimetre observatory. It employs two Ge:Ga photoconductor arrays (stressed and unstressed) with 16x25 pixels, each, and two filled silicon bolometer arrays with 16x32 and 32x64 pixels, respectively, to perform integral-field spectroscopy and imaging photometry in the 60-210\mu\ m wavelength regime. In photometry mode, it simultaneously images two bands, 60-85\mu\ m or 85-125\mu\m and 125-210\mu\ m, over a field of view of ~1.75'x3.5', with close to Nyquist beam sampling in each band. In spectroscopy mode, it images a field of 47"x47", resolved into 5x5 pixels, with an instantaneous spectral coverage of ~1500km/s and a spectral resolution of ~175km/s. We summarise the design of the instrument, describe observing modes, calibration, and data analysis methods, and present our current assessment of the in-orbit performance of the instrument based on the Performance Verification tests. PACS is fully operational, and the achieved performance is close to or better than the pre-launch predictions.

We present a power spectrum analysis of the Herschel-SPIRE observations of the Polaris flare, a high Galactic latitude cirrus cloud midway between the diffuse and molecular phases. The SPIRE images of the Polaris flare reveal for the first time the structure of the diffuse interstellar medium down to 0.01 parsec over a 10 square degrees region. These exceptional observations highlight the highly filamentary and clumpy structure of the interstellar medium even in diffuse regions of the map. The power spectrum analysis shows that the structure of the interstellar medium is well described by a single power law with an exponent of -2.7 +- 0.1 at all scales from 30" to 8 degrees. That the power spectrum slope of the dust emission is constant down to the SPIRE angular resolution is an indication that the inertial range of turbulence extends down to the 0.01 pc scale. The power spectrum analysis also allows the identification of a Poissonian component at sub-arcminute scales in agreement with predictions of the cosmic infrared background level at SPIRE wavelengths. Finally, the comparison of the SPIRE and IRAS 100 micron data of the Polaris flare clearly assesses the capability of SPIRE in maping diffuse emission over large areas.

We measure the angular correlation function, w(theta), from 0.5 to 30 arcminutes of detected sources in two wide fields of the Herschel Multi-tiered Extragalactic Survey (HerMES). Our measurements are consistent with the expected clustering shape from a population of sources that trace the dark matter density field, including non-linear clustering at arcminute angular scales arising from multiple sources that occupy the same dark matter halos. By making use of the halo model to connect the spatial clustering of sources to the dark matter halo distribution, we estimate source bias and halo occupation number for dusty sub-mm galaxies at z ~ 2. We find that sub-mm galaxies with 250 micron flux densities above 30 mJy reside in dark matter halos with mass above (5\pm4) x 10^12 M_sun, while (14\pm8)% of such sources appear as satellites in more massive halos.

We present clumps of dust emission from Herschel observations of the Large Magellanic Cloud (LMC) and their physical and statistical properties. We catalog cloud features seen in the dust emission from Herschel observations of the LMC, the Magellanic type irregular galaxy closest to the Milky Way, and compare these features with HI catalogs from the ATCA+Parkes HI survey. Using an automated cloud-finding algorithm, we identify clouds and clumps of dust emission and examine the cumulative mass distribution of the detected dust clouds. The mass of cold dust is determined from physical parameters that we derive by performing spectral energy distribution fits to 250, 350, and 500 micronm emission from SPIRE observations using DUSTY and GRASIL radiative transfer calculation with dust grain size distributions for graphite/silicate in low-metallicity extragalactic environments. The dust cloud mass spectrum follows a power law distribution with an exponent of gamma=-1.8 for clumps larger than 400 solar mass and is similar to the HI mass distribution. This is expected from the theory of ISM structure in the vicinity of star formation.

We use the Science Demonstration Phase data of the Herschel Virgo Cluster Survey to search for dust emission of early-type dwarf galaxies in the central regions of the Virgo Cluster as an alternative way of identifying the interstellar medium.We present the first possible far-infrared detection of cluster early-type dwarf galaxies: VCC781 and VCC951 are detected at the 10 sigma level in the SPIRE 250 micron image. Both detected galaxies have dust masses of the order of 10^5 Msun and average dust temperatures ~20K. The detection rate (less than 1%) is quite high compared to the 1.7% detection rate for Hi emission, considering that dwarfs in the central regions are more Hi deficient. We conclude that the removal of interstellar dust from dwarf galaxies resulting from ram pressure stripping, harassment, or tidal effects must be as efficient as the removal of interstellar gas.

Passive early-type galaxies (ETGs) provide an ideal laboratory for studying the interplay between dust formation around evolved stars and its subsequent destruction in a hot gas. Using Spitzer-IRS and Herschel data we compare the dust production rate in the envelopes of evolved AGB stars with a constraint on the total dust mass. Early-type galaxies which appear to be truly passively evolving are not detected by Herschel. We thus derive a distance independent upper limit to the dust grain survival time in the hostile environment of ETGs of < 46 +/- 25 Myr for amorphous silicate grains. This implies that ETGs which are detected at far-infrared wavelengths have acquired a cool dusty medium via interaction. Given likely time-scales for ram-pressure stripping, this also implies that only galaxies with dust in a cool (atomic) medium can release dust into the intra-cluster medium.

We present a resolved dust analysis of three of the largest angular size spiral galaxies, NGC 4501 and NGC 4567/8, in the Herschel Virgo Cluster Survey (HeViCS) Science Demonstration field. Herschel has unprecedented spatial resolution at far-infrared wavelengths and with the PACS and SPIRE instruments samples both sides of the peak in the far infrared spectral energy distribution (SED).We present maps of dust temperature, dust mass, and gas-to-dust ratio, produced by fitting modified black bodies to the SED for each pixel. We find that the distribution of dust temperature in both systems is in the range ~19 - 22 K and peaks away from the centres of the galaxies. The distribution of dust mass in both systems is symmetrical and exhibits a single peak coincident with the galaxy centres. This Letter provides a first insight into the future analysis possible with a large sample of resolved galaxies to be observed by Herschel.

We present the dust properties of a small sample of Virgo cluster dwarf galaxies drawn from the science demonstration phase data set of the Herschel Virgo Cluster Survey. These galaxies have low metallicities (7.8 < 12 + log(O/H) < 8.3) and star-formation rates < 10^{-1} M_{sun}/yr. We measure the spectral energy distribution (SED) from 100 to 500 um and derive dust temperatures and dust masses. The SEDs are fitted by a cool component of temperature T < 20 K, implying dust masses around 10^{5} M_{sun} and dust-to-gas ratios D within the range 10^{-3}-10^{-2}. The completion of the full survey will yield a larger set of galaxies, which will provide more stringent constraints on the dust content of star-forming dwarf galaxies.

The origin of the far-infrared emission from the nearby radio galaxy M87 remains a matter of debate. Some studies find evidence of a far-infrared excess due to thermal dust emission, whereas others propose that the far-infrared emission can be explained by synchrotron emission without the need for an additional dust emission component. We present Herschel PACS and SPIRE observations of M87, taken as part of the science demonstration phase observations of the Herschel Virgo Cluster Survey. We compare these data with a synchrotron model based on mid-infrared, far-infrared, submm and radio data from the literature to investigate the origin of the far-infrared emission. Both the integrated SED and the Herschel surface brightness maps are adequately explained by synchrotron emission. At odds with previous claims, we find no evidence of a diffuse dust component in M87, which is not unexpected in the harsh X-ray environment of this radio galaxy sitting at the core of the Virgo Cluster.

Context: Sh2-104 is a Galactic H ii region with a bubble morphology, detected at optical and radio wavelengths. It is considered the first observational confirmation of the collect-and-collapse model of triggered star-formation. Aims: We aim to analyze the dust and gas properties of the Sh2-104 region to better constrain its effect on local future generations of stars. In addition, we investigate the relationship between the dust emissivity index {\beta} and the dust temperature, T_dust. Methods: Using Herschel PACS and SPIRE images at 100, 160, 250, 350 and 500 {\mu}m we determine T_dust and {\beta} throughout Sh2-104, fitting the spectral energy distributions (SEDs) obtained from aperture photometry. With the SPIRE Fourier transform spectrometer (FTS) we obtained spectra at different positions in the Sh2-104 region. We detect J-ladders of CO and 13CO, with which we derive the gas temperature and column density. We also detect proxies of ionizing flux as the [NII] 3P1-3P0 and [CI] 3P2-3P1 transitions. Results: We find an average value of {\beta} ~ 1.5 throughout Sh2-104, as well as a T dust difference between the photodissociation region (PDR, ~ 25 K) and the interior (~ 40 K) of the bubble. We recover the anti-correlation between {\beta} and dust temperature reported numerous times in the literature. The relative isotopologue abundances of CO appear to be enhanced above the standard ISM values, but the obtained value is very preliminary and is still affected by large uncertainties.

The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies.

We present far-infrared spectra and maps of the DR21 molecular cloud core between 196 and 671 microns, using the Herschel-SPIRE spectrometer. Nineteen molecular lines originating from CO, 13CO, HCO+ and H2O, plus lines of [N II] and [CI] were recorded, including several transitions not previously detected. The CO lines are excited in warm gas with Tkin ~ 125 K and nH2 ~ 7 x 10^4 cm-3, CO column density N(CO) ~ 3.5 x 10^18 cm^-2 and a filling factor of ~ 12%, and appear to trace gas associated with an outflow. The rotational temperature analysis incorporating observations from ground-based telescopes reveals an additional lower excitation CO compoment which has a temperature ~ 78 K and N(CO) ~ 4.5 x 10^21 cm^-2. Astronomy & Astrophysics HERSCHEL special Issue, in press.

We describe the Herschel Virgo Cluster Survey (HeViCS) and the first data obtained as part of the Science Demonstration Phase (SDP). The data cover a central 4x4 sq deg region of the cluster. We use SPIRE and PACS photometry data to produce 100, 160, 250, 350 and 500 micron luminosity functions (LFs) for optically bright galaxies that are selected at 500 micron and detected in all bands. We compare these LFs with those previously derived using IRAS, BLAST and Herschel-ATLAS data. The Virgo Cluster LFs do not have the large numbers of faint galaxies or examples of very luminous galaxies seen previously in surveys covering less dense environments.

By combining Herschel-SPIRE observations obtained as part of the Herschel Virgo Cluster Survey with 21 cm HI data from the literature, we investigate the role of the cluster environment on the dust content of Virgo spiral galaxies.We show for the first time that the extent of the dust disk is significantly reduced in HI-deficient galaxies, following remarkably well the observed 'truncation' of the HI disk. The ratio of the submillimetre-to- optical diameter correlates with the HI-deficiency, suggesting that the cluster environment is able to strip dust as well as gas. These results provide important insights not only into the evolution of cluster galaxies but also into the metal enrichment of the intra-cluster medium.

The properties of the dust grains (e.g., temperature and mass) can be derived from fitting far-IR SEDs (>100 micron). Only with SPIRE on Herschel has it been possible to get high spatial resolution at 200 to 500 micron that is beyond the peak (~160 micron) of dust emission in most galaxies. We investigate the differences in the fitted dust temperatures and masses determined using only <200 micron data and then also including >200 micron data (new SPIRE observations) to determine how important having >200 micron data is for deriving these dust properties. We fit the 100 to 350 micron observations of the Large Magellanic Cloud (LMC) point-by-point with a model that consists of a single temperature and fixed emissivity law. The data used are existing observations at 100 and 160 micron (from IRAS and Spitzer) and new SPIRE observations of 1/4 of the LMC observed for the HERITAGE Key Project as part of the Herschel Science Demonstration phase. The dust temperatures and masses computed using only 100 and 160 micron data can differ by up to 10% and 36%, respectively, from those that also include the SPIRE 250 & 350 micron data. We find that an emissivity law proportional to lambda^-1.5 minimizes the 100-350 micron fractional residuals. We find that the emission at 500 micron is ~10% higher than expected from extrapolating the fits made at shorter wavelengths. We find the fractional 500 micron excess is weakly anti-correlated with MIPS 24 micron flux and the total gas surface density. This argues against a flux calibration error as the origin of the 500 micron excess. Our results do not allow us to distinguish between a systematic variation in the wavelength dependent emissivity law or a population of very cold dust only detectable at lambda > 500 micron for the origin of the 500 micron excess.