No Arabic abstract
We present the first radiative transfer (RT) model of a non-edge-on disk galaxy in which the large-scale geometry of stars and dust is self-consistently derived through fitting of multiwavelength imaging observations from the UV to the submm. To this end we used the axi-symmetric RT model of Popescu et al. and a new methodology for deriving geometrical parameters, and applied this to decode the{spectral energy distribution (SED) of M33. We successfully account for both the spatial and spectral energy distribution, with residuals typically within $7%$ in the profiles of surface brightness and within $8%$ in the spatially-integrated SED. We predict well the energy balance between absorption and re-emission by dust, with no need to invoke modified grain properties, and we find no submm emission that is in excess of our model predictions. We calculate that $80pm8%$ of the dust heating is powered by the young stellar populations. We identify several morphological components in M33, a nuclear, an inner, a main and an outer disc, showing a monotonic trend in decreasing star-formation surface-density ($Sigma_{rm SFR}$) from the nuclear to the outer disc. In relation to surface density of stellar mass, the $Sigma_{rm SFR}$ of these components define a steeper relation than the main sequence of star-forming galaxies, which we call a structurally resolved main sequence. Either environmental or stellar feedback mechanisms could explain the slope of the newly defined sequence. We find the star-formation rate to be ${rm SFR}=0.28^{+0.02}_{-0.01}{rm M}_{odot}{rm yr}^{-1}$.
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 present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time flexibility in the assumed stellar physics. We include a model for nebular line emission that can employ either precomputed Cloudy lookup tables (for efficiency), or direct photoionization calculations for all young stars (for flexibility). The dust content follows either observationally-motivated prescriptions, direct modeling from galaxy formation simulations, or a novel approach that includes the dust content via learning-based algorithms from the SIMBA cosmological galaxy formation simulation. AGN can additionally be included via a range of prescriptions. The output of these models are broadband SEDs, as well as filter-convolved images. Powderday is designed to eliminate last-mile efforts by researchers that employ different hydrodynamic galaxy formation models, and seamlessly interfaces with GIZMO, AREPO, GASOLINE, CHANGA, and ENZO. We demonstrate the capabilities of the code via three applications: a model for the star formation rate (SFR) - infrared luminosity relation in galaxies (including the impact of AGN); the impact of circumstellar dust around AGB stars on the mid-infrared emission from galaxy SEDs; and the impact of galaxy inclination angle on dust attenuation laws.
We present a Giant Molecular Cloud (GMC) catalog toward M33, containing 71 GMCs in total, based on wide field and high sensitivity CO(J=3-2) observations with a spatial resolution of 100 pc using the ASTE 10 m telescope. Employing archival optical data, we identify 75 young stellar groups (YSGs) from the excess of the surface stellar density, and estimate their ages by comparing with stellar evolution models. A spatial comparison among the GMCs, YSGs, and HII regions enable us to classify GMCs into four categories: Type A showing no sign of massive star formation (SF), Type B being associated only with HII regions, Type C with both HII regions and <10 Myr-old YSGs and Type-D with both HII regions and 10--30 Myr YSGs. Out of 65 GMCs (discarding those at the edges of the observed fields), 1 (1%), 13 (20%), 29 (45%), and 22 (34%) are Types A, B, C, and D, respectively. We interpret these categories as stages in a GMC evolutionary sequence. Assuming that the timescale for each evolutionary stage is proportional to the number of GMCs, the lifetime of a GMC with a mass >10^5 Mo is estimated to be 20--40 Myr. In addition, we find that the dense gas fraction as traced by the CO(J=3-2)/CO(J=1-0) ratio is enhanced around SF regions. This confirms a scenario where dense gas is preferentially formed around previously generated stars, and will be the fuel for the next stellar generation. In this way, massive SF gradually propagates in a GMC until gas is exhausted.
We perform simulations of isolated galaxies in order to investigate the likely origin of the spiral structure in M33. In our models, we find that gravitational instabilities in the stars and gas are able to reproduce the observed spiral pattern and velocity field of M33, as seen in HI, and no interaction is required. We also find that the optimum models have high levels of stellar feedback which create large holes similar to those observed in M33, whilst lower levels of feedback tend to produce a large amount of small scale structure, and undisturbed long filaments of high surface density gas, hardly detected in the M33 disc. The gas component appears to have a significant role in producing the structure, so if there is little feedback, both the gas and stars organise into clear spiral arms, likely due to a lower combined $Q$ (using gas and stars), and the ready ability of cold gas to undergo spiral shocks. By contrast models with higher feedback have weaker spiral structure, especially in the stellar component, compared to grand design galaxies. We did not see a large difference in the behaviour of $Q_{stars}$ with most of these models, however, because $Q_{stars}$ stayed relatively constant unless the disc was more strongly unstable. Our models suggest that although the stars produce some underlying spiral structure, this is relatively weak, and the gas physics has a considerable role in producing the large scale structure of the ISM in flocculent spirals.
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.