ABSTRACT
ALMA observations have revealed the presence of dust in the first generations of galaxies in the Universe. However, the dust temperature Td remains mostly unconstrained due to the few ...available FIR continuum data at redshift $z$ > 5. This introduces large uncertainties in several properties of high-$z$ galaxies, namely their dust masses, infrared luminosities, and obscured fraction of star formation. Using a new method based on simultaneous C $\scriptstyle \rm II$ 158-μm line and underlying dust continuum measurements, we derive Td in the continuum and C $\scriptstyle \rm II$ detected $z$ ≈ 7 galaxies in the ALMA Large Project REBELS sample. We find 39 < Td < 58 K, and dust masses in the narrow range Md = (0.9−3.6) × 107 M⊙. These results allow us to extend for the first time the reported Td($z$) relation into the Epoch of Reionization. We produce a new physical model that explains the increasing Td($z$) trend with the decrease of gas depletion time, tdep = Mg/SFR, induced by the higher cosmological accretion rate at early times; this hypothesis yields Td ∝ (1 + $z$)0.4. The model also explains the observed Td scatter at a fixed redshift. We find that dust is warmer in obscured sources, as a larger obscuration results in more efficient dust heating. For UV-transparent (obscured) galaxies, Td only depends on the gas column density (metallicity), $T_{\rm d} \propto N_{\rm H}^{1/6}$ (Td ∝ Z−1/6). REBELS galaxies are on average relatively transparent, with effective gas column densities around NH ≃ (0.03−1) × 1021 cm−2. We predict that other high-$z$ galaxies (e.g. MACS0416-Y1, A2744-YD4), with estimated Td ≫ 60 K, are significantly obscured, low-metallicity systems. In fact, Td is higher in metal-poor systems due to their smaller dust content, which for fixed LIR results in warmer temperatures.
Interstellar dust and starlight are modeled for the galaxies of the project "Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel." The galaxies were observed by the Infrared Array ...Camera and the Multiband Imaging Photometer for Spitzer on Spitzer Space Telescope, and the Photodetector Array Camera and Spectrometer and the Spectral and Photometric Imaging Receiver on Herschel Space Observatory. With data from 3.6 to 500 m, dust models are strongly constrained. Using a physical dust model, for each pixel in each galaxy we estimate (1) dust surface density, (2) dust mass fraction in polycyclic aromatic hydrocarbons (PAHs), (3) distribution of starlight intensities heating the dust, (4) total infrared (IR) luminosity emitted by the dust, and (5) IR luminosity originating in subregions with high starlight intensity. The dust models successfully reproduce the observed global and resolved spectral energy distributions. With the angular resolution of Herschel, we obtain well-resolved maps (available online) for the dust properties. As in previous studies, we find the PAH fraction to be an increasing function of metallicity, with a threshold oxygen abundance Z/Z 0.1, but we find the data to be fitted best with increasing linearly with above a threshold value of 0.15(O/H) . We obtain total dust masses for each galaxy by summing the dust mass over the individual map pixels; these "resolved" dust masses are consistent with the masses inferred from a model fit to the global photometry. The global dust-to-gas ratios obtained from this study are found to correlate with galaxy metallicities. Systems with Z/Z 0.5 have most of their refractory elements locked up in dust, whereas in systems with Z/Z 0.3 most of these elements tend to remain in the gas phase. Within galaxies, we find that is suppressed in regions with unusually warm dust with . With knowledge of one long-wavelength flux density ratio (e.g., f160/f500), the minimum starlight intensity heating the dust ( ) can be estimated to within ∼50%, despite a variation in of more than two orders of magnitude. For the adopted dust model, dust masses can be estimated to within ∼0.2 dex accuracy using the f160/f500 flux ratio and the integrated dust luminosity, and to ∼0.07 dex accuracy using the 500 m luminosity alone. There are additional systematic errors arising from the choice of dust model, but these are hard to estimate. These calibrated prescriptions for estimating starlight heating intensity and dust mass may be useful for studies of high-redshift galaxies.
ABSTRACT
We analyse FIR dust continuum measurements for 14 galaxies (redshift z ≈ 7) in the ALMA Reionization Era Bright Emission Line Survey (REBELS) Large Program to derive their physical ...properties. Our model uses three input data, i.e. (a) the UV spectral slope, β, (b) the observed UV continuum flux at 1500 Å, F1500, (c) the observed continuum flux at $\approx 158\, \mu$m, F158, and considers Milky Way (MW) and SMC extinction curves, along with different dust geometries. We find that REBELS galaxies have 28−90.5 per cent of their star formation obscured; the total (UV+IR) star formation rates are in the range $31.5 \lt {\rm SFR}/({\rm M}_\odot \, {\rm yr}^{-1}) \lt 129.5$. The sample-averaged dust mass and temperature are $(1.3\pm 1.1)\times 10^7 \, \mathrm{M}_\odot$ and 52 ± 11 K, respectively. However, in some galaxies dust is particularly abundant (REBELS-14, $M^{\prime }_{\rm d} \approx 3.4 \times 10^7 \, \mathrm{M}_\odot$), or hot (REBELS-18, $T^{\prime }_{\rm d} \approx 67$ K). The dust distribution is compact (<0.3 kpc for 70 per cent of the galaxies). The inferred dust yield per supernova is $0.1 \le y_{\rm d}/\, \mathrm{M}_\odot \le 3.3$, with 70 per cent of the galaxies requiring $y_{\rm d} \lt 0.25 \, \mathrm{M}_\odot$. Three galaxies (REBELS-12, 14, 39) require $y_{\rm d} \gt 1 \, \mathrm{M}_\odot$, which is likely inconsistent with pure SN production, and might require dust growth via accretion of heavy elements from the interstellar medium. With the SFR predicted by the model and a MW extinction curve, REBELS galaxies detected in C ii nicely follow the local LCII−SFR relation, and are approximately located on the Kennicutt–Schmidt relation. The sample-averaged gas depletion time is $0.11\, y_{\rm P}^{-2}$ Gyr, where yP is the ratio of the gas-to-stellar distribution radius. For some systems, a solution simultaneously matching the observed (β, F1500, F158) values cannot be found. This occurs when the index Im = (F158/F1500)/(β − βint), where βint is the intrinsic UV slope, exceeds $I_m^{*}\approx 1120$ for an MW curve. For these objects, we argue that the FIR and UV emitting regions are not co-spatial, questioning the use of the IRX–β relation.
Aims. The purpose of this work is the characterization of the radial distribution of dust, stars, gas, and star-formation rate (SFR) in a sub-sample of 18 face-on spiral galaxies extracted from the ...DustPedia sample. Methods. This study is performed by exploiting the multi-wavelength DustPedia database, from ultraviolet (UV) to sub-millimeter bands, in addition to molecular (12CO) and atomic (Hi) gas maps and metallicity abundance information available in the literature. We fitted the surface-brightness profiles of the tracers of dust and stars, the mass surface-density profiles of dust, stars, molecular gas, and total gas, and the SFR surface-density profiles with an exponential curve and derived their scale-lengths. We also developed a method to solve for the CO-to-H2 conversion factor (αCO) per galaxy by using dust- and gas-mass profiles. Results. Although each galaxy has its own peculiar behavior, we identified a common trend of the exponential scale-lengths versus wavelength. On average, the scale-lengths normalized to the B-band 25 mag/arcsec2 radius decrease from UV to 70 μm, from 0.4 to 0.2, and then increase back up to ~0.3 at 500 microns. The main result is that, on average, the dust-mass surface-density scale-length is about 1.8 times the stellar one derived from IRAC data and the 3.6 μm surface brightness, and close to that in the UV. We found a mild dependence of the scale-lengths on the Hubble stage T: the scale-lengths of the Herschel bands and the 3.6 μm scale-length tend to increase from earlier to later types, the scale-length at 70 μm tends to be smaller than that at longer sub-mm wavelength with ratios between longer sub-mm wavelengths and 70 μm that decrease with increasing T. The scale-length ratio of SFR and stars shows a weak increasing trend towards later types. Our αCO determinations are in the range (0.3−9) M⊙ pc-2 (K km s-1)-1, almost invariant by using a fixed dust-to-gas ratio mass (DGR) or a DGR depending on metallicity gradient.
ABSTRACT
Core-collapse supernovae (CCSNe) produce large ($\gtrsim0.1\,{\rm M}_\odot$) masses of dust, and are potentially the primary source of dust in the Universe, but much of this dust may be ...destroyed before reaching the interstellar medium. Cassiopeia A (Cas A) is the only supernova remnant where an observational measurement of the dust destruction efficiency in the reverse shock is possible at present. We determine the pre- and post-shock dust masses in Cas A using a substantially improved dust emission model. In our preferred models, the unshocked ejecta contains $0.6\!-\!0.8\,{\rm M}_\odot$ of $0.1\,{\rm \mu m}$ silicate grains, while the post-shock ejecta has $0.02\!-\!0.09\,{\rm M}_\odot$ of $5\!-\!10 \, {\rm nm}$ grains in dense clumps, and $2 \times 10^{-3}\,{\rm M}_\odot$ of $0.1 \, {\rm \mu m}$ grains in the diffuse X-ray emitting shocked ejecta. The implied dust destruction efficiency is $74\!-\!94\,{\rm per\,cent}$ in the clumps and $92\!-\!98\,{\rm per\,cent}$ overall, giving Cas A a final dust yield of $0.05\!-\!0.30\,{\rm M}_\odot$. If the unshocked ejecta grains are larger than $0.1\,{\rm \mu m}$, the dust masses are higher, the destruction efficiencies are lower, and the final yield may exceed $0.5\,{\rm M}_\odot$. As Cas A has a dense circumstellar environment and thus a much stronger reverse shock than is typical, the average dust destruction efficiency across all CCSNe is likely to be lower, and the average dust yield higher. This supports a mostly stellar origin for the cosmic dust budget.
Aims. The DustPedia project is capitalising on the legacy of the Herschel Space Observatory, using cutting-edge modelling techniques to study dust in the 875 DustPedia galaxies – representing the ...vast majority of extended galaxies within 3000 km s-1 that were observed by Herschel. This work requires a database of multiwavelength imagery and photometry that greatly exceeds the scope (in terms of wavelength coverage and number of galaxies) of any previous local-Universe survey. Methods. We constructed a database containing our own custom Herschel reductions, along with standardised archival observations from GALEX, SDSS, DSS, 2MASS, WISE, Spitzer, and Planck. Using these data, we performed consistent aperture-matched photometry, which we combined with external supplementary photometry from IRAS and Planck. Results. We present our multiwavelength imagery and photometry across 42 UV-microwave bands for the 875 DustPedia galaxies. Our aperture-matched photometry, combined with the external supplementary photometry, represents a total of 21 857 photometric measurements. A typical DustPedia galaxy has multiwavelength photometry spanning 25 bands. We also present the Comprehensive & Adaptable Aperture Photometry Routine (CAAPR), the pipeline we developed to carry out our aperture-matched photometry. CAAPR is designed to produce consistent photometry for the enormous range of galaxy and observation types in our data. In particular, CAAPR is able to determine robust cross-compatible uncertainties, thanks to a novel method for reliably extrapolating the aperture noise for observations that cover a very limited amount of background. Our rich database of imagery and photometry is being made available to the community.
The luminous Type IIn SN 2010jl shows strong signs of interaction between the SN ejecta and dense circumstellar material. Dust may be present in the unshocked ejecta; the cool, dense shell (CDS) ...between the shocks in the interaction region; or in the circumstellar medium (CSM). We present and model new optical and infrared photometry and spectroscopy of SN 2010jl from 82 to 1367 days since explosion. We evaluate the photometric and spectroscopic evolution using the radiative transfer codes mocassin and damocles, respectively. We propose an interaction scenario and investigate the resulting dust formation scenarios and dust masses. We find that SN 2010jl has been continuously forming dust based on the evolution of its infrared emission and optical spectra. There is evidence for preexisting dust in the CSM as well as new dust formation in the CDS and/or ejecta. We estimate that 0.005-0.01 M of predominantly carbon dust grains has formed in SN 2010jl by ∼1400 days post-outburst.
ABSTRACT
Grain growth by accretion of gas-phase metals is a common assumption in models of dust evolution, but in dense gas, where the time-scale is short enough for accretion to be effective, ...material is accreted in the form of ice mantles rather than adding to the refractory grain mass. It has been suggested that negatively charged small grains in the diffuse interstellar medium (ISM) can accrete efficiently due to the Coulomb attraction of positively-charged ions, avoiding this issue. We show that this inevitably results in the growth of the small-grain radii until they become positively charged, at which point further growth is effectively halted. The resulting gas-phase depletions under diffuse ISM conditions are significantly overestimated when a constant grain size distribution is assumed. While observed depletions can be reproduced by changing the initial size distribution or assuming highly efficient grain shattering, both options result in unrealistic levels of far-ultraviolet extinction. We suggest that the observed elemental depletions in the diffuse ISM are better explained by higher initial depletions, combined with inefficient dust destruction by supernovae at moderate ($n_{\rm H}\sim 30 \, {\rm cm}^{-3}$) densities, rather than by higher accretion efficiences.
ABSTRACT
Dust destruction by supernovae is one of the main processes removing dust from the interstellar medium (ISM). Estimates of the efficiency of this process, both theoretical and observational, ...typically assume a shock propagating into a homogeneous medium, whereas the ISM possesses significant substructure in reality. We self-consistently model the dust and gas properties of the shocked ISM in three supernova remnants (SNRs), using X-ray and infrared (IR) data combined with corresponding emission models. Collisional heating by gas with properties derived from X-ray observations produces dust temperatures too high to fit the far-IR fluxes from each SNR. An additional colder dust component is required, which has a minimum mass several orders of magnitude larger than that of the warm dust heated by the X-ray emitting gas. Dust-to-gas mass ratios indicate that the majority of the dust in the X-ray emitting material has been destroyed, while the fraction of surviving dust in the cold component is plausibly close to unity. As the cold component makes up virtually all the total dust mass, destruction time-scales based on homogeneous models, which cannot account for multiple phases of shocked gas and dust, may be significantly overestimating actual dust destruction efficiencies, and subsequently underestimating grain lifetimes.
In high density environments, the gas content of galaxies is stripped, leading to a rapid quenching of their star formation activity. This dramatic environmental effect, which is not related to ...typical passive evolution, is generally not taken into account in the star formation histories (SFHs) usually assumed to perform spectral energy distribution (SED) fitting of these galaxies, yielding a poor fit of their stellar emission and, consequently, biased estimate of the star formation rate (SFR). In this work, we aim at reproducing this rapid quenching using a truncated delayed SFH that we implemented in the SED fitting code CIGALE. We show that the ratio between the instantaneous SFR and the SFR just before the quenching (rSFR) is well constrained as long as rest-frame UV data are available. This SED modeling is applied to the Herschel Reference Survey (HRS) containing isolated galaxies and sources falling in the dense environment of the Virgo cluster. The latter are Hi-deficient because of ram pressure stripping. We show that the truncated delayed SFH successfully reproduces their SED, while typical SFH assumptions fail. A good correlation is found between rSFR and Hi−def, the parameter that quantifies the gas deficiency of cluster galaxies, meaning that SED fitting results can be used to provide a tentative estimate of the gas deficiency of galaxies for which Hi observations are not available. The HRS galaxies are placed on the SFR-M∗ diagram showing that the Hi-deficient sources lie in the quiescent region, thus confirming previous studies. Using the rSFR parameter, we derive the SFR of these sources before quenching and show that they were previously on the main sequence relation. We show that the rSFR parameter is also recovered well for deeply obscured high redshift sources, as well as in the absence of IR data. SED fitting is thus a powerful tool for identifying galaxies that underwent a rapid star formation quenching.