ABSTRACT
Shockwaves driven by supernovae both destroy dust and reprocess the surviving grains, greatly affecting the resulting dust properties of the interstellar medium (ISM). While these processes ...have been extensively studied theoretically, observational constraints are limited. We use physically motivated models of dust emission to fit the infrared (IR) spectral energy distributions of seven Galactic supernova remnants, allowing us to determine the distribution of dust mass between diffuse and dense gas phases, and between large and small grain sizes. We find that the dense ($\sim \! 10^3\ {\rm cm}^{-3}$), relatively cool ($\sim \! 10^3\ {\rm K}$) gas phase contains $\gt 90{{\ \rm per\ cent}}$ of the dust mass, making the warm dust located in the X-ray emitting plasma ($\sim \! 1\ {\rm cm}^{-3}$/$10^6\ {\rm K}$) a negligible fraction of the total, despite dominating the mid-IR emission. The ratio of small ($\lesssim\!{10}\ {\rm nm}$) to large ($\gtrsim \! 0.1\ {\rm \mu m}$) grains in the cold component is consistent with that in the ISM, and possibly even higher, whereas the hot phase is almost entirely devoid of small grains. This suggests that grain shattering, which processes large grains into smaller ones, is ineffective in the low-density gas, contrary to model predictions. Single-phase models of dust destruction in the ISM, which do not account for the existence of the cold swept-up material containing most of the dust mass, are likely to greatly overestimate the rate of dust destruction by supernovae.
Dynamic and thermal processes regulate the structure of the multiphase interstellar medium (ISM), and ultimately establish how galaxies evolve through star formation. Thus, to constrain ISM models ...and better understand the interplay of these processes, it is of great interest to measure the thermal pressure ( ) of the diffuse, neutral gas. By combining C ii 158 m, H i, and CO data from 31 galaxies selected from the Herschel KINGFISH sample, we have measured thermal pressures in 534 predominantly atomic regions with typical sizes of ∼1 kiloparsec. We find a distribution of thermal pressures in the K cm−3 range. For a sub-sample of regions with conditions similar to those of the diffuse, neutral gas in the Galactic plane, we find thermal pressures that follow a log-normal distribution with a median value of Pth/k 3600 K cm−3. These results are consistent with thermal pressure measurements using other observational methods. We find that increases with radiation field strength and star formation activity, as expected from the close link between the heating of the gas and the star formation rate. Our thermal pressure measurements fall in the regime where a two-phase ISM with cold and warm neutral media could exist in pressure equilibrium. Finally, we find that the midplane thermal pressure of the diffuse gas is about ∼30% of the vertical weight of the overlying ISM, consistent with results from hydrodynamical simulations of self-regulated star formation in galactic disks.
ABSTRACT
The dust mass absorption coefficient, κd is the conversion function used to infer physical dust masses from observations of dust emission. However, it is notoriously poorly constrained, and ...it is highly uncertain how it varies, either between or within galaxies. Here we present the results of a proof-of-concept study, using the DustPedia data for two nearby face-on spiral galaxies M 74 (NGC 628) and M 83 (NGC 5236), to create the first ever maps of κd in galaxies. We determine κd using an empirical method that exploits the fact that the dust-to-metals ratio of the interstellar medium is constrained by direct measurements of the depletion of gas-phase metals. We apply this method pixel-by-pixel within M 74 and M 83, to create maps of κd. We also demonstrate a novel method of producing metallicity maps for galaxies with irregularly sampled measurements, using the machine learning technique of Gaussian process regression. We find strong evidence for significant variation in κd. We find values of κd at 500 $\mu$m spanning the range 0.11–0.25 ${\rm m^{2}\, kg^{-1}}$ in M 74, and 0.15–0.80 ${\rm m^{2}\, kg^{-1}}$ in M 83. Surprisingly, we find that κd shows a distinct inverse correlation with the local density of the interstellar medium. This inverse correlation is the opposite of what is predicted by standard dust models. However, we find this relationship to be robust against a large range of changes to our method – only the adoption of unphysical or highly unusual assumptions would be able to suppress it.
ABSTRACT
The N
ii
122 and 205
μ
m transitions are powerful tracers of the ionized gas in the interstellar medium. By combining data from 21 galaxies selected from the
Herschel
KINGFISH and Beyond ...the Peak surveys, we have compiled 141 spatially resolved regions with a typical size of ∼1 kpc, with observations of both N
ii
far-infrared lines. We measure N
ii
122/205 line ratios in the ∼0.6–6 range, which corresponds to electron gas densities of
n
e
∼ 1–300 cm
−3
, with a median value of
n
e
= 30 cm
−3
. Variations in the electron density within individual galaxies can be as high as a factor of ∼50, frequently with strong radial gradients. We find that
n
e
increases as a function of infrared color, dust-weighted mean starlight intensity, and star-formation rate (SFR) surface density (Σ
SFR
). As the intensity of the N
ii
transitions is related to the ionizing photon flux, we investigate their reliability as tracers of the SFR. We derive relations between the N
ii
emission and SFR in the low-density limit and in the case of a log-normal distribution of densities. The scatter in the correlation between N
ii
surface brightness and Σ
SFR
can be understood as a property of the
n
e
distribution. For regions with
n
e
close to or higher than the N
ii
line critical densities, the low-density limit N
ii
-based SFR calibration systematically underestimates the SFR because the N
ii
emission is collisionally quenched. Finally, we investigate the relation between N
ii
emission, SFR, and
n
e
by comparing our observations to predictions from the MAPPINGS-III code.
Context. To compute the star formation rate (SFR) of galaxies from the rest-frame ultraviolet (UV), it is essential to take the obscuration by dust into account. To do so, one of the most popular ...methods consists in combining the UV with the emission from the dust itself in the infrared (IR). Yet, different studies have derived different estimators, showing that no such hybrid estimator is truly universal. Aims. In this paper we aim at understanding and quantifying what physical processes fundamentally drive the variations between different hybrid estimators. In so doing, we aim at deriving new universal UV+IR hybrid estimators to correct the UV for dust attenuation at local and global scales, taking the intrinsic physical properties of galaxies into account. Methods. We use the CIGALE code to model the spatially resolved far-UV to far-IR spectral energy distributions of eight nearby star-forming galaxies drawn from the KINGFISH sample. This allows us to determine their local physical properties, and in particular their UV attenuation, average SFR, average specific SFR (sSFR), and their stellar mass. We then examine how hybrid estimators depend on said properties. Results. We find that hybrid UV+IR estimators strongly depend on the stellar mass surface density (in particular at 70 μm and 100 μm) and on the sSFR (in particular at 24 μm and the total infrared). Consequently, the IR scaling coefficients for UV obscuration can vary by almost an order of magnitude: from 1.55 to 13.45 at 24 μm for instance. This result contrasts with other groups who found relatively constant coefficients with small deviations. We exploit these variations to construct a new class of adaptative hybrid estimators based on observed UV to near-IR colours and near-IR luminosity densities per unit area. We find that they can reliably be extended to entire galaxies. Conclusions. The new estimators provide better estimates of attenuation-corrected UV emission than classical hybrid estimators published in the literature. Taking naturally variable impact of dust heated by old stellar populations into account, they constitute an important step towards universal estimators.
Context. The relative abundance of the dust grain types in the interstellar medium is directly linked to physical quantities that trace the evolution of galaxies. Because of the poor spatial ...resolution of the infrared and submillimetre data, we are able to study the dependence of the resolved infrared spectral energy distribution (SED) across regions of the interstellar medium (ISM) with different physical properties in just a few objects. Aims. We aim to study the dust properties of the whole disc of M 33 at spatial scales of ~170 pc. This analysis allows us to infer how the relative dust grain abundance changes with the conditions of the ISM, study the existence of a submillimetre excess and look for trends of the gas-to-dust mass ratio (GDR) with other physical properties of the galaxy. Methods. For each pixel in the disc of M 33 we have fitted the infrared SED using a physically motivated dust model that assumes an emissivity index β close to two. We applied a Bayesian statistical method to fit the individual SEDs and derived the best output values from the study of the probability density function of each parameter. We derived the relative amount of the different dust grains in the model, the total dust mass, and the strength of the interstellar radiation field (ISRF) heating the dust at each spatial location. Results. The relative abundance of very small grains tends to increase, and for big grains to decrease, at high values of Hα luminosity. This shows that the dust grains are modified inside the star-forming regions, in agreement with a theoretical framework of dust evolution under different physical conditions. The radial dependence of the GDR is consistent with the shallow metallicity gradient observed in this galaxy. The strength of the ISRF derived in our model correlates with the star formation rate in the galaxy in a pixel by pixel basis. Although this is expected, it is the first time that a correlation between the two quantities has been reported. We have produced a map of submillimetre excess in the 500 μm SPIRE band for the disc of M 33. The excess can be as high as 50% and increases at large galactocentric distances. We further studied the relation of the excess with other physical properties of the galaxy and find that the excess is prominent in zones of diffuse ISM outside the main star-forming regions, where the molecular gas and dust surface density are low.
Aims. We aim to study the fraction of stellar radiation absorbed by dust, fabs, in 814 galaxies of different morphological types. The targets constitute the vast majority (93%) of the DustPedia ...sample, including almost all large (optical diameter larger than 1′), nearby (v ≤ 3000 km s−1) galaxies observed with the Herschel Space Observatory. Methods. For each object, we modelled the spectral energy distribution from the ultraviolet to the sub-millimetre using the dedicated, aperture-matched DustPedia photometry and the Code Investigating GALaxy Evolution (CIGALE). The value of fabs was obtained from the total luminosity emitted by dust and from the bolometric luminosity, which are estimated by the fit. Results. On average, 19% of the stellar radiation is absorbed by dust in DustPedia galaxies. The fraction rises to 25% if only late-type galaxies are considered. The dependence of fabs on morphology, showing a peak for Sb-Sc galaxies, is weak; it reflects a stronger, yet broad, positive correlation with the bolometric luminosity, which is identified for late-type, disk-dominated, high-specific-star-formation rate, gas-rich objects. We find no variation of fabs with inclination, at odds with radiative transfer models of edge-on galaxies. These results call for a self-consistent modelling of the evolution of the dust mass and geometry along the build-up of the stellar content. We also provide template spectral energy distributions in bins of morphology and luminosity and study the variation of fabs with stellar mass and specific star-formation rate. We confirm that the local Universe is missing the high fabs, luminous and actively star-forming objects necessary to explain the energy budget in observations of the extragalactic background light.
ABSTRACT
We present the first observational infrared luminosity function (IRLF) measurement in the Epoch of Reionization (EoR) based on a ultraviolet (UV)-selected galaxy sample with the Atacama ...Large Millimeter Array (ALMA) spectroscopic observations. Our analysis is based on the ALMA large program Reionization Era Bright Emission Line Survey (REBELS), which targets 42 galaxies at z = 6.4–7.7 with C ii 158 $\rm{\mu m}$ line scans. 16 sources exhibit dust detection, 15 of which are also spectroscopically confirmed through the C ii line. The infrared (IR) luminosities of the sample range from log LIR/L⊙ = 11.4 to 12.2. Using the UV luminosity function as a proxy to derive the effective volume for each of our target sources, we derive IRLF estimates, both for detections and for the full sample including IR luminosity upper limits. The resulting IRLFs are well reproduced by a Schechter function with the characteristic luminosity of $\log L_{*}/\mathrm{ L}_\odot =11.6^{+0.2}_{-0.1}$ . Our observational results are in broad agreement with the average of predicted IRLFs from simulations at z ∼ 7. Conversely, our IRLFs lie significantly below lower redshift estimates, suggesting a rapid evolution from z ∼ 4 to z ∼ 7, into the reionization epoch. The IR obscured contribution to the cosmic star formation rate density at z ∼ 7 amounts to $\mathrm{log(SFRD/{\rm M}_{\odot }\,yr^{-1}\,Mpc^{-3}) = -2.66^{+0.17}_{-0.14} }$ that is at least ∼10 per cent of UV-based estimates. We conclude that the presence of dust is already abundant in the EoR and discuss the possibility of unveiling larger samples of dusty galaxies with future ALMA and JWST observations.
We have fit the far-ultraviolet (FUV) to sub-millimeter (850 μm) spectral energy distributions (SEDs) of the 61 galaxies from the Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel ...(KINGFISH). The fitting has been performed using three models: the Code for Investigating GALaxy Evolution (CIGALE), the GRAphite-SILicate approach (GRASIL), and the Multiwavelength Analysis of Galaxy PHYSical properties (MAGPHYS). We have analyzed the results of the three codes in terms of the SED shapes, and by comparing the derived quantities with simple “recipes” for stellar mass (Mstar), star-formation rate (SFR), dust mass (Mdust), and monochromatic luminosities. Although the algorithms rely on different assumptions for star-formation history, dust attenuation and dust reprocessing, they all well approximate the observed SEDs and are in generally good agreement for the associated quantities. However, the three codes show very different behavior in the mid-infrared regime: in the 5–10 μm region dominated by PAH emission, and also between 25 and 70 μm where there are no observational constraints for the KINGFISH sample. We find that different algorithms give discordant SFR estimates for galaxies with low specific SFR, and that the standard recipes for calculating FUV absorption overestimate the extinction compared to the SED-fitting results. Results also suggest that assuming a “standard” constant stellar mass-to-light ratio overestimates Mstar relative to the SED fitting, and we provide new SED-based formulations for estimating Mstar from WISE W1 (3.4 μm) luminosities and colors. From a principal component analysis of Mstar, SFR, Mdust, and O/H, we reproduce previous scaling relations among Mstar, SFR, and O/H, and find that Mdust can be predicted to within ∼0.3 dex using only Mstar and SFR.