High-resolution, 3D radiative transfer modelling

Abstract

Context. Interstellar dust absorbs stellar light very efficiently, thus shaping the energy output of galaxies. Studying the impact of different stellar populations on the dust heating continues to be a challenge because it requires decoupling the relative geometry of stars and dust and also involves complex processes such as scattering and non-local dust heating. Aims. We aim to constrain the relative distribution of dust and stellar populations in the spiral galaxy M 81 and create a realistic model of the radiation field that adequately describes the observations. By 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 these can be used for comparable modelling of other nearby galaxies. Methods. We present a semi-automated radiative transfer modelling pipeline that implements necessary steps such as the geometric model construction and the normalisation of the components through an optimisation routine. We used 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 were able to quantify to what extent different stellar age populations contribute to the heating of dust. Our method takes into account the effects of non-local heating. Results. We obtained a realistic 3D radiative transfer model of the face-on galaxy M 81. We find that only 50.2% of the dust heating can be attributed to young stellar populations (≲100 Myr). We confirm that there is a tight correlation between the specific star formation rate and the heating fraction by young stellar populations, both in sky projections and in 3D, which is also found for radiative transfer models of M 31 and M 51. Conclusions. We conclude that old stellar populations can be a major contributor to the heating of dust. In M 81, old stellar populations are the dominant heating agent in the central regions, contributing to half of the absorbed radiation. Regions of higher star formation do not correspond to the highest dust temperatures. On the contrary, it is the dominant bulge which is most efficient in heating the dust. The approach we present here can immediately be applied to other galaxies. It does contain a number of caveats, which we discuss in detail.