Search

In the framework of the IMEGIN Large Program, we used the NIKA2 camera on the IRAM 30-m telescope to observe the edge-on galaxy NGC 891 at 1.15 mm and 2 mm and at a FWHM of 11.1" and 17.6", respectively. Multiwavelength data enriched with the new NIKA2 observations fitted by the HerBIE SED code (coupled with the THEMIS dust model) were used to constrain the physical properties of the ISM. Emission originating from the diffuse dust disk is detected at all wavelengths from mid-IR to mm, while mid-IR observations reveal warm dust emission from compact HII regions. Indications of mm excess emission have also been found in the outer parts of the galactic disk. Furthermore, our SED fitting analysis constrained the mass fraction of the small (< 15 Angstrom) dust grains. We found that small grains constitute 9.5% of the total dust mass in the galactic plane, but this fraction increases up to ~ 20% at large distances (|z| > 3 kpc) from the galactic plane., Comment: To appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

As the millimeter wavelength range remains a largely unexplored spectral region for galaxies, the IMEGIN large program aims to map the millimeter continuum emission of 22 nearby galaxies at 1.15 and 2 mm. Using the high-resolution maps produced by the NIKA2 camera, we explore the existence of very cold dust and take possible contamination by free-free and synchrotron emission into account. We study the IR-to-radio emission coming from different regions along the galactic plane and at large vertical distances. New observations of NGC 891, using the NIKA2 camera on the IRAM 30m telescope, along with a suite of observations at other wavelengths were used to perform a multiwavelength study of the spectral energy distribution in the interstellar medium in this galaxy. This analysis was performed globally and locally, using the advanced hierarchical Bayesian fitting code, HerBIE, coupled with the THEMIS dust model. Our dust modeling is able to reproduce the near-IR to millimeter emission of NGC 891, with the exception of an excess at a level of 25% obtained by the NIKA2 observations in the outermost parts of the disk. The radio continuum and thermal dust emission are distributed differently in the disk and galaxy halo. Different dusty environments are also revealed by a multiwavelength investigation of the emission features. Our detailed decomposition at millimeter and centimeter wavelengths shows that emission at 1 mm is purely originated by dust. Radio components become progressively important with increasing wavelengths. Finally, we find that emission arising from small dust grains accounts for ~ 9.5% of the total dust mass, reaching up to 20% at large galactic latitudes. Shock waves in the outflows that shatter the dust grains might explain this higher fraction of small grains in the halo., Comment: 15 pages, 10 figures, Accepted for publication in A&A

Context: The far-infrared (FIR) and sub-millimeter (submm) emissivity of the Milky Way (MW) cirrus is an important benchmark for dust grain models. Dust masses in other galaxies are generally derived from the FIR/submm using the emission properties of these MW-calibrated models. Aims: We seek to derive the FIR/submm emissivity in nine nearby spiral galaxies to check its compatibility with MW cirrus measurements. Methods: We obtained values of the emissivity at 70 to 500 um, using maps of dust emission from the Herschel satellite and of gas surface density from the THINGS and HERACLES surveys on a scale generally corresponding to 440 pc. We studied the variation of the emissivity with the surface brightness ratio I(250um)/I(500um), a proxy for the intensity of the interstellar radiation field heating the dust. Results: We find that the average value of the emissivity agrees with MW estimates for pixels sharing the same color as the cirrus, namely, for I(250um)/I(500um) = 4.5. For I(250um)/I(500um) > 5, the measured emissivity is instead up to a factor ~2 lower than predicted from MW dust models heated by stronger radiation fields. Regions with higher I(250um)/I(500um) are preferentially closer to the galactic center and have a higher overall (stellar+gas) surface density and molecular fraction. The results do not depend strongly on the adopted CO-to-molecular conversion factor and do not appear to be affected by the mixing of heating conditions. Conclusions: Our results confirm the validity of MW dust models at low density, but are at odds with predictions for grain evolution in higher density environments. If the lower-than-expected emissivity at high I(250um)/I(500um) is the result of intrinsic variations in the dust properties, it would imply an underestimation of the dust mass surface density of up to a factor ~2 when using current dust models., Comment: 16 pages, 11 figures, A&A accepted

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

The Nearby Evolved Stars Survey (NESS) is a volume-complete sample of $\sim$850 Galactic evolved stars within 3\,kpc at (sub-)mm wavelengths, observed in the CO $J = $ (2$-$1) and (3$-$2) rotational lines, and the sub-mm continuum, using the James Clark Maxwell Telescope and Atacama Pathfinder Experiment. NESS consists of five tiers, based on distances and dust-production rate (DPR). We define a new metric for estimating the distances to evolved stars and compare its results to \emph{Gaia} EDR3. Replicating other studies, the most-evolved, highly enshrouded objects in the Galactic Plane dominate the dust returned by our sources, and we initially estimate a total DPR of $4.7\times 10^{-5}$ M$_\odot$ yr$^{-1}$ from our sample. Our sub-mm fluxes are systematically higher and spectral indices are typically shallower than dust models typically predict. The 450/850 $\mu$m spectral indices are consistent with the blackbody Rayleigh--Jeans regime, suggesting a large fraction of evolved stars have unexpectedly large envelopes of cold dust., Comment: 20 pages, 13 figures, 5 tables, accepted for publications in MNRAS

Scaling laws of dust, HI gas and metal mass with stellar mass, specific star formation rate and metallicity are crucial to our understanding of the buildup of galaxies through their enrichment with metals and dust. In this work, we analyse how the dust and metal content varies with specific gas mass ($M_{\text{HI}}$/$M_{\star}$) across a diverse sample of 423 nearby galaxies. The observed trends are interpreted with a set of Dust and Element evolUtion modelS (DEUS) - incluidng stellar dust production, grain growth, and dust destruction - within a Bayesian framework to enable a rigorous search of the multi-dimensional parameter space. We find that these scaling laws for galaxies with $-1.0\lesssim \log M_{\text{HI}}$/$M_{\star}\lesssim0$ can be reproduced using closed-box models with high fractions (37-89$\%$) of supernova dust surviving a reverse shock, relatively low grain growth efficiencies ($\epsilon$=30-40), and long dus lifetimes (1-2\,Gyr). The models have present-day dust masses with similar contributions from stellar sources (50-80\,$\%$) and grain growth (20-50\,$\%$). Over the entire lifetime of these galaxies, the contribution from stardust ($>$90\,$\%$) outweighs the fraction of dust grown in the interstellar medium ($<$10$\%$). Our results provide an alternative for the chemical evolution models that require extremely low supernova dust production efficiencies and short grain growth timescales to reproduce local scaling laws, and could help solving the conundrum on whether or not grains can grow efficiently in the interstellar medium., Comment: 56 pages, 30 figures, accepted for publication in MNRAS

Dust plays an important role in shaping a galaxy's spectral energy distribution (SED). It absorbs ultraviolet (UV) to near-infrared (NIR) radiation and re-emits this energy in the far-infrared (FIR). The FIR is essential to understand dust in galaxies. However, deep FIR observations require a space mission, none of which are still active today. We aim to infer the FIR emission across six Herschel bands, along with dust luminosity, mass, and effective temperature, based on the available UV to mid-infrared (MIR) observations. We also want to estimate the uncertainties of these predictions, compare our method to energy balance SED fitting, and determine possible limitations of the model. We propose a machine learning framework to predict the FIR fluxes from 14 UV-MIR broadband fluxes. We used a low redshift sample by combining DustPedia and H-ATLAS, and extracted Bayesian flux posteriors through SED fitting. We trained shallow neural networks to predict the far-infrared fluxes, uncertainties, and dust properties. We evaluated them on a test set using a root mean square error (RMSE) in log-space. Our results (RMSE = 0.19 dex) significantly outperform UV-MIR energy balance SED fitting (RMSE = 0.38 dex), and are inherently unbiased. We can identify when the predictions are off, for example when the input has large uncertainties on WISE 22, or when the input does not resemble the training set. The galaxies for which we have UV-FIR observations can be used as a blueprint for galaxies that lack FIR data. This results in a 'virtual FIR telescope', which can be applied to large optical-MIR galaxy samples. This helps bridge the gap until the next FIR mission., Comment: 25 pages, 16 figures, accepted by A&A. Interactive figures + jupyter notebooks at https://wdobbels.github.io/FIREnet/

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

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., Comment: 13 pages, 6 figures, 2 tables, A&A accepted

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., Comment: Accepted for publication in Astronomy and Astrophysics

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., Comment: 22 pages, 18 figures

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., Comment: 24 pages, 15 figures, published in A&A, bibref: 2019A&A...623A...5D, This Version 2 has increased accuracy in Table B1 in appendix B

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.751.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., Comment: 24 pages, 17 figures, Accepted for publication in Astronomy & Astrophysics

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., Comment: 21 pages, 14 figures, A&A accepted. Template SEDs available at the DustPedia Archive: http://dustpedia.astro.noa.gr/

The thermal emission of dust grains is a powerful tool for probing cold, dense regions of molecular gas in the ISM, and so constraining dust properties is key to obtaining accurate measurements of dust mass and temperature. By placing constraints on the dust emissivity spectral index, beta, towards two star-forming infrared dark clouds, SDC18.888 and SDC24.489, we evaluate the role of mass concentration in the associated star-formation activity. We exploit the simultaneous 1.2mm and 2.0mm imaging capability of NIKA on the IRAM 30m telescope to construct maps of beta for both clouds, and by incorporating Herschel observations, we create H2 column density maps with 13" resolution. While we find no significant systematic radial variations around the most massive clumps in either cloud on >0.1 pc scales, their mean beta values are significantly different, with beta = 2.07 +/- 0.09 (rand) +/- 0.25 (syst) for SDC18.888 and beta = 1.71 +/- 0.09 (rand) +/- 0.25 (syst) for SDC24.489. These differences could be a consequence of the very different environments in which both clouds lie, and we suggest that the proximity of SDC18.888 to the W39 HII region may raise beta on scales of 1 pc. We also find that the mass in SDC24.489 is more centrally concentrated and circularly symmetric than in SDC18.888, and is consistent with a scenario in which spherical globally-collapsing clouds concentrate a higher fraction of their mass into a single core than elongated clouds that will more easily fragment, distributing their mass into many cores. We demonstrate that beta variations towards interstellar clouds can be robustly constrained with high-SNR NIKA observations, providing more accurate estimates of their masses. The methods presented here will be applied to the Galactic Star Formation with NIKA2 (GASTON) large programme, extending our analysis to a statistically significant sample of star-forming clouds., Comment: Abstract abridged. 18 pages, 16 figures, 2 tables. Accepted for publication in Astronomy & Astrophysics

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

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

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 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 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., Comment: 25 pages, 19 figures, A&A accepted

Using results from the Herschel Astrophysical Terrahertz Large-Area Survey and the Galaxy and Mass Assembly project, we show that, for galaxy masses above approximately 1.0e8 solar masses, 51% of the stellar mass-density in the local Universe is in early-type galaxies (ETGs: Sersic n > 2.5) while 89% of the rate of production of stellar mass-density is occurring in late-type galaxies (LTGs: Sersic n < 2.5). From this zero-redshift benchmark, we have used a calorimetric technique to quantify the importance of the morphological transformation of galaxies over the history of the Universe. The extragalactic background radiation contains all the energy generated by nuclear fusion in stars since the Big Bang. By resolving this background radiation into individual galaxies using the deepest far-infrared survey with the Herschel Space Observatory and a deep near-infrared/optical survey with the Hubble Space Telescope (HST), and using measurements of the Sersic index of these galaxies derived from the HST images, we estimate that approximately 83% of the stellar mass-density formed over the history of the Universe occurred in LTGs. The difference between this and the fraction of the stellar mass-density that is in LTGs today implies there must have been a major transformation of LTGs into ETGs after the formation of most of the stars., Comment: Accepted for publication in MNRAS

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

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

We present an analysis of CO molecular gas tracers in a sample of 500{\mu}m-selected Herschel-ATLAS galaxies at z<0.05 (cz<14990km/s). Using 22-500{\mu}m photometry from WISE, IRAS and Herschel, with HI data from the literature, we investigate correlations between warm and cold dust, and tracers of the gas in different phases. The correlation between global CO(3-2) line fluxes and FIR-submillimetre fluxes weakens with increasing IR wavelength ({\lambda}>60{\mu}m), as a result of colder dust being less strongly associated with dense gas. Conversely, CO(2-1) and HI line fluxes both appear to be better correlated with longer wavelengths, suggesting that cold dust is more strongly associated with diffuse atomic and molecular gas phases, consistent with it being at least partially heated by radiation from old stellar populations. The increased scatter at long wavelengths implies that submillimetre fluxes are a poorer tracer of SFR. Fluxes at 22 and 60{\mu}m are also better correlated with diffuse gas tracers than dense CO(3-2), probably due to very-small-grain emission in the diffuse interstellar medium, which is not correlated with SFR. The FIR/CO luminosity ratio and the dust mass/CO luminosity ratio both decrease with increasing luminosity, as a result of either correlations between mass and metallicity (changing CO/H2) or between CO luminosity and excitation [changing CO(3-2)/CO(1-0)]., Comment: 27 pages, 11 figures. MNRAS accepted

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

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

The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) provides an unprecedented opportunity to search for blazars at sub-mm wavelengths. We cross-matched the FIRST radio source catalogue with the 11655 sources brighter than 35 mJy at 500{\mu}m in the \sim 135 square degrees of the sky covered by the H-ATLAS equatorial fields at 9 h and 15 h, plus half of the field at 12 h. We found that 379 of the H-ATLAS sources have a FIRST counterpart within 10 arcsec, including 8 catalogued blazars (plus one known blazar that was found at the edge of one the H-ATLAS maps). To search for additional blazar candidates we have devised new diagnostic diagrams and found that known blazars occupy a region of the log(S500{\mu}m/S350{\mu}m) vs. log(S500{\mu}m/S1.4GHz) plane separated from that of the other sub-mm sources with radio counterparts. Using this diagnostic we have selected 12 further candidates that turn out to be scattered in the (r-z) vs. (u-r) plane or in the WISE colour-colour diagram proposed by Massaro et al. (2012), where known blazars are concentrated in well defined strips. This suggests that the majority of them either are not blazars or have spectral energy distributions contaminated by their host galaxies. A significant fraction of true blazars are found to be hosted by star-forming galaxies. This finding, supported by an analysis of blazars detected in Planck 545 and 857 GHz bands, is at odds with the notion that blazar hosts are passive ellipticals and indicates that the sub-mm selection is providing a novel prospect on blazar properties. Based on an inspection of the available photometric data, including the WISE all-sky survey, the unpublished VIKING survey and new radio observations, we tentatively estimate that there are 11 blazars with synchrotron flux density S500{\mu}m > 35mJy over the considered area. This result already allows us to constrain blazar evolution models., Comment: 13 pages, 5 figures, 5 tables. Submitted to MNRAS

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