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

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

Here we introduce the interstellar dust modelling framework THEMIS (The Heterogeneous dust Evolution Model for Interstellar Solids), which takes a global view of dust and its evolution in response to the local conditions in interstellar media. This approach is built upon a core model that was developed to explain the dust extinction and emission in the diffuse interstellar medium. The model was then further developed to self-consistently include the effects of dust evolution in the transition to denser regions. The THEMIS approach is under continuous development and currently we are extending the framework to explore the implications of dust evolution in HII regions and the photon-dominated regions associated with star formation. We provide links to the THEMIS, DustEM and DustPedia websites where more information about the model, its input data and applications can be found.

Dust modeling is crucial to infer dust properties and budget for galaxy studies. However, there are systematic disparities between dust grain models that result in corresponding systematic differences in the inferred dust properties of galaxies. Quantifying these systematics requires a consistent fitting analysis. We compare the output dust parameters and assess the differences between two dust grain models, Compiegne et al (2011), and THEMIS (Jones et al 2013, Kohler et al 2015). In this study, we use a single fitting method applied to all the models to extract a coherent and unique statistical analysis. We fit the models to the dust emission seen by Spitzer and Herschel in the Small and Large Magellanic Clouds (SMC and LMC). The observations cover the infrared (IR) spectrum from a few microns to the sub-millimeter range. For each fitted pixel, we calculate the full n-D likelihood, based on the method described in Gordon et al (2014). The free parameters are both environmental ($U$, the interstellar radiation field strength; $\alpha_\mathrm{ISRF}$, power-law coefficient for a multi-U environment; $\Omega^*$, the starlight strength) and intrinsic to the model ($Y_\mathrm{i}$: abundances of the grain species $i$; $\alpha_\mathrm{sCM20}$, coefficient in the small carbon grain size distribution). Fractional residuals of 5 different sets of parameters show that fitting THEMIS brings a more accurate reproduction of the observations than the Compiegne model. However, independent variations of the dust species show strong model-dependencies. We find that the abundance of silicates can only be constrained to an upper-limit and the silicate/carbon ratio is different than that seen in our Galaxy. In the LMC, our fits result in dust masses slightly lower than those found in literature , by a factor lower than 2. In the SMC, we find dust masses in agreement with previous studies.

Dust grains are classically thought to form in the winds of asymptotic giant branch (AGB) stars. However, there is increasing evidence today for dust formation in supernovae (SNe). To establish the relative importance of these two classes of stellar sources of dust, it is important to know the fraction of freshly formed dust in SN ejecta that is able to survive the passage of the reverse shock and be injected in the interstellar medium. We have developed a new code (GRASH\_Rev) which follows the newly-formed dust evolution throughout the supernova explosion until the merging of the forward shock with the circumstellar ISM. We have considered four well studied SNe in the Milky Way and Large Magellanic Cloud: SN1987A, CasA, the Crab Nebula, and N49. For all the simulated models, we find good agreement with observations and estimate that between 1 and 8$\%$ of the observed mass will survive, leading to a SN dust production rate of $(3.9 \pm 3.7) \times 10^{-4}$ M$_{\odot}$yr$^{-1}$ in the Milky Way. This value is one order of magnitude larger than the dust production rate by AGB stars but insufficient to counterbalance the dust destruction by SNe, therefore requiring dust accretion in the gas phase.

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

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

Astrophysical images issued from different instruments and/or spectral bands often require to be processed together, either for fitting or comparison purposes. However each image is affected by an instrumental response, also known as PSF, that depends on the characteristics of the instrument as well as the wavelength and the observing strategy. Given the knowledge of the PSF in each band, a straightforward way of processing images is to homogenise them all to a target PSF using convolution kernels, so that they appear as if they had been acquired by the same instrument. We propose an algorithm that generates such PSF-matching kernels, based on Wiener filtering with a tunable regularisation parameter. This method ensures all anisotropic features in the PSFs to be taken into account. We compare our method to existing procedures using measured Herschel/PACS and SPIRE PSFs and simulated JWST/MIRI PSFs. Significant gains up to two orders of magnitude are obtained with respect to the use of kernels computed assuming Gaussian or circularised PSFs.

The knowledge of the point spread function (PSF) of imaging instruments represents a fundamental requirement for astronomical observations. The Herschel PACS PSFs delivered by the instrument control centre are obtained from observations of the Vesta asteroid, providing a characterisation of the central part therefore excluding fainter features. However, in many cases information on both the core and the wings of the PSFs is needed. With this aim, we combine Vesta and Mars dedicated observations and obtain PACS PSFs with an unprecedented dynamic range $(\sim 10^6)$, at slow and fast scan speeds and for the three photometric bands.

Dust grains are classically thought to form in the winds of AGB stars. However, nowadays there is increasing evidence for dust formation in SNe. In order to establish the relative importance of these two classes of stellar sources of dust it is important to know what is the fraction of freshly formed dust in SN ejecta that is able to survive the passage of the reverse shock and be injected in the interstellar medium. With this aim, we have developed a new code, GRASH_Rev, that allows to follow the dynamics of dust grains in the shocked SN ejecta and to compute the time evolution of the mass, composition and size distribution of the grains. We consider four well studied SNe in the Milky Way and LMC: SN 1987a, Cas A, the Crab Nebula, and N49. For all the simulated models, we find good agreement with observations. Our study suggests that SN 1987A is too young for the reverse shock to have affected the dust mass. Conversely, in the other three SNe, the reverse shock has already destroyed between 10 and 40% of the initial dust mass. However, the largest dust mass destruction is predicted to occur between 10^3 and 10^5 yr after the explosions. Since the oldest SN in the sample has an estimated age of 4800 yr, current observations can only provide an upper limit to the mass of SN dust that will enrich the interstellar medium, the so-called effective dust yields. We find that only between 1 and 8% of the currently observed mass will survive. This is in good agreement with the values adopted in chemical evolution models which consider the effect of the SN reverse shock. We discuss the astrophysical implications of our results for dust enrichment in local galaxies and at high redshift.

We quantify the role of Population (Pop) III core-collapse supernovae (SNe) as the first cosmic dust polluters. Starting from a homogeneous set of stellar progenitors with masses in the range [13 - 80] Msun, we find that the mass and composition of newly formed dust depend on the mixing efficiency of the ejecta and the degree of fallback experienced during the explosion. For standard Pop III SNe, whose explosions are calibrated to reproduce the average elemental abundances of Galactic halo stars with [Fe/H] < -2.5, between 0.18 and 3.1 Msun (0.39 - 1.76 Msun) of dust can form in uniformly mixed (unmixed) ejecta, and the dominant grain species are silicates. We also investigate dust formation in the ejecta of faint Pop III SN, where the ejecta experience a strong fallback. By examining a set of models, tailored to minimize the scatter with the abundances of carbon-enhanced Galactic halo stars with [Fe/H ] < -4, we find that amorphous carbon is the only grain species that forms, with masses in the range 2.7 10^{-3} - 0.27 Msun (7.5 10^{-4} - 0.11 Msun) for uniformly mixed (unmixed) ejecta models. Finally, for all the models we estimate the amount and composition of dust that survives the passage of the reverse shock, and find that, depending on circumstellar medium densities, between 3 and 50% (10 - 80%) of dust produced by standard (faint) Pop III SNe can contribute to early dust enrichment.

Observations of edge-on galaxies allow us to investigate the vertical extent and properties of dust, gas and stellar distributions. NGC 891 has been studied for decades and represents one of the best studied cases of an edge-on galaxy. We use deep PACS data together with IRAC, MIPS and SPIRE data to study the vertical extent of dust emission around NGC 891. We also test the presence of a more extended, thick dust component. By performing a convolution of an intrinsic vertical profile emission with each instrument PSF and comparing it with observations we derived the scaleheight of a thin and thick dust disc component. For all wavelengths considered the emission is best fit with the sum of a thin and a thick dust component. The scaleheight of both dust components shows a gradient passing from 70 $\mu$m to 250 $\mu$m. This could be due to a drop in dust heating (and thus dust temperature) with the distance from the plane, or to a sizable contribution ($\sim 15 - 80%$) of an unresolved thin disc of hotter dust to the observed surface brightness at shorter wavelengths. The scaleheight of the thick dust component, using observations from 70 $\mu$m to 250 $\mu$m has been estimated to be $(1.44\pm 0.12)$ kpc, consistent with previous estimates (extinction and scattering in optical bands and MIR emission). The amount of dust mass at distances larger than $\sim 2$ kpc from the midplane represents $2 - 3.3$ % of the total galactic dust mass and the relative abundance of small grains with respect to large grains is almost halved comparing to that in the midplane. The paucity of small grains high above the midplane might indicate that dust is hit by interstellar shocks or galactic fountains and entrained together with gas. The halo dust component is likely to be embedded in an atomic / molecular gas and heated by a thick stellar disc.

Context. The evolution of amorphous hydrocarbon materials, a-C(:H), principally resulting from ultraviolet (UV) photon absorption- induced processing, are likely at the heart of the variations in the observed properties of dust in the interstellar medium. Aims. The consequences of the size-dependent and compositional variations in a-C(:H), from aliphatic-rich a-C:H to aromatic-rich a-C, are studied within the context of the interstellar dust extinction and emission. Methods. Newly-derived optical property data for a-C(:H) materials, combined with that for an amorphous forsterite-type silicate with iron nano-particle inclusions, a-SilFe, are used to explore dust evolution in the interstellar medium. Results. We present a new dust model that consists of a power-law distribution of small a-C grains and log-normal distributions of large a-SilFe and a-C(:H) grains. The model, which is firmly anchored by laboratory-data, is shown to quite naturally explain the variations in the infrared (IR) to far-ultraviolet (FUV) extinction, the 217 nm UV bump, the IR absorption and emission bands and the IR-mm dust emission. Conclusions. The major strengths of the new model are its inherent simplicity and built-in capacity to follow dust evolution in interstellar media. We show that mantle accretion in molecular clouds and UV photo-processing in photo-dominated regions are likely the major drivers of dust evolution.

Observational evidence seems to indicate that the depletion of interstellar carbon into dust shows rather wide variations and that carbon undergoes rather rapid recycling in the interstellar medium (ISM). Small hydrocarbon grains are processed in photo-dissociation regions by UV photons, by ion and electron collisions in interstellar shock waves and by cosmic rays. A significant fraction of hydrocarbon dust must therefore be re-formed by accretion in the dense, molecular ISM. A new dust model (Jones et al., Astron. Astrophys., 2013, 558, A62) shows that variations in the dust observables in the diffuse interstellar medium (nH = 1000 cm^3), can be explained by systematic and environmentally-driven changes in the small hydrocarbon grain population. Here we explore the consequences of gas-phase carbon accretion onto the surfaces of grains in the transition regions between the diffuse ISM and molecular clouds (e.g., Jones, Astron. Astrophys., 2013, 555, A39). We find that significant carbonaceous dust re-processing and/or mantle accretion can occur in the outer regions of molecular clouds and that this dust will have significantly different optical properties from the dust in the adjacent diffuse ISM. We conclude that the (re-)processing and cycling of carbon into and out of dust is perhaps the key to advancing our understanding of dust evolution in the ISM.

The $Herschel$ Fornax Cluster Survey (HeFoCS) is a deep, far-infrared (FIR) survey of the Fornax cluster. The survey is in 5 $Herschel$ bands (100 - 500 $\mu$m) and covers an area of 16 deg$^2$ centred on NGC1399. This paper presents photometry, detection rates, dust masses and temperatures using an optically selected sample from the Fornax Cluster Catalogue (FCC). Our results are compared with those previously obtained using data from the $Herschel$ Virgo Cluster Survey (HeViCS). In Fornax, we detect 30 of the 237 (13%) optically selected galaxies in at least one $Herschel$ band. The global detection rates are significantly lower than Virgo, reflecting the morphological make up of each cluster - Fornax has a lower fraction of late-type galaxies. For galaxies detected in at least 3 bands we fit a modified blackbody with a $\beta = 2$ emissivity. Detected early-type galaxies (E/S0) have a mean dust mass, temperature, and dust-to-stars ratio of $\log_{10}(<M_{dust}>/\mathrm{M_{\odot}}) = 5.82 \pm 0.20$, $<T_{dust}> = 20.82 \pm 1.77$K, and $\log_{10}(M_{dust}/M_{stars}) = -3.87 \pm 0.28$, respectively. Late-type galaxies (Sa to Sd) have a mean dust mass, temperature, and dust-to-stars ratio of $\log_{10}(<M_{dust}>/\mathrm{M_{\odot}}) = 6.54 \pm 0.19$, $<T_{dust}> = 17.47 \pm 0.97$K, and $\log_{10}(M_{dust}/M_{stars}) = -2.93 \pm 0.09$, respectively. The different cluster environments seem to have had little effect on the FIR properties of the galaxies and so we conclude that any environment dependent evolution, has taken place before the cluster was assembled.

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