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.
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.
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.
Abstract
With its relatively low ionization potential, C
+
can be found throughout the interstellar medium (ISM) and provides one of the main cooling channels of the ISM via the C
ii
157
μ
m ...emission. While the strength of the C
ii
line correlates with the star formation rate, the contributions of the various gas phases to the C
ii
emission on galactic scales are not well established. In this study we establish an empirical multi-component model of the ISM, including dense H
ii
regions, dense photon dissociation regions (PDRs), the warm ionized medium (WIM), low density and
surfaces of molecular clouds (SfMCs), and the cold neutral medium (CNM). We test our model on ten luminous regions within the two nearby galaxies NGC 3184 and NGC 628 on angular scales of 500–600 pc. Both galaxies are part of the
Herschel
key program KINGFISH, and are complemented by a large set of ancillary ground- and space-based data. The five modeled phases together reproduce the observed C
ii
emission quite well, overpredicting the total flux slightly (about 45%) averaged over all regions. We find that dense PDRs are the dominating component, contributing 68% of the C
ii
flux on average, followed by the WIM and the SfMCs, with mean contributions of about half of the contribution from dense PDRs, each. CNM and dense H
ii
regions are only minor contributors with less than 5% each. These estimates are averaged over the selected regions, but the relative contributions of the various phases to the C
ii
flux vary significantly between these regions.
ABSTRACT
In order to understand the contribution of core-collapse supernovae to the dust budget of the early Universe, it is important to understand not only the mass of dust that can form in ...core-collapse supernovae but also the location and rate of dust formation. SN 2005ip is of particular interest since dust has been inferred to have formed in both the ejecta and the post-shock region behind the radiative reverse shock. We have collated eight optical archival spectra that span the lifetime of SN 2005ip and we additionally present a new X-shooter optical–near-IR spectrum of SN 2005ip at 4075 d post-discovery. Using the Monte Carlo line transfer code damocles, we have modelled the blueshifted broad and intermediate-width H α, H β, and He i lines from 48 to 4075 d post-discovery using an ejecta dust model. We find that dust in the ejecta can account for the asymmetries observed in the broad and intermediate-width H α, H β, and He i line profiles at all epochs and that it is not necessary to invoke post-shock dust formation to explain the blueshifting observed in the intermediate-width post-shock lines. Using a Bayesian approach, we have determined the evolution of the ejecta dust mass in SN 2005ip over 10 yr presuming an ejecta dust model, with an increasing dust mass from ∼10−8 M⊙ at 48 d to a current dust mass of ∼0.1 M⊙.