Search

In the framework of the IMEGIN Large Program, we used the NIKA2 camera on the IRAM 30-m telescope to observe the edge-on galaxy NGC 891 at 1.15 mm and 2 mm and at a FWHM of 11.1" and 17.6", respectively. Multiwavelength data enriched with the new NIKA2 observations fitted by the HerBIE SED code (coupled with the THEMIS dust model) were used to constrain the physical properties of the ISM. Emission originating from the diffuse dust disk is detected at all wavelengths from mid-IR to mm, while mid-IR observations reveal warm dust emission from compact HII regions. Indications of mm excess emission have also been found in the outer parts of the galactic disk. Furthermore, our SED fitting analysis constrained the mass fraction of the small (< 15 Angstrom) dust grains. We found that small grains constitute 9.5% of the total dust mass in the galactic plane, but this fraction increases up to ~ 20% at large distances (|z| > 3 kpc) from the galactic plane., Comment: To appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

As the millimeter wavelength range remains a largely unexplored spectral region for galaxies, the IMEGIN large program aims to map the millimeter continuum emission of 22 nearby galaxies at 1.15 and 2 mm. Using the high-resolution maps produced by the NIKA2 camera, we explore the existence of very cold dust and take possible contamination by free-free and synchrotron emission into account. We study the IR-to-radio emission coming from different regions along the galactic plane and at large vertical distances. New observations of NGC 891, using the NIKA2 camera on the IRAM 30m telescope, along with a suite of observations at other wavelengths were used to perform a multiwavelength study of the spectral energy distribution in the interstellar medium in this galaxy. This analysis was performed globally and locally, using the advanced hierarchical Bayesian fitting code, HerBIE, coupled with the THEMIS dust model. Our dust modeling is able to reproduce the near-IR to millimeter emission of NGC 891, with the exception of an excess at a level of 25% obtained by the NIKA2 observations in the outermost parts of the disk. The radio continuum and thermal dust emission are distributed differently in the disk and galaxy halo. Different dusty environments are also revealed by a multiwavelength investigation of the emission features. Our detailed decomposition at millimeter and centimeter wavelengths shows that emission at 1 mm is purely originated by dust. Radio components become progressively important with increasing wavelengths. Finally, we find that emission arising from small dust grains accounts for ~ 9.5% of the total dust mass, reaching up to 20% at large galactic latitudes. Shock waves in the outflows that shatter the dust grains might explain this higher fraction of small grains in the halo., Comment: 15 pages, 10 figures, Accepted for publication in A&A

Interstellar dust plays an important role in the formation of molecular gas and the heating and cooling of the interstellar medium. The spatial distribution of the mm-wavelength dust emission from galaxies is largely unexplored. The NIKA2 Guaranteed Time Project IMEGIN (Interpreting the Millimeter Emission of Galaxies with IRAM and NIKA2) has recently mapped the mm emission in the grand design spiral galaxy NGC6946. By subtracting the contributions from the free-free, synchrotron, and CO line emission, we map the distribution of the pure dust emission at 1:15mm and 2mm. Separating the arm/interarm regions, we find a dominant 2mm emission from interarms indicating the significant role of the general interstellar radiation field in heating the cold dust. Finally, we present maps of the dust mass, temperature, and emissivity index using the Bayesian MCMC modeling of the spectral energy distribution in NGC6946., Comment: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

The mm-to-cm range of the Spectral Energy Distribution of spiral galaxies remains largely unexplored. Its coverage is required to disentangle the contribution of dust emission, free-free and synchrotron radiation and can provide constraints on dust models, star-formation rates and ISM properties. We present the case for a synergy between NIKA2 observations of nearby spirals and those from planned and current instrumentation at the Sardinia Radio Telescope, and report on a pilot K-band program to search for Anomalous Microwave Emission, an elusive emission component which is presumably related to dust., Comment: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

New observations of the edge-on galaxy NGC 891, at 1.15 and 2 mm obtained with the IRAM 30-m telescope and the NIKA2 camera, within the framework of the IMEGIN (Interpreting the Millimetre Emission of Galaxies with IRAM and NIKA2) Large Program, are presented in this work. By using multiwavelength maps (from the mid-IR to the cm wavelengths) we perform SED fitting in order to extract the physical properties of the galaxy on both global and local ($\sim$kpc) scales. For the interpretation of the observations we make use of a state-of-the-art SED fitting code, HerBIE (HiERarchical Bayesian Inference for dust Emission). The observations indicate a galaxy morphology, at mm wavelengths, similar to that of the cold dust emission traced by sub-mm observations and to that of the molecular gas. The contribution of the radio emission at the NIKA2 bands is very small (negligible at 1.15 mm and $\sim10\%$ at 2 mm) while it dominates the total energy budget at longer wavelengths (beyond 5 mm). On local scales, the distribution of the free-free emission resembles that of the dust thermal emission while the distribution of the synchrotron emission shows a deficiency along the major axis of the disc of the galaxy., Comment: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

We present the results of high spatial resolution HCO$^{+}$($1-0$) and HCN($1-0$) observations of N55 south region (N55-S) in the Large Magellanic Cloud (LMC), obtained with the Atacama Large Millimeter/submillimeter Array (ALMA). N55-S is a relatively less extreme star-forming region of the LMC characterized by a low radiation field. We carried out a detailed analysis of the molecular emission to investigate the relation between dense molecular clumps and star formation in the quiescent environment of N55-S. We detect ten molecular clumps with significant HCO$^{+}(1-0)$ emission and eight with significant HCN($1-0$) emission, and estimate the molecular clump masses by virial and local thermodynamic equilibrium analysis. All identified young stellar objects (YSOs) in the N55-S are found to be near the HCO$^{+}$ and HCN emission peaks showing the association of these clumps with recent star formation activity. The molecular clumps that have associated YSOs show relatively larger linewidths and masses than those without YSOs. We compare the clump properties of the N55-S with those of other giant molecular clouds (GMCs) in the LMC and find that N55-S clumps possess similar size but relatively lower linewidth and larger HCN/HCO$^{+}$(1$-$0) flux ratio. These results can be attributed to the low radiation field in N55-S resulted by relatively low star formation activity compared to other active star-forming regions like 30Doradus-10 and N159. The dense gas fraction of N55-S is $\sim$ 0.025, lower compared to other GMCs of the LMC supporting the low star formation efficiency of this region., Comment: Accepted for publication in ApJ, 13 pages, 4 figures

While the CO(1-0) transition is often used to deduce the total molecular hydrogen in galaxies, it is challenging to detect in low metallicity galaxies, in spite of the star formation taking place. In contrast, the [CII] 158 micron line is relatively bright, highlighting a potentially important reservoir of H2 that is not traced by CO(1-0), but residing in the C+ - emitting regions. We explore a method to quantify the total H2 mass (MH2) in galaxies and learn what parameters control the CO-dark gas reservoir. We present Cloudy grids of density, radiation field and metallicity in terms of observed quantities, such as [OI], [CI], CO(1-0), [CII], total infrared luminosity and the total MH2 and provide recipes based on these models to derive total MH2 mass estimates from observations. The models are applied to the Herschel Dwarf Galaxy Survey, extracting the total MH2 for each galaxy which is compared to the H2 determined from the observed CO(1-0) line. While the H2 traced by CO(1-0) can be negligible, the [CII] 158 micron line can trace the total H2. 70% to 100% of the total H2 mass is not traced by CO(1-0) in the dwarf galaxies, but is well-traced by [CII] 158 micron line. The CO-dark gas mass fraction correlates with the observed L[CII]/LCO(1-0) ratio. A conversion factor for [CII] luminosity to total H2 and a new CO-to-total-MH2 conversion factor, as a function of metallicity, is presented. A recipe is provided to quantify the total mass of H2 in galaxies, taking into account the CO and [CII] observations. Accounting for this CO-dark H2 gas, we find that the star forming dwarf galaxies now fall on the Schmidt-Kennicutt relation. Their star-forming efficiency is rather normal, since the reservoir from which they form stars is now more massive when introducing the [CII] measures of the total H2, compared to the little amount of H2 in the CO-emitting region., Comment: Accepted for publication in A&A

Investigating star formation requires precise knowledge of the properties of the dense molecular gas. The low metallicity and wide range of star formation activity of the Large and Small Magellanic Clouds make them prime laboratories to study how local physical conditions impact the dense gas reservoirs. The aim of the Dense Gas Survey for the Magellanic Clouds (DeGaS-MC) project is to expand our knowledge of the relation between dense gas properties and star formation activity by targeting the LMC and SMC observed in the HCO+(2-1) and HCN(2-1) transitions. We carried out a pointing survey toward 30 LMC and SMC molecular clouds using the SEPIA180 instrument installed on the APEX telescope and a follow-up mapping campaign in 13 star-forming regions. This first paper provides line characteristic catalogs and integrated line-intensity maps of the sources. HCO+(2-1) is detected in 20 and HCN(2-1) in 8 of the 29 pointings observed. The dense gas velocity pattern follows the line-of-sight velocity field derived from the stellar population. The HCN emission is less extended than the HCO+ emission. The HCO+(2-1)/HCN(2-1) brightness temperature ratios range from 1 to 7, which is consistent with the large ratios commonly observed in low-metallicity environments. A larger number of young stellar objects are found at high HCO+ intensities and lower HCO+/HCN flux ratios, and thus toward denser lines of sight. The dense gas luminosities correlate with the star formation rate traced by the total infrared luminosity over the two orders of magnitude covered by our observations, although substantial region-to-region variations are observed., Comment: Accepted for publication in A&A; 25 pages, 13 figures, 5 tables

The H$\alpha$ and H$\beta$ emission line luminosities measured in a single integrated spectrum are affected in non-trivial ways by point-to-point variations in dust attenuation in a galaxy. This work investigates the impact of this variation when estimating global H$\alpha$ luminosities corrected for the presence of dust by a global Balmer decrement. Analytical arguments show that the dust-corrected H$\alpha$ luminosity is always underestimated when using the global H$\alpha$/H$\beta$ flux ratio to correct for dust attenuation. We measure this effect on 156 face-on star-forming galaxies from the Mapping Nearby Galaxies at APO (MaNGA) survey. At 1-2 kpc spatial resolution, the effect is small but systematic, with the integrated dust-corrected H$\alpha$ luminosity underestimated by $2$-$4$ per cent (and typically not more than by $10$ per cent), and depends on the specific star formation rate of the galaxy. Given the spatial resolution of MaNGA, these are lower limits for the effect. From Multi Unit Spectroscopic Explorer (MUSE) observations of NGC 628 with a resolution of 36 pc we find the discrepancy between the globally and the point-by-point dust-corrected H$\alpha$ luminosity to be $14 \pm 1$ per cent, which may still underestimate the true effect. We use toy models and simulations to show that the true difference depends strongly on the spatial variance of the H$\alpha$/H$\beta$ flux ratio, and on the slope of the relation between H$\alpha$ luminosity and dust attenuation within a galaxy. Larger samples of higher spatial resolution observations are required to quantify the dependence of this effect as a function of galaxy properties., Comment: MNRAS accepted

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 have 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. In this particular study, we focus on NGC 1068, which contains an active galactic nucleus (AGN). The effect of diffuse dust heating by AGN (beyond the torus) was so far only investigated for quasars. This additional dust heating source further contaminates the broadband fluxes on which classic galaxy modelling tools rely 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 submm. We find a strong wavelength dependency of AGN contamination to the broadband fluxes. It peaks in the MIR, drops in the FIR, but 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 100 pc, 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 can be stronger in galaxies with a more luminous AGN. This could significantly impact the derived star formation rates and ISM masses for such systems., Comment: 18 pages, 11 figures, accepted for publication in A&A

[abridged] Understanding how the infalling gas redistribute most of its initial angular momentum inherited from prestellar cores before reaching the stellar embryo is a key question. Disk formation has been naturally considered as a possible solution to this "angular momentum problem". However, how the initial angular momentum of protostellar cores is distributed and evolves during the main accretion phase and the beginning of disk formation has largely remained unconstrained up to now. In the framework of the IRAM CALYPSO survey, we used high dynamic range C$^{18}$O (2-1) and N$_2$H$^+$ (1-0) observations to quantify the distribution of specific angular momentum along the equatorial axis in a sample of 12 Class 0 protostellar envelopes from scales ~50 to 10000 au. The radial distributions of specific angular momentum in the CALYPSO sample suggest two distinct regimes within protostellar envelopes: the specific angular momentum decreases as $j \propto r^{1.6 \pm 0.2}$ down to ~1600 au and then tends to become relatively constant around 6 $\times$ 10$^{-4}$ km s$^{-1}$ pc down to ~50 au. The values of specific angular momentum measured in the inner Class 0 envelopes, namely that of the material directly involved in the star formation process ($<$1600 au), is on the same order of magnitude as what is inferred in small T-Tauri disks. Thus, disk formation appears to be a direct consequence of angular momentum conservation during the collapse. Our analysis reveals a dispersion of the directions of velocity gradients at envelope scales $>$1600 au, suggesting that they may not be related to rotational motions of the envelopes. We conclude that the specific angular momentum observed at these scales could find its origin in core-forming motions (infall, turbulence) or trace an imprint of the initial conditions for the formation of protostellar cores., Comment: 62 pages, 66 figures

Dust and starlight are modeled for the KINGFISH project galaxies. With data from 3.6 micron to 500 micron, models are strongly constrained. For each pixel in each galaxy we estimate (1) dust surface density; (2) q_PAH, the dust mass fraction in PAHs; (3) distribution of starlight intensities heating the dust; (4) luminosity emitted by the dust; and (5) dust luminosity from regions with high starlight intensity. The models successfully reproduce both global and resolved spectral energy distributions. We provide well-resolved maps for the dust properties. As in previous studies, we find q_PAH to be an increasing function of metallicity, above a threshold Z/Z_sol approx 0.15. Dust masses are obtained by summing the dust mass over the map pixels; these "resolved" dust masses are consistent with the masses inferred from model fits to the global photometry. The global dust-to-gas ratios obtained from this study correlate with galaxy metallicities. Systems with Z/Z_sol > 0.5 have most of their refractory elements locked up in dust, whereas when Z/Z_sol < 0.3 most of these elements tend to remain in the gas phase. Within galaxies, we find that q_PAH is suppressed in regions with unusually warm dust with nu L_nu(70 um) > 0.4L_dust. With knowledge of one long-wavelength flux density ratio (e.g., f_{160}/f_{500}), the minimum starlight intensity heating the dust (U_min) can be estimated to within ~50%. For the adopted dust model, dust masses can be estimated to within ~0.07 dex accuracy using the 500 micron luminosity nu L_nu(500) alone. There are additional systematic errors arising from the choice of dust model, but these are hard to estimate. These calibrated prescriptions may be useful for studies of high-redshift galaxies., Comment: ApJ, accepted

(abridged) The ambiguous origin of [CII] 158um in the interstellar medium complicates its use for diagnostics concerning the star-formation rate and physical conditions in photodissociation regions (PDRs). We observed the giant HII region N11 in the Large Magellanic Cloud with SOFIA/GREAT in order to investigate the origin of [CII] to obtain the total H2 gas content, the fraction of CO-dark H2 gas, and the influence of environmental effects such as stellar feedback. We present an innovative spectral decomposition method that allows statistical trends to be derived. The [CII] line is resolved in velocity and compared to HI and CO, using a Bayesian approach to decompose the profiles. A simple model accounting for collisions in the neutral atomic and molecular gas was used in order to derive the H2 column density traced by C+. The profile of [CII] most closely resembles that of CO, but the integrated [CII] line width lies between that of CO and that of HI. Using various methods, we find that [CII] mostly originates from the neutral gas. We show that [CII] mostly traces the CO-dark H2 gas but there is evidence of a weak contribution from neutral atomic gas preferentially in the faintest components. Most of the molecular gas is CO-dark. The fraction of CO-dark H2 gas decreases with increasing CO column density, with a slope that seems to depend on the impinging radiation field from nearby massive stars. Finally we extend previous measurements of the photoelectric-effect heating efficiency, which we find is constant across regions probed with Herschel, with [CII] and [OI] being the main coolants in faint and diffuse, and bright and compact regions, respectively, and with PAH emission tracing the CO-dark H2 gas heating where [CII] and [OI] emit. Our study highlights the importance of velocity-resolved PDR diagnostics and higher spatial resolution for HI observations., Comment: Accepted for publication in A&A

Dust plays an important role in shaping a galaxy's spectral energy distribution (SED). It absorbs ultraviolet (UV) to near-infrared (NIR) radiation and re-emits this energy in the far-infrared (FIR). The FIR is essential to understand dust in galaxies. However, deep FIR observations require a space mission, none of which are still active today. We aim to infer the FIR emission across six Herschel bands, along with dust luminosity, mass, and effective temperature, based on the available UV to mid-infrared (MIR) observations. We also want to estimate the uncertainties of these predictions, compare our method to energy balance SED fitting, and determine possible limitations of the model. We propose a machine learning framework to predict the FIR fluxes from 14 UV-MIR broadband fluxes. We used a low redshift sample by combining DustPedia and H-ATLAS, and extracted Bayesian flux posteriors through SED fitting. We trained shallow neural networks to predict the far-infrared fluxes, uncertainties, and dust properties. We evaluated them on a test set using a root mean square error (RMSE) in log-space. Our results (RMSE = 0.19 dex) significantly outperform UV-MIR energy balance SED fitting (RMSE = 0.38 dex), and are inherently unbiased. We can identify when the predictions are off, for example when the input has large uncertainties on WISE 22, or when the input does not resemble the training set. The galaxies for which we have UV-FIR observations can be used as a blueprint for galaxies that lack FIR data. This results in a 'virtual FIR telescope', which can be applied to large optical-MIR galaxy samples. This helps bridge the gap until the next FIR mission., Comment: 25 pages, 16 figures, accepted by A&A. Interactive figures + jupyter notebooks at https://wdobbels.github.io/FIREnet/

Aims: We compare the far-infrared to sub-millimetre dust emission properties measured in high Galactic latitude cirrus with those determined in a sample of 204 late-type DustPedia galaxies. The aim is to verify if it is appropriate to use Milky Way dust properties to derive dust masses in external galaxies. Methods: We used Herschel observations and atomic and molecular gas masses to estimate the disc-averaged dust emissivity at 250 micrometres, and from this, the absorption cross section per H atom and per dust mass. The emissivity requires one assumption, which is the CO-to-H_2 conversion factor, and the dust temperature is additionally required for the absorption cross section per H atom; yet another constraint on the dust-to-hydrogen ratio D/H, depending on metallicity, is required for the absorption cross section dust mass. Results: We find epsilon(250) = 0.82 +/- 0.07 MJy sr^-1 (1E20 H cm^-2)^-1 for galaxies with 4 < F(250)/F(500) < 5. This depends only weakly on the adopted CO-to-H_2 conversion factor. The value is almost the same as that for the Milky Way at the same colour ratio. Instead, for F(250)/F(500) > 6, epsilon(250) is lower than predicted by its dependence on the heating conditions. The reduction suggests a variation in dust emission properties for spirals of earlier type, higher metallicity, and with a higher fraction of molecular gas. When the standard emission properties of Galactic cirrus are used for these galaxies, their dust masses might be underestimated by up to a factor of two. Values for the absorption cross sections at the Milky Way metallicity are also close to those of the cirrus. Mild trends of the absorption cross sections with metallicity are found, although the results depend on the assumptions made., Comment: 13 pages, 6 figures, 2 tables, A&A accepted

SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars., Comment: Accepted for publication in PASA; 37 pages, 16 figures

We use a sub-set of the DustPedia galaxy sample (461 galaxies) to investigate the effect the environment has had on galaxies. We consider Virgo cluster and field samples and also assign a density contrast parameter to each galaxy, as defined by the local density of SDSS galaxies. We consider their chemical evolution (using M_{Dust}/M_{Baryon} and M_{Gas}/M_{Baryon}), their specific star formation rate (SFR/M_{Stars}), star formation efficiency (SFR/M_{Gas}), stars-to-dust mass ratio (M_{Stars}/M_{Dust}), gas-to-dust mass ratio (M_{Gas}/M_{Dust}) and the relationship between star formation rate per unit mass of dust and dust temperature (SFR/M_{Dust} and T_{Dust}). Late type galaxies (later than Sc) in all of the environments can be modelled using simple closed box chemical evolution and a simple star formation history (SFR(t) \propto t\exp{-t/\tau}). For earlier type galaxies the physical mechanisms that give rise to their properties are clearly much more varied and require a more complicated model (mergers, gas in or outflow). However, we find little or no difference in the properties of galaxies of the same morphological type within the cluster, field or with different density contrasts. It appears that it is morphology, how and whenever this is laid down, and consistent internal physical processes that primarily determine the derived properties of galaxies in the DustPedia sample and not processes related to differences in the local environment., Comment: Accepted for publication in Astronomy and Astrophysics

The sensitive infrared telescopes, Spitzer and Herschel, have been used to target low-metallicity star-forming galaxies, allowing us to investigate the properties of their interstellar medium (ISM) in unprecedented detail. Interpretation of the observations in physical terms relies on careful modeling of those properties. We have employed a multiphase approach to model the ISM phases (HII region and photodissociation region) with the spectral synthesis code Cloudy. Our goal is to characterize the physical conditions (gas densities, radiation fields, etc.) in the ISM of the galaxies from the Herschel Dwarf Galaxy Survey. We are particularly interested in correlations between those physical conditions and metallicity or star-formation rate. Other key issues we have addressed are the contribution of different ISM phases to the total line emission, especially of the [CII]157um line, and the characterization of the porosity of the ISM. We find that the lower-metallicity galaxies of our sample tend to have higher ionization parameters and galaxies with higher specific star-formation rates have higher gas densities. The [CII] emission arises mainly from PDRs and the contribution from the ionized gas phases is small, typically less than 30% of the observed emission. We also find correlation - though with scatter - between metallicity and both the PDR covering factor and the fraction of [CII] from the ionized gas. Overall, the low metal abundances appear to be driving most of the changes in the ISM structure and conditions of these galaxies, and not the high specific star-formation rates. These results demonstrate in a quantitative way the increase of ISM porosity at low metallicity. Such porosity may be typical of galaxies in the young Universe., Comment: Accepted in A&A. 20 pages, 12 figures, 5 tables plus appendices

Within the framework of the DustPedia project we investigate the properties of cosmic dust and its interaction with the stellar radiation (originating from different stellar populations) for 814 galaxies in the nearby Universe, all observed by the Herschel Space Observatory. We take advantage of the widely used galaxy SED fitting code CIGALE, properly adapted to include the state-of-the-art dust model THEMIS. 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 as well as the unattenuated and the absorbed by dust stellar light, for both the old (> 200 Myr) and young (<= 200 Myr) stellar populations. 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 different 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 efficient 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, donating up to ~60% of their luminosity to this purpose. Finally, we find a strong correlation of the dust heating fraction by young stars with morphology and the specific star-formation rate., Comment: 22 pages, 18 figures

Observations of evolution in the dust-to-metal ratio allow us to constrain the dominant dust processing mechanisms. In this work, we present a study of the dust-to-metal and dust-to-gas ratios in a sub-sample of ~500 DustPedia galaxies. Using literature and MUSE emission line fluxes, we derived gas-phase metallicities (oxygen abundances) for over 10000 individual regions and determine characteristic metallicities for each galaxy. We study how the relative dust, gas, and metal contents of galaxies evolve by using metallicity and gas fraction as proxies for evolutionary state. The global oxygen abundance and nitrogen-to-oxygen ratio are found to increase monotonically as galaxies evolve. Additionally, unevolved galaxies (gas fraction > 60%, metallicity 12 + log(O/H) < 8.2) have dust-to-metal ratios that are about a factor of 2.1 lower (a factor of six lower for galaxies with gas fraction > 80%) than the typical dust-to-metal ratio (Md/MZ ~ 0.214) for more evolved sources. However, for high gas fractions, the scatter is larger due to larger observational uncertainties as well as a potential dependence of the dust grain growth timescale and supernova dust yield on local conditions and star formation histories. We find chemical evolution models with a strong contribution from dust grain growth describe these observations reasonably well. The dust-to-metal ratio is also found to be lower for low stellar masses and high specific star formation rates (with the exception of some sources undergoing a starburst). Finally, the metallicity gradient correlates weakly with the HI-to-stellar mass ratio, the effective radius and the dust-to-stellar mass ratio, but not with stellar mass., Comment: 24 pages, 15 figures, published in A&A, bibref: 2019A&A...623A...5D, This Version 2 has increased accuracy in Table B1 in appendix B

We present a spatially resolved study of the relation between dust and metallicity in the nearby spiral galaxies M101 (NGC 5457) and NGC 628 (M74). We explore the relation between the chemical abundances of their gas and stars with their dust content and their chemical evolution. The empirical spatially resolved oxygen effective yield and the gas to dust mass ratio (GDR) across both disc galaxies are derived, sampling one dex in oxygen abundance. We find that the metal budget of the NGC 628 disc and most of the M101 disc appears consistent with the predictions of the simple model of chemical evolution for an oxygen yield between half and one solar, whereas the outermost region (R<0.8R25) of M101 presents deviations suggesting the presence of gas flows. The GDR-metallicity relation shows a two slopes behaviour, with a break at 12+log(O/H)~8.4, a critical metallicity predicted by theoretical dust models when stardust production equals grain growth. A relation between GDR and the fraction of molecular to total gas, Sigma(H2)/Sigma(gas) is also found. We suggest an empirical relationship between GDR and the combination of 12+log(O/H), for metallicity, and Sigma(H2)/Sigma(gas), a proxy for the molecular clouds fraction. The GDR is closely related with metallicity at low abundance and with Sigma(H2)/Sigma(gas) for higher metallicities suggesting ISM dust growth. The ratio Sigma(dust)/Sigma(star) correlates well with 12 + log(O/H) and strongly with log(N/O) in both galaxies. For abundances below the critical one, the 'stardust' production gives us a constant value suggesting a stellar dust yield similar to the oxygen yield., Comment: 16 pages, 14 Figures, accepted in MNRAS

Most radiative transfer models assume that dust in spiral galaxies is distributed exponentially. In this paper our goal is to verify this assumption by analysing the two-dimensional large-scale distribution of dust in galaxies from the DustPedia sample. For this purpose, we make use of Herschel imaging in five bands, from 100 to 500{\mu}m, in which the cold dust constituent is primarily traced and makes up the bulk of the dust mass in spiral galaxies. For a subsample of 320 disc galaxies, we successfully perform a simultaneous fitting with a single S\'ersic model of the Herschel images in all five bands using the multiband modelling code GALFITM. We report that the S\'ersic index $n$, which characterises the shape of the S\'ersic profile, lies systematically below 1 in all Herschel bands and is almost constant with wavelength. The average value at 250{\mu}m is $0.67\pm0.37$ (187 galaxies are fitted with $n_{250}\leq0.75$, 87 galaxies have $0.751.25$). Most observed profiles exhibit a depletion in the inner region (at $r<0.3-0.4$ of the optical radius $r_{25}$ ) and are more or less exponential in the outer part. We also find breaks in the dust emission profiles at longer distances $(0.5-0.6)r_{25}$ which are associated with the breaks in the optical and near-infrared. We assume that the observed deficit of dust emission in the inner galaxy region is related to the depression in the radial profile of the HI surface density in the same region because the atomic gas reaches high enough surface densities there to be transformed into molecular gas. If a galaxy has a triggered star formation in the inner region (for example, because of a strong bar instability, which transfers the gas inwards to the centre, or a pseudobulge formation), no depletion or even an excess of dust emission in the centre is observed., Comment: 24 pages, 17 figures, Accepted for publication in Astronomy & Astrophysics

Our objective is to investigate the physical properties of the ionised gas of the low-metallicity dwarf galaxy, IC 10, at various spatial scales: from individual HII regions to the entire galaxy scale and examine whether diagnostics for integrated measurements introduce bias in the results. We modeled the ionised gas combining the mid- and far-infrared fine-structure cooling lines observed with Spitzer/IRS and Herschel/PACS, with the photoionisation code Cloudy. The free parameters of the models are the age of the stellar cluster, the density and the ionisation parameter of the ionised gas as well as the depth of the cloud. The latter is used to investigate the leakage of the ionising photons from the analysed regions of IC 10. We investigate HII regions in the main star-forming body, on scales of ~25 pc, three in the main star-forming region in the center of the galaxy and two on the first arc. We then consider larger sizes on the scale of ~200 pc. We find that most clumps have nearly identical properties, density ~10$^{2.} $ - 10$^{2.6}$ cm$^{-3}$, ionisation parameter between 10$^{-2.2}$ and 10$^{-1.6}$ and age of the stellar cluster ~5.5 Myr. All of them are matter-bounded regions, allowing ionising photons to leak. The relatively uniform physical properties of the clumps suggest a common origin for their star formation activity, which could be related to the feedback from stellar winds or supernovae of a previous generation of stars. The properties derived for ~200 pc size "zones" have similar properties as the HII regions they encompass, but with the larger regions tending to be more radiation-bounded. Finally, we investigate the fraction of [CII] 157.7 {\mu}m, [SiII] 34.8 {\mu}m, and [FeII] 25.9 {\mu}m, emission arising from the ionised gas phase and we find that most of the emission originates from the neutral gas, not from the ionised gas., Comment: accepted for publication in A&A

We study the fraction of stellar radiation absorbed by dust, f_abs, 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) galaxies observed with the Herschel Space Observatory. For each object, we model the spectral energy distribution from the ultraviolet to the sub-millimetre using the dedicated, aperture-matched DustPedia photometry and the fitting code CIGALE. The value of f_abs is obtained from the total luminosity emitted by dust and from the bolometric luminosity, which are estimated by the fit. 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 f_abs 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 f_abs 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 f_abs with stellar mass and specific star formation rate. We confirm that the local Universe is missing the high f_abs}, luminous and actively star-forming objects necessary to explain the energy budget in observations of the extragalactic background light., Comment: 21 pages, 14 figures, A&A accepted. Template SEDs available at the DustPedia Archive: http://dustpedia.astro.noa.gr/

We have fit the far-ultraviolet (FUV) to sub-millimeter (850 micron) 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 Multi-wavelength 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 particular between 25 and 70 micron 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 micron) 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 roughly 0.3 dex using only Mstar and SFR., Comment: 42 pages, 23 figures, 8 tables, 2 appendices: updated to v2, A&A, accepted

This article gives an overview of the constitution, physical conditions and observables of dust in the interstellar medium of nearby galaxies. We first review the macroscopic, spatial distribution of dust in these objects, and its consequence on our ability to study grain physics. We also discuss the possibility to use dust tracers as diagnostic tools. We then survey our current understanding of the microscopic, intrinsic properties of dust in different environments, derived from different observables: emission, extinction, polarization, depletions, over the whole electromagnetic spectrum. Finally, we summarize the clues of grain evolution, evidenced either on local scales or over cosmic time. We put in perspective the different evolution scenarios. We attempt a comprehensive presentation of the main observational constraints, analysis methods and modelling frameworks of the distinct processes. We do not cover the dust properties of the Milky Way and distant galaxies, nor circumstellar or active galactic nucleus torus dust., Comment: 41 pages, 9 figures, accepted for publication by Annual Review of Astronomy & Astrophysics. Version 2: only the acknowledgements have been edited. Version 3: the last sentence of the first paragraph in Sect. 3.3.2 (p. 21) has been edited. Version 4: Figure 1 has been corrected; the orientation of the polarization relative to the magnetic field was wrong; an erratum has been submitted to ARAA

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. 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. 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, Comment: Accepted for publication in Astronomy & Astrophysics

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. This study is performed by exploiting the multi-wavelength, from UV to sub-mm bands, DustPedia database, 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 (alpha_CO) per galaxy by using dust and gas mass profiles. Although each galaxy has its own peculiar behaviour, we identified a common trend of the exponential scale-lengths vs. wavelength. On average, the scale-lengths normalized to the B-band 25 mag/arcsec^2 radius decrease from UV to 70 micron, 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 micron 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 micron scale-length tend to increase from earlier to later types, the scale-length at 70 micron tends to be smaller than that at longer sub-mm wavelength with ratios between longer sub-mm wavelengths and 70 micron that decrease with increasing T. The scale-length ratio of SFR and stars shows a weak increasing trend towards later types., Comment: 33 pages, 9 figures (at low resolution here), 8 tables, Accepted for publication in A&A

We present an update to the ultraviolet-to-radio database of global broadband photometry for the 79 nearby galaxies that comprise the union of the KINGFISH (Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel) and SINGS (Spitzer Infrared Nearby Galaxies Survey) samples. The 34-band dataset presented here includes contributions from observational work carried out with a variety of facilities including GALEX, SDSS, PS, NOAO, 2MASS, WISE, Spitzer, Herschel, Planck, JCMT, and the VLA. Improvements of note include recalibrations of previously-published SINGS BVRcIc and KINGFISH far-infrared/submillimeter photometry. Similar to previous results in the literature, an excess of submillimeter emission above model predictions is seen primarily for low-metallicity dwarf/irregular galaxies. This 34-band photometric dataset for the combined KINGFISH$+$SINGS sample serves as an important multi-wavelength reference for the variety of galaxies observed at low redshift. A thorough analysis of the observed spectral energy distributions is carried out in a companion paper., Comment: Accepted for publication in the Astrophysical Journal; object, sky, and contaminant apertures are available from the first author

Dynamic and thermal processes regulate the structure of the multi-phase 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 ($P_{\rm th}$) of the diffuse, neutral gas. By combining [C II] 158 $\mu$m, HI, 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 $\sim$1 kiloparsec. We find a distribution of thermal pressures in the $P_{\rm th}/k\sim10^3-10^5$ 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 $P_{\rm th}/k\approx3600$ K cm$^{-3}$. These results are consistent with thermal pressure measurements using other observational methods. We find that $P_{\rm th}$ 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 medium could exist in pressure equilibrium. Finally, we find that the midplane thermal pressure of the diffuse gas is about $\sim30$% of the vertical weight of the overlying ISM, consistent with results from hydrodynamical simulations of self-regulated star formation in galactic disks., Comment: Accepted for publication in ApJ on Dec. 20th, 2016; 17 pages; 7 figures and 1 table

Extragalactic observations of water emission can provide valuable insights into the excitation of the interstellar medium. In addition, extragalactic megamasers are powerful probes of kinematics close to active nuclei. Therefore, it is paramount to determine the true origin of the water emission, whether it is excited by processes close to an AGN or in star-forming regions. We use ALMA Band 5 science verification observations to analyse the emission of the 183 GHz water line in Arp 220 on sub-arcsecond scales, in conjunction with new ALMA Band 7 data at 325 GHz. Specifically, the nature of the process leading to the excitation of emission at these water lines is studied in this context. Supplementary 22 GHz VLA observations are used to better constrain the parameter space in the excitation modelling of the water lines. We detect 183 GHz H2O and 325 GHz water emission towards the two compact nuclei at the center of Arp 220, being brighter in Arp 220 West. The emission at these two frequencies is compared to previous single-dish data and does not show evidence of variability. The 183 and 325 GHz lines show similar spectra and kinematics, but the 22 GHz profile is significantly different in both nuclei due to a blend with an NH3 absorption line. Our findings suggest that the most likely scenario to cause the observed water emission in Arp 220 is a large number of independent masers originating from numerous star-forming regions., Comment: 10 pages, 7 figures, 3 tables, accepted for publication in A&A

We study the spectral energy distribution (SED) of the radio continuum emission from the KINGFISH sample of nearby galaxies to understand the energetics and origin of this emission. Effelsberg multi-wavelength observations at 1.4GHz, 4.8GHz, 8.5GHz, and 10.5GHz combined with archive data allow us, for the first time, to determine the mid-radio continuum (1-10 GHz, MRC) bolometric luminosities and further present calibration relations vs. the monochromatic radio luminosities. The 1-10 GHz radio SED is fitted using a Bayesian Markov Chain Monte Carlo (MCMC) technique leading to measurements for the nonthermal spectral index and the thermal fraction f_th with mean values of alpha_nt=0.97+-0.16 (0.79+-0.15 for the total spectral index) and f_th= 10% +- 9% at 1.4 GHz. The MRC luminosity changes over ~3 orders of magnitude in the sample. The thermal emission is responsible for ~23% of the MRC on average. We also compare the extinction-corrected diagnostics of star formation rate with the thermal and nonthermal radio tracers and derive the first star formation calibration relations using the MRC radio luminosity. The nonthermal spectral index flattens with increasing star formation rate surface density, indicating the effect of the star formation feedback on the cosmic ray electron population in galaxies. Comparing the radio and IR SEDs, we find that the FIR-to-MRC ratio could decrease with star formation rate, due to the amplification of the magnetic fields in star forming regions. This particularly implies a decrease in the ratio at high redshifts, where mostly luminous/star forming galaxies are detected., Comment: 29 pages, 15 figures, accepted for publication in Astrophysical Journal

The radiation of stars heats dust grains in the diffuse interstellar medium and in star-forming regions in galaxies. Modelling this interaction provides information on dust in galaxies, a vital ingredient for their evolution. It is not straightforward to identify the stellar populations heating the dust, and to link attenuation to emission on a sub-galactic scale. Radiative transfer models are able to simulate this dust-starlight interaction in a realistic, three-dimensional setting. We investigate the dust heating mechanisms on a local and global galactic scale, using the Andromeda galaxy (M31) as our laboratory. We perform a series of panchromatic radiative transfer simulations of Andromeda with our code SKIRT. The high inclination angle of M31 complicates the 3D modelling and causes projection effects. However, the observed morphology and flux density are reproduced fairly well from UV to sub-millimeter wavelengths. Our model reveals a realistic attenuation curve, compatible with previous, observational estimates. We find that the dust in M31 is mainly (91 % of the absorbed luminosity) heated by the evolved stellar populations. The bright bulge produces a strong radiation field and induces non-local heating up to the main star-forming ring at 10 kpc. The relative contribution of unevolved stellar populations to the dust heating varies strongly with wavelength and with galactocentric distance.The dust heating fraction of unevolved stellar populations correlates strongly with NUV-r colour and specific star formation rate. These two related parameters are promising probes for the dust heating sources at a local scale., Comment: 16 pages + appendix, 13 figures, accepted for publication in A&A

The European Space Agency has invested heavily in two cornerstones missions; Herschel and Planck. The legacy data from these missions provides us with an unprecedented opportunity to study cosmic dust in galaxies so that we can answer fundamental questions about, for example: the origin of the chemical elements, physical processes in the interstellar medium (ISM), its effect on stellar radiation, its relation to star formation and how this relates to the cosmic far infrared background. In this paper we describe the DustPedia project, which is enabling us to develop tools and computer models that will help us relate observed cosmic dust emission to its physical properties (chemical composition, size distribution, temperature), to its origins (evolved stars, super novae, growth in the ISM) and the processes that destroy it (high energy collisions and shock heated gas). To carry out this research we will combine the Herschel/Planck data with that from other sources of data, providing observations at numerous wavelengths (< 41) across the spectral energy distribution, thus creating the DustPedia database. To maximise our spatial resolution and sensitivity to cosmic dust we limit our analysis to 4231 local galaxies (v < 3000 km/s) selected via their near infrared luminosity (stellar mass). To help us interpret the data we have developed a new physical model for dust (THEMIS), a new Bayesian method of fitting and interpreting spectral energy distributions (HerBIE) and a state-of-the-art Monte Carlo photon tracing radiative transfer model (SKIRT). In this the first of the DustPedia papers we describe the project objectives, data sets used and provide an insight into the new scientific methods we plan to implement., Comment: 66 pages, accepted for publication in PASP

We present a line survey of the ultra-luminous infrared galaxy Arp 220, taken with the newly installed SEPIA Band 5 instrument on APEX. We illustrate the capacity of SEPIA to detect the 183.3 GHz H2O 31,3-22,0 line against the atmospheric H2O absorption feature. We confirm the previous detection of the HCN(2-1) line, and detect new transitions of standard dense gas tracers such as HNC(2-1), HCO+(2-1), CS(4-3), C34S(4-3), HC3N(20-19). We also detect HCN(2-1) v2=1 and the 193.5 GHz methanol (4-3) group for the first time. The absence of time variations in the megamaser water line compared to previous observations seems to rule out an AGN nuclear origin for the line. It could, on the contrary, favor a thermal origin instead, but also possibly be a sign that the megamaser emission is associated with star-forming cores washed-out in the beam. We finally discuss how the new transitions of HCN, HNC, HCO+ refine our knowledge of the ISM physical conditions in Arp 220., Comment: 5 pages, 2 figures, 1 table, accepted for publication in MNRAS

The [NII] 122 and 205 \mu 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 kiloparsec, with observations of both [NII] far-infrared lines. We measure [NII] 122/205 line ratios in the ~0.6-6 range, which corresponds to electron gas densities $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 a 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 surface density ($\Sigma_{SFR}$). As the intensity of the [NII] transitions is related to the ionizing photon flux, we investigate their reliability as tracers of the star formation rate (SFR). We derive relations between the [NII] 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 [NII] surface brightness and $\Sigma_{SFR}$ can be understood as a property of the $n_e$ distribution. For regions with $n_e$ close to or higher than the [NII] line critical densities, the low-density limit [NII]-based SFR calibration systematically underestimates the SFR since [NII] emission is collisionally quenched. Finally, we investigate the relation between [NII] emission, SFR, and $n_e$ by comparing our observations to predictions from the MAPPINGS-III code., Comment: 18 pages, 9 figures, accepted for publication in The Astrophysical Journal

To compute the SFR of galaxies from the rest-frame UV it is essential to take into account the obscuration by dust. To do so, one of the most popular methods consists in combining the UV with the emission from the dust itself in the IR. Yet, different studies have derived different estimators, showing that no such hybrid estimator is truly universal. In this paper we aim at understanding and quantifying what physical processes drive the variations between different hybrid estimators. Doing so, we aim at deriving new universal UV+IR hybrid estimators to correct the UV for dust attenuation, taking into account the intrinsic physical properties of galaxies. We use the CIGALE code to model the spatially-resolved FUV to FIR SED 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. We find that hybrid UV+IR estimators strongly depend on the stellar mass surface density (in particular at 70 and 100 micron) and on the sSFR (in particular at 24 micron and the TIR). Consequently, the IR scaling coefficients for UV obscuration can vary by almost an order of magnitude. This result contrasts with other groups who found relatively constant coefficients with small deviations. We exploit these variations to construct a new class of 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. The new estimators provide better estimates of attenuation-corrected UV emission than classical hybrid estimators. Naturally taking into account the variable impact of dust heated by old stellar populations, they constitute a step towards universal estimators., Comment: 25 pages, 12 figures, 8 tables. Accepted for publication in A&A

We combine Spitzer and Herschel data of the star-forming region N11 in the Large Magellanic Cloud to produce detailed maps of the dust properties in the complex and study their variations with the ISM conditions. We also compare APEX/LABOCA 870um observations with our model predictions in order to decompose the 870um emission into dust and non-dust (free-free emission and CO(3-2) line) contributions. We find that in N11, the 870um can be fully accounted for by these 3 components. The dust surface density map of N11 is combined with HI and CO observations to study local variations in the gas-to-dust mass ratios. Our analysis leads to values lower than those expected from the LMC low-metallicity as well as to a decrease of the gas-to-dust mass ratio with the dust surface density. We explore potential hypotheses that could explain the low observed gas-to-dust mass ratios (variations in the XCO factor, presence of CO-dark gas or of optically thick HI or variations in the dust abundance in the dense regions). We finally decompose the local SEDs using a Principal Component Analysis (i.e. with no a priori assumption on the dust composition in the complex). Our results lead to a promising decomposition of the local SEDs in various dust components (hot, warm, cold) coherent with that expected for the region. Further analysis on a larger sample of galaxies will follow in order to understand how unique this decomposition is or how it evolves from one environment to another., Comment: 24 pages, 16 figures, accepted for publication in MNRAS

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&A

The far-infrared (FIR) lines are key tracers of the physical conditions of the interstellar medium (ISM) and are becoming workhorse diagnostics for galaxies throughout the universe. Our goal is to explain the differences and trends observed in the FIR line emission of dwarf galaxies compared to more metal-rich galaxies. We present Herschel PACS spectroscopic observations of the CII157um, OI63 and 145um, OIII88um, NII122 and 205um, and NIII57um 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. 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 OIII/NII122 and NIII/NII122 are high, indicative of hard radiation fields. In the photodissociation region (PDR), the CII/OI63 ratio is slightly higher than in metal-rich galaxies, with a small increase with metallicity, and the OI145/OI63 ratio is generally lower than 0.1, demonstrating that optical depth effects should be small on the scales probed. The OIII/OI63 ratio can be used as an indicator of the ionized gas/PDR filling factor, and is found ~4 times higher in the dwarfs than in metal-rich galaxies. The high CII/LTIR, OI/LTIR, and OIII/LTIR ratios, which decrease with increasing LTIR and LTIR/LB, are interpreted as a combination of moderate FUV fields and low PDR covering factor. Harboring compact phases of 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 in metal-rich galaxies., Comment: Accepted for publication in A&A

The manner in which gas accretes and orbits within circumnuclear rings has direct implications for the star formation process. In particular, gas may be compressed and shocked at the inflow points, resulting in bursts of star formation at these locations. Afterwards the gas and young stars move together through the ring. In addition, star formation may occur throughout the ring, if and when the gas reaches sufficient density to collapse under gravity. These two scenarios for star formation in rings are often referred to as the `pearls on a string' and `popcorn' paradigms. In this paper, we use new Herschel PACS observations, obtained as part of the KINGFISH Open Time Key Program, along with archival Spitzer and ground-based observations from the SINGS Legacy project, to investigate the heating and cooling of the interstellar medium in the nearby star-forming ring galaxy, NGC4736. By comparing spatially resolved estimates of the stellar FUV flux available for heating, with the gas and dust cooling derived from the FIR continuum and line emission, we show that while star formation is indeed dominant at the inflow points in NGC 4736, additional star formation is needed to balance the gas heating and cooling throughout the ring. This additional component most likely arises from the general increase in gas density in the ring over its lifetime. Our data provide strong evidence, therefore, for a combination of the two paradigms for star formation in the ring in NGC4736., Comment: accepted for publication in A&A

The rate at which interstellar gas is converted into stars, and its dependence on environment, is one of the pillars on which our understanding of the visible Universe is build. We present a comparison of the surface density of young stars (Sigma_*) and dust surface density (Sigma_d) across NGC346 (N66) in 115 independent pixels of 6x6 pc^2. We find a correlation between Sigma_* and Sigma_d with a considerable scatter. A power law fit to the data yields a steep relation with an exponent of 2.6+-0.2. We convert Sigma_d to gas surface density (Sigma_g) and Sigma_* to star formation rate (SFR) surface densities (Sigma_SFR), using simple assumptions for the gas-to-dust mass ratio and the duration of star formation. The derived total SFR (4+-1 10^-3 M_sun/yr) is consistent with SFR estimated from the Ha emission integrated over the Ha nebula. On small scales the Sigma_SFR derived using Ha systematically underestimates the count-based Sigma_SFR, by up to a factor of 10. This is due to ionizing photons escaping the area, where the stars are counted. We find that individual 36 pc^2 pixels fall systematically above integrated disc-galaxies in the Schmidt-Kennicutt diagram by on average a factor of ~7. The NGC346 average SFR over a larger area (90 pc radius) lies closer to the relation but remains high by a factor of ~3. The fraction of the total mass (gas plus young stars) locked in young stars is systematically high (~10 per cent) within the central 15 pc and systematically lower outside (2 per cent), which we interpret as variations in star formation efficiency. The inner 15 pc is dominated by young stars belonging to a centrally condensed cluster, while the outer parts are dominated by a dispersed population. Therefore, the observed trend could reflect a change of star formation efficiency between clustered and non-clustered star-formation., Comment: Accepted for publication in MNRAS

With the newly available SPIRE images at 250 and 500 micron from Herschel Space Observatory, we study quantitative correlations over a sub-kpc scale among three distinct emission components in the interstellar medium of the nearby spiral galaxy M 81 (NGC 3031): (a) $I_{8}$ or $I_{24}$, the surface brightness of the mid-infrared emission observed in the Spitzer IRAC 8 or MIPS 24 micron band, with $I_8$ and $I_{24}$ being dominated by the emissions from Polycyclic Aromatic Hydrocarbons (PAHs) and very small grains (VSGs) of dust, respectively; (b) $I_{500}$, that of the cold dust continuum emission in the Herschel SPIRE 500 micron band, dominated by the emission from large dust grains heated by evolved stars, and (c) $I_{{\rm H}\alpha}$, a nominal surface brightness of the H$\alpha$ 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 micron (24 micron) emission derives its heating from evolved stars, with the remainder attributed to radiation heating associated with ionizing stars., Comment: 11 pages, 7 figures; accepted for publication in the Astrophysical Journal

The [CII] 157.74 $\mu$m transition is the dominant coolant of the neutral interstellar gas, and has great potential as a star formation rate (SFR) tracer. Using the Herschel KINGFISH sample of 46 nearby galaxies, we investigate the relation of [CII] surface brightness and luminosity with SFR. We conclude that [CII] can be used for measurements of SFR on both global and kiloparsec scales in normal star-forming galaxies in the absence of strong active galactic nuclei (AGN). The uncertainty of the $\Sigma_{\rm [CII]}-\Sigma_{\rm SFR}$ calibration is $\pm$0.21 dex. The main source of scatter in the correlation is associated with regions that exhibit warm IR colors, and we provide an adjustment based on IR color that reduces the scatter. We show that the color-adjusted $\Sigma_{\rm[CII]}-\Sigma_{\rm SFR}$ correlation is valid over almost 5 orders of magnitude in $\Sigma_{\rm SFR}$, holding for both normal star-forming galaxies and non-AGN luminous infrared galaxies. Using [CII] luminosity instead of surface brightness to estimate SFR suffers from worse systematics, frequently underpredicting SFR in luminous infrared galaxies even after IR color adjustment (although this depends on the SFR measure employed). We suspect that surface brightness relations are better behaved than the luminosity relations because the former are more closely related to the local far-UV field strength, most likely the main parameter controlling the efficiency of the conversion of far-UV radiation into gas heating. A simple model based on Starburst99 population-synthesis code to connect SFR to [CII] finds that heating efficiencies are $1\%-3\%$ in normal galaxies., Comment: 25 pages, 14 figures, accepted for publication in ApJ

Physical conditions of the interstellar medium in galaxies are closely linked to the ambient radiation field and the heating of dust grains. In order to characterize dust properties in galaxies over a wide range of physical conditions, we present here the radial surface brightness profiles of the entire sample of 61 galaxies from Key Insights into Nearby Galaxies: Far-Infrared Survey with Herschel (KINGFISH). The main goal of our work is the characterization of the grain emissivities, dust temperatures, and interstellar radiation fields responsible for heating the dust. After fitting the dust and stellar radial profiles with exponential functions, we fit the far-infrared spectral energy distribution (SED) in each annular region with single-temperature modified black bodies using both variable (MBBV) and fixed (MBBF) emissivity indices beta, as well as with physically motivated dust models. Results show that while most SED parameters decrease with radius, the emissivity index beta also decreases with radius in some galaxies, but in others is increasing, or rising in the inner regions and falling in the outer ones. Despite the fixed grain emissivity (average beta~ 2.1) of the physically-motivated models, they are well able to accommodate flat spectral slopes with beta<= 1. We find that flatter slopes (beta<= 1.5) are associated with cooler temperatures, contrary to what would be expected from the usual Tdust-beta degeneracy. This trend is related to variations in Umin since beta and Umin are very closely linked over the entire range in Umin sampled by the KINGFISH galaxies: low Umin is associated with flat beta<=1. Both these results strongly suggest that the low apparent \beta values (flat slopes) in MBBV fits are caused by temperature mixing along the line-of-sight, rather than by intrinsic variations in grain properties. Abstract truncated for arXiv., Comment: 28 pages, 20 figures, accepted for publication in A&A

We search for variations in the disk of Centaurus A of the emission from atomic fine structure lines using Herschel PACS and SPIRE spectroscopy. In particular we observe the [C II](158 $\mu$m), [N II](122 and 205 $\mu$m), [O I](63 and 145 $\mu$m) and [O III](88 $\mu$m) lines, which all play an important role in cooling the gas in photo-ionized and photodissociation regions. We determine that the ([C II]+[O I]$_{63}$)/$F_{TIR}$ line ratio, a proxy for the heating efficiency of the gas, shows no significant radial trend across the observed region, in contrast to observations of other nearby galaxies. We determine that 10 - 20% of the observed [C II] emission originates in ionized gas. Comparison between our observations and a PDR model shows that the strength of the far-ultraviolet radiation field, $G_0$, varies between $10^{1.75}$ and $10^{2.75}$ and the hydrogen nucleus density varies between $10^{2.75}$ and $10^{3.75}$ cm$^{-3}$, with no significant radial trend in either property. In the context of the emission line properties of the grand-design spiral galaxy M51 and the elliptical galaxy NGC 4125, the gas in Cen A appears more characteristic of that in typical disk galaxies rather than elliptical galaxies., Comment: Accepted for publication in the Astrophysical Journal. 22 pages, 10 figures, 5 tables

We investigate the connection between dust and gas in the nearby edge-on spiral galaxy NGC 891. High resolution Herschel PACS and SPIRE 70, 100, 160, 250, 350, and 500 $\mu$m images are combined with JCMT SCUBA 850 $\mu$m observations to trace the far-infrared/submillimetre spectral energy distribution (SED). Maps of the HI 21 cm line and CO(J=3-2) emission trace the atomic and molecular hydrogen gas, respectively. We fit one-component modified blackbody models to the integrated SED, finding a global dust mass of 8.5$\times$10$^{7}$ M$_{\odot}$ and an average temperature of 23$\pm$2 K. We also fit the pixel-by-pixel SEDs to produce maps of the dust mass and temperature. The dust mass distribution correlates with the total stellar population as traced by the 3.6 $\mu$m emission. The derived dust temperature, which ranges from approximately 17 to 24 K, is found to correlate with the 24 $\mu$m emission. Allowing the dust emissivity index to vary, we find an average value of $\beta$ = 1.9$\pm$0.3. We confirm an inverse relation between the dust emissivity spectral index and dust temperature, but do not observe any variation of this relationship with vertical height from the mid-plane of the disk. A comparison of the dust properties with the gaseous components of the ISM reveals strong spatial correlations between the surface mass densities of dust and the molecular hydrogen and total gas surface densities. Observed asymmetries in the dust temperature, and the H$_{2}$-to-dust and total gas-to-dust ratios hint that an enhancement in the star formation rate may be the result of larger quantities of molecular gas available to fuel star formation in the NE compared to the SW. Whilst the asymmetry likely arises from dust obscuration due to the geometry of the line-of-sight projection of the spiral arms, we cannot exclude an enhancement in the star formation rate in the NE side of the disk., Comment: Accepted for publication in A&A. 21 pages, including 13 figures and 4 tables

We gather infrared (IR) photometric data from 8 to 500 microns (Spitzer, WISE, IRAS and Herschel) for all of the HRS galaxies. Draine & Li (2007) 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 models parameters and that a strong constraint on the 20-60 microns range is mandatory to estimate the relative contribution of the photo-dissociation regions to the IR SED. The SED models tend to systematically under-estimate the observed 500 microns flux densities, especially for low mass systems. We provide the output parameters for all of the galaxies: the minimum intensity of the interstellar radiation field (ISRF), the fraction of PAH, the relative contribution of PDR and evolved stellar population to the dust heating, the $M_{dust}$ and the $L_{IR}$. For a subsample of gas-rich galaxies, we analyze the relations between these parameters and the integrated properties of galaxies, such as $M_*$, SFR, metallicity, H$\alpha$ and H-band surface brightness, and the FUV 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. The intensity of the IRSF and the H-band and H$\alpha$ 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 IR 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 put our sample on the SFR-$M_*$ diagram. The templates are compared to the most popular IR SED libraries, enlightening a large discrepancy between all of them in the 20-100 microns range., Comment: 37 pages, 17 figures, 11 tables, accepted for publication in A&A

We model the infrared to submillimeter spectral energy distribution of 11 nearby galaxies of the KINGFISH sample using Spitzer and Herschel data and compare model extrapolations at 870um (using different fitting techniques) with LABOCA 870um observations. We investigate how the differences between predictions and observations vary with model assumptions or environment. At global scales, we find that modified blackbody models using realistic cold emissivity indices (beta_c=2 or 1.5) are able to reproduce the 870um observed emission within the uncertainties for most of the sample. Low values (beta_c<1.3) would be required in NGC0337, NGC1512 and NGC7793. At local scales, we observe a systematic 870um excess when using beta_=2.0. The beta_c=1.5 or the Draine and Li (2007) models can reconcile predictions with observations in part of the disks. Some of the remaining excesses occur towards the centres and can be partly or fully accounted for by non-dust contributions such as CO(3-2) or, to a lesser extent, free-free or synchrotron emission. In three non-barred galaxies, the remaining excesses rather occur in the disk outskirts. This could be a sign of a flattening of the submm slope (and decrease of the effective emissivity index) with radius in these objects., Comment: 31 pages (including appendix), 7 figures, accepted for publication in MNRAS