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

We combine data from ALMA and MUSE to study the resolved (~300 pc scale) star formation relation (star formation rate vs. molecular gas surface density) in cluster galaxies. Our sample consists of 9 Fornax cluster galaxies, including spirals, ellipticals, and dwarfs, covering a stellar mass range of ~10^8.8 - 10^11 M_Sun. CO(1-0) and extinction corrected Halpha were used as tracers for the molecular gas mass and star formation rate, respectively. We compare our results with Kennicutt (1998) and Bigiel et al. (2008). Furthermore, we create depletion time maps to reveal small-scale variations in individual galaxies. We explore these further in FCC290, using the 'uncertainty principle for star formation' (Kruijssen & Longmore, 2014a) to estimate molecular cloud lifetimes, which we find to be short (<10 Myr) in this galaxy. Galaxy-averaged depletion times are compared with other parameters such as stellar mass and cluster-centric distance. We find that the star formation relation in the Fornax cluster is close to those from Kennicutt (1998) and Bigiel et al. (2008}), but overlaps mostly with the shortest depletion times predicted by Bigiel et al. (2008). This slight decrease in depletion time is mostly driven by dwarf galaxies with disturbed molecular gas reservoirs close to the virial radius. In FCC90, a dwarf galaxy with a molecular gas tail, we find that depletion times are a factor >~10 higher in its tail than in its stellar body.

Globular clusters (GCs) are found ubiquitously in massive galaxies and due to their old ages, they are regarded as fossil records of galaxy evolution. Spectroscopic studies of GC systems are often limited to the outskirts of galaxies, where GCs stand out against the galaxy background and serve as bright tracers of galaxy assembly. In this work, we use the capabilities of the Multi Unit Explorer Spectrograph (MUSE) to extract a spectroscopic sample of 722 GCs in the inner regions ($\lesssim 3 R_\text{eff}$) of 32 galaxies in the Fornax cluster. These galaxies were observed as part of the Fornax 3D project, a MUSE survey that targets early and late-type galaxies within the virial radius of Fornax. After accounting for the galaxy background in the GC spectra, we extracted line-of-sight velocities and determined metallicities of a sub-sample of 238 GCs. We found signatures of rotation within GC systems, and comparing the GC kinematics and that of the stellar body shows that the GCs trace the spheroid of the galaxies. While the red GCs prove to closely follow the metallicity profile of the host galaxy, the blue GCs show a large spread of metallicities but they are generally more metal-poor than the host.

Globular cluster (GC) systems of massive galaxies often show a bimodal colour distribution. This has been interpreted as a metallicity bimodality, created by a two-stage galaxy formation where the red, metal-rich GCs were formed in the parent halo and the blue metal-poor GCs were accreted. This interpretation, however, crucially depends on the assumption that GCs are exclusively old stellar systems with a linear colour-metallicity relation (CZR). The shape of the CZR and range of GC ages are currently under debate, because their study requires high quality spectra to derive reliable stellar population properties. We determined metallicities with full spectral fitting from a sample of 187 GCs with high spectral signal-to-noise ratio in 23 galaxies of the Fornax cluster that were observed as part of the Fornax 3D project. The derived CZR from this sample is non-linear and can be described by a piecewise linear function with a break point at ($g - z$) $\sim$ 1.1 mag. The less massive galaxies in our sample ($M_\ast < 10^{10} M_\odot$) appear to have slightly younger GCs, but the shape of the CZR is insensitive to the GC ages. Although the least massive galaxies lack red, metal-rich GCs, a non-linear CZR is found irrespective of the galaxy mass, even in the most massive galaxies ($M_\ast \geq 10^{11} M_\odot$). Our CZR predicts narrow unimodal GC metallicity distributions for low mass and broad unimodal distributions for very massive galaxies, dominated by a metal-poor and metal-rich peak, respectively, and bimodal distributions for galaxies with intermediate masses (10$^{10}$ $\leq$ $M_\ast < 10^{11} M_\odot$) as a consequence of the relative fraction of red and blue GCs. The diverse metallicity distributions challenge the simple differentiation of GC populations solely based on their colour.

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

Extragalactic planetary nebulae (PNe) are detectable through relatively strong nebulous [OIII] emission and act as direct probes into the local stellar population. Because they have an apparently universal invariant magnitude cut-off, PNe are also considered to be a remarkable standard candle for distance estimation. Through detecting PNe within the galaxies, we aim to connect the relative abundances of PNe to the properties of their host galaxy stellar population. By removing the stellar background components from FCC 167 and FCC 219, we aim to produce PN luminosity functions (PNLF) of these galaxies, and thereby also estimate the distance modulus to these two systems. Finally, we test the reliability and robustness of our novel detection and analysis method. It detects unresolved point sources by their [OIII] 5007{\AA} emission within regions that have previously been unexplored. We model the [OIII] emissions in the spatial and spectral dimensions together, as afforded to us by the Multi Unit Spectroscopic Explorer (MUSE), and we draw on data gathered as part of the Fornax3D survey. For each source, we inspect the properties of the nebular emission lines to remove other sources that might hinder the safe construction of the PNLF, such as supernova remnants and HII regions. As a further step, we characterise any potential limitations and draw conclusions about the reliability of our modelling approach through a set of simulations. By applying this novel detection and modelling approach to integral field unit observations, we report for the distance estimates and luminosity-specific PNe frequency values for the two galaxies. Furthermore, we include an overview of source contamination, galaxy differences, and possible effects on the PNe populations in the dense stellar environments.

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

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

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

In our local Galactic neighborhood, molecular clouds are best studied using a combination of dust measurements, to determine robust masses, sizes and internal structures of the clouds, and molecular-line observations to determine cloud kinematics and chemistry. We present here the first results of a program designed to extend such studies to nearby galaxies beyond the Magellanic Clouds. Utilizing the wideband upgrade of the Submillimeter Array (SMA) at 230 GHz we have obtained the first continuum detections of the thermal dust emission on sub-GMC scales ($\sim$ 15 pc) within the Andromeda galaxy (M31). These include the first resolved continuum detections of dust emission from individual GMCs beyond the Magellanic Clouds. Utilizing a powerful capability of the SMA, we simultaneously recorded CO(2-1) emission with identical $(u,\,v)$ coverage, astrometry and calibration, enabling the first measurements of the CO conversion factor, $\alpha_{\rm\,CO(2-1)}$, toward individual GMCs across an external galaxy. Our direct measurement yields an average CO--to--dust mass conversion factor of $\alpha^\prime_{\rm CO-dust} = 0.042\pm0.018$ $M_\odot$ (K km s$^{-1}$ pc$^2$)$^{-1}$ for the $J= 2-1$ transition. This value does not appear to vary with galactocentric radius. Assuming a constant gas-to-dust ratio of 136, the resulting $\alpha_{\rm CO}$ $=$ 5.7 $\pm$ 2.4 $M_\odot$ (K km s$^{-1}$ pc$^2$)$^{-1}$ for the 2-1 transition is in excellent agreement with that of Milky Way GMCs, given the uncertainties. Finally, using the same analysis techniques, we compare our results with observations of the local Orion molecular clouds, placed at the distance of M31 and simulated to appear as they would if observed by the SMA.

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

We present a CO(3-2) survey of selected regions in the M31 disc as part of the JCMT large programme, HARP and SCUBA-2 High-Resolution Terahertz Andromeda Galaxy Survey (HASHTAG). The 12 CO(3-2) fields in this survey cover a total area of 60 square arcminutes, spanning a deprojected radial range of 2 - 14 kpc across the M31 disc. Combining these observations with existing IRAM 30m CO(1-0) observations and JCMT CO(3-2) maps of the nuclear region of M31, as well as dust temperature and star formation rate surface density maps, we are able to explore the radial distribution of the CO(3-2)/CO(1-0) integrated intensity ratio (R31) and its relationship with dust temperature and star formation. We find that the value of R31 between 2 - 9 kpc galactocentric radius is 0.14, significantly lower than what is seen in the nuclear ring at ~1 kpc (R31 ~ 0.8), only to rise again to 0.27 for the fields centred on the 10 kpc star forming ring. We also found that R31 is positively correlated with dust temperature, with Spearman's rank correlation coefficient $\rho$ = 0.55. The correlation between star formation rate surface density and CO(3--2) intensity is much stronger than with CO(1-0), with $\rho$ = 0.54 compared to -0.05, suggesting that the CO(3-2) line traces warmer and denser star forming gas better. We also find that R31 correlates well with star formation rate surface density, with $\rho$ = 0.69.

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.

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.

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

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.

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

The 31 brightest galaxies (m_B < 15 mag) inside the virial radius of the Fornax cluster were observed from the centres to the outskirts with the Multi Unit Spectroscopic Explorer on the Very Large Telescope. These observations provide detailed high-resolution maps of the line-of-sight kinematics and line strengths of the stars and ionised gas reaching 2-3 Re for 21 early-type galaxies and 1-2 Re for 10 late-type galaxies. The majority of the galaxies are regular rotators, with eight hosting a kinematically distinct core. Only two galaxies are slow rotators. The mean age, total metallicity, and [Mg/Fe] abundance ratio in the bright central region inside 0.5 Re and in the galaxy outskirts are presented. Extended emission-line gas is detected in 13 galaxies, most of them are late-type objects with wide-spread star formation. The measured structural properties are analysed in relation to the galaxies' position in the projected phase space of the cluster. This shows that the Fornax cluster appears to consist of three main groups of galaxies inside the virial radius: the old core; a clump of galaxies, which is aligned with the local large-scale structure and was accreted soon after the formation of the core; and a group of galaxies that fell in more recently.

In this work, we characterise the contributions from both ongoing star formation and the ambient radiation field in Local Group galaxy M33, as well as estimate the scale of the local dust-energy balance (i.e. the scale at which the dust is re-emitting starlight generated in that same region) in this galaxy through high-resolution radiative transfer (RT) modelling, with defined stellar and dust geometries. We have characterised the spectral energy distribution (SED) of M33 from UV to sub-mm wavelengths, at a spatial scale of 100 pc. We constructed input maps of the various stellar and dust geometries for use in the RT modelling. By modifying our dust mix (fewer very small carbon grains and a lower silicate-to-carbon ratio as compared to the Milky Way), we can much better fit the sub-mm dust continuum. Using this new dust composition, we find that we are able to well reproduce the observed SED of M33 using our adopted model. In terms of stellar attenuation by dust, we find a reasonably strong, broad UV bump, as well as significant systematic differences in the amount of dust attenuation when compared to standard SED modelling. We also find discrepancies in the residuals of the spiral arms versus the diffuse interstellar medium (ISM), indicating a difference in properties between these two regimes. The dust emission is dominated by heating due to the young stellar populations at all wavelengths ($\sim$80% at 10 $\mu$m to $\sim$50% at 1 mm). We find that the local dust-energy balance is restored at spatial scales greater than around 1.5 kpc.

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.

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.

We used deep MUSE observations to perform a stellar-kinematic and population analysis of FCC 153 and FCC 177, two edge-on S0 galaxies in the Fornax cluster. The geometrical definition of the different structural components of these two galaxies allows us to describe the nature of their thick disks. These are both old, relatively metal poor and [Mg/Fe]-enhanced, and their star formation history (SFH) reveals a minor younger component whose chemical properties suggest its later accretion. Moreover, the outer regions of these geometrically defined thick disks show higher values of metallicity and lower values of [Mg/Fe]. These stars probably formed in the thin-disk region and they were dynamically heated to form the flares present in these two galaxies. We propose different formation scenarios for the three populations of these thick disks: in-situ formation, accretion and disk heating. A clear distinction in age is found between the metal poor and [Mg/Fe]-enhanced thick disks (old, $\sim 12-13$ Gyr), and the metal rich and less [Mg/Fe]-enhanced thin disks (young, $\sim 4-5$ Gyr). These two galaxies show signs of relatively recent star formation in their thin disks and nuclear regions. While the thin disks show more continuous SFHs, the nuclei display a rather bursty SFH. These two galaxies are located outside of the densest region of the Fornax cluster where FCC 170 resides. This other edge-on S0 galaxy was studied by \citet{Pinna2019}. We compare and discuss our results with this previous study. The differences between these three galaxies, at different distances from the cluster center, suggest that the environment can have a strong effect on the galaxy evolutionary path.

The stellar initial mass function (IMF) regulates the baryonic cycle within galaxies, and is a key ingredient to translate observations into physical quantities. Although for decades it was assumed to be universal, there is now growing observational evidence showing that the center of massive early-type galaxies host an enhanced population of low-mass stars compared to the expectations from the Milky Way. Moreover, these variations in the IMF have been found to be related to the radial metallicity variations in massive galaxies. We present here a two-dimensional stellar population analysis of the massive lenticular galaxy FCC 167 (NGC 1380) as part of the Fornax3D project. Using a newly developed stellar population fitting scheme, we derive a full two-dimensional IMF map of an early-type galaxy. This two-dimensional analysis allows us go further than a radial analysis, showing how the metallicity changes along a disc-like structure while the IMF follows a distinct, less disky distribution. Thus, our findings indicate that metallicity cannot be the sole driver of the observed radial IMF variations. In addition, a comparison with the orbital decomposition shows suggestive evidence of a coupling between stellar population properties and the internal dynamical structure of FCC 167, where metallicity and IMF maps seem to track the distribution of cold and warm orbits, respectively.

We present the severe, systematic uncertainties currently facing our understanding of dust emission, which stymie our ability to truly exploit dust as a tool for studying galaxy evolution. We propose a program of study to tackle these uncertainties, describe the necessary facilities, and discuss the potential science gains that will result. This white paper was submitted to the US National Academies' Astro2020 Decadal Survey on Astronomy and Astrophysics.

One of the most important properties of a galaxy is the total stellar mass, or equivalently the stellar mass-to-light ratio (M/L). It is not directly observable, but can be estimated from stellar population synthesis. Currently, a galaxy's M/L is typically estimated from global fluxes. For example, a single global g - i colour correlates well with the stellar M/L. Spectral energy distribution (SED) fitting can make use of all available fluxes and their errors to make a Bayesian estimate of the M/L. We want to investigate the possibility of using morphology information to assist predictions of M/L. Our first goal is to develop and train a method that only requires a g-band image and redshift as input. This will allows us to study the correlation between M/L and morphology. Next, we can also include the i-band flux, and determine if morphology provides additional constraints compared to a method that only uses g- and i-band fluxes. We used a machine learning pipeline that can be split in two steps. First, we detected morphology features with a convolutional neural network. These are then combined with redshift, pixel size and g-band luminosity features in a gradient boosting machine. Our training target was the M/L acquired from the GALEX-SDSS-WISE Legacy Catalog, which uses global SED fitting and contains galaxies with z ~ 0.1. Morphology is a useful attribute when no colour information is available, but can not outperform colour methods on its own. When we combine the morphology features with global g- and i-band luminosities, we find an improved estimate compared to a model which does not make use of morphology. While our method was trained to reproduce global SED fitted M/L, galaxy morphology gives us an important additional constraint when using one or two bands. Our framework can be extended to other problems to make use of morphological information.

We present and discuss the stellar kinematics and populations of the S0 galaxy FCC 170 (NGC 1381) in the Fornax cluster, using deep MUSE data from the Fornax 3D survey. We show the maps of the first four moments of the stellar line-of-sight velocity distribution and of the mass-weighted mean stellar age, metallicity and [Mg/Fe] abundance ratio. The high-quality MUSE stellar kinematic measurements unveil the structure of this massive galaxy: a nuclear disk, a bar seen as a boxy bulge with a clear higher-velocity-dispersion X shape, a fast-rotating and flaring thin disk and a slower rotating thick disk. Whereas their overall old age makes it difficult to discuss differences in the formation epoch between these components, we find a clear-cut distinction between metal-rich and less [Mg/Fe]-enhanced populations in the thin-disk, boxy-bulge and nuclear disk, and more metal-poor and [Mg/Fe]-enhanced stars in the thick disk. Located in the densest region of the Fornax cluster, where signs of tidal stripping have been recently found, the evolution of FCC 170 might have been seriously affected by its environment. We discuss the possibility of its "pre-processing" in a subgroup before falling into the present-day cluster, which would have shaped this galaxy a long time ago. The thick disk displays a composite star formation history, as a significant fraction of younger stars co-exist with the main older thick-disk population. The former sub-population is characterized by even lower-metallicity and higher-[Mg/Fe] values, suggesting that these stars formed later and faster in a less chemically evolved satellite, which was subsequently accreted. Finally, we discuss evidence that metal-rich and less [Mg/Fe]-enhanced stars were brought in the outer parts of the thick disk by the flaring of the thin disk.

Galaxies continuously reprocess their interstellar material. One can therefore expect changing dust grain properties in galaxies which have followed different evolutionary pathways. Determining the intrinsic dust grain mix of a galaxy helps in reconstructing its evolutionary history. Early-type galaxies occasionally display regular dust lanes in their central regions. Due to the relatively simple geometry and composition of their stellar bodies, these galaxies are ideal to disentangle dust mix variations from geometric effects. We therefore model the various components of such a galaxy (FCC 167). We reconstruct its recent history, and investigate the possible fate of the dust lane. MUSE and ALMA observations reveal a nested ISM structure. An ionised-gas disk pervades the central regions of FCC 167, including those occupied by the main dust lane. Inward of the dust lane, we also find a disk/ring of cold molecular gas where stars are forming and HII regions contribute to the ionised-gas emission. Further in, the gas ionisation points towards an active galactic nucleus and the fuelling of a central supermassive black hole from its surrounding ionised and molecular reservoir. Observational constraints and radiative transfer models suggest the dust and gas are distributed in a ring-like geometry and the dust mix lacks small grains. The derived dust destruction timescales from sputtering in hot gas are short and we conclude that the dust must be strongly self-shielding and clumpy, or will quickly be eroded and disappear. Our findings show how detailed analysis of individual systems can complement statistical studies of dust-lane ETGs.

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

We have fit the far-ultraviolet (FUV) to sub-millimeter (850 micron) spectral energy distributions (SEDs) of the 61 galaxies from the "Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel" (KINGFISH). The fitting has been performed using three models: the Code for Investigating GALaxy Evolution (CIGALE), the GRAphite-SILicate approach (GRASIL), and the Multi-wavelength Analysis of Galaxy PHYSical properties (MAGPHYS). We have analyzed the results of the three codes in terms of the SED shapes, and by comparing the derived quantities with simple "recipes" for stellar mass (Mstar), star-formation rate (SFR), dust mass (Mdust), and monochromatic luminosities. Although the algorithms rely on different assumptions for star-formation history, dust attenuation and dust reprocessing, they all well approximate the observed SEDs and are in generally good agreement for the associated quantities. However, the three codes show very different behavior in the mid-infrared regime, in particular between 25 and 70 micron where there are no observational constraints for the KINGFISH sample. We find that different algorithms give discordant SFR estimates for galaxies with low specific SFR, and that the standard "recipes" for calculating FUV absorption overestimate the extinction compared to the SED-fitting results. Results also suggest that assuming a "standard" constant stellar mass-to-light ratio overestimates Mstar relative to the SED fitting, and we provide new SED-based formulations for estimating Mstar from WISE W1 (3.4 micron) luminosities and colors. From a Principal Component Analysis of Mstar, SFR, Mdust, and O/H, we reproduce previous scaling relations among Mstar, SFR, and O/H, and find that Mdust can be predicted to within roughly 0.3 dex using only Mstar and SFR.

Unresolved gas and dust observations show a surprising diversity in the amount of interstellar matter in early-type galaxies. Using ALMA observations we resolve the ISM in z$\sim$0.05 early-type galaxies. From a large sample of early-type galaxies detected in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) we selected five of the dustiest cases, with dust masses M$_d\sim$several$\times10^7$M$_\odot$, with the aim of mapping their submillimetre continuum and $^{12}$CO(2-1) line emission distributions. These observations reveal molecular gas disks. There is a lack of associated, extended continuum emission in these ALMA observations, most likely because it is resolved out or surface brightness limited, if the dust distribution is as extended as the CO gas. However, two galaxies have central continuum ALMA detections. An additional, slightly offset, continuum source is revealed in one case, which may have contributed to confusion in the Herschel fluxes. Serendipitous continuum detections further away in the ALMA field are found in another case. Large and massive rotating molecular gas disks are mapped in three of our targets, reaching a few$\times10^{9}$M$_\odot$. One of these shows evidence of kinematic deviations from a pure rotating disc. The fields of our two remaining targets contain only smaller, weak CO sources, slightly offset from the optical galaxy centres. These may be companion galaxies seen in ALMA observations, or background objects. These heterogeneous findings in a small sample of dusty early-type galaxies reveal the need for more such high spatial resolution studies, to understand statistically how dust and gas are related in early-type galaxies.

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.

The radio and far-infrared luminosities of star-forming galaxies are tightly correlated over several orders of magnitude; this is known as the far-infrared radio correlation (FIRC). Previous studies have shown that a host of factors conspire to maintain a tight and linear FIRC, despite many models predicting deviation. This discrepancy between expectations and observations is concerning since a linear FIRC underpins the use of radio luminosity as a star-formation rate indicator. Using LOFAR 150MHz, FIRST 1.4 GHz, and Herschel infrared luminosities derived from the new LOFAR/H-ATLAS catalogue, we investigate possible variation in the monochromatic (250$\mathrm{\mu m}$) FIRC at low and high radio frequencies. We use statistical techniques to probe the FIRC for an optically-selected sample of 4,082 emission-line classified star-forming galaxies as a function of redshift, effective dust temperature, stellar mass, specific star formation rate, and mid-infrared colour (an empirical proxy for specific star formation rate). Although the average FIRC at high radio frequency is consistent with expectations based on a standard power-law radio spectrum, the average correlation at 150MHz is not. We see evidence for redshift evolution of the FIRC at 150MHz, and find that the FIRC varies with stellar mass, dust temperature and specific star formation rate, whether the latter is probed using MAGPHYS fitting, or using mid-infrared colour as a proxy. We can explain the variation, to within 1$\sigma$, seen in the FIRC over mid-infrared colour by a combination of dust temperature, redshift, and stellar mass using a Bayesian partial correlation technique.

We present results of the detailed dust energy balance study for the seven large edge-on galaxies in the HEROES sample using 3D radiative transfer (RT) modelling. Based on available optical and near-infrared observations of the HEROES galaxies, we derive the 3D distribution of stars and dust in these galaxies. For the sake of uniformity, we apply the same technique to retrieve galaxy properties for the entire sample: we use a stellar model consisting of a S\'ersic bulge and three double-exponential discs (a superthin disc for a young stellar population and thin and thick discs for old populations). For the dust component, we adopt a double-exponential disc with the new THEMIS dust-grain model. We fit oligochromatic radiative transfer (RT) models to the optical and near-infrared images with the fitting algorithm FitSKIRT and do panchromatic simulations with the SKIRT code at wavelengths ranging from ultraviolet to submillimeter. We confirm the previously stated dust energy balance problem in galaxies: for the HEROES galaxies, the dust emission derived from our RT calculations underestimates the real observations by a factor 1.5-4 for all galaxies except NGC 973 and NGC 5907 (apparently, the latter galaxy has a more complex geometry than we used). The comparison between our RT simulations and the observations at mid-infrared-submillimeter wavelengths shows that most of our galaxies exhibit complex dust morphologies (possible spiral arms, star-forming regions, more extended dust structure in the radial and vertical directions). We suggest that, in agreement with the results from Saftly et al. (2015), the large- and small-scale structure is the most probable explanation for the dust energy balance problem.

Studies both of entire galaxies and of local Galactic star formation indicate a dependency of a molecular cloud's star formation rate (SFR) on its dense gas mass. In external galaxies, such measurements are derived from HCN(1-0) observations, usually encompassing many Giant Molecular Clouds (GMCs) at once. The Andromeda galaxy (M31) is a unique laboratory to study the relation of the SFR and HCN emission down to GMC scales at solar-like metallicities. In this work, we correlate our composite SFR determinations with archival HCN, HCO$^+$, and CO observations, resulting in a sample of nine reasonably representative GMCs. We find that, at the scale of individual clouds, it is important to take into account both obscured and unobscured star formation to determine the SFR. When correlated against the dense-gas mass from HCN, we find that the SFR is low, in spite of these refinements. We nevertheless retrieve an SFR - dense-gas mass correlation, confirming that these SFR tracers are still meaningful on GMC scales. The correlation improves markedly when we consider the HCN/CO ratio instead of HCN by itself. This nominally indicates a dependency of the SFR on the dense-gas fraction, in contradiction to local studies. However, we hypothesize that this partly reflects the limited dynamic range in dense-gas mass, and partly that the ratio of single-pointing HCN and CO measurements may be less prone to systematics like sidelobes. In this case, the HCN/CO ratio would importantly be a better empirical measure of the dense-gas content itself.

We present the largest submillimeter images that have been made of the extragalactic sky. The 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{\mu}m. In this paper we present the images from our two largest fields which account for ~75% of the survey. The first field is 180.1 deg$^2$ in size centered on the North Galactic Pole (NGP) and the second field is 317.6 deg$^2$ in size centered on the South Galactic Pole. The NGP field serendipitously contains the Coma cluster. Over most (~80%) of the images, the pixel noise, including both instrumental noise and confusion noise, is approximately 3.6, and 3.5 mJy/pix at 100 and 160{\mu}m, and 11.0, 11.1 and 12.3 mJy/beam at 250, 350 and 500{\mu}m, respectively, but reaches lower values in some parts of the images. If a matched filter is applied to optimize point-source detection, our total 1{\sigma} map sensitivity is 5.7, 6.0, and 7.3 mJy at 250, 350, and 500{\mu}m, respectively. We describe the results of an investigation of the noise properties of the images. We make the most precise estimate of confusion in SPIRE maps to date finding values of 3.12+/-0.07, 4.13+/-0.02 and 4.45+/-0.04 mJy/beam at 250, 350, and 500{\mu}m in our un-convolved maps. For PACS we find an estimate of the confusion noise in our fast-parallel observations of 4.23 and 4.62 mJy/beam at 100 and 160{\mu}m. Finally, we give recipes for using these images to carry out photometry, both for unresolved and extended sources.

We present a kinematic analysis of the dust-lane elliptical NGC 5626 based on MUSE observations. These data allow to robustly classify this galaxy as a fast rotator and to infer a virial mass of $10^{11.7} M_\odot$, making it one of the most massive fast rotators known. In addition, the depth and extent of the MUSE data reveal a strong kinematic twist in the stellar velocity field (by up to $45$ degrees beyond $1.5R_e$). A comparison with the ATLAS$^\mathrm{3D}$ sample underlines the rareness of this system, although we show that such a large-scale kinematic twist could have been missed by the ATLAS$^\mathrm{3D}$ data due to the limited spatial sampling of this survey (typically extending to $0.6R_e$ for massive ETGs). MUSE thus has the potential to unveil more examples of this type of galaxies. We discuss the environment and possible formation history of NGC 5626 and finally argue how a merger between the Milky Way and Andromeda could produce a galaxy of the same class as NGC 5626.

The Herschel Space Observatory has revealed a very different galaxyscape from that shown by optical surveys which presents a challenge for galaxy-evolution models. The Herschel surveys reveal (1) that there was rapid galaxy evolution in the very recent past and (2) that galaxies lie on a a single Galaxy Sequence (GS) rather than a star-forming `main sequence' and a separate region of `passive' or `red-and-dead' galaxies. The form of the GS is now clearer because far-infrared surveys such as the Herschel ATLAS pick up a population of optically-red star-forming galaxies that would have been classified as passive using most optical criteria. The space-density of this population is at least as high as the traditional star-forming population. By stacking spectra of H-ATLAS galaxies over the redshift range 0.001 < z < 0.4, we show that the galaxies responsible for the rapid low-redshift evolution have high stellar masses, high star-formation rates but, even several billion years in the past, old stellar populations - they are thus likely to be relatively recent ancestors of early-type galaxies in the Universe today. The form of the GS is inconsistent with rapid quenching models and neither the analytic bathtub model nor the hydrodynamical EAGLE simulation can reproduce the rapid cosmic evolution. We propose a new gentler model of galaxy evolution that can explain the new Herschel results and other key properties of the galaxy population.

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

The newest generation of integral field unit spectrographs brings three-dimensional mapping of nearby galaxies one step closer. While the focus up to this point was mostly on stars and ionised gas, it is also possible to look at dust in a new, more complete way. Using MUSE science verification observations of NGC 5626, we map the interstellar matter in this dusty lenticular. We use the resolving power of MUSE to measure the optical attenuation with a spectral resolution of 6.25 \AA, at physical scales of 0.1-1 kpc. The integrated attenuation curve of NGC 5626 shows a smooth, slightly steeper than Milky Way and SMC attenuation curves. Several sharp features are superimposed: we measure lower attenuation at spectral emission lines and higher attenuation for the sodium line doublet. No correlation was observed between sodium line strength and reddening by dust on spatially resolved scales. Additionally, the continuum attenuation was found to be independent from the Balmer decrement (tracing ionised gas attenuation). We model and interpret the variations in the attenuation curves of each spatial resolution element of NGC 5626. We find that the amount and distribution of dust along the line-of-sight is highly degenerate with any variation in the intrinsic extinction law. Our analysis shows that the interstellar matter in NGC 5626 resides in a regular and well-settled disk. Our results preach caution in the application of simple recipes to de-redden global galaxy spectra and underlines the need for more realistic dust geometries when constructing such correction formulas.

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

We present the calibration between the dust continuum luminosity and interstellar gas content obtained from the Valpara\'{i}so ALMA Line Emission Survey (VALES) sample of 67 main-sequence star-forming galaxies at 0.02<$z$<0.35. We use CO(1-0) observations from the Atacama Large Millimetre/submillimetre Array (ALMA) to trace the molecular gas mass, $M_{\mathrm{H}_{2}}$, and estimate the rest-frame monochromatic luminosity at 850 $\mu$m, $L_{\nu_{850}}$, by extrapolating the dust continuum from MAGPHYS modelling of the far-ultraviolet to submillimetre spectral energy distribution sampled by the Galaxy And Mass Assembly (GAMA) survey. Adopting $\alpha_{\rm CO}$ = 6.5 (K km s$^{-1}$ pc$^{2}$)$^{-1}$, the average ratio of $L_{\nu_{850}}/M_{\mathrm{H}_{2}}$ = (6.4$\pm$1.4)$\times10^{19}$ erg s$^{-1}$ Hz$^{-1}$ $\mathrm{M}_{\odot}^{-1}$, in excellent agreement with literature values. We obtain a linear fit of $\log_{10}$ ($M_{\mathrm{H}_{2}}/\mathrm{M}_{\odot}$) = (0.92$\pm$0.02) $\log_{10}$ ($L_{\nu_{850}}$/erg s$^{-1}$ Hz$^{-1}$)-(17.31$\pm$0.59). We provide relations between $L_{\nu_{850}}$, $M_{\mathrm{H}_{2}}$ and $M_{\mathrm{ISM}}$ when combining the VALES and literature samples, and adopting a Galactic $\alpha_{\rm CO}$ value.

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

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

This paper is the second in a pair of articles presenting data release 1 (DR1) of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), the largest single open-time key project carried out with the Herschel Space Observatory. The H-ATLAS is a wide-area imaging survey carried out in five photometric bands at 100, 160, 250, 350 and 500$\mu$m covering a total area of 600deg$^2$. In this paper we describe the identification of optical counterparts to submillimetre sources in DR1, comprising an area of 161 deg$^2$ over three equatorial fields of roughly 12$^\circ$x4.5$^\circ$ centred at 9$^h$, 12$^h$ and 14.5$^h$ respectively. Of all the H-ATLAS fields, the equatorial regions benefit from the greatest overlap with current multi-wavelength surveys spanning ultraviolet (UV) to mid-infrared regimes, as well as extensive spectroscopic coverage. We use a likelihood-ratio technique to identify SDSS counterparts at r<22.4 for 250-$\mu$m-selected sources detected at $\geq$ 4$\sigma$ ($\approx$28mJy). We find `reliable' counterparts (reliability R$\geq$0.8) for 44,835 sources (39 per cent), with an estimated completeness of 73.0 per cent and contamination rate of 4.7 per cent. Using redshifts and multi-wavelength photometry from GAMA and other public catalogues, we show that H-ATLAS-selected galaxies at $z<0.5$ span a wide range of optical colours, total infrared (IR) luminosities, and IR/UV ratios, with no strong disposition towards mid-IR-classified AGN in comparison with optical selection. The data described herein, together with all maps and catalogues described in the companion paper (Valiante et al. 2016), are available from the H-ATLAS website at www.h-atlas.org.

In most studies of dust in galaxies, dust is only detected from its emission to approximately the optical radius of the galaxy. By combining the signal of 110 spiral galaxies observed as part of the Herschel Reference Survey, we are able to improve our sensitivity by an order-of-magnitude over that for a single object. Here we report the direct detection of dust from its emission that extends out to at least twice the optical radius. We find that the distribution of dust is consistent with an exponential at all radii with a gradient of ~-1.7 dex R$_{25}^{-1}$. Our dust temperature declines linearly from ~25 K in the centre to 15 K at R$_{25}$ from where it remains constant out to ~2.0 R$_{25}$. The surface-density of dust declines with radius at a similar rate to the surface-density of stars but more slowly than the surface-density of the star-formation rate. Studies based on dust extinction and reddening of high-redshift quasars have concluded that there are substantial amounts of dust in intergalactic space. By combining our results with the number counts and angular correlation function from the SDSS, we show that with Milky Way type dust we can explain the reddening of the quasars by the dust within galactic disks alone. Given the uncertainties in the properties of any intergalactic dust, we cannot rule out its existence, but our results show that statistical investigations of the dust in galactic halos that use the reddening of high-redshift objects must take account of the dust in galactic disks.

We investigate the dust energy balance for the edge-on galaxy IC 2531, one of the seven galaxies in the HEROES sample. We perform a state-of-the-art radiative transfer modelling based, for the first time, on a set of optical and near-infrared galaxy images. We show that taking into account near-infrared imaging in the modelling significantly improves the constraints on the retrieved parameters of the dust content. We confirm the result from previous studies that including a young stellar population in the modelling is important for explaining the observed stellar energy distribution. However, the discrepancy between the observed and modelled thermal emission at far-infrared wavelengths, the so-called dust energy balance problem, is still present: the model underestimates the observed fluxes by a factor of about two. We compare two different dust models, and find that dust parameters and thus the spectral energy distribution in the infrared domain are sensitive to the adopted dust model. In general, the THEMIS model reproduces the observed emission in the infrared wavelength domain better than the popular Zubko et al. BARE-GR-S model. Our study of IC 2531 is a pilot case for detailed and uniform radiative transfer modelling of the entire HEROES sample, which will shed more light on the strength and origins of the dust energy balance problem.

We combine dust, atomic (HI) and molecular (H$_{2}$) hydrogen mass measurements for 176 galaxies in the Herschel Reference Survey to investigate the effect of environment on the gas-to-dust mass ($M_{\rm gas}/M_{\rm dust}$) ratio of nearby galaxies. We find that, at fixed stellar mass, the average $M_{\rm gas}/M_{\rm dust}$ ratio varies by no more than a factor of $\sim$2 when moving from field to cluster galaxies, with Virgo galaxies being slightly more dust rich (per unit of gas) than isolated systems. Remarkably, once the molecular and atomic hydrogen phases are investigated separately, we find that \hi-deficient galaxies have at the same time lower $M_{\rm HI}/M_{\rm dust}$ ratio but higher $M_{\rm H_{2}}/M_{\rm dust}$ ratio than \hi-normal systems. In other words, they are poorer in atomic but richer in molecular hydrogen if normalized to their dust content. By comparing our findings with the predictions of theoretical models, we show that the opposite behavior observed in the $M_{\rm HI}/M_{\rm dust}$ and $M_{\rm H_{2}}/M_{\rm dust}$ ratios is fully consistent with outside-in stripping of the interstellar medium (ISM), and is simply a consequence of the different distribution of dust, \hi\ and H$_{2}$ across the disk. Our results demonstrate that the small environmental variations in the total $M_{\rm gas}/M_{\rm dust}$ ratio, as well as in the gas-phase metallicity, do not automatically imply that environmental mechanisms are not able to affect the dust and metal content of the ISM in galaxies.

The Star formation rate (SFR) is a crucial parameter to investigate galaxy evolution. At low redshift the cosmic SFR density declines smoothly, and massive active galaxies become passive, reducing their star formation activity. This implies that the bulk of the SFR density at low redshift is mainly driven by low mass objects. We investigate the properties of a sample of low luminosity Far-Infrared (FIR) sources selected at 250 microns from Pappalardo et al. (2015). We have collected data from Ultraviolet to FIR to perform a multi-wavelengths analysis. The main goal is to investigate the correlation between SFR, stellar mass, and dust mass for a galaxy population with a wide range in dust content and stellar mass, including the low mass regime that most probably dominates the SFR density at low z. We define a main sample of ~800 sources with full Spectral Energy Distribution (SED) coverage between 0.15 < lambda < 500 microns and an extended sample with ~5000 sources in which we remove the constraints on the Ultraviolet and Near-Infrared bands. We analyze both samples with two different SED fitting methods: MAGPHYS and CIGALE. In the SFR versus stellar mass plane our samples occupy a region included between local spirals and higher redshift star forming galaxies. The galaxies subsample with the higher masses (M* > 3e10 Msol) does not lie on the main sequence, but shows a small offset, as a consequence of the decreased star formation. Low mass galaxies (M* < 1e10 Msol) settle in the main sequence with SFR and stellar mass consistent with local spirals. Deep Herschel data allow the identification of a mixed galaxy population, with galaxies still in an assembly phase, or galaxies at the beginning of their passive evolution. We find that the dust luminosity is the parameter that discriminates these two galaxy populations.

Thanks to deep UV observations with GALEX and Swift, diffuse UV haloes have recently been discovered around galaxies. Based on UV-optical colours, it has been advocated that the UV haloes around spiral galaxies are due to UV radiation emitted from the disc and scattered off dust grains at high latitudes. Detailed UV radiative transfer models that take into account scattering and absorption can explain the morphology of the UV haloes, and they require the presence of an additional thick dust disc next the to traditional thin disc for half of the galaxies in their sample. We test whether such an additional thick dust disc agrees with the observed infrared emission in NGC 3628, an edge-on galaxy with a clear signature of a thick dust disc. We extend the far-ultraviolet radiative transfer models to full-scale panchromatic models. Our model, which contains no fine-tuning, can almost perfectly reproduce the observed spectral energy distribution from UV to mm wavelengths. These results corroborate the interpretation of the extended UV emission in NGC 3628 as scattering off dust grains, and hence of the presence of a substantial amount of diffuse extra-planar dust. A significant caveat, however, is the geometrical simplicity and non-uniqueness of our model: other models with a different geometrical setting could lead to a similar spectral energy distribution. More detailed radiative transfer simulations that compare the model results to images from UV to submm wavelengths are a way to break this degeneracy, as are UV polarisation measurements.

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

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

[Abridged] In high density environment, the gas content of galaxies is stripped, leading to a rapid quenching of their star formation activity. This dramatic environmental effect is generally not taken into account in the SFHs usually assumed to perform spectral energy distribution (SED) fitting of these galaxies, yielding to a poor fit of their stellar emission and, consequently, a biased estimate of the SFR. We aim at reproducing the SFH of galaxies that underwent a rapid star formation quenching using a truncated delayed SFH that we implemented in the SED fitting code CIGALE. We show that the ratio between the instantaneous SFR and the SFR just before the quenching ($r_{SFR}$) is well constrained as long as rest frame UV data are available. This SED modelling is applied to the Herschel Reference Survey (HRS) containing isolated galaxies and sources falling in the dense environment of the Virgo cluster. The latter are HI-deficient due to ram pressure stripping. We show that the truncated delayed SFH successfully reproduces their SED while typical SFH assumptions fail. A good correlation is found between $r_{SFR}$ and HI-def, the parameter quantifying the gas deficiency of cluster galaxies, meaning that SED fitting results can be used to provide a tentative estimate of the gas deficiency of galaxies for which HI observations are not available. The HRS galaxies are placed on the SFR-$M_*$ diagram showing that the HI-deficient sources lie in the quiescent region confirming previous studies. Using the $r_{SFR}$ parameter, we derive the SFR of these sources before quenching and show that they were previously on the main sequence relation. We show that the $r_{SFR}$ parameter is also well recovered for deeply obscured high redshift sources, as well as in absence of IR data. SED fitting is thus a powerful tool to identify galaxies that underwent a rapid star formation quenching.

Context. Edge-on galaxies can offer important insights in galaxy evolution as they are the only systems where the distribution of the different components can be studied both radially and vertically. The HEROES project was designed to investigate the interplay between the gas, dust, stars and dark matter (DM) in a sample of 7 massive edge-on spiral galaxies. Aims. In this second HEROES paper we present an analysis of the atomic gas content of 6 out of 7 galaxies in our sample. The remaining galaxy was recently analysed according to the same strategy. The primary aim of this work is to constrain the surface density distribution, the rotation curve and the geometry of the gas disks in a homogeneous way. In addition we identify peculiar features and signs of recent interactions. Methods. We construct detailed tilted-ring models of the atomic gas disks based on new GMRT 21-cm observations of NGC 973 and UGC 4277 and re-reduced archival HI data of NGC 5907, NGC 5529, IC 2531 and NGC 4217. Potential degeneracies between different models are resolved by requiring a good agreement with the data in various representations of the data cubes. Results. From our modelling we find that all but one galaxy are warped along the major axis. In addition, we identify warps along the line of sight in three galaxies. A flaring gas layer is required to reproduce the data only for one galaxy, but (moderate) flares cannot be ruled for the other galaxies either. A coplanar ring-like structure is detected outside the main disk of NGC 4217, which we suggest could be the remnant of a recent minor merger event. We also find evidence for a radial inflow of 15 +- 5 km/s in the disk of NGC 5529, which might be related to the ongoing interaction with two nearby companions. (Abridged)

A fraction of the early-type galaxy population hosts a prominent dust lane. Methods to quantify the dust content of these systems based on optical imaging data usually yield dust masses which are an order of magnitude lower than dust masses derived from the observed FIR emission. High-quality optical data from the Next Generation Virgo cluster Survey (NGVS) and FIR/submm observations from the Herschel Virgo Cluster Survey (HeViCS) allow us to revisit previous methods to determine the dust content in galaxies and explore new ones. We aim to derive the dust mass in NGC 4370 from both optical and FIR data, and investigate the need to invoke a putative diffuse dust component. We create color and attenuation maps, which are converted to approximate dust mass maps based on simple dust geometries. Dust masses are also derived from SED fits to FIR/submm observations. Finally, inverse radiative transfer fitting is performed to investigate more complex dust geometries. The empirical methods applied to the optical data yield lower limits of 3.4e5 solar masses, an order of magnitude below the total dust masses derived from SED fitting. In contrast, radiative transfer models yield dust masses which are slightly lower, but fully consistent with the FIR-derived mass. Dust is more likely to be distributed in a ring around the centre of NGC 4370 as opposed to an exponential disc or a simple foreground screen. Moreover, using inverse radiative transfer fitting, we are able to constrain most of the parameters describing these geometries. The resulting dust masses are high enough to account for the dust observed at FIR/submm wavelengths, so that no diffuse dust component needs to be invoked. We furthermore caution for the interpretation of dust masses and optical depths based on optical data alone, using overly simplistic star-dust geometries and the neglect of scattering effects. [ABRIDGED]

The Herschel Space Observatory has had a tremendous impact on the study of extragalactic dust. Specifically, early-type galaxies (ETG) have been the focus of several studies. In this paper we combine results from two Herschel studies - a Virgo cluster study HeViCS and a broader, low-redshift H-ATLAS/GAMA study - and contrast the dust and associated properties for similar mass galaxies. This comparison is motivated by differences in results exhibited between multiple Herschel studies of early-type galaxies. A comparison between consistent modified blackbody derived dust mass is carried out, revealing strong differences between the two samples in both dust mass and dust-to-stellar mass ratio. In particular, the HeViCS sample lacks massive ETG with as high a specific dust content as found in H-ATLAS. This is most likely connected with the difference in environment for the two samples. We calculate nearest neighbour environment densities in a consistent way, showing that H-ATLAS ETG occupy sparser regions of the local Universe, whereas HeViCS ETG occupy dense regions. This is also true for ETG that are not Herschel-detected but are in the Virgo and GAMA parent samples. Spectral energy distributions are fit to the panchromatic data. From these we find that in H-ATLAS the specific star formation rate anticorrelates with stellar mass and reaches values as high as in our Galaxy. On the other hand HeViCS ETG appear to have little star formation. Based on the trends found here, H-ATLAS ETG are thought to have more extended star formation histories and a younger stellar population than HeViCS ETG.

Interstellar dust in galaxies can be traced either through its extinction effects on the star light, or through its thermal emission at infrared wavelengths. Recent radiative transfer studies of several nearby edge-on galaxies have found an apparent inconsistency in the dust energy balance: the radiative transfer models that successfully explain the optical extinction underestimate the observed fluxes by an average factor of three. We investigate the dust energy balance for IC4225 and NGC5166, two edge-on spiral galaxies observed by the Herschel Space Observatory in the frame of the H-ATLAS survey. We start from models which were constrained from optical data and extend them to construct the entire spectral energy distribution of our galaxies. These predicted values are subsequently compared to the observed far-infrared fluxes. We find that including a young stellar population in the modelling is necessary as it plays a non-negligible part in the heating of the dust grains. While the modelling approach for both galaxies is nearly identical, we find two very different results. As is often seen in other edge-on spiral galaxies, the far-infrared emission of our radiative transfer model of IC4225 underestimates the observed fluxes by a factor of about three. For NGC5166 on the other hand, we find that both the predicted spectral energy distribution as well as the simulated images match the observations particularly well. We explore possible reasons for this difference and conclude that it is unlikely that one single mechanism is the cause of the dust energy balance problem in spiral galaxies. We discuss the different approaches that can be considered in order to get a conclusive answer on the origin this discrepancy.

We present the rest-frame optical spectral energy distribution and stellar masses of six Herschel- selected gravitationally lensed dusty, star-forming galaxies (DSFGs) at 1 < z < 3. These galaxies were first identified with Herschel/SPIRE imaging data from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). The targets were observed with Spitzer/IRAC at 3.6 and 4.5um. Due to the spatial resolution of the IRAC observations at the level of 2 arcseconds, the lensing features of a background DSFG in the near-infrared are blended with the flux from the foreground lensing galaxy in the IRAC imaging data. We make use of higher resolution Hubble/WFC3 or Keck/NIRC2 Adaptive Optics imaging data to fit light profiles of the foreground lensing galaxy (or galaxies) as a way to model the foreground components, in order to successfully disentangle the foreground lens and background source flux densities in the IRAC images. The flux density measurements at 3.6 and 4.5um, once combined with Hubble/WFC3 and Keck/NIRC2 data, provide important constraints on the rest-frame optical spectral energy distribution of the Herschel-selected lensed DSFGs. We model the combined UV- to millimeter-wavelength SEDs to establish the stellar mass, dust mass, star-formation rate, visual extinction, and other parameters for each of these Herschel-selected DSFGs. These systems have inferred stellar masses in the range 8 x 10^10 to 4 x 10^11 Msun and star-formation rates of around 100 Msun yr-1. This puts these lensed sub-millimeter systems well above the SFR-M* relation observed for normal star-forming galaxies at similar redshifts. The high values of SFR inferred for these systems are consistent with a major merger-driven scenario for star formation.

We present the stellar and gaseous kinematics of an Sb galaxy, NGC 3223, with the aim of determining the vertical and radial stellar velocity dispersion as a function of radius, which can help to constrain disk heating theories. Together with the observed NIR photometry, the vertical velocity dispersion is also used to determine the stellar mass-to-light (M/L) ratio, typically one of the largest uncertainties when deriving the dark matter distribution from the observed rotation curve. We find a vertical-to-radial velocity dispersion ratio of sigma_z/sigma_R=1.21+-0.14, significantly higher than expectations from known correlations, and a weakly-constrained Ks-band stellar M/L ratio in the range 0.5-1.7, at the high end of (but consistent with) the predictions of stellar population synthesis models. Such a weak constraint on the stellar M/L ratio, however, does not allow us to securely determine the dark matter density distribution. To achieve this, either a statistical approach or additional data (e.g. integral-field unit) are needed.

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

We present Herschel PACS and SPIRE spectroscopy of the most important far-infrared cooling lines in the nearby edge-on spiral galaxy, NGC 891: [CII] 158 $\mu$m, [NII] 122, 205 $\mu$m, [OI] 63, 145 $\mu$m, and [OIII] 88 $\mu$m. We find that the photoelectric heating efficiency of the gas, traced via the ([CII]+[OII]63)/$F_{\mathrm{TIR}}$ ratio, varies from a mean of 3.5$\times$10$^{-3}$ in the centre up to 8$\times$10$^{-3}$ at increasing radial and vertical distances in the disc. A decrease in ([CII]+[OII]63)/$F_{\mathrm{TIR}}$ but constant ([CII]+[OI]63)/$F_{\mathrm{PAH}}$ with increasing FIR colour suggests that polycyclic aromatic hydrocarbons (PAHs) may become important for gas heating in the central regions. We compare the observed flux of the FIR cooling lines and total IR emission with the predicted flux from a PDR model to determine the gas density, surface temperature and the strength of the incident far-ultraviolet (FUV) radiation field, $G_{0}$. Resolving details on physical scales of ~0.6 kpc, a pixel-by-pixel analysis reveals that the majority of the PDRs in NGC 891's disc have hydrogen densities of 1 < log ($n$/cm$^{-3}$) < 3.5 experiencing an incident FUV radiation field with strengths of 1.7 < log $G_0$ < 3. Although these values we derive for most of the disc are consistent with the gas properties found in PDRs in the spiral arms and inter-arm regions of M51, observed radial trends in $n$ and $G_0$ are shown to be sensitive to varying optical thickness in the lines, demonstrating the importance of accurately accounting for optical depth effects when interpreting observations of high inclination systems. With an empirical relationship between the MIPS 24 $\mu$m and [NII] 205 $\mu$m emission, we estimate an enhancement of the FUV radiation field strength in the far north-eastern side of the disc.

We present three independent catalogs of point-sources extracted from SPIRE images at 250, 350, and 500 micron as a part of the Herschel Virgo Cluster Survey (HeViCS). The source positions are determined by estimating the likelihood to be a real source for each peak on the maps and the flux densities are estimated using the sourceExtractorTimeline, a timeline-based point source fitter. Afterwards, each source is subtracted from the maps, removing a Gaussian function in every position with the full width half maximum equal to that estimated in sourceExtractorTimeline. This procedure improves the robustness of our algorithm in terms of source identification. The HeViCS catalogs contain about 52000, 42200, and 18700 sources selected at 250, 350, and 500 micron above 3sigma and are ~ 75%, 62%, and 50% complete at flux densities of 20 mJy at 250, 350, 500 micron, respectively. We then measured source number counts at 250, 350, and 500 micron and we also cross-correlated the catalogs with the Sloan Digital Sky Survey to investigate the redshift distribution of the nearby sources. From this cross-correlation, we select ~2000 sources with reliable fluxes and a high signal-to-noise ratio, finding an average redshift z~0.3+/-0.22. The number counts at 250, 350, and 500 micron show an increase in the slope below 200 mJy, indicating a strong evolution in number of density for galaxies at these fluxes. In general, models tend to overpredict the counts at brighter flux densities, underlying the importance of studying the Rayleigh-Jeans part of the spectral energy distribution to refine the theoretical recipes of the models. Our iterative method for source identification allowed the detection of a family of 500 micron sources that are not foreground objects belonging to Virgo and not found in other catalogs.

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.

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

We examined variations in the 160/250 and 250/350 micron surface brightness ratios within 24 nearby (<30 Mpc) face-on spiral galaxies observed with the Herschel Space Observatory to identify the heating mechanisms for dust emitting at these wavelengths. The analysis consisted of both qualitative and quantitative comparisons of the 160/250 and 250/350 micron ratios to H alpha and 24 micron surface brightnesses, which trace the light from star forming regions, and 3.6 micron emission, which traces the light from the older stellar populations of the galaxies. We find broad variations in the heating mechanisms for the dust. In one subset of galaxies, we found evidence that emission at <=160 microns (and in rare cases potentially at <=350 microns) originates from dust heated by star forming regions. In another subset, we found that the emission at >=250 microns (and sometimes at >=160 microns) originates from dust heated by the older stellar population. In the rest of the sample, either the results are indeterminate or both of these stellar populations may contribute equally to the global dust heating. The observed variations in dust heating mechanisms does not necessarily match what has been predicted by dust emission and radiative transfer models, which could lead to overestimated dust temperatures, underestimated dust masses, false detections of variability in dust emissivity, and inaccurate star formation rate measurements.

The imprints of dust-starlight interactions are visible in scaling relations between stellar characteristics, star formation parameters and dust properties. We aim to examine dust scaling relations on a sub-kpc resolution in the Andromeda galaxy (M31) by comparing the properties on a local and global scale to other galaxies of the local universe. New Herschel observations are combined with available data from GALEX, SDSS, WISE and Spitzer to construct a dataset covering UV to submm wavelengths. We work at the resolution of the SPIRE $500\; \mu$m beam, with pixels corresponding to physical regions of 137 x 608 pc in the galaxy's disk. A panchromatic spectral energy distribution was modelled for each pixel and several dust scaling relations are investigated. We find, on a sub-kpc scale, strong correlations between $M_d/M_\star$ and NUV-r, and between $M_d/M_\star$ and $\mu_\star$ (the stellar mass surface density). Striking similarities with corresponding relations based on integrated galaxies are found. We decompose M31 in four macro-regions based on their FIR morphology; the bulge, inner disk, star forming ring and the outer disk. All regions closely follow the galaxy-scale average trends. The specific star formation characteristics we derive for these macro-regions give strong hints of an inside-out formation of the bulge-disk geometry, as well as an internal downsizing process. However, within each macro-region, a great diversity in individual micro-regions is found. Furthermore, we confirm that dust in the bulge of M31 is heated only by the old stellar populations. In general, the local dust scaling relations indicate that the dust content in M31 is maintained by a subtle interplay of past and present star formation. The similarity with galaxy-based relations strongly suggests that they are in situ correlations, with underlying processes that must be local in nature. (Abriged)

Evidence is mounting that a significant fraction of the early-type galaxy population contains substantial reservoirs of cold interstellar gas and dust. We investigate the gas and dust in NGC 5485, an early-type galaxy with a prominent minor-axis dust lane. Using new Herschel PACS and SPIRE imaging data, we detect 3.8 x 10^6 Msun of cool interstellar dust in NGC 5485, which is in stark contrast with the non-detection of the galaxy in sensitive HI and CO observations from the ATLAS3D consortium. The resulting gas-to-dust ratio upper limit is Mgas/Md < 14.5, almost an order of magnitude lower than the canonical value for the Milky Way. We scrutinize the reliability of the dust, atomic gas and molecular gas mass estimates, but these do not show systematic uncertainties that can explain the extreme gas-to-dust ratio. Also a warm or hot ionized gas medium does not offer an explanation. A possible scenario could be that NGC 5485 merged with an SMC-type metal-poor galaxy with a substantial CO-dark molecular gas component and that the bulk of atomic gas was lost during the interaction, but it remains to be investigated whether such a scenario is possible.

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

We investigate the amount and spatial distribution of interstellar dust in edge-on spiral galaxies, using detailed radiative transfer modeling of a homogeneous sample of 12 galaxies selected from the CALIFA survey. Our automated fitting routine, FitSKIRT, was first validated against artificial data. This is done by simultaneously reproducing the SDSS $g$-, $r$-, $i$- and $z$-band observations of a toy model in order to combine the information present in the different bands. We show that this combined, oligochromatic fitting, has clear advantages over standard monochromatic fitting especially regarding constraints on the dust properties. We model all galaxies in our sample using a three-component model, consisting of a double exponential disc to describe the stellar and dust discs and using a S\'ersic profile to describe the central bulge. The full model contains 19 free parameters, and we are able to constrain all these parameters to a satisfactory level of accuracy without human intervention or strong boundary conditions. Apart from two galaxies, the entire sample can be accurately reproduced by our model. We find that the dust disc is about 75% more extended but only half as high as the stellar disc. The average face-on optical depth in the V-band is $0.76$ and the spread of $0.60$ within our sample is quite substantial, which indicates that some spiral galaxies are relatively opaque even when seen face-on.

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.

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

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

We cross-correlate the Planck Catalogue of Compact Sources (PCCS) with the fully sampled 84 deg2 Herschel Virgo Cluster Survey (HeViCS) fields. We search for and identify the 857 and 545 GHz PCCS sources in the HeViCS fields by studying their FIR/submm and optical counterparts. We find 84 and 48 compact Planck sources in the HeViCS fields at 857 and 545 GHz, respectively. Almost all sources correspond to individual bright Virgo Cluster galaxies. The vast majority of the Planck detected galaxies are late-type spirals, with the Sc class dominating the numbers, while early-type galaxies are virtually absent from the sample, especially at 545 GHz. We compare the HeViCS SPIRE flux densities for the detected galaxies with the four different PCCS flux density estimators and find an excellent correlation with the aperture photometry flux densities, even at the highest flux density levels. We find only seven PCCS sources in the HeViCS fields without a nearby galaxy as obvious counterpart, and conclude that all of these are dominated by Galactic cirrus features or are spurious detections. No Planck sources in the HeViCS fields seem to be associated to high-redshift proto-clusters of dusty galaxies or strongly lensed submm sources. Finally, our study is the first empirical confirmation of the simulation-based estimated completeness of the PCCS, and provides a strong support of the internal PCCS validation procedure.

Context. Edge-on spiral galaxies with prominent dust lanes provide us with an excellent opportunity to study the distribution and properties of the dust within them. The HEROES project was set up to observe a sample of seven large edge-on galaxies across various wavelengths for this investigation. Aims. Within this first paper, we present the Herschel observations and perform a qualitative and quantitative analysis on them, and we derive some global properties of the far infrared and submillimetre emission. Methods. We determine horizontal and vertical profiles from the Herschel observations of the galaxies in the sample and describe the morphology. Modified black-body fits to the global fluxes, measured using aperture photometry, result in dust temperatures and dust masses. The latter values are compared to those that are derived from radiative transfer models taken from the literature. Results. On the whole, our Herschel flux measurements agree well with archival values. We find that the exponential horizontal dust distribution model often used in the literature generally provides a good description of the observed horizontal profiles. Three out of the seven galaxies show signatures of extended vertical emission at 100 and 160 {\mu}m at the 5{\sigma} level, but in two of these it is probably due to deviations from an exactly edge-on orientation. Only for NGC 4013, a galaxy in which vertically extended dust has already been detected in optical images, we can detect vertically extended dust, and the derived scaleheight agrees with the value estimated through radiative transfer modelling. Our analysis hints at a correlation between the dust scaleheight and its degree of clumpiness, which we infer from the difference between the dust masses as calculated from modelling of optical data and from fitting the spectral energy distribution of Herschel datapoints.