Search

Low-metallicity dwarf galaxies often show no or little CO emission, despite the intense star formation observed in local samples. Both simulations and resolved observations indicate that molecular gas in low-metallicity galaxies may reside in small dense clumps, surrounded by a substantial amount of more diffuse gas, not traced by CO. Constraining the relative importance of CO-bright versus CO-dark H2 star-forming reservoirs is crucial to understand how star formation proceeds at low metallicity. We put to the test classically used single component radiative transfer models and compare their results to those obtained assuming an increasingly complex structure of the interstellar gas, mimicking an inhomogeneous distribution of clouds with various physical properties. We compute representative models of the interstellar medium as combinations of several gas components, each with a specific set of physical parameters. We introduce physically-motivated models assuming power-law distributions for the density, ionization parameter, and the depth of molecular clouds. We confirm the presence of a predominantly CO-dark molecular reservoir in low-metallicity galaxies. The predicted total H2 mass is best traced by [C II]158um and, to a lesser extent, by [CI] 609um, rather than by CO(1-0). We examine the CO-to-H2 conversion factor vs. metallicity relation and find that its dispersion increases significantly when different geometries of the gas are considered. We define a clumpiness parameter that anti-correlates with [CII]/CO and explains the dispersion of the CO-to-H2 conversion factor vs. metallicity relation. We find that low-metallicity galaxies with high clumpiness may have CO-to-H2 conversion factor as low as the Galactic value. We identify the clumpiness of molecular gas as a key parameter to understand variations of geometry-sensitive quantities, such as CO-to-H2 conversion factor., Comment: Submitted to A&A, comments welcome

Context. Molecular filaments and hubs have received special attention recently thanks to new studies showing their key role in star formation. While the (column) density and velocity structures of both filaments and hubs have been carefully studied, their magnetic field (B-field) properties have yet to be characterized. Consequently, the role of B-fields in the formation and evolution of hub-filament systems is not well constrained. Aims. We aim to understand the role of the B-field and its interplay with turbulence and gravity in the dynamical evolution of the NGC 6334 filament network that harbours cluster-forming hubs and high-mass star formation. Methods. We present new observations of the dust polarized emission at 850 μm toward the 2 pc × 10 pc map of NGC 6334 at a spatial resolution of 0.09 pc obtained with the James Clerk Maxwell Telescope (JCMT) as part of the B-field In STar-forming Region Observations (BISTRO) survey. We study the distribution and dispersion of the polarized intensity (PI), the polarization fraction (PF), and the plane-of-the-sky B-field angle (χB_POS) toward the whole region, along the 10 pc-long ridge and along the sub-filaments connected to the ridge and the hubs. We derived the power spectra of the intensity and χBPOS along the ridge crest and compared them with the results obtained from simulated filaments. Results. The observations span ~3 orders of magnitude in Stokes I and PI and ~2 orders of magnitude in PF (from ~0.2 to ~ 20%). A large scatter in PI and PF is observed for a given value of I. Our analyses show a complex B-field structure when observed over the whole region (~ 10 pc); however, at smaller scales (~1 pc), χBPOS varies coherently along the crests of the filament network. The observed power spectrum of χBPOS can be well represented with a power law function with a slope of − 1.33 ± 0.23, which is ~20% shallower than that of I. We find that this result is compatible with the properties of simulated filaments and may indicate the physical processes at play in the formation and evolution of star-forming filaments. Along the sub-filaments, χBPOS rotates frombeing mostly perpendicular or randomly oriented with respect to the crests to mostly parallel as the sub-filaments merge with the ridge and hubs. This variation of the B-field structure along the sub-filaments may be tracing local velocity flows of infalling matter in the ridge and hubs. Our analysis also suggests a variation in the energy balance along the crests of these sub-filaments, from magnetically critical or supercritical at their far ends to magnetically subcritical near the ridge and hubs. We also detect an increase in PF toward the high-column density (NH2 ≳ 1023 cm−2) star cluster-forming hubs. These latter large PF values may be explained by the increase in grain alignment efficiency due to stellar radiation from the newborn stars, combined with an ordered B-field structure. Conclusions. These observational results reveal for the first time the characteristics of the small-scale (down to ~ 0.1 pc) B-field structure of a 10 pc-long hub-filament system. Our analyses show variations in the polarization properties along the sub-filaments that may be tracing the evolution of their physical properties during their interaction with the ridge and hubs. We also detect an impact of feedback from young high-mass stars on the local B-field structure and the polarization properties, which could put constraints on possible models for dust grain alignment and provide important hints as to the interplay between the star formation activity and interstellar B-fields.

Aims. The purpose of this work is the characterization of the main scaling relations between all of the interstellar medium (ISM) components, namely dust, atomic, molecular, and total gas, and gas-phase metallicity, as well as other galaxy properties, such as stellar mass (Mstar) and galaxy morphology, for late-type galaxies in the Local Universe. Methods. This study was 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 spanning from T = 1-10, Mstar from 6 × 107 to 3 × 1011 M⊙ , star formation rate from 6 × 10-4 to 60 M⊙ yr-1, and oxygen abundance from 12 + log(O/H) = 8-9.5. Molecular and atomic gas data were collected from the literature and properly homogenized. All the masses involved in our analysis refer to the values within the optical disks of galaxies. The scaling relations involving the molecular gas are studied by assuming both a constant and a metallicity-dependent CO-to-H2 conversion factor (XCO). The analysis was performed by means of the survival analysis technique, in order to properly take into account the presence of both detection and nondetection in the data. Results. 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 with the molecular one. We characterize important mass ratios such as the gas fraction, the molecular-to-atomic gas mass ratio, the dust-to-total gas mass ratio (DGR), and the dust-to-stellar mass ratio, and study how they relate to each other, to galaxy morphology, and to gas-phase metallicity. 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. The DGR, the gas-phase metallicity, and the dust-to-stellar mass ratio are, for our galaxy sample, directly linked to galaxy morphology. The molecular-to-atomic gas mass ratio and the DGR show a positive correlation for low molecular gas fractions, but for galaxies rich in molecular gas this trend breaks down. To our knowledge, this trend has never been found before, and provides new constraints for theoretical models of galaxy evolution and a reference for high-redshift studies. We discuss several scenarios related to this finding. Conclusions. The DustPedia database of late-type galaxies is an extraordinary tool for the study of the ISM scaling relations, thanks to its homogeneous collection of data for the different ISM components. The database is made publicly available to the whole community. © ESO 2020.

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 ultraviolet and optical photons emitted by the youngest (100 Myr) stars. Consequently, thermal re-emission by ust at far-infrared 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 stellar populations might be more significant than previously thought. Advances in radiation transfer simulations finally allow us to actually quantify the heating mechanisms of difuse dust by the stellar radiation field. Aims. As one of the main goals in the DustPedia project, we have developed a framework to construct detailed 3D stellar and dust radiative transfer models for nearby galaxies. In this study, we analyse the contribution of the di erent stellar populations to the dust heating in four nearby face-on barred galaxies: NGC 1365, M 83, M 95, and M 100. We aim to quantify the fraction directly related to young stellar populations, both globally and on local scales, and to assess the influence of the bar on the heating fraction. Methods. From 2D images we derive the 3D distributions of stars and dust. To model the complex geometries, we used skirt, a state-of-the-art 3D Monte Carlo radiative transfer code designed to self-consistently simulate the absorption, scattering, and thermal re-emission by the dust for arbitrary 3D distributions. Results. We derive global attenuation laws for each galaxy and confirm that galaxies of high specific star-formation rate have shallower attenuation curves and weaker UV bumps. On average, 36.5% of the bolometric luminosity is absorbed by dust in our galaxy sample. We report a clear e ect of the bar structure on the radial profiles of the dust-heating fraction by the young stellar populations, and the dust temperature.We find that the young stellar populations 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 stellar populations and the specific star-formation rate. © 2020 ESO.

The star formation rate and the mass of interstellar medium (ISM) have a high predictive power for the future evolution of a galaxy. Nevertheless, deriving such properties is not straightforward. 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 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. We take special care in modelling the morphological features present for each source of radiation. In this particular study, we focus on NGC 1068, which in addition contains an active galactic nucleus (AGN). The effect of diffuse dust heating by an AGN (beyond the torus) has so far only been investigated for quasars. This additional dust heating source further contaminates the broadband fluxes that are used by classic galaxy modelling tools 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 submillimetre (submm). This is partly due to beam smearing effects, but also because the input dust distribution is not sufficiently peaked in the centre. We find that AGN contamination of the broadband fluxes has a strong dependency on wavelength. It peaks in the mid-infrared, drops in the far-infrared, and then 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 hundred parsecs, 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 could be stronger in galaxies with a more luminous AGN. This could significantly impact the derived star formation rates and ISM masses for such systems. © ESO 2020.

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. © ESO 2020.

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. © 2020 The Author(s)

Aims. Within the framework of the DustPedia project we investigate the properties of cosmic dust and its interaction with stellar radiation (originating from di erent stellar populations) for 814 galaxies in the nearby Universe, all observed by the Herschel Space Observatory. Methods. We take advantage of the widely used fitting code CIGALE, properly adapted to include the state-of-the-art dust model THEMIS. For comparison purposes, an estimation of the dust properties is provided by approximating the emission at far-infrared and sub-millimeter wavelengths with a modified blackbody. 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, the unattenuated and the absorbed by dust stellar light, for both the old (>200 Myr) and young (≤200 Myr) stellar populations. Results. 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 di erent 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 ecient 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, with ~60% of their luminosity contributing to that purpose. Finally, we find a strong correlation of the dust heating fraction by young stars with morphology and the specific star-formation rate. © ESO 2019.