S. Viaene

Dust spectral energy distributions of nearby galaxies: an insight from the Herschel Reference Survey

Although it accounts only for a small fraction of the baryonic mass, dust has a profound impact on the physical processes at play in galaxies. Thus, to understand the evolution of galaxies, it is essential not only to characterize dust properties per se, but also in relation to global galaxy properties. To do so, we derive the dust properties of galaxies in a volume limited, K-band selected sample, the Herschel Reference Survey (HRS). We gather infrared photometric data from 8 μm to 500 μm from Spitzer, WISE, IRAS, and Herschel for all of the HRS galaxies. Draine & Li (2007, ApJ, 663, 866) 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 parameters of the Draine & Li models and that a strong constraint on the 20−60 μm range is mandatory to estimate the relative contribution of the photo-dissociation regions to the infrared spectral energy distribution (SED). The SED models tend to systematically underestimate the observed 500 μm flux densities, especially for low-mass systems. We provide the output parameters for all of the galaxies, i.e., the minimum intensity of the interstellar radiation field, the fraction of polycyclic aromatic hydrocarbon (PAH), the relative contribution of PDR and evolved stellar population to the dust heating, the dust mass, and the infrared luminosity. For a subsample of gas-rich galaxies, we analyze the relations between these parameters and the main integrated properties of galaxies, such as stellar mass, star formation rate, infraredluminosity, metallicity, Hα and H-band surface brightness, and the far-ultraviolet 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 and, thus, low stellar masses. The intensity of the diffuse interstellar radiation field and the H-band and Hα 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 infrared 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 place our sample on the SFR − M∗ diagram. The templates are compared to the most popular infrared SED libraries, enlightening a large discrepancy between all of them in the 20−100 μm range.

The distribution of interstellar dust in CALIFA edge-on galaxies via oligochromatic radiative transfer fitting

We investigate the amount and spatial distribution of interstellar dust in edge-on spiral galaxies, using detailed radiative transfer modelling of a homogeneous sample of 12 galaxies selected from the Calar Alto Legacy Integral Field Area survey. Our automated fitting routine, FITSKIRT, was first validated against artificial data. This is done by simultaneously reproducing the Sloan Digital Sky Survey 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 Sersic 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 per cent 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 identification of dust heating mechanisms in nearby galaxies using Herschel 160/250 and 250/350μm surface brightness ratios

We examined variations in the 160/250 and 250/350 mu m surface brightness ratios within 24 nearby (\textless30 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 mu m ratios to H alpha and 24 mu m emission, which trace the light from star-forming regions, and 3.6 mu m 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 \textless= 160 mu m (and in rare cases potentially at \textless= 350 mu m) originates from dust heated by star-forming regions. In another subset, we found that the emission at \textgreater= 250 mu m (and sometimes at \textgreater= 160 mu m) 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 do 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 Herschel Exploitation of Local Galaxy Andromeda (HELGA). IV. Dust scaling relations at sub-kpc resolution

Context. Dust and stars play a complex game of interactions in the interstellar medium and around young stars. The imprints of these processes are visible in scaling relations between stellar characteristics, star formation parameters, and dust properties. Aims. In the present work, we aim to examine dust scaling relations on a sub-kpc resolution in the Andromeda galaxy (M 31). The goal is to investigate the properties of M 31 on both a global and local scale and compare them to other galaxies of the local universe. Methods. New Herschel observations are combined with available data from GALEX, SDSS, WISE, and Spitzer to construct a dataset covering UV to submm wavelengths. All images were brought to the beam size and pixel grid of the SPIRE 500 μm frame. This divides M 31 in 22 437 pixels of 36 arcseconds in size on the sky, corresponding to physical regions of 137 × 608 pc in the galaxy’s disk. A panchromatic spectral energy distribution was modelled for each pixel and maps of the physical quantities were constructed. Several scaling relations were investigated, focussing on the interactions of dust with starlight. Results. We find, on a sub-kpc scale, strong correlations between M_dust/M_⋆ and NUV-r, and between M_dust/M_⋆ and μ_⋆ (the stellar mass surface density). Striking similarities with corresponding relations based on integrated galaxies are found. We decompose M 31 in four macro-regions based on their far-infrared morphology; the bulge, inner disk, star forming ring, and the outer disk region. In the scaling relations, all regions closely follow the galaxy-scale average trends and behave like galaxies of different morphological types. 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. Within each macro-region, however, a great diversity in individual micro-regions is found, regardless of the properties of the macro-regions. Furthermore, we confirm that dust in the bulge of M 31 is heated only by the old stellar populations. Conclusions. In general, the local dust scaling relations indicate that the dust content in M 31 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.

High-resolution, 3D radiative transfer modeling - I. The grand-design spiral galaxy M 51

Context. Dust reprocesses about half of the stellar radiation in galaxies. The thermal re-emission by dust of absorbed energy is considered to be driven merely by young stars so is often applied to tracing the star formation rate in galaxies. Recent studies have argued that the old stellar population might be responsible for 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 M 51 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 to 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 to account for the absorption, scattering, and non-local thermal equilibrium (NLTE) emission of dust in M 51. The 3D distribution of stars is derived from the 2D morphology observed in the IRAC 3.6 μm, GALEX FUV, Hα, and MIPS 24 μ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. Results. The dust emission derived from RT calculations is consistent with far-infrared and submillimeter observations of M 51, implying that the absorbed stellar energy is balanced by the thermal re-emission of dust. The young stars provide 63% of the energy for heating the dust responsible for the total infrared emission (8−1000 μm), while 37% of the dust emission is governed through heating by the evolved stellar population. In individual wavebands, the contribution from young stars to the dust heating dominates at all infrared wavebands but gradually decreases towards longer infrared and submillimeter wavebands for which the old stellar population becomes a non-negligible source of heating. Upon extrapolation of the results for M 51, we present prescriptions for estimating the contribution of young stars to the global dust heating based on a tight correlation between the dust heating fraction and specific star formation rate.

A resolved analysis of cold dust and gas in the nearby edge-on spiral NGC 891

We investigate the connection between dust and gas in the nearby edge-on spiral galaxy NGC 891, a target of the Very Nearby Galaxies Survey. High resolution Herschel PACS and SPIRE 70, 100, 160, 250, 350, and 500 μm images are combined with JCMT SCUBA 850 μ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 ± 2.0) × 107M⊙ and an average temperature of 23 ± 2 K, consistent with results from previous far-infrared experiments. We also fit one-component modified blackbody models to 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 μm emission. The derived dust temperature, which ranges from approximately 17 to 24 K, is found to correlate with the 24 μm emission. Allowing the dust emissivity index to vary, we find an average value of β = 1.9 ± 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 disc. A comparison of the dust properties with the gaseous components of the ISM reveals strong spatial correlations between the surface mass densities of dust (Σdust) and the molecular hydrogen (ΣH2) and total gas surface densities (Σgas). These observations reveal the presence of regions of dense, cold dust that are coincident with peaks in the gas distribution and are associated with a molecular ring. Furthermore, the observed asymmetries in the dust temperature, the H2-to-dust ratio and the total gas-to-dust ratio 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 that there is also an enhancement in the star formation rate in the NE part of the disc.

HERschel Observations of Edge-on Spirals (HEROES) - I. Far-infrared morphology and dust mass determination

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 μm at the 5σ 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.

The Herschel-ATLAS Data Release 1 – II. Multi-wavelength counterparts to submillimetre sources

This paper is the second in a pair of papers 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 HerschelSpace Observatory. The H-ATLAS is a wide-area imaging survey carried out in five photometric bands at 100, 160, 250, 350 and 500 μm covering a total area of 600 deg2. In this paper, we describe the identification of optical counterparts to submillimetre sources in DR1, comprising an area of 161 deg2 over three equatorial fields of roughly 12 × 4.5 deg centred at 9h, 12h and 14 h.5, 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 Sloan Digital Sky Survey counterparts at r < 22.4 for 250-μm-selected sources detected at ≥4σ (≈28 mJy). We find ‘reliable’ counterparts (reliability R ≥ 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 active galactic nuclei in comparison with optical selection. The data described herein, together with all maps and catalogues described in the companion paper, are available from the H-ATLAS website at www.h-atlas.org.

The Herschel exploitation of local galaxy Andromeda (HELGA) V: Strengthening the case for substantial interstellar grain growth

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 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 ejecta and/or stellar atmospheres, and grain destruction in the interstellar medium 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.

The imprint of rapid star formation quenching on the spectral energy distributions of galaxies

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