Search

We present JWST NIRCam (F356W and F444W filters) and MIRI (F770W) images and NIRSpec- IFU spectroscopy of the young supernova remnant Cassiopeia A (Cas A). We obtained the data as part of a JWST survey of Cas A. The NIRCam and MIRI images map the spatial distributions of synchrotron radiation, Ar-rich ejecta, and CO on both large and small scales, revealing remarkably complex structures. The CO emission is stronger at the outer layers than the Ar ejecta, which indicates the reformation of CO molecules behind the reverse shock. NIRSpec-IFU spectra (3 - 5.5 microns) were obtained toward two representative knots in the NE and S fields. Both regions are dominated by the bright fundamental rovibrational band of CO in the two R and P branches, with strong [Ar VI] and relatively weaker, variable strength ejecta lines of [Si IX], [Ca IV], [Ca V] and [Mg IV]. The NIRSpec-IFU data resolve individual ejecta knots and filaments spatially and in velocity space. The fundamental CO band in the JWST spectra reveals unique shapes of CO, showing a few tens of sinusoidal patterns of rovibrational lines with pseudo-continuum underneath, which is attributed to the high-velocity widths of CO lines. The CO also shows high J lines at different vibrational transitions. Our results with LTE modeling of CO emission indicate a temperature of 1080 K and provide unique insight into the correlations between dust, molecules, and highly ionized ejecta in supernovae, and have strong ramifications for modeling dust formation that is led by CO cooling in the early Universe., Comment: accepted for the ApJ letter (17 pages and 10 figures)

JWST/NIRCam obtained high angular-resolution (0.05-0.1''), deep near-infrared 1--5 micron imaging of Supernova (SN) 1987A taken 35 years after the explosion. In the NIRCam images, we identify: 1) faint H2 crescents, which are emissions located between the ejecta and the equatorial ring, 2) a bar, which is a substructure of the ejecta, and 3) the bright 3-5 micron continuum emission exterior to the equatorial ring. The emission of the remnant in the NIRCam 1-2.3 micron images is mostly due to line emission, which is mostly emitted in the ejecta and in the hot spots within the equatorial ring. In contrast, the NIRCam 3-5 micron images are dominated by continuum emission. In the ejecta, the continuum is due to dust, obscuring the centre of the ejecta. In contrast, in the ring and exterior to the ring, synchrotron emission contributes a substantial fraction to the continuum. Dust emission contributes to the continuum at outer spots and diffuse emission exterior to the ring, but little within the ring. This shows that dust cooling and destruction time scales are shorter than the synchrotron cooling time scale, and the time scale of hydrogen recombination in the ring is even longer than the synchrotron cooling time scale. With the advent of high sensitivity and high angular resolution images provided by JWST/NIRCam, our observations of SN 1987A demonstrate that NIRCam opens up a window to study particle-acceleration and shock physics in unprecedented details, probed by near-infrared synchrotron emission, building a precise picture of how a SN evolves., Comment: Accepted for publication in MNRAS; 18 pages

Shockwaves driven by supernovae both destroy dust and reprocess the surviving grains, greatly affecting the resulting dust properties of the interstellar medium (ISM). While these processes have been extensively studied theoretically, observational constraints are limited. We use physically-motivated models of dust emission to fit the infrared (IR) spectral energy distributions of seven Galactic supernova remnants, allowing us to determine the distribution of dust mass between diffuse and dense gas phases, and between large and small grain sizes. We find that the dense ($\sim 10^3 \,{\rm cm}^{-3}$), relatively cool ($\sim 10^3 \, {\rm K}$) gas phase contains $>90\%$ of the dust mass, making the warm dust located in the X-ray emitting plasma ($\sim 1 \,{\rm cm}^{-3}$/$10^6 \, {\rm K}$) a negligible fraction of the total, despite dominating the mid-IR emission. The ratio of small ($\lesssim 10 \, {\rm nm}$) to large ($\gtrsim 0.1 \, {\rm \mu m}$) grains in the cold component is consistent with that in the ISM, and possibly even higher, whereas the hot phase is almost entirely devoid of small grains. This suggests that grain shattering, which processes large grains into smaller ones, is ineffective in the low-density gas, contrary to model predictions. Single-phase models of dust destruction in the ISM, which do not account for the existence of the cold swept-up material containing most of the dust mass, are likely to greatly overestimate the rate of dust destruction by supernovae., Comment: 13 pages, 12 figures. MNRAS accepted

We present a self-consistent cross-calibration of the three main molecular gas mass tracers in galaxies, the $\rm ^{12}CO$(1-0), [CI]($^3P_1$-$^3P_0$) lines, and the submm dust continuum emission, using a sample of 407 galaxies, ranging from local disks to submillimetre-selected galaxies (SMGs) up to $z \approx 6$. A Bayesian method is used to produce galaxy-scale universal calibrations of these molecular gas indicators, that hold over 3-4 orders of magnitude in infrared luminosity, $L_{\rm IR}$. Regarding the dust continuum, we use a mass-weighted dust temperature, $T_{\rm mw}$, determined using new empirical relations between temperature and luminosity. We find the average $L/M_{\rm mol}$ gas mass conversion factors to be $\alpha_{850}= 6.9\times10^{12}\,\rm W\,Hz^{-1}\,M_{\odot}^{-1}$, $\alpha_{\rm CO} = \rm 4\,M_{\odot} (K\,km\,s^{-1}\,pc^2)^{-1}$ and $\alpha_{\rm CI} = \rm 17.0 \,M_{\odot} (K\,km\,s^{-1}\,pc^2)^{-1}$, based on the assumption that the mean dust properties of the sample ($\kappa_H$ = gas-to-dust ratio/dust emissivity) will be similar to those of local metal rich galaxies and the MW. The tracer with the least intrinsic scatter is [CI](1-0), while CO(1-0) has the highest. The conversion factors show a weak but significant correlation with $L_{\rm IR}$. Assuming dust properties typical of metal-rich galaxies, we infer a neutral carbon abundance $X_{\rm CI} = [C^0/\rm mol]=1.6\times 10^{-5}$, similar to that in the MW. We find no evidence for bimodality of $\alpha_{\rm CO}$ between main-sequence (MS) galaxies and those with extreme star-formation intensity, i.e. ULIRGs and SMGs. The means of the three conversion factors are found to be similar between MS galaxies and ULIRGs/SMGs, to within 10-20%. We show that for metal-rich galaxies, near-universal average values for $\alpha_{\rm CO}$, $X_{\rm CI}$ and $\kappa_H$ are adequate for global molecular gas estimates., Comment: 42 pages, accepted for publication in MNRAS. Abstract abridged. Data tables will be available on CDS imminently. Minor changes from proof stage including swapping order of Tables 8 and 9 and small corrections to the bibliography

We have carried out the first spatially-resolved investigation of the multi-phase interstellar medium (ISM) at high redshift, using the z=4.24 strongly-lensed sub-millimetre galaxy H-ATLASJ142413.9+022303 (ID141). We present high-resolution (down to ~350 pc) ALMA observations in dust continuum emission and in the CO(7-6), H_2O (2_{1,1} - 2_{0,2}), CI(1-0) and CI(2-1) lines, the latter two allowing us to spatially resolve the cool phase of the ISM for the first time. Our modelling of the kinematics reveals that the system appears to be dominated by a rotationally-supported gas disk with evidence of a nearby perturber. We find that the CI(1-0) line has a very different distribution to the other lines, showing the existence of a reservoir of cool gas that might have been missed in studies of other galaxies. We have estimated the mass of the ISM using four different tracers, always obtaining an estimate in the range (3.2-3.8) x 10^{11} M_sol, significantly higher than our dynamical mass estimate of (0.8-1.3) x 10^{11} M_sol. We suggest that this conflict and other similar conflicts reported in the literature is because the gas-to-tracer ratios are ~4 times lower than the Galactic values used to calibrate the ISM in high-redshift galaxies. We demonstrate that this could result from a top-heavy initial mass function and strong chemical evolution. Using a variety of quantitative indicators, we show that, extreme though it is at z=4.24, ID141 will likely join the population of quiescent galaxies that appears in the Universe at z~3., Comment: This article has been accepted for publication in MNRAS, published by Oxford University Press on behalf of the Royal Astronomical Society. 24 pages

We present ALMA observations of a small but statistically complete sample of twelve 250 micron selected galaxies at $z=0.35$ designed to measure their dust submillimeter continuum emission as well as their CO(1-0) and atomic carbon [CI](3P1-3P0) spectral lines. This is the first sample of galaxies with global measures of all three $H_2$-mass tracers and which show star formation rates (4-26 Msun yr$^{-1}$) and infra-red luminosities ($1-6\times10^{11}$ Lsun) typical of star forming galaxies in their era. We find a surprising diversity of morphology and kinematic structure; one-third of the sample have evidence for interaction with nearby smaller galaxies, several sources have disjoint dust and gas morphology. Moreover two galaxies have very high $L_{CI}/L_{CO}$ ratios for their global molecular gas reservoirs; if confirmed, such extreme intensity ratios in a sample of dust selected, massive star forming galaxies presents a challenge to our understanding of ISM. Finally, we use the emission of the three molecular gas tracers, to determine the carbon abundance, $X_{ci}$, and CO-$\rm{H_2}$ conversion $\alpha_{co}$ in our sample, using a weak prior that the gas-to-dust ratio is similar to that of the Milky Way for these massive and metal rich galaxies. Using a likelihood method which simultaneously uses all three gas tracer measurements, we find mean values and errors on the mean of $\alpha_{co}=3.0\pm0.5\,\rm{Msun\,(K\,kms^{-1}\,pc^2)^{-1}}$ and $X_{ci}=1.6\pm0.1\times 10^{-5}$ (or $\alpha_{ci}=18.8\,K kms^{-1}\,pc^2 (Msun)^{-1}$) and $\delta_{GDR}=128\pm16$ (or $\alpha_{850}=5.9\times10^{12}\,\rm{W\,Hz^{-1}\, Msun^{-1}}$), where our starting assumption is that these metal rich galaxies have an average gas-to-dust ratio similar to that of the Milky Way centered on $\delta_{GDR}=135$., Comment: 34 pages, published in MNRAS

We build a rigorous statistical framework to provide constraints on the chemical and dust evolution parameters for nearby late-type galaxies with a wide range of gas fractions ($3\%

Dust destruction by supernovae is one of the main processes removing dust from the interstellar medium (ISM). Estimates of the efficiency of this process, both theoretical and observational, typically assume a shock propagating into a homogeneous medium, whereas the ISM possesses significant substructure in reality. We self-consistently model the dust and gas properties of the shocked ISM in three supernova remnants (SNRs), using X-ray and infrared (IR) data combined with corresponding emission models. Collisional heating by gas with properties derived from X-ray observations produces dust temperatures too high to fit the far-IR fluxes from each SNR. An additional colder dust component is required, which has a minimum mass several orders of magnitude larger than that of the warm dust heated by the X-ray emitting gas. Dust-to-gas mass ratios indicate that the majority of the dust in the X-ray emitting material has been destroyed, while the fraction of surviving dust in the cold component is plausibly close to unity. As the cold component makes up virtually all the total dust mass, destruction timescales based on homogeneous models, which cannot account for multiple phases of shocked gas and dust, may be significantly overestimating actual dust destruction efficiencies, and subsequently underestimating grain lifetimes., Comment: 12 pages, 9 figures. MNRAS accepted 31/10/2020

We investigate the evolution of the gas mass fraction for galaxies in the COSMOS field using submillimetre emission from dust at 850$\mu$m. We use stacking methodologies on the 850$\mu$m S2COSMOS map to derive the gas mass fraction of galaxies out to high redshifts, 0 <= $z$ <= 5, for galaxies with stellar masses of $10^{9.5} < M_* (\rm M_{\odot}) < 10^{11.75}$. In comparison to previous literature studies we extend to higher redshifts, include more normal star-forming galaxies (on the main sequence), and also investigate the evolution of the gas mass fraction split by star-forming and passive galaxy populations. We find our stacking results broadly agree with scaling relations in the literature. We find tentative evidence for a peak in the gas mass fraction of galaxies at around $z$ ~ 2.5-3, just before the peak of the star formation history of the Universe. We find that passive galaxies are particularly devoid of gas, compared to the star-forming population. We find that even at high redshifts, high stellar mass galaxies still contain significant amounts of gas., Comment: 17 pages (+ 9 pages of appendices), 8 figures (+ 7 figures in appendices), accepted for publication in MNRAS

The Herschel-ATLAS unbiased survey of cold dust in the local Universe is dominated by a surprising population of very blue (FUV-K < 3.5), dust-rich galaxies with high gas fractions (f_HI = M_HI/(M*+M_HI)>0.5)). Dubbed `Blue and Dusty Gas Rich Sources' (BADGRS) they have cold diffuse dust temperatures, and the highest dust-to-stellar mass ratios of any galaxies in the local Universe. Here, we explore the molecular ISM in a representative sample of BADGRS, using very deep CO(J_up=1,2,3) observations across the central and outer disk regions. We find very low CO brightnesses (Tp=15-30 mK), despite the bright far-infrared emission and metallicities in the range 0.5

We have modelled the near-infrared to radio images of the Crab Nebula with a Bayesian SED model to simultaneously fit its synchrotron, interstellar and supernova dust emission. We infer an interstellar dust extinction map with an average $A_{\text{V}}$=1.08$\pm$0.38 mag, consistent with a small contribution (<22%) to the Crab's overall infrared emission. The Crab's supernova dust mass is estimated to be between 0.032 and 0.049 M$_{\odot}$ (for amorphous carbon grains) with an average dust temperature $T_{\text{dust}}$=41$\pm$3K, corresponding to a dust condensation efficiency of 8-12%. This revised dust mass is up to an order of magnitude lower than some previous estimates, which can be attributed to our different interstellar dust corrections, lower SPIRE flux densities, and higher dust temperature than were used in previous studies. The dust within the Crab is predominantly found in dense filaments south of the pulsar, with an average V-band dust extinction of $A_{\text{V}}$=0.20-0.39 mag, consistent with recent optical dust extinction studies. The modelled synchrotron power-law spectrum is consistent with a radio spectral index $\alpha_{\text{radio}}$=0.297$\pm$0.009 and an infrared spectral index $\alpha_{\text{IR}}$=0.429$\pm$0.021. We have identified a millimetre excess emission in the Crab's central regions, and argue that it most likely results from two distinct populations of synchrotron emitting particles. We conclude that the Crab's efficient dust condensation (8-12%) provides further evidence for a scenario where supernovae can provide substantial contributions to the interstellar dust budgets in galaxies., Comment: 33 pages, 17 figures, Manuscript accepted for publication in MNRAS. Accepted on 31/05/2019. Deposited on 31/05/2019

The {\it Herschel} Astrophysical Terahertz Large Area Survey (H-ATLAS) is a survey of 660 deg$^2$ with the PACS and SPIRE cameras in five photometric bands: 100, 160, 250, 350 and 500\mic. This is the second of three papers describing the data release for the large fields at the south and north Galactic poles (NGP and SGP). In this paper we describe the catalogues of far-infrared and submillimetre sources for the NGP and SGP, which cover 177 deg$^2$ and 303 deg$^2$, respectively. The catalogues contain 153,367 sources for the NGP field and 193,527 sources for the SGP field detected at more than 4$\sigma$ significance in any of the 250, 350 or 500\mic\ bands. The source detection is based on the 250\mic\ map, and we present photometry in all five bands for each source, including aperture photometry for sources known to be extended. The rms positional accuracy for the faintest sources is about 2.4 arc seconds in both right ascension and declination. We present a statistical analysis of the catalogues and discuss the practical issues -- completeness, reliability, flux boosting, accuracy of positions, accuracy of flux measurements -- necessary to use the catalogues for astronomical projects.

We present the dust mass function (DMF) of 15,750 galaxies with redshift $z< 0.1$, drawn from the overlapping area of the GAMA and {\it H-}ATLAS surveys. The DMF is derived using the density corrected $V_{\rm max}$ method, where we estimate $V_{\rm max}$ using: (i) the normal photometric selection limit ($pV_{\rm max}$) and (ii) a bivariate brightness distribution (BBD) technique, which accounts for two selection effects. We fit the data with a Schechter function, and find $M^{*}=(4.65\pm0.18)\times 10^{7}\,h^2_{70}\, M_{\odot}$, $\alpha=(1.22\pm 0.01)$, $\phi^{*}=(6.26\pm 0.28)\times 10^{-3}\,h^3_{70}\,\rm Mpc^{-3}\,dex^{-1}$. The resulting dust mass density parameter integrated down to $10^4\,M_{\odot}$ is $\Omega_{\rm d}=(1.11 \pm0.02)\times 10^{-6}$ which implies the mass fraction of baryons in dust is $f_{m_b}=(2.40\pm0.04)\times 10^{-5}$; cosmic variance adds an extra 7-17\,per\,cent uncertainty to the quoted statistical errors. Our measurements have fewer galaxies with high dust mass than predicted by semi-analytic models. This is because the models include too much dust in high stellar mass galaxies. Conversely, our measurements find more galaxies with high dust mass than predicted by hydrodynamical cosmological simulations. This is likely to be from the long timescales for grain growth assumed in the models. We calculate DMFs split by galaxy type and find dust mass densities of $\Omega_{\rm d}=(0.88\pm0.03)\times 10^{-6}$ and $\Omega_{\rm d}=(0.060\pm0.005)\times 10^{-6}$ for late-types and early-types respectively. Comparing to the equivalent galaxy stellar mass functions (GSMF) we find that the DMF for late-types is well matched by the GMSF scaled by $(8.07\pm0.35) \times 10^{-4}$., Comment: 25 Pages, 18 Figures. Submitted December 2017

The process of planet formation offers a rich source of dust production via grain growth in protostellar discs, and via grinding of larger bodies in debris disc systems. Chemical evolution models, designed to follow the build up of metals and dust in galaxies, do not currently account for planet formation. We consider the possibility that the apparent under-prediction of dust mass in high redshift galaxies by chemical evolution models could be in part, due to these models neglecting this process, specifically due to their assumption that a large fraction of the dust mass is removed from the interstellar medium during star formation (so-called astration). By adding a planet formation phase into galaxy chemical evolution, we demonstrate that the dust budget crisis can be partially ameliorated by a factor of 1.3-1.5 only if a) circumstellar discs prevent a large fraction of the dust mass entering the star during its birth, and b) that dust mass is preferentially liberated via jets, winds and outflows rather than accreted into planetary-mass bodies., Comment: 9 pages, 11 figures, accepted for publication in MNRAS

We present infrared (IR) and submillimeter observations of the Crab-like supernova remnant (SNR) G54.1+0.3 including 350 micron (SHARC-II), 870 micron (LABOCA), 70, 100, 160, 250, 350, 500 micron (Herschel) and 3-40 micron (Spitzer). We detect dust features at 9, 11 and 21 micron and a long wavelength continuum dust component. The 21 micron dust coincides with [Ar II] ejecta emission, and the feature is remarkably similar to that in Cas A. The IRAC 8 micron image including Ar ejecta is distributed in a shell-like morphology which is coincident with dust features, suggesting that dust has formed in the ejecta. We create a cold dust map that shows excess emission in the northwestern shell. We fit the spectral energy distribution of the SNR using the continuous distributions of ellipsoidal (CDE) grain model of pre-solar grain SiO2 that reproduces the 21 and 9 micron dust features and discuss grains of SiC and PAH that may be responsible for the 10-13 micron dust features. To reproduce the long-wavelength continuum, we explore models consisting of different grains including Mg2SiO4, MgSiO3, Al2O3, FeS, carbon, and Fe3O4. We tested a model with a temperature-dependent silicate absorption coefficient. We detect cold dust (27-44 K) in the remnant, making this the fourth such SNR with freshly-formed dust. The total dust mass in the SNR ranges from 0.08-0.9 Msun depending on the grain composition, which is comparable to predicted masses from theoretical models. Our estimated dust masses are consistent with the idea that SNe are a significant source of dust in the early Universe., Comment: MNRAS: accepted on June 28, 2018 and published on July 4, 2018

Most massive stars end their lives in core-collapse supernova explosions and enrich the interstellar medium with explosively nucleosynthesized elements. Following core collapse, the explosion is subject to instabilities as the shock propagates outwards through the progenitor star. Observations of the composition and structure of the innermost regions of a core-collapse supernova provide a direct probe of the instabilities and nucleosynthetic products. SN 1987A in the Large Magellanic Cloud (LMC) is one of very few supernovae for which the inner ejecta can be spatially resolved but are not yet strongly affected by interaction with the surroundings. Our observations of SN 1987A with the Atacama Large Millimeter/submillimeter Array (ALMA) are of the highest resolution to date and reveal the detailed morphology of cold molecular gas in the innermost regions of the remnant. The 3D distributions of carbon and silicon monoxide (CO and SiO) emission differ, but both have a central deficit, or torus-like distribution, possibly a result of radioactive heating during the first weeks ("nickel heating"). The size scales of the clumpy distribution are compared quantitatively to models, demonstrating how progenitor and explosion physics can be constrained., Comment: ApJL accepted, 9 pages, 5 figures

We combine samples of nearby galaxies with Herschel photometry selected on their dust, metal, HI, and stellar mass content, and compare these to chemical evolution models in order to discriminate between different dust sources. In a companion paper, we used a HI-selected sample of nearby galaxies to reveal a sub-sample of very gas rich (gas fraction > 80 per cent) sources with dust masses significantly below predictions from simple chemical evolution models, and well below $M_d/M_*$ and $M_d/M_{gas}$ scaling relations seen in dust and stellar-selected samples of local galaxies. We use a chemical evolution model to explain these dust-poor, but gas-rich, sources as well as the observed star formation rates (SFRs) and dust-to-gas ratios. We find that (i) a delayed star formation history is required to model the observed SFRs; (ii) inflows and outflows are required to model the observed metallicities at low gas fractions; (iii) a reduced contribution of dust from supernovae (SNe) is needed to explain the dust-poor sources with high gas fractions. These dust-poor, low stellar mass galaxies require a typical core-collapse SN to produce 0.01 - 0.16 $M_{\odot}$ of dust. To match the observed dust masses at lower gas fractions, significant grain growth is required to counteract the reduced contribution from dust in SNe and dust destruction from SN shocks. These findings are statistically robust, though due to intrinsic scatter it is not always possible to find one single model that successfully describes all the data. We also show that the dust-to-metals ratio decreases towards lower metallicity., Comment: 15 pages (+10 pages appendix), 8 figures, Accepted for publication in MNRAS

We report the first molecular line survey of Supernova 1987A in the millimetre wavelength range. In the ALMA 210--300 and 340--360 GHz spectra, we detected cold (20--170 K) CO, 28SiO, HCO+ and SO, with weaker lines of 29SiO from ejecta. This is the first identification of HCO+ and SO in a young supernova remnant. We find a dip in the J=6--5 and 5--4 SiO line profiles, suggesting that the ejecta morphology is likely elongated. The difference of the CO and SiO line profiles is consistent with hydrodynamic simulations, which show that Rayleigh-Taylor instabilities cause mixing of gas, with heavier elements much more disturbed, making more elongated structure. We obtained isotopologue ratios of 28SiO/29SiO>13, 28SiO/30SiO>14, and 12CO/13CO>21, with the most likely limits of 28SiO/29SiO>128, 28SiO/30SiO>189. Low 29Si and 30Si abundances in SN 1987A are consistent with nucleosynthesis models that show inefficient formation of neutron-rich isotopes in a low metallicity environment, such as the Large Magellanic Cloud. The deduced large mass of HCO+ (~5x10^-6 Msun) and small SiS mass (<6x10^-5 Msun) might be explained by some mixing of elements immediately after the explosion. The mixing might have caused some hydrogen from the envelope to sink into carbon and oxygen-rich zones after the explosion, enabling the formation of a substantial mass of HCO+. Oxygen atoms may have penetrated into silicon and sulphur zones, suppressing formation of SiS. Our ALMA observations open up a new window to investigate chemistry, dynamics and explosive-nucleosynthesis in supernovae., Comment: Accepted by MNRAS, 16 pages

Theoretical models predict that core-collapse supernovae (CCSNe) can be efficient dust producers (0.1-1.0 Msun), potentially accounting for most of the dust production in the early Universe. Observational evidence for this dust production efficiency is however currently limited to only a few CCSN remnants (e.g., SN1987A, Crab Nebula). In this paper, we revisit the dust mass produced in Cassiopeia A (Cas A), a ~330-year old O-rich Galactic supernova remnant (SNR) embedded in a dense interstellar foreground and background. We present the first spatially resolved analysis of Cas A based on Spitzer and Herschel infrared and submillimetre data at a common resolution of ~0.6 arcmin for this 5 arcmin diameter remnant following a careful removal of contaminating line emission and synchrotron radiation. We fit the dust continuum from 17 to 500 micron with a four-component interstellar medium (ISM) and supernova (SN) dust model. We find a concentration of cold dust in the unshocked ejecta of Cas A and derive a mass of 0.3-0.5 Msun of silicate grains freshly produced in the SNR, with a lower limit of >=0.1-0.2 Msun. For a mixture of 50% of silicate-type grains and 50% of carbonaceous grains, we derive a total SN dust mass between 0.4 Msun and 0.6 Msun. These dust mass estimates are higher than from most previous studies of Cas A and support the scenario of supernova dominated dust production at high redshifts. We furthermore derive an interstellar extinction map for the field around Cas A which towards Cas A gives average values of A_V=6-8 mag, up to a maximum of A_V=15 mag., Comment: 37 pages, 24 figures, Manuscript published in MNRAS, including minor corrections. Accepted on 01/11/2016. Deposited on 08/01/2017

We present a study of the dust, stars and atomic gas (HI) in an HI-selected sample of local galaxies (z<0.035) in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) fields. This HI-selected sample reveals a population of very high gas fraction (>80 per cent), low stellar mass sources that appear to be in the earliest stages of their evolution. We compare this sample with dust and stellar mass selected samples to study the dust and gas scaling relations over a wide range of gas fraction (proxy for evolutionary state of a galaxy). The most robust scaling relations for gas and dust are those linked to NUV-r (SSFR) and gas fraction, these do not depend on sample selection or environment. At the highest gas fractions, our additional sample shows the dust content is well below expectations from extrapolating scaling relations for more evolved sources, and dust is not a good tracer of the gas content. The specific dust mass for local galaxies peaks at a gas fraction of ~75 per cent. The atomic gas depletion time is also longer for high gas fraction galaxies, opposite to the trend found for molecular gas depletion timescale. We link this trend to the changing efficiency of conversion of HI to H2 as galaxies increase in stellar mass surface density as they evolve. Finally, we show that galaxies start out barely obscured and increase in obscuration as they evolve, yet there is no clear and simple link between obscuration and global galaxy properties., Comment: 16 pages (+11 pages appendix), 12 figures, Published in 2017 in MNRAS, 464, 4680; This version 2 has corrected a small error in the table heading of Table 1, but is otherwise unchanged

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

In our series of papers presenting the Herschel imaging of evolved planetary nebulae, we present images of the dust distribution in the Helix nebula (NGC 7293). Images at 70, 160, 250, 350, and 500 micron were obtained with the PACS and SPIRE instruments on board the Herschel satellite. The broadband maps show the dust distribution over the main Helix nebula to be clumpy and predominantly present in the barrel wall. We determined the spectral energy distribution of the main nebula in a consistent way using Herschel, IRAS, and Planck flux values. The emissivity index of 0.99 +/- 0.09, in combination with the carbon rich molecular chemistry of the nebula, indicates that the dust consists mainly of amorphous carbon. The dust excess emission from the central star disk is detected at 70 micron and the flux measurement agree with previous measurement. We present the temperature and dust column density maps. The total dust mass across the Helix nebula (without its halo) is determined to be 0.0035 solar mass at a distance of 216 pc. The temperature map shows dust temperatures between 22 and 42 K, which is similar to the kinetic temperature of the molecular gas, strengthening the fact that the dust and gas co-exist in high density clumps. Archived images are used to compare the location of the dust emission in the far infrared (Herschel) with the ionized (GALEX, Hbeta) and molecular hydrogen component. The different emission components are consistent with the Helix consisting of a thick walled barrel-like structure inclined to the line of sight. The radiation field decreases rapidly through the barrel wall., Comment: 8 pages, 9 figures, revised version A&A in press

In this paper we consider the implications of the distributions of dust and metals in the disc of M31. We derive mean radial dust distributions using a dust map created from Herschel images of M31 sampling the entire far-infrared (FIR) peak. Modified blackbodies are fit to approximately 4000 pixels with a varying, as well as a fixed, dust emissivity index (beta). An overall metal distribution is also derived using data collected from the literature. We use a simple analytical model of the evolution of the dust in a galaxy with dust contributed by stellar sources and interstellar grain growth, and fit this model to the radial dust-to-metals distribution across the galaxy. Our analysis shows that the dust-to-gas gradient in M31 is steeper than the metallicity gradient, suggesting interstellar dust growth is (or has been) important in M31. We argue that M31 helps build a case for cosmic dust in galaxies being the result of substantial interstellar grain growth, while the net dust production from stars may be limited. We note, however, that the efficiency of dust production in stars, e.g., in supernovae (SNe) ejecta and/or stellar atmospheres, and grain destruction in the interstellar medium (ISM) may be degenerate in our simple model. We can conclude that interstellar grain growth by accretion is likely at least as important as stellar dust production channels in building the cosmic dust component in M31., Comment: 12 pages, 7 figures. Published in MNRAS 444, 797. This version is updated to match the published version

We present a comparison of the physical properties of a rest-frame $250\mu$m selected sample of massive, dusty galaxies from $01$ SMGs have an average SFR of $390^{+80}_{-70}\,$M$_\odot$yr$^{-1}$ which is 120 times that of the low-redshift sample matched in stellar mass to the SMGs (SFR$=3.3\pm{0.2}$ M$_\odot$yr$^{-1}$). The SMGs harbour a substantial mass of dust ($1.2^{+0.3}_{-0.2}\times{10}^9\,$M$_\odot$), compared to $(1.6\pm0.1)\times{10}^8\,$M$_\odot$ for low-redshift dusty galaxies. At low redshifts the dust luminosity is dominated by the diffuse ISM, whereas a large fraction of the dust luminosity in SMGs originates from star-forming regions. At the same dust mass SMGs are offset towards a higher SFR compared to the low-redshift H-ATLAS galaxies. This is not only due to the higher gas fraction in SMGs but also because they are undergoing a more efficient mode of star formation, which is consistent with their bursty star-formation histories. The offset in SFR between SMGs and low-redshift galaxies is similar to that found in CO studies, suggesting that dust mass is as good a tracer of molecular gas as CO., Comment: 17 pages, 8 pages in Appendix, 12 figures, 2 tables. Accepted for publication in MNRAS. Minor typos corrected to match published version

We apply a chemical evolution model to investigate the sources and evolution of dust in a sample of 26 high-redshift ($z>1$) submillimetre galaxies (SMGs) from the literature, with complete photometry from ultraviolet to the submillimetre. We show that dust produced only by low-intermediate mass stars falls a factor 240 short of the observed dust masses of SMGs, the well-known `dust-budget crisis'. Adding an extra source of dust from supernovae can account for the dust mass in 19 per cent of the SMG sample. Even after accounting for dust produced by supernovae the remaining deficit in the dust mass budget provides support for higher supernova yields, substantial grain growth in the interstellar medium or a top-heavy IMF. Including efficient destruction of dust by supernova shocks increases the tension between our model and observed SMG dust masses. The models which best reproduce the physical properties of SMGs have a rapid build-up of dust from both stellar and interstellar sources and minimal dust destruction. Alternatively, invoking a top-heavy IMF or significant changes in the dust grain properties can solve the dust budget crisis only if dust is produced by both low mass stars and supernovae and is not efficiently destroyed by supernova shocks., Comment: 21 pages including 1 page appendix, 5 figures, 6 tables. Accepted for publication in MNRAS

Noble gas molecules have not hitherto been detected in space. From spectra obtained with the Herschel Space Observatory, we report the detection of emission in the 617.5 GHz and 1234.6 GHz J = 1-0 and 2-1 rotational lines of {36}ArH^+ at several positions in the Crab Nebula, a supernova remnant known to contain both H2 molecules and regions of enhanced ionized argon emission. {36}Ar is believed to have originated from explosive nucleosynthesis in massive stars during core-collapse supernova events. Its detection in the Crab Nebula, the product of such a supernova event, confirms this expectation. The likely excitation mechanism for the observed {36}ArH^+ emission lines is electron collisions in partially ionized regions with electron densities of a few hundred per centimeter cubed., Comment: Published in Science on 13th December 2013; 14 pages, 2 figures, 1 table

Upon its completion the Herschel ATLAS (H-ATLAS) will be the largest submillimetre survey to date, detecting close to half-a-million sources. It will only be possible to measure spectroscopic redshifts for a small fraction of these sources. However, if the rest-frame spectral energy distribution (SED) of a typical H-ATLAS source is known, this SED and the observed Herschel fluxes can be used to estimate the redshifts of the H-ATLAS sources without spectroscopic redshifts. In this paper, we use a subset of 40 H-ATLAS sources with previously measured redshifts in the range 0.5 1 we estimate that the rms of Delta z/(1+z) = 0.12. We have used this template to estimate the redshift distribution for the sources detected in the H-ATLAS equatorial fields, finding a bimodal distribution with a mean redshift of 1.2, 1.9 and 2.5 for 250, 350 and 500 um selected sources respectively. \end{abstract}, Comment: 13 pages, 11 figures

As part of the Herschel Guaranteed Time Key Project MESS (Mass loss of Evolved StarS) we have imaged a sample of planetary nebulae. In this paper we present the PACS and SPIRE images of the classical bipolar planetary nebula NGC 650. We used these images to derive a temperature map of the dust. We also constructed a photoionization and dust radiative transfer model using the spectral synthesis code Cloudy. To constrain this model, we used the PACS and SPIRE fluxes and combined these with hitherto unpublished IUE and Spitzer IRS spectra as well as various other data from the literature. The temperature map combined with the photoionization model were used to study various aspects of the central star, the nebula, and in particular the dust grains in the nebula. The central star parameters are determined to be T_eff = 208 kK and L = 261 L_sol assuming a distance of 1200 pc. The stellar temperature is much higher than previously published values. We confirm that the nebula is carbon-rich with a C/O ratio of 2.1. The nebular abundances are typical for a type IIa planetary nebula. With the photoionization model we determined that the grains in the ionized nebula are large (assuming single-sized grains, they would have a radius of 0.15 micron. Most likely these large grains were inherited from the asymptotic giant branch phase. The PACS 70/160 micron temperature map shows evidence for two radiation components heating the grains. The first component is direct emission from the central star, while the second component is diffuse emission from the ionized gas (mainly Ly alpha). We show that previous suggestions that there is a photo-dissociation region surrounding the ionized region are incorrect. The neutral material resides in dense clumps inside the ionized region. These may also harbor stochastically heated very small grains in addition to the large grains., Comment: 18 pages, 15 figures, 7 tables

We report spectroscopic and imaging observations of rotational transitions of cold CO and SiO in the ejecta of SN1987A, the first such emission detected in a supernova remnant. In addition to line luminosities for the CO J=1-0, 2-1, 6-5, and 7-6 transitions, we present upper limits for all other transitions up to J=13-12, collectively measured from the Atacama Large Millimeter Array (ALMA), the Atacama Pathfinder EXperiment (APEX), and the Herschel Spectral and Photometric Imaging REceiver (SPIRE). Simple models show the lines are emitted from at least 0.01 solar masses of CO at a temperature > 14 K, confined within at most 35% of a spherical volume expanding at ~ 2000 km/s. Moreover, we locate the emission within 1'' of the central debris. These observations, along with a partial observation of SiO, confirm the presence of cold molecular gas within supernova remnants and provide insight into the physical conditions and chemical processes in the ejecta. Furthermore, we demonstrate the powerful new window into supernova ejecta offered by submillimeter observations., Comment: Accepted to the Astrophysical Journal Letters, 6 pages, 3 figures

We present two large, nearby (0.013$\le$z$\le$0.06) samples of Early-Type Galaxies (ETGs): a visually classified sample of 220 ETGs, created using source-matched data from the Galaxy and Mass Assembly (GAMA) database with FIR/sub-mm detections from $Herschel$-ATLAS; and a visually classified sample of 551 ETGs which are undetected with $Herschel$-ATLAS. Active galactic nuclei (AGN) are removed from our samples using optical emission line diagnostics. These samples are scrutinised to determine characteristics of sub-mm detected versus undetected ETGs. We find similarities in the stellar mass distributions of the two ETG samples but testing other properties uncovers significant differences. The sub-mm detected sample is shown to have lower concentration and S\'ersic indices than those in the undetected sample - a result which may be linked to the presence of dust in the former. Optical and UV-optical colours are also shown to be much bluer, indicating that the dust is linked with recent star formation. The intrinsic effective radii are on average 1.5 times larger for the sub-mm detected ETGs. Surface densities and groups data from the GAMA database are examined for the two samples, leading to the conclusion that dusty ETGs inhabit sparser environments than non-dusty ETGs in the nearby universe, although environments of the brightest ETGs are shown to differ the least. Modified Planck functions are fit to the H-ATLAS detected PACS and SPIRE fluxes for ETGs with sub-mm flux densities of at least 3$\sigma$ in the 350$\mu$m SPIRE band, giving a resultant mean cold dust temperature of T$_{d}$=22.1K, with a range of 9-30K. The corresponding mean dust mass is 1.8$\times10^{7}$M$_{\odot}$, with a range of (0.08-35.0)$\times10^{7}$M$_{\odot}$. The dust masses calculated from these fits, normalised by stellar mass, are shown to increase with decreasing stellar mass and bluer colours., Comment: 20 pages, 15 figures, accepted for publication by MNRAS

Context. The interaction between stellar winds and the interstellar medium (ISM) can create complex bow shocks. The photometers on board the Herschel Space Observatory are ideally suited to studying the morphologies of these bow shocks. Aims. We aim to study the circumstellar environment and wind-ISM interaction of the nearest red supergiant, Betelgeuse. Methods. Herschel PACS images at 70, 100, and 160 micron and SPIRE images at 250, 350, and 500 micron were obtained by scanning the region around Betelgeuse. These data were complemented with ultraviolet GALEX data, near-infrared WISE data, and radio 21 cm GALFA-HI data. The observational properties of the bow shock structure were deduced from the data and compared with hydrodynamical simulations. Results. The infrared Herschel images of the environment around Betelgeuse are spectacular, showing the occurrence of multiple arcs at 6-7 arcmin from the central target and the presence of a linear bar at 9 arcmin. Remarkably, no large-scale instabilities are seen in the outer arcs and linear bar. The dust temperature in the outer arcs varies between 40 and 140 K, with the linear bar having the same colour temperature as the arcs. The inner envelope shows clear evidence of a non-homogeneous clumpy structure (beyond 15 arcsec), probably related to the giant convection cells of the outer atmosphere. The non-homogeneous distribution of the material even persists until the collision with the ISM. A strong variation in brightness of the inner clumps at a radius of 2 arcmin suggests a drastic change in mean gas and dust density some 32 000 yr ago. Using hydrodynamical simulations, we try to explain the observed morphology of the bow shock around Betelgeuse. Conclusions: [abbreviated], Comment: 26 pages

Whether supernovae are major sources of dust in galaxies is a long-standing debate. We present infrared and submillimeter photometry and spectroscopy from the Herschel Space Observatory of the Crab Nebula between 51 and 670 micron as part of the Mass Loss from Evolved StarS program. We compare the emission detected with Herschel with multiwavelength data including millimeter, radio, mid-infrared and archive optical images. We carefully remove the synchrotron component using the Herschel and Planck fluxes measured in the same epoch. The contribution from line emission is removed using Herschel spectroscopy combined with Infrared Space Observatory archive data. Several forbidden lines of carbon, oxygen and nitrogen are detected where multiple velocity components are resolved, deduced to be from the nitrogen-depleted, carbon-rich ejecta. No spectral lines are detected in the SPIRE wavebands; in the PACS bands, the line contribution is 5% and 10% at 70 and 100 micron and negligible at 160 micron. After subtracting the synchrotron and line emission, the remaining far-infrared continuum can be fit with two dust components. Assuming standard interstellar silicates, the mass of the cooler component is 0.24(+0.32)(-0.08) Msolar for T = 28.1(+5.5)(-3.2)K. Amorphous carbon grains require 0.11 +/- 0.01 Msolar of dust with T = 33.8(+2.3)(-1.8) K. A single temperature modified-blackbody with 0.14Msolar and 0.08Msolar for silicate and carbon dust respectively, provides an adequate fit to the far- infrared region of the spectral energy distribution but is a poor fit at 24-500 micron. The Crab Nebula has condensed most of the relevant refractory elements into dust, suggesting the formation of dust in core-collapse supernova ejecta is efficient., Comment: 12 pages, 6 figures, ApJ in press, final version after proofs corrected

We present a pan-chromatic analysis of an unprecedented sample of 1402 250 micron-selected galaxies at z < 0.5 (mean z = 0.24) from the Herschel-ATLAS survey. We complement our Herschel 100-500 micron data with UV-K-band photometry from the Galaxy And Mass Assembly (GAMA) survey and apply the MAGPHYS energy-balance technique to produce pan-chromatic SEDs for a representative sample of 250 micron selected galaxies spanning the most recent 5 Gyr of cosmic history. We derive estimates of physical parameters, including star formation rates, stellar masses, dust masses and infrared luminosities. The typical H-ATLAS galaxy at z < 0.5 has a far-infrared luminosity in the range 10^10 - 10^12 Lsolar (SFR: 1-50 Msolar/yr) thus is broadly representative of normal star forming galaxies over this redshift range. We show that 250 micron-selected galaxies contain a larger mass of dust at a given infra-red luminosity or star formation rate than previous samples selected at 60 micron from IRAS. We derive typical SEDs for H-ATLAS galaxies, and show that the emergent SED shape is most sensitive to specific star formation rate. The optical-UV SEDs also become more reddened due to dust at higher redshifts. Our template SEDs are significantly cooler than existing infra-red templates. They may therefore be most appropriate for inferring total IR luminosities from moderate redshift submillimetre selected samples and for inclusion in models of the lower redshift submillimetre galaxy populations., Comment: 26 pages, 24 figures, Accepted by MNRAS

CW Leo has been observed six times between October 2009 and June 2012 with the SPIRE instrument on board the Herschel satellite. Variability has been detected in the flux emitted by the central star with a period of 639 \pm 4 days, in good agreement with determinations in the literature. Variability is also detected in the bow shock around CW Leo that had previously been detected in the ultraviolet and Herschel PACS/SPIRE data. Although difficult to prove directly, our working hypothesis is that this variability is directly related to that of the central star. In this case, fitting a sine curve with the period fixed to 639 days results in a time-lag in the variability between bow shock and the central star of 402 \pm 37 days. The orientation of the bow shock relative to the plane of the sky is unknown (but see below). For an inclination angle of zero degrees, the observed time-lag translates into a distance to CW Leo of 130 \pm 13 pc, and for non-zero inclination angles the distance is smaller. Fitting the shape of the bow shock with an analytical model (Wilkin 1996), the effect of the inclination angle on the distance may be estimated. Making the additional assumption that the relative peculiar velocity between the interstellar medium (ISM) and CW Leo is determined entirely by the star space velocity with respect to the local standard of rest (i.e. a stationary ISM), the inclination angle is found to be (-33.3 \pm 0.8) degrees based on the observed proper motion and radial velocity. Using the Wilkin model, our current best estimate of the distance to CW Leo is 123 \pm 14 pc. For a distance of 123 pc, we derive a mean luminosity of 7790 \pm 150 Lsol (internal error)., Comment: Accepted A&A Letters

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., Comment: 26 pages, 9 figures, accepted for publication in A&A. The complete paper containing the atlas images of the 323 galaxies in the three SPIRE bands and in the SDSS r' band is available on http://hedam.oamp.fr/HRS/publication_download.php

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., Comment: 19 pages, 18 figures. Submitted to ApJ April 2012; Accepted July 2012

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., Comment: Accepted for publication in MNRAS; 11 pages, 8 figures, 1 table

UV and optical surveys are essential to gain insight into the processes driving galaxy formation and evolution. The rest-frame UV emission is key to measure the cosmic SFR. However, UV light is strongly reddened by dust. In starburst galaxies, the UV colour and the attenuation are linked, allowing to correct for dust extinction. Unfortunately, evidence has been accumulating that the relation between UV colour and attenuation is different for normal star-forming galaxies when compared to starburst galaxies. It is still not understood why star-forming galaxies deviate from the UV colour-attenuation relation of starburst galaxies. Previous work and models hint that the role of the shape of the attenuation curve and the age of stellar populations have an important role. In this paper we aim at understanding the fundamental reasons to explain this deviation. We have used the CIGALE SED fitting code to model the far UV to the far IR emission of a set of 7 reasonably face-on spiral galaxies from the HRS. We have explored the influence of a wide range of physical parameters to quantify their influence and impact on the accurate determination of the attenuation from the UV colour, and why normal galaxies do not follow the same relation as starburst galaxies. We have found that the deviation can be best explained by intrinsic UV colour differences between different regions in galaxies. Variations in the shape of the attenuation curve can also play a secondary role. Standard age estimators of the stellar populations prove to be poor predictors of the intrinsic UV colour. These results are also retrieved on a sample of 58 galaxies when considering their integrated fluxes. When correcting the emission of normal star-forming galaxies for the attenuation, it is crucial to take into account possible variations in the intrinsic UV colour as well as variations of the shape of the attenuation curve., Comment: Accepted for publication in A&A, 18 pages, 14 figures. The paper with high resolution figures can be downloaded at http://www.oamp.fr/people/mboquien/HRS/boquien_IRX_beta.pdf

We use the Herschel-ATLAS survey to conduct the first large-scale statistical study of the submm properties of optically selected galaxies. Using ~80,000 r-band selected galaxies from 126 deg^2 of the GAMA survey, we stack into submm imaging at 250, 350 and 500{\mu}m to gain unprecedented statistics on the dust emission from galaxies at z < 0.35. We find that low redshift galaxies account for 5% of the cosmic 250{\mu}m background (4% at 350{\mu}m; 3% at 500{\mu}m), of which approximately 60% comes from 'blue' and 20% from 'red' galaxies (rest-frame g - r). We compare the dust properties of different galaxy populations by dividing the sample into bins of optical luminosity, stellar mass, colour and redshift. In blue galaxies we find that dust temperature and luminosity correlate strongly with stellar mass at a fixed redshift, but red galaxies do not follow these correlations and overall have lower luminosities and temperatures. We make reasonable assumptions to account for the contaminating flux from lensing by red sequence galaxies and conclude that galaxies with different optical colours have fundamentally different dust emission properties. Results indicate that while blue galaxies are more luminous than red galaxies due to higher temperatures, the dust masses of the two samples are relatively similar. Dust mass is shown to correlate with stellar mass, although the dust/stellar mass ratio is much higher for low stellar mass galaxies, consistent with the lowest mass galaxies having the highest specific star formation rates. We stack the 250{\mu}m/NUV luminosity ratio, finding results consistent with greater obscuration of star formation at lower stellar mass and higher redshift. Submm luminosities and dust masses of all galaxies are shown to evolve strongly with redshift, indicating a fall in the amount of obscured star formation in ordinary galaxies over the last four billion years., Comment: 38 pages, 25 figures, MNRAS Accepted

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], Comment: 15 Pages, 9 Figures. Accepted for publication in Astronomy and Astrophysics. A high resolution version of the paper can be found at http://wazn.ugent.be/jfritz/HelgaI_final.pdf

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., Comment: 28 Pages, 12 Figures. Submitted to ApJ December 2011; accepted January 2012

The origin of interstellar dust in galaxies is poorly understood, particularly the relative contributions from supernovae and the cool stellar winds of low-intermediate mass stars. Here, we present Herschel PACS and SPIRE photometry at 70-500um of the historical young supernova remnants: Kepler and Tycho; both thought to be the remnants of Type Ia explosion events. We detect a warm dust component in Kepler's remnant with T = 82K and mass 0.0031Msun; this is spatially coincident with thermal X-ray emission optical knots and filaments, consistent with the warm dust originating in the circumstellar material swept up by the primary blast wave of the remnant. Similarly for Tycho's remnant, we detect warm dust at 90K with mass 0.0086Msun. Comparing the spatial distribution of the warm dust with X-rays from the ejecta and swept-up medium, and Ha emission arising from the post-shock edge, we show that the warm dust is swept up interstellar material. We find no evidence of a cool (25-50 K) component of dust with mass >0.07Msun as observed in core-collapse remnants of massive stars. Neither the warm or cold dust components detected here are spatially coincident with supernova ejecta material. We compare the lack of observed supernova dust with a theoretical model of dust formation in Type Ia remnants which predicts dust masses of 0.088(0.017)Msun for ejecta expanding into surrounding densities of 1(5)cm-3. The model predicts that silicon- and carbon-rich dust grains will encounter the interior edge of the observed dust emission at 400 years confirming that the majority of the warm dust originates from swept up circumstellar or interstellar grains (for Kepler and Tycho respectively). The lack of cold dust grains in the ejecta suggests that Type Ia remnants do not produce substantial quantities of iron-rich dust grains and has important consequences for the 'missing' iron mass observed in ejecta., Comment: 17 pages, 14 figures, accepted for publication in MNRAS, final version including corrected typos and references

In 2007, R Coronae Borealis (R CrB) went into an historically deep and long decline. In this state, the dust acts like a natural coronagraph at visible wavelengths, allowing faint nebulosity around the star to be seen. Imaging has been obtained from 0.5 to 500 micron with Gemini/GMOS, HST/WFPC2, Spitzer/MIPS, and Herschel/SPIRE. Several of the structures around R CrB are cometary globules caused by wind from the star streaming past dense blobs. The estimated dust mass of the knots is consistent with their being responsible for the R CrB declines if they form along the line of sight to the star. In addition, there is a large diffuse shell extending up to 4 pc away from the star containing cool 25 K dust that is detected all the way out to 500 micron. The SED of R CrB can be well fit by a 150 AU disk surrounded by a very large diffuse envelope which corresponds to the size of the observed nebulosity. The total masses of the disk and envelope are 10^-4 and 2 M(Sun), respectively, assuming a gas-to-dust ratio of 100. The evidence pointing toward a white-dwarf merger or a final-helium-shell flash origin for R CrB is contradictory. The shell and the cometary knots are consistent with a fossil planetary nebula. Along with the fact that R CrB shows significant Lithium in its atmosphere, this supports the final-helium-shell flash. However, the relatively high inferred mass of R CrB and its high fluorine abundance support a white-dwarf merger., Comment: 30 pages, 11 figures, ApJ in press

We present the dust properties and star-formation histories of local submillimetre-selected galaxies in Herschel-ATLAS, classified by optical morphology. The early-type galaxies (ETGs) that are detected contain as much dust as typical spirals, and form a unique sample that has been blindly selected at submillimetre wavelengths. Comparing H-ATLAS galaxies to a control sample of optically selected galaxies, we find 5.5% of luminous ETGs are detected in H-ATLAS. The H-ATLAS ETGs contain a significant mass of cold dust: the mean dust mass is 5.5x10^7 Msun, with individual galaxies ranging from 9x10^5-4x10^8 Msun. This is comparable to that of spirals in our sample, and is an order of magnitude more dust than that found for the control ETGs, which have a median dust mass inferred from stacking of (0.8-4.0)x10^6 Msun. The ETGs detected in H-ATLAS have bluer NUV-r colours, higher specific star-formation rates and younger stellar populations than ETGs which are optically selected, and may be transitioning from the blue cloud to the red sequence. We also find that H-ATLAS and control ETGs inhabit similar low-density environments. We conclude that the dust in H-ATLAS and control ETGs cannot be solely from stellar sources, and a large contribution from dust formed in the ISM or external sources is required. Alternatively, dust destruction may not be as efficient as predicted. We also explore the properties of the most passive spiral galaxies in our sample with SSFR<10^-11/yr. We find these passive spirals have lower dust-to-stellar mass ratios, higher stellar masses and older stellar population ages than normal spirals. The passive spirals inhabit low density environments similar to those of the normal spiral galaxies in our sample. This shows that the processes which turn spirals passive do not occur solely in the intermediate density environments of group and cluster outskirts. (Abridged), Comment: 40 pages, 17 figures, 8 tables, accepted by MNRAS (27/09/2011). A high resolution version of the paper can be found at http://www.nottingham.ac.uk/~ppxkr/ETGs_spirals_Rowlands_accepted.pdf

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., Comment: 28 pages, 17 figures, accepted for publication in MNRAS

We investigate the ultraviolet and optical properties and environment of low redshift galaxies detected in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) science demonstration data. We use the Sloan Digital Sky Survey seventh release and the Galaxy And Mass Assembly database to select galaxies with r_Petro < 19.0 mag in the redshift range 0.02 < z < 0.2 and look for their submillimeter counterparts in H-ATLAS. Our results show that at low redshift, H-ATLAS detects mainly blue/star-forming galaxies with a minor contribution from red systems which are highly obscured by dust. In addition we find that the colour of a galaxy rather than the local density of its environment determines whether it is detectable by H-ATLAS. The average dust temperature of galaxies that are simultaneously detected by both PACS and SPIRE is 25K \pm 4K, independent of environment. This analysis provides a glimpse of the potential of the H-ATLAS data to investigate the submillimeter properties of galaxies in the local universe., Comment: 10 pages, 7 figures (accepted for publication in MNRAS)

In 2001, the discovery of circumstellar water vapour around the ageing carbon star IRC+10216 was announced. This detection challenged the current understanding of chemistry in old stars, since water vapour was predicted to be absent in carbon-rich stars. Several explanations for the occurrence of water vapour were postulated, including the vaporization of icy bodies (comets or dwarf planets) in orbit around the star, grain surface reactions, and photochemistry in the outer circumstellar envelope. However, the only water line detected so far from one carbon-rich evolved star can not discriminate, by itself, between the different mechanisms proposed. Here we report on the detection by the Herschel satellite of dozens of water vapour lines in the far-infrared and sub-millimetre spectrum of IRC+10216, including some high-excitation lines with energies corresponding to ~1000 K. The emission of these high-excitation water lines can only be explained if water vapour is present in the warm inner region of the envelope. A plausible explanation for the formation of warm water vapour appears to be the penetration of ultraviolet (UV) photons deep into a clumpy circumstellar envelope. This mechanism triggers also the formation of other molecules such as ammonia, whose observed abundances are much higher than hitherto predicted., Comment: Main article: 18 pages, 4 figures Supplementary information: 12 pages, 6 figures

MESS (Mass-loss of Evolved StarS) is a Guaranteed Time Key Program that uses the PACS and SPIRE instruments on board the Herschel Space Observatory to observe a representative sample of evolved stars, that include asymptotic giant branch (AGB) and post-AGB stars, planetary nebulae and red supergiants, as well as luminous blue variables, Wolf-Rayet stars and supernova remnants. In total, of order 150 objects are observed in imaging and about 50 objects in spectroscopy. This paper describes the target selection and target list, and the observing strategy. Key science projects are described, and illustrated using results obtained during Herschel's science demonstration phase. Aperture photometry is given for the 70 AGB and post-AGB stars observed up to October 17, 2010, which constitutes the largest single uniform database of far-IR and sub-mm fluxes for late-type stars., Comment: A&A accepted

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., Comment: 6 pages, 4 figures, accepted for publication in MNRAS