Search

The ionising radiation of young and massive stars is a crucial form of stellar feedback. Most ionising (Lyman-continuum; LyC, $\lambda < 912A$) photons are absorbed close to the stars that produce them, forming compact HII regions, but some escape into the wider galaxy. Quantifying the fraction of LyC photons that escape is an open problem. In this work, we present a semi-novel method to estimate the escape fraction by combining broadband photometry of star clusters from the Legacy ExtraGalactic UV Survey (LEGUS) with HII regions observed by the Star formation, Ionized gas, and Nebular Abundances Legacy Survey (SIGNALS) in the nearby spiral galaxy NGC 628. We first assess the completeness of the combined catalogue, and find that 49\% of HII regions lack corresponding star clusters as a result of a difference in the sensitivities of the LEGUS and SIGNALS surveys. For HII regions that do have matching clusters, we infer the escape fraction from the difference between the ionising power required to produce the observed HII luminosity and the predicted ionising photon output of their host star clusters; the latter is computed using a combination of LEGUS photometric observations and a stochastic stellar population synthesis code SLUG (Stochastically Lighting Up Galaxies). Overall, we find an escape fraction of $f_{esc} = 0.09^{+0.06}_{-0.06}$ across our sample of 42 HII regions; in particular, we find HII regions with high $f_{esc}$ are predominantly regions with low H$\alpha$-luminosity. We also report possible correlation between $f_{esc}$ and the emission lines [O ii]/[N ii] and [O ii]/H$\beta$., Comment: Accepted for publication at MNRAS

We study the young star cluster populations in 23 dwarf and irregular galaxies observed by the HST Legacy ExtraGalactic Ultraviolet Survey (LEGUS), and examine relationships between the ensemble properties of the cluster populations and those of their host galaxies: star formation rate (SFR) density ($\Sigma_{\rm SFR}$). A strength of this analysis is the availability of SFRs measured from temporally resolved star formation histories which provide the means to match cluster and host-galaxy properties on several timescales (1-10, 1-100, and 10-100~Myr). Nevertheless, studies of this kind are challenging for dwarf galaxies due to the small numbers of clusters in each system. We mitigate these issues by combining the clusters across different galaxies with similar $\Sigma_{\rm SFR}$ properties. We find good agreement with a well-established relationship ($M_{V}^{brightest}$-SFR), but find no significant correlations between $\Sigma_{\rm SFR}$ and the slopes of the cluster luminosity function, mass function, nor the age distribution. We also find no significant trend between the the fraction of stars in bound clusters at different age ranges ($\Gamma_{1-10}$, $\Gamma_{10-100}$, and $\Gamma_{1-100}$) and $\Sigma_{\rm SFR}$ of the host galaxy. Our data show a decrease in $\Gamma$ over time (from 1-10 to 10-100~Myr) suggesting early cluster dissolution, though the presence of unbound clusters in the youngest time bin makes it difficult to quantify the degree of dissolution. While our data do not exhibit strong correlations between $\Sigma_{\rm SFR}$ and ensemble cluster properties, we cannot rule out that a weak trend might exist given the relatively large uncertainties due to low number statistics and the limited $\Sigma_{\rm SFR}$ range probed., Comment: 30 pages, 15 figures, 7 tables, Accepted for publication in MNRAS

We present the second data release of high-resolution ($\leq0.2$ arcsec) 1.5-GHz radio images of 177 nearby galaxies from the Palomar sample, observed with the e-MERLIN array, as part of the LeMMINGs (Legacy e-MERLIN Multi-band Imaging of Nearby Galaxy Sample) survey. Together with the 103 targets of the first LeMMINGs data release, this represents a complete sample of 280 local active (LINER and Seyfert) and inactive galaxies HII galaxies and Absorption Line Galaxies, ALG). This large program is the deepest radio survey of the local Universe, $\gtrsim$10$^{17.6}$ W Hz$^{-1}$, regardless of the host and nuclear type: we detect radio emission $\gtrsim$0.25 mJy beam$^{-1}$ for 125/280 galaxies (44.6 per cent) with sizes of typically $\lesssim$100 pc. Of those 125, 106 targets show a core which coincides within 1.2 arcsec with the optical nucleus. Although we observed mostly cores, around one third of the detected galaxies features jetted morphologies. The detected radio core luminosities of the sample range between $\sim$10$^{34}$ and 10$^{40}$ erg s$^{-1}$. LINERs and Seyferts are the most luminous sources, whereas HII galaxies are the least. LINERs show FRI-like core-brightened radio structures, while Seyferts reveal the highest fraction of symmetric morphologies. The majority of HII galaxies have single radio core or complex extended structures, which probably conceal a nuclear starburst and/or a weak active nucleus (seven of them show clear jets). ALGs, which are typically found in evolved ellipticals, although the least numerous, exhibit on average the most luminous radio structures, similar to LINERs., Comment: Accepted for publication on MNRAS (48 pages, 7 figures, 10 tables, Appendix)

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

Scaling laws of dust, HI gas and metal mass with stellar mass, specific star formation rate and metallicity are crucial to our understanding of the buildup of galaxies through their enrichment with metals and dust. In this work, we analyse how the dust and metal content varies with specific gas mass ($M_{\text{HI}}$/$M_{\star}$) across a diverse sample of 423 nearby galaxies. The observed trends are interpreted with a set of Dust and Element evolUtion modelS (DEUS) - incluidng stellar dust production, grain growth, and dust destruction - within a Bayesian framework to enable a rigorous search of the multi-dimensional parameter space. We find that these scaling laws for galaxies with $-1.0\lesssim \log M_{\text{HI}}$/$M_{\star}\lesssim0$ can be reproduced using closed-box models with high fractions (37-89$\%$) of supernova dust surviving a reverse shock, relatively low grain growth efficiencies ($\epsilon$=30-40), and long dus lifetimes (1-2\,Gyr). The models have present-day dust masses with similar contributions from stellar sources (50-80\,$\%$) and grain growth (20-50\,$\%$). Over the entire lifetime of these galaxies, the contribution from stardust ($>$90\,$\%$) outweighs the fraction of dust grown in the interstellar medium ($<$10$\%$). Our results provide an alternative for the chemical evolution models that require extremely low supernova dust production efficiencies and short grain growth timescales to reproduce local scaling laws, and could help solving the conundrum on whether or not grains can grow efficiently in the interstellar medium., Comment: 56 pages, 30 figures, accepted for publication in MNRAS

Dust is formed out of stellar material and is constantly affected by different mechanisms occurring in the ISM. Dust grains behave differently under these mechanisms depending on their sizes, and therefore the dust grain size distribution also evolves as part of the dust evolution itself. Following how the grain size distribution evolves is a difficult computing task that is just recently being overtaking. Smoothed particle hydrodynamic (SPH) simulations of a single galaxy as well as cosmological simulations are producing the first predictions of the evolution of the dust grain size distribution. We compare for the first time the evolution of the dust grain size distribution predicted by the SPH simulations with the results provided by the observations. We analyse how the radial distribution of the small to large grain mass ratio (D(S)/D(L)) changes over the whole discs in three galaxies: M 101, NGC 628 and M 33. We find good agreement between the observed radial distribution of D(S)/D(L) and what is obtained from the SPH simulations of a single galaxy. The central parts of NGC 628, at high metallicity and with a high molecular gas fraction, are mainly affected not only by accretion but also by coagulation of dust grains. The centre of M 33, having lower metallicity and lower molecular gas fraction, presents an increase of D(S)/D(L), showing that shattering is very effective in creating a large fraction of small grains. Observational results provided by our galaxies confirm the general relations predicted by the cosmological simulations based on the two grain size approximation. However, we present evidence that the simulations could be overestimating the amount of large grains in high massive galaxies., Comment: 11 pages, 6 Figures, to be published in A&A

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

The morphology of HII regions around young star clusters provides insight into the timescales and physical processes that clear a cluster's natal gas. We study ~700 young clusters (<10Myr) in three nearby spiral galaxies (NGC 7793, NGC 4395, and NGC 1313) using Hubble Space Telescope (HST) imaging from LEGUS (Legacy ExtraGalactic Ultraviolet Survey). Clusters are classified by their H$\alpha$ morphology (concentrated, partially exposed, no-emission) and whether they have neighboring clusters (which could affect the clearing timescales). Through visual inspection of the HST images, and analysis of ages, reddenings, and stellar masses from spectral energy distributions fitting, together with the (U-B), (V-I) colors, we find: 1) the median ages indicate a progression from concentrated (~3 Myr), to partially exposed (~4 Myr), to no H$\alpha$ emission (>5Myr), consistent with the expected temporal evolution of HII regions and previous results. However, 2) similarities in the age distributions for clusters with concentrated and partially exposed H$\alpha$ morphologies imply a short timescale for gas clearing (<1Myr). 3) our cluster sample's median mass is ~1000 M, and a significant fraction (~20%) contain one or more bright red sources (presumably supergiants), which can mimic reddening effects. Finally, 4) the median E(B-V) values for clusters with concentrated H$\alpha$ and those without H$\alpha$ emission appear to be more similar than expected (~0.18 vs. ~0.14, respectively), but when accounting for stochastic effects, clusters without H$\alpha$ emission are less reddened. To mitigate stochastic effects, we experiment with synthesizing more massive clusters by stacking fluxes of clusters within each H$\alpha$ morphological class. Composite isolated clusters also reveal a color and age progression for H$\alpha$ morphological classes, consistent with analysis of the individual clusters., Comment: 20 pages, 16 figures, 7 tables; accepted for publication in MNRAS

SIGNALS, the Star formation, Ionized Gas, and Nebular Abundances Legacy Survey, is a large observing program designed to investigate massive star formation and HII regions in a sample of local extended galaxies. The program will use the imaging Fourier transform spectrograph SITELLE at the Canada-France-Hawaii Telescope. Over 355 hours (54.7 nights) have been allocated beginning in fall 2018 for eight consecutive semesters. Once completed, SIGNALS will provide a statistically reliable laboratory to investigate massive star formation, including over 50 000 resolved HII regions : the largest, most complete, and homogeneous database of spectroscopically and spatially resolved extragalactic HII regions ever assembled. For each field observed, three datacubes covering the spectral bands of the filters SN1 (363 -386 nm), SN2 (482 - 513 nm), and SN3 (647 - 685 nm) are gathered. The spectral resolution selected for each spectral band is 1000, 1000, and 5000, respectively. As defined, the project sample will facilitate the study of small-scale nebular physics and many other phenomena linked to star formation at a mean spatial resolution of 20 pc. This survey also has considerable legacy value for additional topics including planetary nebulae, diffuse ionized gas, andsupernova remnants. The purpose of this paper is to present a general outlook of the survey, notably the observing strategy, galaxy sample, and science requirements., Comment: 19 pages, 14 figures, submitted to MNRAS

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations., Comment: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

We present the star cluster catalogs for 17 dwarf and irregular galaxies in the $HST$ Treasury Program "Legacy ExtraGalactic UV Survey" (LEGUS). Cluster identification and photometry in this subsample are similar to that of the entire LEGUS sample, but special methods were developed to provide robust catalogs with accurate fluxes due to low cluster statistics. The colors and ages are largely consistent for two widely used aperture corrections, but a significant fraction of the clusters are more compact than the average training cluster. However, the ensemble luminosity, mass, and age distributions are consistent suggesting that the systematics between the two methods are less than the random errors. When compared with the clusters from previous dwarf galaxy samples, we find that the LEGUS catalogs are more complete and provide more accurate total fluxes. Combining all clusters into a composite dwarf galaxy, we find that the luminosity and mass functions can be described by a power law with the canonical index of $-2$ independent of age and global SFR binning. The age distribution declines as a power law, with an index of $\approx-0.80\pm0.15$, independent of cluster mass and global SFR binning. This decline of clusters is dominated by cluster disruption since the combined star formation histories and integrated-light SFRs are both approximately constant over the last few hundred Myr. Finally, we find little evidence for an upper-mass cutoff ($<2\sigma$) in the composite cluster mass function, and can rule out a truncation mass below $\approx10^{4.5}$M$_{\odot}$ but cannot rule out the existence of a truncation at higher masses., Comment: 25 pages; 28 figures (Accepted to MNRAS)

We present a study correlating the spatial locations of young star clusters with those of molecular clouds in NGC~5194, in order to investigate the timescale over which clusters separate from their birth clouds. The star cluster catalogues are from the Legacy ExtraGalactic UV Survey (LEGUS) and the molecular clouds from the Plateau de Bure Interefrometer Arcsecond Whirpool Survey (PAWS). We find that younger star clusters are spatially closer to molecular clouds than older star clusters. The median ages for clusters associated with clouds is 4~Myr whereas it is 50~Myr for clusters that are sufficiently separated from a molecular cloud to be considered unassociated. After $\sim$6~Myr, the majority of the star clusters lose association with their molecular gas. Younger star clusters are also preferentially located in stellar spiral arms where they are hierarchically distributed in kpc-size regions for 50-100~Myr before dispersing. The youngest star clusters are more strongly clustered, yielding a two-point correlation function with $\alpha=-0.28\pm0.04$, than the GMCs ($\alpha=-0.09\pm0.03$) within the same PAWS field. However, the clustering strength of the most massive GMCs, supposedly the progenitors of the young clusters for a star formation efficiency of a few percent, is comparable ($\alpha=-0.35\pm0.05$) to that of the clusters. We find a galactocentric-dependence for the coherence of star formation, in which clusters located in the inner region of the galaxy reside in smaller star-forming complexes and display more homogeneous distributions than clusters further from the centre. This result suggests a correlation between the survival of a cluster complex and its environment., Comment: 17 pages, 12 figures. Accepted for publication in MNRAS

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

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

We present an investigation of the relationship between giant molecular cloud (GMC) properties and the associated stellar clusters in the nearby flocculent galaxy NGC 7793. We combine the star cluster catalog from the HST LEGUS (Legacy ExtraGalactic UV Survey) program with the 15 parsec resolution ALMA CO(2-1) observations. We find a strong spatial correlation between young star clusters and GMCs such that all clusters still associated with a GMC are younger than 11 Myr and display a median age of 2 Myr. The age distribution increases gradually as the cluster-GMC distance increases, with star clusters that are spatially unassociated with molecular gas exhibiting a median age of 7 Myr. Thus, star clusters are able to emerge from their natal clouds long before the timescale required for clouds to disperse. To investigate if the hierarchy observed in the stellar components is inherited from the GMCs, we quantify the amount of clustering in the spatial distributions of the components and find that the star clusters have a fractal dimension slope of $-0.35 \pm 0.03$, significantly more clustered than the molecular cloud hierarchy with slope of $-0.18 \pm 0.04$ over the range 40-800 pc. We find, however, that the spatial clustering becomes comparable in strength for GMCs and star clusters with slopes of $-0.44\pm0.03$ and $-0.45\pm0.06$ respectively, when we compare massive ($>$10$^5$ M$_{\odot}$) GMCs to massive and young star clusters. This shows that massive star clusters trace the same hierarchy as their parent GMCs, under the assumption that the star formation efficiency is a few percent., Comment: 12 page, 8 figures. Accepted for publication in MNRAS

We present a study of the dust-to-gas ratios in five nearby galaxies NGC 628 (M74), NGC 6503, NGC 7793, UGC 5139 (Holmberg I), and UGC 4305 (Holmberg II). Using Hubble Space Telescope broad band WFC3/UVIS UV and optical images from the Treasury program LEGUS (Legacy ExtraGalactic UV Survey) combined with archival HST/ACS data, we correct thousands of individual stars for extinction across these five galaxies using an isochrone-matching (reddening-free Q) method. We generate extinction maps for each galaxy from the individual stellar extinctions using both adaptive and fixed resolution techniques, and correlate these maps with neutral HI and CO gas maps from literature, including The HI Nearby Galaxy Survey (THINGS) and the HERA CO-Line Extragalactic Survey (HERACLES). We calculate dust-to-gas ratios and investigate variations in the dust-to-gas ratio with galaxy metallicity. We find a power law relationship between dust-to-gas ratio and metallicity, consistent with other studies of dust-to-gas ratio compared to metallicity. We find a change in the relation when H$_2$ is not included. This implies that underestimation of $N_{H_2}$ in low-metallicity dwarfs from a too-low CO-to-H$_2$ conversion factor $X_{CO}$ could have produced too low a slope in the derived relationship between dust-to-gas ratio and metallicity. We also compare our extinctions to those derived from fitting the spectral energy distribution (SED) using the Bayesian Extinction and Stellar Tool (BEAST) for NGC 7793 and find systematically lower extinctions from SED-fitting as compared to isochrone matching., Comment: 22 pages, 14 figures, accepted by ApJ

Two dwarf irregular galaxies DDO 187 and NGC 3738 exhibit a striking pattern of star formation: intense star formation is taking place in a large region occupying roughly half of the inner part of the optical galaxy. We use data on the HI distribution and kinematics and stellar images and colors to examine the properties of the environment in the high star formation rate (HSF) halves of the galaxies in comparison with the low star formation rate (LSF) halves. We find that the pressure and gas density are higher on the HSF sides by 30-70%. In addition we find in both galaxies that the HI velocity fields exhibit significant deviations from ordered rotation and there are large regions of high velocity dispersion and multiple velocity components in the gas beyond the inner regions of the galaxies. The conditions in the HSF regions are likely the result of large-scale external processes affecting the internal environment of the galaxies and enabling the current star formation there., Comment: In press in ApJ

The Legacy ExtraGalactic UV Survey (LEGUS) is a multiwavelength Cycle 21 Treasury program on the Hubble Space Telescope. It studied 50 nearby star-forming galaxies in five bands from the near UV to the I-band, combining new Wide Field Camera 3 observations with archival Advanced Camera for Surveys data. LEGUS was designed to investigate how star formation occurs and develops on both small and large scales, and how it relates to the galactic environments. In this paper we present the photometric catalogs for all the apparently single stars identified in the 50 LEGUS galaxies. Photometric catalogs and mosaicked images for all filters are available for download. We present optical and near UV color-magnitude diagrams for all the galaxies. For each galaxy we derived the distance from the tip of the red giant branch. We then used the NUV color-magnitude diagrams to identify stars more massive than 14 Mo, and compared their number with the number of massive stars expected from the GALEX FUV luminosity. Our analysis shows that the fraction of massive stars forming in star clusters and stellar associations is about constant with the star formation rate. This lack of a relation suggests that the time scale for evaporation of unbound structures is comparable or longer than 10 Myr. At low star formation rates this translates to an excess of mass in clustered environments as compared to model predictions of cluster evolution, suggesting that a significant fraction of stars form in unbound systems., Comment: 71 pages, 15 figures, 4 tables, Accpeted for publication on ApJS

The relative abundance of the dust grain types in the interstellar medium (ISM) is directly linked to physical quantities that trace the evolution of galaxies. We study the dust properties of the whole disc of M33 at spatial scales of ~170 pc. This analysis allows us to infer how the relative dust grain abundance changes with the conditions of the ISM, study the existence of a submillimetre excess and look for trends of the gas-to-dust mass ratio (GDR) with other physical properties of the galaxy. For each pixel in the disc of M33 we fit the infrared SED using a physically motivated dust model that assumes an emissivity index beta close to 2. We derive the relative amount of the different dust grains in the model, the total dust mass, and the strength of the interstellar radiation field (ISRF) heating the dust at each spatial location. The relative abundance of very small grains tends to increase, and for big grains to decrease, at high values of Halpha luminosity. This shows that the dust grains are modified inside the star-forming regions, in agreement with a theoretical framework of dust evolution under different physical conditions. The radial dependence of the GDR is consistent with the shallow metallicity gradient observed in this galaxy. The strength of the ISRF derived in our model correlates with the star formation rate in the galaxy in a pixel by pixel basis. Although this is expected it is the first time that a correlation between both quantities is reported. We produce a map of submillimetre excess in the 500 microns SPIRE band for the disc of M33. The excess can be as high as 50% and increases at large galactocentric distances. We further study the relation of the excess with other physical properties of the galaxy and find that the excess is prominent in zones of diffuse ISM outside the main star-forming regions, where the molecular gas and dust surface density are low., Comment: Accepted to Astronomy and Astrophysics, 19 pages, 18 figures plus appendix

We estimate the current extinction-corrected H$\alpha$ star formation rate (SFR) of the different morphological components that shape galaxies (bulges, bars, and disks). We use a multi-component photometric decomposition based on SDSS imaging to CALIFA Integral Field Spectroscopy datacubes for a sample of 219 galaxies. This analysis reveals an enhancement of the central SFR and specific SFR (sSFR $=$ SFR/$M_{\star}$) in barred galaxies. Along the Main Sequence, we find more massive galaxies in total have undergone efficient suppression (quenching) of their star formation, in agreement with many studies. We discover that more massive disks have had their star formation quenched as well. We evaluate which mechanisms might be responsible for this quenching process. The presence of type-2 AGNs plays a role at damping the sSFR in bulges and less efficiently in disks. Also, the decrease in the sSFR of the disk component becomes more noticeable for stellar masses around 10$^{10.5}$ M$_{\odot}$; for bulges, it is already present at $\sim$10$^{9.5}$ M$_{\odot}$. The analysis of the line-of-sight stellar velocity dispersions ($\sigma$) for the bulge component and of the corresponding Faber-Jackson relation shows that AGNs tend to have slightly higher $\sigma$ values than star-forming galaxies for the same mass. Finally, the impact of environment is evaluated by means of the projected galaxy density, $\Sigma$$_{5}$. We find that the SFR of both bulges and disks decreases in intermediate-to-high density environments. This work reflects the potential of combining IFS data with 2D multi-component decompositions to shed light on the processes that regulate the SFR., Comment: 22 pages, 8 figures, 1 table (Accepted for publication in ApJ)

We report the large effort which is producing comprehensive high-level young star cluster (YSC) catalogues for a significant fraction of galaxies observed with the Legacy ExtraGalactic UV Survey (LEGUS) Hubble treasury program. We present the methodology developed to extract cluster positions, verify their genuine nature, produce multiband photometry (from NUV to NIR), and derive their physical properties via spectral energy distribution fitting analyses. We use the nearby spiral galaxy NGC628 as a test case for demonstrating the impact that LEGUS will have on our understanding of the formation and evolution of YSCs and compact stellar associations within their host galaxy. Our analysis of the cluster luminosity function from the UV to the NIR finds a steepening at the bright end and at all wavelengths suggesting a dearth of luminous clusters. The cluster mass function of NGC628 is consistent with a power-law distribution of slopes $\sim -2$ and a truncation of a few times $10^5$ M$_\odot$. After their formation YSCs and compact associations follow different evolutionary paths. YSCs survive for a longer timeframe, confirming their being potentially bound systems. Associations disappear on time scales comparable to hierarchically organized star-forming regions, suggesting that they are expanding systems. We find mass-independent cluster disruption in the inner region of NGC628, while in the outer part of the galaxy there is little or no disruption. We observe faster disruption rates for low mass ($\leq$ $10^4$ M$_\odot$) clusters suggesting that a mass-dependent component is necessary to fully describe the YSC disruption process in NGC628., Comment: accepted for publication in ApJ; 27 pages; complete stellar cluster catalogues will be released in June 2017 at https://legus.stsci.edu

We present a study of the hierarchical clustering of the young stellar clusters in six local (3--15 Mpc) star-forming galaxies using Hubble Space Telescope broad band WFC3/UVIS UV and optical images from the Treasury Program LEGUS (Legacy ExtraGalactic UV Survey). We have identified 3685 likely clusters and associations, each visually classified by their morphology, and we use the angular two-point correlation function to study the clustering of these stellar systems. We find that the spatial distribution of the young clusters and associations are clustered with respect to each other, forming large, unbound hierarchical star-forming complexes that are in general very young. The strength of the clustering decreases with increasing age of the star clusters and stellar associations, becoming more homogeneously distributed after ~40--60 Myr and on scales larger than a few hundred parsecs. In all galaxies, the associations exhibit a global behavior that is distinct and more strongly correlated from compact clusters. Thus, populations of clusters are more evolved than associations in terms of their spatial distribution, traveling significantly from their birth site within a few tens of Myr whereas associations show evidence of disruption occurring very quickly after their formation. The clustering of the stellar systems resembles that of a turbulent interstellar medium that drives the star formation process, correlating the components in unbound star-forming complexes in a hierarchical manner, dispersing shortly after formation, suggestive of a single, continuous mode of star formation across all galaxies., Comment: 21 pages, 15 figures, 4 tables, accepted for publication in ApJ

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

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 conversion of the IR emission into star formation rate can be strongly dependent on the physical properties of the dust, which are affected by the environmental conditions where the dust is embedded. We study here the dust properties of a set of HII regions in the Local Group Galaxy M33 presenting different spatial configurations between the stars, gas and dust to understand the dust evolution under different environments. We model the SED of each region using the DustEM tool and obtain the mass relative to hydrogen for Very Small Grains (YVSG), Polycyclic Aromatic Hydrocarbons (YPAH) and Big Grains (YBG). The relative mass of the VSGs (YVSG/YTOT) is a factor of 1.7 higher for HII regions classified as filled and mixed than for regions presenting a shell structure. The enhancement of VSGs within NGC 604 and NGC 595 is correlated to expansive gas structures with velocities greater than 50 km/s. The gas-to-dust ratio derived for the HII regions in our sample exhibits two regimes related to the HI-H2 transition of the ISM. Regions corresponding to the HI diffuse regime present a gas-to-dust ratio compatible with the expected value if we assume that the gas-to-dust ratio scales linearly with metallicity, while regions corresponding to a H2 molecular phase present a flatter dust-gas surface density distribution. The fraction of VSGs can be affected by the conditions of the interstellar environment: strong shocks of 50-90 km/s existing in the interior of the most luminous HII regions can lead to fragmentation of BGs into smaller ones, while the more evolved shell and clear shell objects provide a more quiescent environment where reformation of dust BG grains might occur. The gas-to-dust variations found in this analysis might imply that grain coagulation and/or gas-phase metals incorporation to the dust mass is occurring in the interior of the HII regions in M33., Comment: Accepted for publication in A&A, 29 pages, 16 figures

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

The HI and CO components of the interstellar medium (ISM) are usually used to derive the dynamical mass M_dyn of nearby galaxies. Both components become too faint to be used as a tracer in observations of high-redshift galaxies. In those cases, the 158 $\mu$m line of atomic carbon [CII] may be the only way to derive M_dyn. As the distribution and kinematics of the ISM tracer affects the determination of M_dyn, it is important to quantify the relative distributions of HI, CO and [CII]. HI and CO are well-characterised observationally, however, for [CII] only very few measurements exist. Here we compare observations of CO, HI, and [CII] emission of a sample of nearby galaxies, drawn from the HERACLES, THINGS and KINGFISH surveys. We find that within R_25, the average [CII] exponential radial profile is slightly shallower than that of the CO, but much steeper than the HI distribution. This is also reflected in the integrated spectrum ("global profile"), where the [CII] spectrum looks more like that of the CO than that of the HI. For one galaxy, a spectrally resolved comparison of integrated spectra was possible; other comparisons were limited by the intrinsic line-widths of the galaxies and the coarse velocity resolution of the [CII] data. Using high-spectral-resolution SOFIA [CII] data of a number of star forming regions in two nearby galaxies, we find that their [CII] linewidths agree better with those of the CO than the HI. As the radial extent of a given ISM tracer is a key input in deriving M_dyn from spatially unresolved data, we conclude that the relevant length-scale to use in determining M_dyn based on [CII] data, is that of the well-characterised CO distribution. This length scale is similar to that of the optical disk., Comment: Accepted for publication in the Astronomical Journal

This paper describes the Third Public Data Release (DR3) of the Calar Alto Legacy Integral Field Area (CALIFA) survey. Science-grade quality data for 667 galaxies are made public, including the 200 galaxies of the Second Public Data Release (DR2). Data were obtained with the integral-field spectrograph PMAS/PPak mounted on the 3.5m telescope at the Calar Alto Observatory. Three different spectral setups are available, i) a low-resolution V500 setup covering the wavelength range 3749-7500 AA (4240-7140 AA unvignetted) with a spectral resolution of 6.0 AA (FWHM), for 646 galaxies, ii) a medium-resolution V1200 setup covering the wavelength range 3650-4840 AA (3650-4620 AA unvignetted) with a spectral resolution of 2.3 AA (FWHM), for 484 galaxies, and iii) the combination of the cubes from both setups (called COMBO), with a spectral resolution of 6.0 AA and a wavelength range between 3700-7500 AA (3700-7140 AA unvignetted), for 446 galaxies. The Main Sample, selected and observed according to the CALIFA survey strategy covers a redshift range between 0.005 and 0.03, spans the color-magnitude diagram and probes a wide range of stellar mass, ionization conditions, and morphological types. The Extension Sample covers several types of galaxies that are rare in the overall galaxy population and therefore not numerous or absent in the CALIFA Main Sample. All the cubes in the data release were processed using the latest pipeline, which includes improved versions of the calibration frames and an even further improved im- age reconstruction quality. In total, the third data release contains 1576 datacubes, including ~1.5 million independent spectra. It is available at http://califa.caha.es/DR3., Comment: 23 pages, 15 figures, accepted to be published in A&A

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

The nearby dwarf starburst galaxy NGC5253 hosts a number of young, massive star clusters, the two youngest of which are centrally concentrated and surrounded by thermal radio emission (the `radio nebula'). To investigate the role of these clusters in the starburst energetics, we combine new and archival Hubble Space Telescope images of NGC5253 with wavelength coverage from 1500 Ang to 1.9 micron in 13 filters. These include H-alpha, P-beta, and P-alpha, and the imaging from the Hubble Treasury Program LEGUS (Legacy Extragalactic UV Survey). The extraordinarily well-sampled spectral energy distributions enable modeling with unprecedented accuracy the ages, masses, and extinctions of the 9 optically brightest clusters (M_V < -8.8) and the two young radio nebula clusters. The clusters have ages ~1-15 Myr and masses ~1x10^4 - 2.5x10^5 M_sun. The clusters' spatial location and ages indicate that star formation has become more concentrated towards the radio nebula over the last ~15 Myr. The most massive cluster is in the radio nebula; with a mass 2.5x10^5 M_sun and an age ~1 Myr, it is 2-4 times less massive and younger than previously estimated. It is within a dust cloud with A_V~50 mag, and shows a clear nearIR excess, likely from hot dust. The second radio nebula cluster is also ~1 Myr old, confirming the extreme youth of the starburst region. These two clusters account for about half of the ionizing photon rate in the radio nebula, and will eventually supply about 2/3 of the mechanical energy in present-day shocks. Additional sources are required to supply the remaining ionizing radiation, and may include very massive stars., Comment: 67 pages; 11 figures; 7 tables. Accepted for publication in the Astrophysical Journal

The Star Formation Rate (SFR) is one of the main parameters used to analyze the evolution of galaxies through time. The need for recovering the light reprocessed by dust commonly requires the use of low spatial resolution far-infrared data. Recombination-line luminosities provide an alternative, although uncertain dust-extinction corrections based on narrow-band imaging or long-slit spectroscopy have traditionally posed a limit to their applicability. Integral Field Spectroscopy (IFS) is clearly the way to overcome such limitation. We obtain integrated H{\alpha}, ultraviolet (UV) and infrared (IR)-based SFR measurements for 272 galaxies from the CALIFA survey at 0.005 < z < 0.03 using single-band and hybrid tracers. We provide updated calibrations, both global and split by properties (including stellar mass and morphological type), referred to H{\alpha}. The extinction-corrected H{\alpha} luminosity agrees with the updated hybrid SFR estimators based on either UV or H{\alpha} plus IR luminosity over the full range of SFRs (0.03-20 M$_{\odot}$ yr$^{-1}$). The coefficient that weights the amount of energy produced by newly-born stars that is reprocessed by dust on the hybrid tracers, a$_{IR}$, shows a large dispersion. However, it does not became increasingly small at high attenuations, as expected if significant highly-obscured H$\alpha$ emission would be missed. Lenticulars, early-type spirals and type-2 AGN host galaxies show smaller coefficients due to the contribution of optical photons and AGN to dust heating. In the Local Universe the H{\alpha} luminosity derived from IFS observations can be used to measure SFR, at least in statistically-significant, optically-selected galaxy samples. The analysis of the SFR calibrations by galaxies properties could be potentially used by other works to study the impact of different selection criteria in the SFR values derived., Comment: 36 pages, 23 Figures, 5 tables. Abridged abstract for arXiv. Accepted for publication in A&A

This paper presents a detailed analysis of two-armed spiral structure in a sample of galax- ies from the Spitzer Infrared Nearby Galaxies Survey (SINGS), with particular focus on the relationships between the properties of the spiral pattern in the stellar disc and the global struc- ture and environment of the parent galaxies. Following Paper I we have used a combination of Spitzer Space Telescope mid-infrared imaging and visible multi-colour imaging to isolate the spiral pattern in the underlying stellar discs, and we examine the systematic behaviours of the observed amplitudes and shapes (pitch angles) of these spirals. In general, spiral morphology is found to correlate only weakly at best with morphological parameters such as stellar mass, gas fraction, disc/bulge ratio, and vflat. In contrast to weak correlations with galaxy structure a strong link is found between the strength of the spiral arms and tidal forcing from nearby companion galaxies. This appears to support the longstanding suggestion that either a tidal interaction or strong bar is a necessary condition for driving grand-design spiral structure. The pitch angles of the stellar arms are only loosely correlated with the pitch angles of the corresponding arms traced in gas and young stars. We find that the strength of the shock in the gas and the contrast in the star formation rate are strongly correlated with the stellar spiral amplitude., Comment: Accepted for publication in MNRAS

The Legacy ExtraGalactic UV Survey (LEGUS) is a Cycle 21 Treasury program on the Hubble Space Telescope, aimed at the investigation of star formation and its relation with galactic environment in nearby galaxies, from the scales of individual stars to those of ~kpc-size clustered structures. Five-band imaging, from the near-ultraviolet to the I-band, with the Wide Field Camera 3, plus parallel optical imaging with the Advanced Camera for Surveys, is being collected for selected pointings of 50 galaxies within the local 12 Mpc. The filters used for the observations with the Wide Field Camera 3 are: F275W(2,704 A), F336W(3,355 A), F438W(4,325 A), F555W(5,308 A), and F814W(8,024 A); the parallel observations with the Advanced Camera for Surveys use the filters: F435W(4,328 A), F606W(5,921 A), and F814W(8,057 A). The multi-band images are yielding accurate recent (<~50 Myr) star formation histories from resolved massive stars and the extinction-corrected ages and masses of star clusters and associations. The extensive inventories of massive stars and clustered systems will be used to investigate the spatial and temporal evolution of star formation within galaxies. This will, in turn, inform theories of galaxy evolution and improve the understanding of the physical underpinning of the gas-star formation relation and the nature of star formation at high redshift. This paper describes the survey, its goals and observational strategy, and the initial science results. Because LEGUS will provide a reference survey and a foundation for future observations with JWST and with ALMA, a large number of data products are planned for delivery to the community., Comment: 67 pages, 9 figures; accepted for publication on the Astronomical 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 have used multi-wavelength Hubble Space Telescope WFC3 data of the starbursting spiral galaxy M83 in order to measure variations in the upper end of the stellar initial mass function (uIMF) using the production rate of ionizing photons in unresolved clusters with ages $\leq$ 8 Myr. As in earlier papers on M51 and NGC 4214, the upper end of the stellar IMF in M83 is consistent with an universal IMF, and stochastic sampling of the stellar populations in the $\lessapprox$ 10$^{3}$ Msun clusters are responsible for any deviations in this universality. The ensemble cluster population, as well as individual clusters, also imply that the most massive star in a cluster does not depend on the cluster mass. In fact, we have found that these small clusters seem to have an over-abundance of ionizing photons when compared to an expected universal or truncated IMF. This also suggests that the presence of massive stars in these clusters does not affect the star formation in a destructive way., Comment: Accepted for publication in The Astrophysical Journal. arXiv admin note: substantial text overlap with arXiv:1302.5006

NGC 2146, a nearby luminous infrared galaxy (LIRG), presents evidence for outflows along the disk minor axis in all gas phases (ionized, neutral atomic and molecular). We present an analysis of the multi-phase starburst driven superwind in the central 5 kpc as traced in spatially resolved spectral line observations, using far-IR Herschel PACS spectroscopy, to probe the effects on the atomic and ionized gas, and optical integral field spectroscopy to examine the ionized gas through diagnostic line ratios. We observe an increased ~250 km/s velocity dispersion in the [OI] 63 micron, [OIII] 88 micron, [NII] 122 micron and [CII] 158 micron fine-structure lines that is spatially coincident with high excitation gas above and below the disk. We model this with a slow ~200 km/s shock and trace the superwind to the edge of our field of view 2.5 kpc above the disk. We present new SOFIA 37 micron observations to explore the warm dust distribution, and detect no clear dust entrainment in the outflow. The stellar kinematics appear decoupled from the regular disk rotation seen in all gas phases, consistent with a recent merger event disrupting the system. We consider the role of the superwind in the evolution of NGC 2146 and speculate on the evolutionary future of the system. Our observations of NGC 2146 in the far-IR allow an unobscured view of the wind, crucial for tracing the superwind to the launching region at the disk center, and provide a local analog for future ALMA observations of outflows in high redshift systems., Comment: 16 pages, 13 figures, accepted for publication in ApJ

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

We investigate the far infrared spectrum of NGC 1266, a S0 galaxy that contains a massive reservoir of highly excited molecular gas. Using the SPIRE-FTS, we detect the $^{12}$CO ladder up to J=(13-12), [C I] and [N II] lines, and also strong water lines more characteristic of UltraLuminous IR Galaxies (ULIRGs). The 12CO line emission is modeled with a combination of a low-velocity C-shock and a PDR. Shocks are required to produce the H2O and most of the high-J 12CO emission. Despite having an infrared luminosity thirty times less than a typical ULIRG, the spectral characteristics and physical conditions of the ISM of NGC 1266 closely resemble those of ULIRGs, which often harbor strong shocks and large-scale outflows., Comment: Accepted for publication in ApJL 11/15/2013

Over the past decade, 3D optical spectroscopy has become the preferred tool for understanding the properties of galaxies and is now increasingly used to carry out galaxy surveys. Low redshift surveys include SAURON, DiskMass, ATLAS3D, PINGS and VENGA. At redshifts above 0.7, surveys such as MASSIV, SINS, GLACE, and IMAGES have targeted the most luminous galaxies to study mainly their kinematic properties. The on-going CALIFA survey ($z\sim0.02$) is the first of a series of upcoming Integral Field Spectroscopy (IFS) surveys with large samples representative of the entire population of galaxies. Others include SAMI and MaNGA at lower redshift and the upcoming KMOS surveys at higher redshift. Given the importance of spatial scales in IFS surveys, the study of the effects of spatial resolution on the recovered parameters becomes important. We explore the capability of the CALIFA survey and a hypothetical higher redshift survey to reproduce the properties of a sample of objects observed with better spatial resolution at lower redshift. Using a sample of PINGS galaxies, we simulate observations at different redshifts. We then study the behaviour of different parameters as the spatial resolution degrades with increasing redshift., Comment: 20 pages, 16 figures. Accepted for publication in A&A

The metal content of a galaxy is a key property for distinguishing between viable galaxy evolutionary scenarios, and it strongly influences many of the physical processes in the interstellar medium. An absolute and robust determination of extragalactic metallicities is essential in constraining models of chemical enrichment and chemical evolution, however, current gas phase abundance determinations from optical fine-structure lines are uncertain to 0.8 dex as conversion of these optical line fluxes to abundances is strongly dependent on the electron temperature of the ionized gas. In contrast, the far-IR emission lines can be used to derive an O++ abundance that is relatively insensitive to temperature, while the ratio of the optical to far-IR lines provides a consistent temperature to be used in the derivation of an O$^+$ abundance. We present observations of the [O III] 88 \micron\ fine-structure line in NGC 628 that were obtained as part of the KINGFISH program. These data are combined with optical IFU data to derive oxygen abundances for seven HII regions. We find the abundance of these regions to all lie between the high and low values of strong line calibrations and in agreement with estimates that assume temperature fluctuations are present in the HII regions., Comment: 12 pages, 7 Figures, 3 Tables, 6 online-only figures, accepted for publication in ApJ

Galaxies, which often contain ionised gas, sometimes also exhibit a so-called low-ionisation nuclear emission line region (LINER). For 30 years this was attributed to a central mass-accreting supermassive black hole (AGN) of low luminosity, making LINER galaxies the largest AGN-sub-population, dominating in numbers over higher luminosity Seyfert galaxies and quasars. This, however, poses a serious problem. While the inferred energy balance is plausible, many LINERs clearly do not contain any other independent signatures of an AGN. Using integral field spectroscopic data from the CALIFA survey, we aim at comparing the observed radial surface brightness profiles with what is expected from illumination by an AGN. Essential for this analysis is a proper extraction of emission-lines, especially weak lines such as the Balmer Hb line which is superposed on an absorption trough. To accomplish this, we use the GANDALF code which simultaneously fits the underlying stellar continuum and emission lines. We show for 48 galaxies with LINER-like emission, that the radial emission-line surface brightness profiles are inconsistent with ionisation by a central point-source and hence cannot be due to an AGN alone. The most probable explanation for the excess LINER-like emission is ionisation by evolved stars during the short but very hot and energetic phase known as post-AGB. This leads us to an entirely new interpretation. Post-AGB stars are ubiquitous and their ionising effect should be potentially observable in every galaxy with gas present and stars older than ~1 Gyr, unless a stronger radiation field from young hot stars or an AGN outshines them. This means that galaxies with LINER-like emission are in fact not a class defined by a property, but rather by the absence of a property. It also explains why LINER emission is observed mostly in massive galaxies with old stars and little star formation., Comment: 8 pages, 7 figures

The use of IFS is since recently allowing to measure the emission line fluxes of an increasingly large number of star-forming galaxies both locally and at high redshift. The main goal of this study is to review the most widely used empirical oxygen calibrations, O3N2 and N2, by using new direct abundance measurements. We pay special attention to the expected uncertainty of these calibrations as a function of the index value or abundance derived and the presence of possible systematic offsets. This is possible thanks to the analysis of the most ambitious compilation of Te-based HII regions to date. This new dataset compiles the Te-based abundances of 603 HII regions extracted from the literature but also includes new measurements from the CALIFA survey. Besides providing new and improved empirical calibrations for the gas abundance, we also present here a comparison between our revisited calibrations with a total of 3423 additional CALIFA HII complexes with abundances derived using the ONS calibration by Pilyugin et al. (2010). The combined analysis of Te-based and ONS abundances allows us to derive their most accurate calibration to date for both the O3N2 and N2 single-ratio indicators, in terms of all statistical significance, quality and coverage of the space of parameters. In particular, we infer that these indicators show shallower abundance dependencies and statistically-significant offsets compared to those of Pettini and Pagel (2004), Nagao et al. (2006) and P\'erez-Montero and Contini (2009). The O3N2 and N2 indicators can be empirically applied to derive oxygen abundances calibrations from either direct abundance determinations with random errors of 0.18 and 0.16, respectively, or from indirect ones (but based on a large amount of data) reaching an average precision of 0.08 and 0.09 dex (random) and 0.02 and 0.08 dex (systematic; compared to the direct estimations),respectively., Comment: 12 pages, 5 figures, accepted for publication in A&A

Chronos is our response to ESA's call for white papers to define the science for the future L2, L3 missions. Chronos targets the formation and evolution of galaxies, by collecting the deepest NIR spectroscopic data, from the formation of the first galaxies at z~10 to the peak of formation activity at z~1-3. The strong emission from the atmospheric background makes this type of survey impossible from a ground-based observatory. The spectra of galaxies represent the equivalent of a DNA fingerprint, containing information about the past history of star formation and chemical enrichment. The proposed survey will allow us to dissect the formation process of galaxies including the timescales of quenching triggered by star formation or AGN activity, the effect of environment, the role of infall/outflow processes, or the connection between the galaxies and their underlying dark matter haloes. To provide these data, the mission requires a 2.5m space telescope optimised for a campaign of very deep NIR spectroscopy. A combination of a high multiplex and very long integration times will result in the deepest, largest, high-quality spectroscopic dataset of galaxies from z=1 to 12, spanning the history of the Universe, from 400 million to 6 billion years after the big bang, i.e. covering the most active half of cosmic history., Comment: White paper for the science definition of ESA's future L2,L3 missions. Updates will be found here: http://www.chronos-mission.eu

Dust has long been identified as a barrier to measuring inherent galaxy properties. However, the link between dust and attenuation is not straightforward and depends on both the amount of dust and its distribution. Herschel imaging of nearby galaxies undertaken as part of the KINGFISH project allows us to map the dust as seen in emission with unprecedented sensitivity and ~1 kpc resolution. We present here new optical integral field unit spectroscopy for eight of these galaxies that provides complementary 100-200 pc scale maps of the dust attenuation through observation of the reddening in both the Balmer decrement and the stellar continuum. The stellar continuum reddening, which is systematically less than that observed in the Balmer decrement, shows no clear correlation with the dust, suggesting that the distribution of stellar reddening acts as a poor tracer of the overall dust content. The brightest HII regions are observed to be preferentially located in dusty regions, and we do find a correlation between the Balmer line reddening and the dust mass surface density for which we provide an empirical relation. Some of the high-inclination systems in our sample exhibit high extinction, but we also find evidence that unresolved variations in the dust distribution on scales smaller than 500 pc may contribute to the scatter in this relation. We caution against the use of integrated A_V measures to infer global dust properties., Comment: 17 pages, 10 figures, accepted for publication in ApJ

Recent results have suggested that the well known mass-metallicity relation has a strong dependence on the star formation rate, to the extent that a three dimensional `fundamental metallicity relation' exists which links the three parameters with minimal scatter. In this work, we use a sample of 4253 local galaxies observed in atomic hydrogen from the ALFALFA survey to demonstrate, for the first time, that a similar fundamental relation (the HI-FMR) also exists between stellar mass, gas-phase metallicity, and HI mass. This latter relation is likely more fundamental, driving the relation between metallicity, SFR and mass. At intermediate masses, the behaviour of the gas fundamental metallicity relation is very similar to that expressed via the star formation rate. However, we find that the dependence of metallicity on HI content persists to the highest stellar masses, in contrast to the `saturation' of metallicity with SFR. It is interesting to note that the dispersion of the relation is very low at intermediate stellar masses (9< log(M*/Msun) <11), suggesting that in this range galaxies evolve smoothy, in an equilibrium between gas inflow, outflow and star formation. At high and low stellar masses, the scatter of the relation is significantly higher, suggesting that merging events and/or stochastic accretion and star formation may drive galaxies outside the relation. We also assemble a sample of galaxies observed in CO. However, due to a small sample size, strong selection bias, and the influence of a metallicity-dependent CO/H2 conversion factor, the data are insufficient to test any influence of molecular gas on metallicity., Comment: 12 pages, 8 figures. Accepted for publication in MNRAS

The production rate of ionizing photons in young < 8 Myr, unresolved stellar clusters in the nearby irregular galaxy NGC 4214 is probed using multi-wavelength Hubble Space Telescope WFC3 data. We normalize the ionizing photon rate by the cluster mass to investigate the upper end of the stellar initial mass function (IMF). We have found that within the uncertainties the upper end of the stellar IMF appears to be universal in this galaxy, and that deviations from a universal IMF can be attributed to stochastic sampling of stars in clusters with masses < 10^3 M_sun. Furthermore, we have found that there does not seem to be a dependence of the maximum stellar mass on the cluster mass. We have also found that for massive clusters, feedback may cause an underrepresentation in Ha luminosities, which needs to be taken into account when conducting this type of analysis., Comment: 11 pages, 7 figures, accepted for publication in ApJ