Dust Energy Balance Research Articles

High-resolution radiative transfer modelling of M33

ABSTRACT In this work, we characterize 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 characterized 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 (∼80 per cent at 10 $\mu$m to ∼50 per cent at 1 mm). We find that the local dust-energy balance is restored at spatial scales greater than around 1.5 kpc.

HERschelObservations of Edge-on Spirals (HEROES)

We present results of the detailed dust energy balance study for the seven large edge-on galaxies in the HEROES sample using three-dimensional (3D) radiative transfer (RT) modelling. Based on available optical and near-infrared (NIR) 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érsic 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 RT models to the optical and NIR images with the fitting algorithmFITSKIRTand run panchromatic simulations with theSKIRTcode 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–submillimetre 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 results from the literature, the large- and small-scale structure is the most probable explanation for the dust energy balance problem.

HERschelObservations of Edge-on Spirals (HEROES)

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.

Dust energy balance study of two edge-on spiral galaxies in the Herschel-ATLAS survey

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 3. We investigate the dust energy balance for IC 4225 and NGC 5166, 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 IC 4225 underestimates the observed fluxes by a factor of about 3. For NGC 5166 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.

Large and small-scale structures and the dust energy balance problem in spiral galaxies

The interstellar dust content in galaxies can be traced in extinction at optical wavelengths, or in emission in the far-infrared. Several studies have found that radiative transfer models that successfully explain the optical extinction in edge-on spiral galaxies generally underestimate the observed FIR/submm fluxes by a factor of about three. In order to investigate this so-called dust energy balance problem, we use two Milky Way-like galaxies produced by high-resolution hydrodynamical simulations. We create mock optical edge-on views of these simulated galaxies (using the radiative transfer code SKIRT), and we then fit the parameters of a basic spiral galaxy model to these images (using the fitting code FitSKIRT). The basic model includes smooth axisymmetric distributions along a S\'ersic bulge and exponential disc for the stars, and a second exponential disc for the dust. We find that the dust mass recovered by the fitted models is about three times smaller than the known dust mass of the hydrodynamical input models. This factor is in agreement with previous energy balance studies of real edge-on spiral galaxies. On the other hand, fitting the same basic model to less complex input models (e.g. a smooth exponential disc with a spiral perturbation or with random clumps), does recover the dust mass of the input model almost perfectly. Thus it seems that the complex asymmetries and the inhomogeneous structure of real and hydrodynamically simulated galaxies are a lot more efficient at hiding dust than the rather contrived geometries in typical quasi-analytical models. This effect may help explain the discrepancy between the dust emission predicted by radiative transfer models and the observed emission in energy balance studies for edge-on spiral galaxies.

The dust energy balance in the edge-on spiral galaxy NGC 4565

We combine new dust continuum observations of the edge-on spiral galaxy NGC 4565 in all Herschel/SPIRE (250, 350, 500 micron) wavebands, obtained as part of the Herschel Reference Survey, and a large set of ancillary data (Spitzer, SDSS, GALEX) to analyze its dust energy balance. We fit a radiative transfer model for the stars and dust to the optical maps with the fitting algorithm FitSKIRT. To account for the observed UV and mid-infrared emission, this initial model was supplemented with both obscured and unobscured star-forming regions. Even though these star-forming complexes provide an additional heating source for the dust, the far-infrared/submillimeter emission long wards of 100 micron is underestimated by a factor of 3-4. This inconsistency in the dust energy budget of NGC 4565 suggests that a sizable fraction (two-thirds) of the total dust reservoir (Mdust ~ 2.9e+8 Msun) consists of a clumpy distribution with no associated young stellar sources. The distribution of those dense dust clouds would be in such a way that they remain unresolved in current far-infrared/submillimeter observations and hardly comtribute to the attenuation at optical wavelengths. More than two-thirds of the dust heating in NGC 4565 is powered by the old stellar population, with localized embedded sources supplying the remaining dust heating in NGC 4565. The results from this detailed dust energy balance study in NGC 4565 is consistent with that of similar analyses of other edge-on spirals.

Panchromatic radiative transfer modelling of stars and dust in the Sombrero galaxy

We present a detailed study of the dust energy balance in the Sombrero galaxy M104. From a full radiative transfer analysis, including scattering, absorption and thermal re-emission, we construct models that can reproduce images at optical/near-infrared wavelengths, the observed stellar SED and the minor axis extinction profiles in the V and R_C band. A standard model, that contains only an old stellar population to heat the dust, underestimates the observations of dust emission at infrared wavelengths by a factor of ~ 3. Supplementing this basic model with a young stellar component of low star formation activity in both the inner disk (SFR ~ 0.21 Msun/yr) and dust ring (SFR ~ 0.05 Msun/yr), we are capable of solving the discrepancy in the dust energy budget of the Sombrero galaxy at wavelengths shortwards of 100 \mum. To account for the increased FIR/submm emission beyond 100 \mum, we propose a additional dust component distributed in quiescent clumps. This model with a clumpy dust structure predicts three-quarters of the total dust content (~ 2.8 x 10^7 Msun) to reside in compact dust clouds with no associated embedded sources. Although the assumption of a clumpy dust structure in the Sombrero galaxy is supported by high-resolution optical data, we cannot rule out the possibility that dust grains with a higher dust emissivity account for part of the discrepancy in the energy budget at submm wavelengths.

A detailed dust energy balance study of the Sombrero galaxy

AbstractThe Sombrero galaxy (M104) is an interesting object for a dust energy balance study due to its very symmetric dust lane, its proximity and its (nearly edge-on) inclination of 84°. From a panchromatic radiative transfer analysis, including scattering, absorption and thermal dust re-emission, we construct a standard model for M104 accounting for observations in the optical wave bands (stellar SED, images and extinction profiles in the V and RC band). This standard model underestimates the observed dust emission at infrared wavelengths by a factor of ~ 3, similar to the discrepancy found in other energy balance studies of edge-on spirals. Supplementing this standard model with a young stellar component of low star formation activity in both the inner disk (SFR ~ 0.21 M⊙ yr−1) and dust ring (SFR ~ 0.05 M⊙ yr−1), we are capable of solving the discrepancy in the dust energy budget of the Sombrero galaxy at wavelengths shortwards of 100 μm. To account for the remaining discrepancy at longer wavelengths, we propose a secondary dust component distributed in quiescent clumps. This model with a clumpy dust structure predicts three-quarters of the total dust content to reside in compact dust clouds with no associated embedded sources.

Herschel-ATLAS: The dust energy balance in the edge-on spiral galaxy UGC 4754

We use Herschel PACS and SPIRE observations of the edge-on spiral galaxy UGC 4754, taken as part of the H-ATLAS SDP observations, to investigate the dust energy balance in this galaxy. We build detailed SKIRT radiative models based on SDSS and UKIDSS maps and use these models to predict the far-infrared emission. We find that our radiative transfer model underestimates the observed FIR emission by a factor of two to three. Similar discrepancies have been found for other edge-on spiral galaxies based on IRAS, ISO, and SCUBA data. Thanks to the good sampling of the SED at FIR wavelengths, we can rule out an underestimation of the FIR emissivity as the cause for this discrepancy. Instead we support highly obscured star formation that contributes little to the optical extinction as a more probable explanation.

IRAS observations of dust heating and energy balance in the Rho Ophiuchi dark cloud

The equilibrium process dust emission in the Rho Ophiuchi dark cloud is studied. The luminosity of the cloud is found to closely match the luminosity of the clouds's known embedded and external radiation sources. There is no evidence for a large population of undetected low-luminosity sources within the cloud and unknown external heating is also only a minor source of energy. Most of the cloud's luminosity is emitted in the mid-to-far-IR. Dust temperature maps indicate that the dust is not hot enough to heat the gas to observed temperatures. A simple cloud model with a radiation field composed of flux HD 147889, S1, and Sco OB2 associations predicts the observed IRAS 60 to 100 micron in-band flux ratios for a mean cloud density n(H2) = 1400. Flattened 12 and 25 micron observations show much extended emission in these bands, suggesting stochastic heating of very small grains or large molecules.