Linking dust emission to fundamental properties in galaxies: The low-metallicity picture

In this work, we aim at providing a consistent analysis of the dust properties from metal-poor to metal-rich environments by linking them to fundamental galactic parameters. We consider two samples of galaxies: the Dwarf Galaxy Survey (DGS) and KINGFISH, totalling 109 galaxies, spanning almost 2 dex in metallicity. We collect infrared (IR) to submillimetre (submm) data for both samples and present the complete data set for the DGS sample. We model the observed spectral energy distributions (SED) with a physically-motivated dust model to access the dust properties. Using a different SED model (modified blackbody), dust composition (amorphous carbon), or wavelength coverage at submm wavelengths results in differences in the dust mass estimate of a factor two to three, showing that this parameter is subject to non-negligible systematic modelling uncertainties. For eight galaxies in our sample, we find a rather small excess at 500 microns (< 1.5 sigma). We find that the dust SED of low-metallicity galaxies is broader and peaks at shorter wavelengths compared to more metal-rich systems, a sign of a clumpier medium in dwarf galaxies. The PAH mass fraction and the dust temperature distribution are found to be driven mostly by the specific star-formation rate, SSFR, with secondary effects from metallicity. The correlations between metallicity and dust mass or total-IR luminosity are direct consequences of the stellar mass-metallicity relation. The dust-to-stellar mass ratios of metal-rich sources follow the well-studied trend of decreasing ratio for decreasing SSFR. The relation is more complex for highly star-forming low-metallicity galaxies and depends on the chemical evolutionary stage of the source (i.e., gas-to-dust mass ratio). Dust growth processes in the ISM play a key role in the dust mass build-up with respect to the stellar content at high SSFR and low metallicity. (abridged)Comment: 44 pages (20 pages main body plus 5 Appendices), 11 figures, 9 tables, accepted for publication in A&

Herschel SPIRE-FTS Observations of Excited CO and [CI] in the Antennae (NGC 4038/39): Warm and Cold Molecular Gas

We present Herschel SPIRE-FTS observations of the Antennae (NGC 4038/39), a well studied, nearby ($22$ Mpc) ongoing merger between two gas rich spiral galaxies. We detect 5 CO transitions ($J=4-3$ to $J=8-7$), both [CI] transitions and the [NII]$205\mu m$ transition across the entire system, which we supplement with ground based observations of the CO $J=1-0$, $J=2-1$ and $J=3-2$ transitions, and Herschel PACS observations of [CII] and [OI]$63\mu m$. Using the CO and [CI] transitions, we perform both a LTE analysis of [CI], and a non-LTE radiative transfer analysis of CO and [CI] using the radiative transfer code RADEX along with a Bayesian likelihood analysis. We find that there are two components to the molecular gas: a cold ($T_{kin}\sim 10-30$ K) and a warm ($T_{kin} \gtrsim 100$ K) component. By comparing the warm gas mass to previously observed values, we determine a CO abundance in the warm gas of $x_{CO} \sim 5\times 10^{-5}$. If the CO abundance is the same in the warm and cold gas phases, this abundance corresponds to a CO $J=1-0$ luminosity-to-mass conversion factor of $\alpha_{CO} \sim 7 \ M_{\odot}{pc^{-2} \ (K \ km \ s^{-1})^{-1}}$in the cold component, similar to the value for normal spiral galaxies. We estimate the cooling from H$_2$, [CII], CO and [OI]$63\mu m$to be$\sim 0.01 L_{\odot}/M_{\odot}$. We compare PDR models to the ratio of the flux of various CO transitions, along with the ratio of the CO flux to the far-infrared flux in NGC 4038, NGC 4039 and the overlap region. We find that the densities recovered from our non-LTE analysis are consistent with a background far-ultraviolet field of strength$G_0\sim 1000$. Finally, we find that a combination of turbulent heating, due to the ongoing merger, and supernova and stellar winds are sufficient to heat the molecular gas.Comment: 50 pages, 15 figures, 8 tables, Accepted for publication in The Astrophysical Journa

Gas-to-Dust mass ratios in local galaxies over a 2 dex metallicity range

This paper analyses the behaviour of the gas-to-dust mass ratio (G/D) of local Universe galaxies over a large metallicity range. We combine three samples: the Dwarf Galaxy Survey, the KINGFISH survey and a subsample from Galametz et al. (2011) totalling 126 galaxies, covering a 2 dex metallicity range, with 30% of the sample with 12+log(O/H) < 8.0. The dust masses are homogeneously determined with a semi-empirical dust model, including submm constraints. The atomic and molecular gas masses are compiled from the literature. Two XCO are used to estimate molecular gas masses: the Galactic XCO, and a XCO depending on the metallicity (as Z^{-2}). Correlations with morphological types, stellar masses, star formation rates and specific star formation rates are discussed. The trend between G/D and metallicity is empirically modelled using power-laws (slope of -1 and free) and a broken power-law. We compare the evolution of the G/D with predictions from chemical evolution models. We find that out of the five tested galactic parameters, metallicity is the galactic property driving the observed G/D. The G/D versus metallicity relation cannot be represented by a power-law with a slope of -1 over the whole metallicity range. The observed trend is steeper for metallicities lower than ~ 8.0. A large scatter is observed in the G/D for a given metallicity, with a dispersion of 0.37 dex in metallicity bins of ~0.1 dex. The broken power-law reproduces best the observed G/D and provides estimates of the G/D that are accurate to a factor of 1.6. The good agreement of the G/D and its scatter with the three tested chemical evolution models shows that the scatter is intrinsic to galactic properties, reflecting the different star formation histories, dust destruction efficiencies, dust grain size distributions and chemical compositions across the sample. (abriged)Comment: 23 pages, 12 figures, accepted in Astronomy & Astrophysic

Quantifying the Heating Sources for Mid-infrared Dust Emissions in Galaxies: The Case of M 81

With the newly available photometric images at 250 and 500 μm from the Herschel Space Observatory, we study quantitative correlations over a sub-kiloparsec scale among three distinct emission components in the interstellar medium of the nearby spiral galaxy M 81 (NGC 3031): (1) I8 or I24, the surface brightness of the mid-infrared emission observed in the Spitzer Space Telescope 8 or 24 μm band, with I8 and I24 being dominated by the emissions from polycyclic aromatic hydrocarbons (PAHs) and very small grains (VSGs) of dust, respectively; (2) I500, that of the cold dust continuum emission in the Herschel Space Observatory 500 μm band, dominated by the emission from large dust grains heated by evolved stars; and (3) IHα, a nominal surface brightness of the Hα line emission, from gas ionized by newly formed massive stars. The results from our correlation study, free from any assumption on or modeling of dust emissivity law or dust temperatures, present solid evidence for significant heating of PAHs and VSGs by evolved stars. In the case of M 81, about 67% (48%) of the 8 μm (24 μm ) emission derives its heating from evolved stars, with the remainder attributed to radiation heating associated with ionizing stars

The Herschel Dwarf Galaxy Survey: I. Properties of the low-metallicity ISM from PACS spectroscopy

International audienceContext. The far-infrared (FIR) lines are important tracers of the cooling and physical conditions of the interstellar medium (ISM) and are rapidly becoming workhorse diagnostics for galaxies throughout the universe. There are clear indications of a different behavior of these lines at low metallicity that needs to be explored. Aims. Our goal is to explain the main differences and trends observed in the FIR line emission of dwarf galaxies compared to more metal-rich galaxies, and how this translates in ISM properties. Methods. We present Herschel/PACS spectroscopic observations of the [C ii] 157 μm, [O i] 63 and 145 μm, [O iii] 88 μm, [N ii] 122 and 205 μm, and [N iii] 57 μm fine-structure cooling lines in a sample of 48 low-metallicity star-forming galaxies of the guaranteed time key program Dwarf Galaxy Survey. We correlate PACS line ratios and line-to-LTIR ratios with LTIR, LTIR/LB, metallicity, and FIR color, and interpret the observed trends in terms of ISM conditions and phase filling factors with Cloudy radiative transfer models. Results. We find that the FIR lines together account for up to 3 percent of LTIR and that star-forming regions dominate the overall emission in dwarf galaxies. Compared to metal-rich galaxies, the ratios of [O iii]88/[N ii]122 and [N iii]57/[N ii]122 are high, indicative of hard radiation fields. In the photodissociation region (PDR), the [C ii]157/[O i]63 ratio is slightly higher than in metal-rich galaxies, with a small increase with metallicity, and the [O i]145/[O i]63 ratio is generally lower than 0.1, demonstrating that optical depth effects should be small on the scales probed. The [O iii]88/[O i]63 ratio can be used as an indicator of the ionized gas/PDR filling factor, and is found to be ~4 times higher in the dwarfs than in metal-rich galaxies. The high [C ii]/LTIR, [O i]/LTIR, and [O iii]/LTIR ratios, which decrease with increasing LTIR and LTIR/LB, are interpreted as a combination of moderate far-UV fields and a low PDR covering factor. Harboring compact phases of a low filling factor and a large volume filling factor of diffuse gas, the ISM of low-metallicity dwarf galaxies has a more porous structure than that of metal-rich galaxies

Gathering dust : A galaxy-wide study of dust emission from cloud complexes in NGC 300

© 2018 ESO. Reproduced with permission from Astronomy & Astrophysics. Content in the UH Research Archive is made available for personal research, educational, and non-commercial purposes only. Unless otherwise stated, all content is protected by copyright, and in the absence of an open license, permissions for further re-use should be sought from the publisher, the author, or other copyright holder.Aims. We use multi-band observations by the Herschel Space Observatory to study the dust emission properties of the nearby spiral galaxy NGC 300. We compile a first catalogue of the population of giant dust clouds (GDCs) in NGC 300, including temperature and mass estimates, and give an estimate of the total dust mass of the galaxy. Methods. We carried out source detection with the multiwavelength source extraction algorithm getsources. We calculated physical properties, including mass and temperature, of the GDCs from five-band Herschel PACS and SPIRE observations from 100 to 500 μm; the final size and mass estimates are based on the observations at 250 μm that have an effective spatial resolution of ~170 pc. We correlated our final catalogue of GDCs to pre-existing catalogues of HII regions to infer the number of GDCs associated with high-mass star formation and determined the Hα emission of the GDCs. Results. Our final catalogue of GDCs includes 146 sources, 90 of which are associated with known HII regions. We find that the dust masses of the GDCs are completely dominated by the cold dust component and range from ~1.1 × 10 3 to 1.4 × 10 4 M. The GDCs have effective temperatures of ~13-23 K and show a distinct cold dust effective temperature gradient from the centre towards the outer parts of the stellar disk. We find that the population of GDCs in our catalogue constitutes ~16% of the total dust mass of NGC 300, which we estimate to be about 5.4 × 10 6 M. At least about 87% of our GDCs have a high enough average dust mass surface density to provide sufficient shielding to harbour molecular clouds. We compare our results to previous pointed molecular gas observations in NGC 300 and results from other nearby galaxies and also conclude that it is very likely that most of our GDCs are associated with complexes of giant molecular clouds.Peer reviewe

Spatially resolved physical conditions of molecular gas and potential star formation tracers in M 83, revealed by the Herschel SPIRE FTS

International audienceWe investigate the physical properties of the molecular and ionized gas, and their relationship to the star formation and dust properties in M 83, based on submillimeter imaging spectroscopy from within the central 3.5′ (~4 kpc in diameter) around the starburst nucleus. The observations use the Fourier Transform Spectrometer (FTS) of the Spectral and Photometric Imaging REceiver (SPIRE) onboard the Herschel Space Observatory. The newly observed spectral lines include [CI] 370 μm, [CI] 609 μm, [NII] 205 μm, and CO transitions from J = 4−3 to J = 13−12. Combined with previously observed J = 1−0 to J = 3−2 transitions, the CO spectral line energy distributions are translated to spatially resolved physical parameters, column density of CO, N(CO), and molecular gas thermal pressure, Pth, with a non-local thermal equilibrium (non-LTE) radiative transfer model, RADEX. Our results show that there is a relationship between the spatially resolved intensities of [NII] 205 μm and the surface density of the star formation rate (SFR), ΣSFR. This relation, when compared to integrated properties of ultra-luminous infrared galaxies (ULIRGs), exhibits a different slope, because the [NII] 205 μm distribution is more extended than the SFR. The spatially resolved [CI] 370 μm, on the other hand, shows a generally linear relationship with ΣSFR and can potentially be a good SFR tracer. Compared with the dust properties derived from broad-band images, we find a positive trend between the emissivity of CO in the J = 1−0 transition with the average intensity of interstellar radiation field (ISRF), ⟨ U ⟩. This trend implies a decrease in the CO-to-H2 conversion factor, XCO, when ⟨ U ⟩ increases. We estimate the gas-to-dust mass ratios to be 77 ± 33 within the central 2 kpc and 93 ± 19 within the central 4 kpc of M 83, which implies a Galactic dust-to-metal mass ratio within the observed region of M 83. The estimated gas-depletion time for the M 83 nucleus is 1.13 ± 0.6 Gyr, which is shorter than the values for nearby spiral galaxies found in the literature (~2.35 Gyr), most likely due to the young nuclear starbursts. A linear relationship between Pth and the radiation pressure generated by ⟨ U ⟩, Prad, is found to be Pth ≈ 30 Prad, which signals that the ISRF alone is insufficient to sustain the observed CO transitions. The spatial distribution of Pth reveals a pressure gradient, which coincides with the observed propagationof starburst activities and the alignment of (possibly background) radio sources. We discover that the off-centered (from the optical nucleus) peak of the molecular gas volume density coincides well with a minimum in the relative aromatic feature strength, indicating a possible destruction of their carriers. We conclude that the observed CO transitions are most likely associated with mechanical heating processes that are directly or indirectly related to very recent nuclear starbursts

High-resolution, 3D radiative transfer modeling: I. The grand-design spiral galaxy M 51

International audienceContext. Dust reprocesses about half of the stellar radiation in galaxies. The thermal re-emission by dust of absorbed energy is considered to be driven merely by young stars so is often applied to tracing the star formation rate in galaxies. Recent studies have argued that the old stellar population might be responsible for 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 M 51 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 to 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 to account for the absorption, scattering, and non-local thermal equilibrium (NLTE) emission of dust in M 51. The 3D distribution of stars is derived from the 2D morphology observed in the IRAC 3.6 μm, GALEX FUV, Hα, and MIPS 24 μ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.Results. The dust emission derived from RT calculations is consistent with far-infrared and submillimeter observations of M 51, implying that the absorbed stellar energy is balanced by the thermal re-emission of dust. The young stars provide 63% of the energy for heating the dust responsible for the total infrared emission (8−1000 μm), while 37% of the dust emission is governed through heating by the evolved stellar population. In individual wavebands, the contribution from young stars to the dust heating dominates at all infrared wavebands but gradually decreases towards longer infrared and submillimeter wavebands for which the old stellar population becomes a non-negligible source of heating. Upon extrapolation of the results for M 51, we present prescriptions for estimating the contribution of young stars to the global dust heating based on a tight correlation between the dust heating fraction and specific star formation rate

Dust spectral energy distributions of nearby galaxies: an insight from the Herschel Reference Survey

International audienceAlthough it accounts only for a small fraction of the baryonic mass, dust has a profound impact on the physical processes at play in galaxies. Thus, to understand the evolution of galaxies, it is essential not only to characterize dust properties per se, but also in relation to global galaxy properties. To do so, we derive the dust properties of galaxies in a volume limited, K-band selected sample, the Herschel Reference Survey (HRS). We gather infrared photometric data from 8 μm to 500 μm from Spitzer, WISE, IRAS, and Herschel for all of the HRS galaxies. Draine & Li (2007, ApJ, 663, 866) models are fit to the data from which the stellar contribution has been carefully removed. We find that our photometric coverage is sufficient to constrain all of the parameters of the Draine & Li models and that a strong constraint on the 20−60 μm range is mandatory to estimate the relative contribution of the photo-dissociation regions to the infrared spectral energy distribution (SED). The SED models tend to systematically underestimate the observed 500 μm flux densities, especially for low-mass systems. We provide the output parameters for all of the galaxies, i.e., the minimum intensity of the interstellar radiation field, the fraction of polycyclic aromatic hydrocarbon (PAH), the relative contribution of PDR and evolved stellar population to the dust heating, the dust mass, and the infrared luminosity. For a subsample of gas-rich galaxies, we analyze the relations between these parameters and the main integrated properties of galaxies, such as stellar mass, star formation rate, infraredluminosity, metallicity, Hα and H-band surface brightness, and the far-ultraviolet attenuation. A good correlation between the fraction of PAH and the metallicity is found, implying a weakening of the PAH emission in galaxies with low metallicities and, thus, low stellar masses. The intensity of the diffuse interstellar radiation field and the H-band and Hα surface brightnesses are correlated, suggesting that the diffuse dust component is heated by both the young stars in star-forming regions and the diffuse evolved population. We use these results to provide a new set of infrared templates calibrated with Herschel observations on nearby galaxies and a mean SED template to provide the z = 0 reference for cosmological studies. For the same purpose, we place our sample on the SFR − M∗ diagram. The templates are compared to the most popular infrared SED libraries, enlightening a large discrepancy between all of them in the 20−100 μm range

The evolution of the dust and gas content in galaxies

We use deep Herschel observations taken with both PACS and SPIRE imaging cameras to estimate the dust mass of a sample of galaxies extracted from the GOODS-S, GOODS-N and the COSMOS fields. We divide the redshift–stellar mass (M star )–star formation rate (SFR) parameter space into small bins and investigate average properties over this grid. In the first part of the work we investigate the scaling relations between dust mass, stellar mass and SFR out to z = 2.5. No clear evolution of the dust mass with redshift is observed at a given SFR and stellar mass. We find a tight correlation between the SFR and the dust mass, which, under reasonable assumptions, is likely a consequence of the Schmidt-Kennicutt (S-K) relation. The previously observed correlation between the stellar content and the dust content flattens or sometimes disappears when considering galaxies with the same SFR. Our finding suggests that most of the correlation between dust mass and stellar mass obtained by previous studies is likely a consequence of the correlation between the dust mass and the SFR combined with the main sequence, i.e., the tight relation observed between the stellar mass and the SFR and followed by the majority of star-forming galaxies. We then investigate the gas content as inferred from dust mass measurements. We convert the dust mass into gas mass by assuming that the dust-to-gas ratio scales linearly with the gas metallicity (as supported by many observations). For normal star-forming galaxies (on the main sequence) the inferred relation between the SFR and the gas mass (integrated S-K relation) broadly agrees with the results of previous studies based on CO measurements, despite the completely different approaches. We observe that all galaxies in the sample follow, within uncertainties, the same S-K relation. However, when investigated in redshift intervals, the S-K relation shows a moderate, but significant redshift evolution. The bulk of the galaxy population at z ∼ 2 converts gas into stars with an efficiency (star formation efficiency, SFE = SFR/M gas , equal to the inverse of the depletion time) about 5 times higher than at z ∼ 0. However, it is not clear what fraction of such variation of the SFE is due to an intrinsic redshift evolution and what fraction is simply a consequence of high-z galaxies having, on average, higher SFR, combined with thesuper-linear slope of the S-K relation (whileother studies finda linear slope). We confirm that the gas fraction (f gas = M gas /(M gas + M star )) decreases with stellar mass and increases with the SFR. We observe no evolution with redshift once M star and SFR are fixed. We explain these trends by introducing a universal relation between gas fraction, stellar mass and SFR that does not evolve with redshift, at least out to z ∼ 2.5. Galaxies move across this relation as their gas content evolves across the cosmic epochs. We use the 3D fundamental f gas –M star –SFR relation, along with the evolution of the main sequence with redshift, to estimate the evolution of the gas fraction in the average population of galaxies as a function of redshift and as a function of stellar mass: we find that M star > ∼ 10 11 M ? galaxies show the strongest evolution at z > ∼ 1.3 and a flatter trend at lower redshift, while f gas decreases more regularly over the entire redshift range probed in M star < ∼ 10 11 Mo galaxies, in agreement with a downsizing scenario