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.

M33 is a gas rich spiral galaxy of the Local Group. We investigate the relationship between the two major gas cooling lines and the total infrared (TIR) dust continuum. We mapped the emission of gas and dust in M33 using the far-infrared lines of [CII] and [OI](63um) and the TIR. The line maps were observed with Herschel/PACS. These maps have 50pc resolution and form a ~370pc wide stripe along its major axis covering the sites of bright HII regions, but also more quiescent arm and inter-arm regions from the southern arm at 2kpc galacto-centric distance to the south out to 5.7kpc distance to the north. Full-galaxy maps of the continuum emission at 24um from Spitzer/MIPS, and at 70um, 100um, and 160um from PACS were combined to obtain a map of the TIR. TIR and [CII] intensities are correlated over more than two orders of magnitude. The range of TIR translates to a range of far ultraviolet (FUV) emission of G0,obs~2 to 200 in units of the average Galactic radiation field. The binned [CII]/TIR ratio drops with rising TIR, with large, but decreasing scatter. Fits of modified black bodies (MBBs) to the continuum emission were used to estimate dust mass surface densities and total gas column densities. A correction for possible foreground absorption by cold gas was applied to the [OI] data before comparing it with models of photon dominated regions (PDRs). Most of the ratios of [CII]/[OI] and ([CII]+[OI])/TIR are consistent with two model solutions. The median ratios are consistent with one solution at n~2x10^2 cm-3, G0~60, and and a second low-FUV solution at n~10^4 cm-3, G0~1.5. The bulk of the gas along the lines-of-sight is represented by a low-density, high-FUV phase with low beam filling factors ~1. A fraction of the gas may, however, be represented by the second solution.

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.

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.

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.

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.

ABRIGED: Quantifying the number, type and distribution of W-R stars is a key component in the context of galaxy evolution, since they put constraints on the age of the star formation bursts. Nearby galaxies (d<5 Mpc) are particularly relevant in this context since they fill the gap between studies in the Local Group, where individual stars can be resolved, and galaxies in the Local Volume and beyond. We intend to characterize the W-R star population in NGC625, a low-metallicity dwarf galaxy suffering a currently declining burst of star formation. Optical IFS data have been obtained with the VIMOS-IFU covering the starburst region. We estimate the number of W-R stars using a linear combination of 3 W-R templates: 1 early-type nitrogen (WN) star, 1 late-type WN star and 1 carbon-type (WC) star (or oxygen-type (WO) star). Fits using several ensembles of templates were tested. Results were confronted with: i) high spatial resolution HST photometry; ii) numbers of W-R stars in nearby galaxies; iii) model predictions. The W-R star population is spread over the main body of the galaxy, not necessarily coincident with the overall stellar distribution. Our best estimation for the number of W-R stars yields a total of 28 W-R stars in the galaxy, out of which 17 are early- type WN, 6 are late-type WN and 5 are WC stars. The width of the stellar features nicely correlates with the dominant W-R type found in each aperture. The distribution of the different types of WR in the galaxy is roughly compatible with the way star formation has propagated in the galaxy, according to previous findings using HST images. Fits using templates at the metallicity of the LMC yield more reasonable number of W-R than those using templates at the metallicity of the SMC. Given the metallicity of NGC 625, this suggests a non-linear relation between the metallicity and the luminosity of the W-R spectral features.

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.

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.

Context. Measuring star formation at a local scale is important to constrain star formation laws. Yet, it is not clear whether and how the measure of star formation is affected by the spatial scale at which a galaxy is observed. Aims. We want to understand the impact of the resolution on the determination of the spatially resolved star formation rate (SFR) and other directly associated physical parameters such as the attenuation. Methods. We have carried out a multi-scale, pixel-by-pixel study of the nearby galaxy M33. Assembling FUV, Halpha, 8, 24, 70, and 100 micron maps, we have systematically compared the emission in individual bands with various SFR estimators from a resolution of 33 pc to 2084 pc. Results. We have found that there are strong, scale-dependent, discrepancies up to a factor 3 between monochromatic SFR estimators and Halpha+24 micron. The scaling factors between individual IR bands and the SFR show a strong dependence on the spatial scale and on the intensity of star formation. Finally, strong variations of the differential reddening between the nebular emission and the stellar continuum are seen, depending on the specific SFR (sSFR) and on the resolution. At the finest spatial scales, there is little differential reddening at high sSFR. The differential reddening increases with decreasing sSFR. At the coarsest spatial scales the differential reddening is compatible with the canonical value found for starburst galaxies. Conclusions. Our results confirm that monochromatic estimators of the SFR are unreliable at scales smaller than 1 kpc. Furthermore, the extension of local calibrations to high redshift galaxies presents non-trivial challenges as the properties of these systems may be poorly known.

We use optical integral field spectroscopy and 8 and 24 micron mid-IR observations of the giant HII region NGC 588 in the disc of M33 as input and constraints for two-dimensional tailor-made photoionisation models. Two different geometrical approaches are followed for the modelling structure: i) Each spatial element of the emitting gas is studied individually using models which assume that the ionisation structure is complete in each element to look for azimuthal variations across gas and dust. ii) A single model is considered, and the two-dimensional structure of the gas and the dust are assumed to be due to the projection of an emitting sphere onto the sky. The models in both assumptions reproduce the radial profiles of Hbeta surface brightness, the observed number of ionising photons, and the strong optical emission-line relative intensities. The first approach produces a constant-density matter-bounded thin shell of variable thickness and dust-to-gas ratio, while the second gives place to a radiation-bounded thick shell sphere of decreasing particle density. However, the radial profile of the 8/24 microns IR ratio, depending on the gas and dust geometry, only fits well when the thick-shell model is used. The resulting dust-to-gas mass ratio, which was obtained empirically from the derived dust mass using data from Spitzer, also has a better fit using the thick-shell solution. In both approaches, models support the chemical homogeneity, and the ionisation-parameter radial decrease, These results must be taken with care in view of the very low extinction values that are derived from the IR, as compared to that derived from the Balmer decrement. Besides, the IR can be possibly contaminated with the emission from a cloud of diffuse gas and dust above the plane of the galaxy detected at 250 micron Herschel image.

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.

We study the wavelength dependence of the dust emission as a function of position and environment across the disk of M33 at a linear resolution of 160 pc using Spitzer and Herschel photometric data. Expressing the emissivity of the dust as a power law, the power-law exponent (beta) is estimated from two independent approaches designed to properly treat the degeneracy between beta and the dust temperature. Both beta and the dust temperature are higher in the inner disk than in the outer disk, contrary to reported beta-T anti-correlations found in other sources. In the cold + warm dust model, the warm component and the ionized gas (Halpha) have a very similar distribution across the galaxy, demonstrating that the model separates the components in an appropriate fashion. The flocculent spiral arms and the dust lanes are evident in the map of the cold component. Both cold and warm dust column densities are high in star forming regions and reach their maxima toward the giant star forming complexes NGC604 and NGC595. beta declines from close to 2 in the center to about 1.3 in the outer disk. beta is positively correlated with star formation and with molecular gas column, as traced by Halpha and CO emission. The lower dust emissivity index in the outer parts of M33 is likely related to the reduced metallicity (different grain composition) and possibly different size distribution. It is not due to the decrease in stellar radiation field or temperature in a simple way because the FIR-bright regions in the outer disk also have a low beta. Like most spirals, M33 has a (decreasing) radial gradient in star formation and molecular-to-atomic gas ratio such that the regions bright in Halpha or CO tend to trace the inner disk, making it difficult to distinguish between their effects on the dust.

We aim to better understand the heating of the gas by observing the prominent gas cooling line [CII] at 158um in the low-metallicity environment of the Local Group spiral galaxy M33 at scales of 280pc. In particular, we aim at describing the variation of the photoelectric heating efficiency with galactic environment. In this unbiased study, we used ISO/LWS [CII] observations along the major axis of M33, in combination with Herschel PACS and SPIRE continuum maps, IRAM 30m CO 2-1 and VLA HI data to study the variation of velocity integrated intensities. The ratio of [CII] emission over the far-infrared continuum is used as a proxy for the heating efficiency, and models of photon-dominated regions are used to study the local physical densities, FUV radiation fields, and average column densities of the molecular clouds. The heating efficiency stays constant at 0.8% in the inner 4.5kpc radius of the galaxy where it starts to increase to reach values of ~3% in the outskirts at about 6kpc radial distance. The rise of efficiency is explained in the framework of PDR models by lowered volume densities and FUV fields, for optical extinctions of only a few magnitudes at constant metallicity. In view of the significant fraction of HI emission stemming from PDRs, and for typical pressures found in the Galactic cold neutral medium (CNM) traced by HI emission, the CNM contributes ~15% to the observed [CII] emission in the inner 2kpc radius of M33. The CNM contribution remains largely undetermined in the south, while positions between 2 and 7.3kpc radial distance in the north of M33 show a contribution of ~40%+-20%.

Within the framework of the Herschel M 33 extended survey HerM33es we study the Spectral Energy Distribution (SED) of a set of HII regions in M 33 as a function of the morphology. We present a catalogue of 119 HII regions morphologically classified: 9 filled, 47 mixed, 36 shell, and 27 clear shell HII regions. For each object we extract the photometry at twelve available wavelength bands (from FUV-1516A to IR-250mi) and obtain the SED. We also obtain emission line profiles across the regions to study the location of the stellar, ionised gas, and dust components. We find trends for the SEDs related to the morphology, showing that the star and gas-dust configuration affects the ratios of the emission in different bands. The mixed and filled regions show higher emission at 24mi than the shells and clear shells, which could be due to the proximity of the dust to the stellar clusters in the case of mixed and filled regions. The FIR peak for shells and clear shells seems to be located towards longer wavelengths, indicating that the dust is colder for this type of objects.The logarithmic 100/70mi ratio for filled and mixed regions remains constant over one order of magnitude in Halpha and FUV surface brightness, while the shells and clear shells exhibit a wider range of values of almost two orders of magnitude. We derive dust masses and temperatures fitting the individual SEDs with dust models proposed in the literature. The derived dust mass range is between 10^2-10^4 Msun and the cold dust temperature spans T(cold)~12-27 K. The spherical geometrical model proposed for the Halpha clear shells is confirmed by the emission profile obtained from the observations and is used to infer the electron density within the envelope: the typical electron density is 0.7+-0.3 cm^-3, while filled regions can reach values two to five times higher.

We study the far-infrared emission from the nearby spiral galaxy M33 in order to investigate the dust physical properties such as the temperature and the luminosity density across the galaxy. Taking advantage of the unique wavelength coverage (100, 160, 250, 350 and 500 micron) of the Herschel Space Observatory and complementing our dataset with Spitzer-IRAC 5.8 and 8 micron and Spitzer-MIPS 24 and 70 micron data, we construct temperature and luminosity density maps by fitting two modified blackbodies of a fixed emissivity index of 1.5. We find that the 'cool' dust grains are heated at temperatures between 11 and 28 K with the lowest temperatures found in the outskirts of the galaxy and the highest ones in the center and in the bright HII regions. The infrared/submillimeter total luminosity (5 - 1000 micron) is estimated to be 1.9x10^9 Lsun. 59% of the total luminosity of the galaxy is produced by the 'cool' dust grains (~15 K) while the rest 41% is produced by 'warm' dust grains (~55 K). The ratio of the cool-to-warm dust luminosity is close to unity (within the computed uncertainties), throughout the galaxy, with the luminosity of the cool dust being slightly enhanced in the center of the galaxy. Decomposing the emission of the dust into two components (one emitted by the diffuse disk of the galaxy and one emitted by the spiral arms) we find that the fraction of the emission in the disk in the mid-infrared (24 micron) is 21%, while it gradually rises up to 57% in the submillimeter (500 micron). We find that the bulk of the luminosity comes from the spiral arm network that produces 70% of the total luminosity of the galaxy with the rest coming from the diffuse dust disk. The 'cool' dust inside the disk is heated at a narrow range of temperatures between 18 and 15 K (going from the center to the outer parts of the galaxy).

We present observational evidence that leakage of ionising photons from star-forming regions can affect the quantification of the star formation rate (SFR) in galaxies. This effect could partially explain the differences between the SFR estimates using the far ultraviolet (FUV) and the Halpha emission. We find that leakage could decrease the SFR(Ha)/SFR(FUV) ratio by up to a 25 per cent. The evidence is based on the observation that the SFR(Ha)/SFR(FUV) ratio is lower for objects showing a shell Halpha structure than for regions exhibiting a much more compact morphology. The study has been performed on three object samples: low luminosity dwarf galaxies from the Local Volume Legacy survey and star-forming regions in the Large Magellanic Cloud and the nearby Local Group galaxy M33. For the three samples we find differences (1.1-1.4sigma) between the SFR(Ha)/SFR(FUV) for compact and shell objects. Although leakage cannot entirely explain the observed trend of SFR(Ha)/SFR(FUV) ratios for systems with low SFR, we show the mechanism can lead to different SFR estimates when using Halpha and FUV luminosities. Therefore, further study is needed to constrain the contribution of leakage to the low SFR(Ha)/SFR(FUV) ratios observed in dwarf galaxies and its impact on the Halpha flux as a SFR indicator in such objects.

Power spectra of de-projected images of late-type galaxies in gas and/or dust emission are very useful diagnostics of the dynamics and stability of their interstellar medium. Previous studies have shown that the power spectra can be approximated as two power-laws, a shallow one at large scales (larger than 500 pc) and a steeper one at small scales, with the break between the two corresponding to the line-of-sight thickness of the galaxy disk. We present a thorough analysis of the power spectra of the dust and gas emission at several wavelengths in the nearby galaxy M33. In particular, we use the recently obtained images at five wavelengths by PACS and SPIRE onboard Herschel. The large dynamical range (2-3 dex in scale) of most images allow us to determine clearly the change in slopes from -1.5 to -4, with some variations with wavelength. The break scale is increasing with wavelength, from 100 pc at 24 and 100micron to 350 pc at 500micron, suggesting that the cool dust lies in a thicker disk than the warm dust, may be due to star formation more confined to the plane. The slope at small scale tends to be steeper at longer wavelength, meaning that the warmer dust is more concentrated in clumps. Numerical simulations of an isolated late-type galaxy, rich in gas and with no bulge, like M33, are carried out, in order to better interpret these observed results. Varying the star formation and feedback parameters, it is possible to obtain a range of power-spectra, with two power-law slopes and breaks, which nicely bracket the data. The small-scale power-law is indeed reflecting the 3D behaviour of the gas layer, steepening strongly while the feedback smoothes the structures, by increasing the gas turbulence. M33 appears to correspond to a fiducial model with an SFR of $\sim$ 0.7 Mo/yr, with 10% supernovae energy coupled to the gas kinematics.

We have carried out a detailed modeling of the dust Spectral Energy Distribution (SED) of the nearby, starbursting dwarf galaxy NGC 4214. A key point of our modeling is that we distinguish the emission from (i) HII regions and their associated photodissociation regions (PDRs) and (ii) diffuse dust. For both components we apply templates from the literature calculated with a realistic geometry and including radiation transfer. The large amount of existing data from the ultraviolet (UV) to the radio allows the direct measurement of most of the input parameters of the models. We achieve a good fit for the total dust SED of NGC 4214. In the present contribution we describe the available data, the data reduction and the determination of the model parameters, whereas a description of the general outline of our work is presented in the contribution of Lisenfeld et al. in this volume.

We have carried out a detailed modelling of the dust heating and emission in the nearby, starbursting dwarf galaxy NGC 4214. Due to its proximity and the great wealth of data from the UV to the millimeter range (from GALEX, HST, {\it Spitzer}, Herschel, Planck and IRAM) it is possible to separately model the emission from HII regions and their associated photodissociation regions (PDRs) and the emission from diffuse dust. Furthermore, most model parameters can be directly determined from the data leaving very few free parameters. We can fit both the emission from HII+PDR regions and the diffuse emission in NGC 4214 with these models with "normal" dust properties and realistic parameters.

Within the framework of the HerM33es Key Project for Herschel and in combination with multi-wavelength data, we study the Spectral Energy Distribution (SED) of a set of HII regions in the Local Group Galaxy M33. Using the Halpha emission, we perform a classification of a selected HII region sample in terms of morphology, separating the objects in filled, mixed, shell and clear shell objects. We obtain the SED for each HII region as well as a representative SED for each class of objects. We also study the emission distribution of each band within the regions. We find different trends in the SEDs for each morphological type that are related to properties of the dust and their associated stellar cluster. The emission distribution of each band within the region is different for each morphological type of object.

Context: The emission line of [CII] at 158 micron is one of the strongest cooling lines of the interstellar medium (ISM) in galaxies. Aims: Disentangling the relative contributions of the different ISM phases to [CII] emission, is a major topic of the HerM33es program, a Herschel key project to study the ISM in the nearby spiral galaxy M33. Methods: Using PACS, we have mapped the emission of [CII] 158 micron, [OI] 63 micron, and other FIR lines in a 2'x2' region of the northern spiral arm of M33, centered on the HII region BCLMP302. At the peak of H-alpha emission, we have observed in addition a velocity resolved [CII] spectrum using HIFI. We use scatterplots to compare these data with PACS 160 micron continuum maps, and with maps of CO and HI data, at a common resolution of 12 arcsec or 50 pc. Maps of H-alpha and 24 micron emission observed with Spitzer are used to estimate the SFR. We have created maps of the [CII] and [OI] 63 micron emission and detected [NII] 122 micron and NIII 57 micron at individual positions. Results: The [CII] line observed with HIFI is significantly broader than that of CO, and slightly blue-shifted. In addition, there is little spatial correlation between [CII] observed with PACS and CO over the mapped region. There is even less spatial correlation between [CII] and the atomic gas traced by HI. Detailed comparison of the observed intensities towards the HII region with models of photo ionization and photon dominated regions, confirms that a significant fraction, 20--30%, of the observed [CII] emission stems from the ionized gas and not from the molecular cloud. The gas heating efficiency, using the ratio between [CII] and the TIR as a proxy, varies between 0.07 and 1.5%, with the largest variations found outside the HII region.

Within the key project "Herschel M33 extended survey" (HerM33es), we are studying the physical and chemical processes driving star formation and galactic evolution in the nearby galaxy M33, combining the study of local conditions affecting individual star formation with properties only becoming apparent on global scales. Here, we present recent results obtained by the HerM33es team. Combining Spitzer and Herschel data ranging from 3.6um to 500um, along with HI, Halpha, and GALEX UV data, we have studied the dust at high spatial resolutions of 150pc, providing estimators of the total infrared (TIR) brightness and of the star formation rate. While the temperature of the warm dust at high brightness is driven by young massive stars, evolved stellar populations appear to drive the temperature of the cold dust. Plane-parallel models of photon dominated regions (PDRs) fail to reproduce fully the [CII], [OI], and CO maps obtained in a first spectroscopic study of one 2'x2' subregion of M33, located on the inner, northern spiral arm and encompassing the HII region BCLMP302.

We present a set of photoionization models that reproduce simultaneously the observed optical and mid-infrared spatial distribution of the HII region NGC595 in the disk of M33 using the code CLOUDY. Both optical (PMAS-Integral Field Spectroscopy) and mid-infrared (8 mi and 24 mi bands from Spitzer) data provide enough spatial resolution to model in a novel approach the inner structure of the HII region. We define a set of elliptical annular regions around the central ionizing cluster with an uniformity in their observed properties and consider each annulus as an independent thin shell structure. For the first time our models fit the relative surface brightness profiles in both the optical (Halpha, [OII], [OIII]) and the mid-infrared emissions (8 mi and 24 mi), under the assumption of a uniform metallicity (12+log(O/H) = 8.45; Esteban et al. 2009) and an age for the stellar cluster of 4.5 Myr (Malumuth et al. 1996). Our models also reproduce the observed uniformity of the R23 parameter and the increase of the [OII]/[OIII] ratio due to the decrease of the ionization parameter. The variation of the Halpha profile is explained in terms of the differences of the occupied volume (the product of filling factor and total volume of the shell) in a matter-bounded geometry, which also allows to reproduce the observed pattern of the extinction. The 8 mi/24 mi ratio is low (ranging between 0.04 and 0.4) because it is dominated by the surviving of small dust grains in the HII region, while the PAHs emit more weakly because they cannot be formed in these thin HII gas shells. The ratio is also well fitted in our models by assuming a dust-to-gas ratio in each annulus compatible with the integrated estimate for the whole HII region after the 70 mi, and 160 mi Spitzer observations.

We study the post-main sequence stars in NGC 604, the most luminous HII region in M33. Previously, a number of Wolf-Rayet (WR) stars and one red supergiant (RSG) have been discovered. Based on broadband photometry of the region, we present evidence that is consistent with the presence of this RSG and with three more RSG candidates. Using SED fitting based on HST UVIJHK photometry we estimate the ages of the WR stars and RSGs finding that the two populations are from distinct formation episodes with ages 3.2$\pm%1.0Myrs and 12.4$\pm$2.1Myrs, respectively. The RSGs have greater extinctions towards their line of sight than the WR stars consistent with the RSGs producing large amount of dust. Using the WR and RSG populations and similar SED fits to the most massive O stars we estimate that the total stellar mass is (3.8 $\pm$ 0.6) x 10^5Msun. We find a large discrepancy between the expected H{\alpha} flux from such a massive cluster and that one observed. This suggests that 49 (+16,-19) percent of the ionizing photons produced by massive stars in NGC 604 is leaking from this HII region. We also suggest that the implications of an old RSG population mean that if NGC 604 was more distant and only observed in the infrared (IR) it would be difficult to study the youngest burst of star formation due to the contamination of RSGs.

We present an analysis of the first space-based far-IR-submm observations of M 33, which measure the emission from the cool dust and resolve the giant molecular cloud complexes. With roughly half-solar abundances, M33 is a first step towards young low-metallicity galaxies where the submm may be able to provide an alternative to CO mapping to measure their H$_2$ content. In this Letter, we measure the dust emission cross-section $\sigma$ using SPIRE and recent CO and \HI\ observations; a variation in $\sigma$ is present from a near-solar neighborhood cross-section to about half-solar with the maximum being south of the nucleus. Calculating the total H column density from the measured dust temperature and cross-section, and then subtracting the \HI\ column, yields a morphology similar to that observed in CO. The H$_2$/\HI\ mass ratio decreases from about unity to well below 10% and is about 15% averaged over the optical disk. The single most important observation to reduce the potentially large systematic errors is to complete the CO mapping of M 33.

Within the framework of the HERM33ES key project, we are studying the star forming interstellar medium in the nearby, metal-poor spiral galaxy M33, exploiting the high resolution and sensitivity of Herschel. We use PACS and SPIRE maps at 100, 160, 250, 350, and 500 micron wavelength, to study the variation of the spectral energy distributions (SEDs) with galacto-centric distance. Detailed SED modeling is performed using azimuthally averaged fluxes in elliptical rings of 2 kpc width, out to 8 kpc galacto-centric distance. Simple isothermal and two-component grey body models, with fixed dust emissivity index, are fitted to the SEDs between 24 and 500 micron using also MIPS/Spitzer data, to derive first estimates of the dust physical conditions. The far-infrared and submillimeter maps reveal the branched, knotted spiral structure of M33. An underlying diffuse disk is seen in all SPIRE maps (250-500 micron). Two component fits to the SEDs agree better than isothermal models with the observed, total and radially averaged flux densities. The two component model, with beta fixed at 1.5, best fits the global and the radial SEDs. The cold dust component clearly dominates; the relative mass of the warm component is less than 0.3% for all the fits. The temperature of the warm component is not well constrained and is found to be about 60K plus/minus 10K. The temperature of the cold component drops significantly from about 24K in the inner 2 kpc radius to 13K beyond 6 kpc radial distance, for the best fitting model. The gas-to-dust ratio for beta=1.5, averaged over the galaxy, is higher than the solar value by a factor of 1.5 and is roughly in agreement with the subsolar metallicity of M33.

We present Integral Field Spectroscopy (IFS) of NGC 595, one of the most luminous HII regions in M33. This type of observations allows studying the variation of the principal emission-line ratios across the surface of the nebula. At each position of the field of view, we fit the main emission-line features of the spectrum within the spectral range 3650-6990A, and create maps of the principal emission-line ratios for the total surface of the region. The extinction map derived from the Balmer decrement and the absorbed H-alpha luminosity show good spatial correlation with the 24 micron emission from Spitzer. We also show here the capability of the IFS to study the existence of Wolf-Rayet (WR) stars, identifying the previously catalogued WR stars and detecting a new candidate towards the north of the region. The ionization structure of the region nicely follows the H-alpha shell morphology and is clearly related to the location of the central ionizing stars. The electron density distribution does not show strong variations within the HII region nor any trend with the H-alpha emission distribution. We study the behaviour within the HII region of several classical emission-line ratios proposed as metallicity calibrators: while [NII]/Ha and [NII]/[OIII] show important variations, the R23 index is substantially constant across the surface of the nebula, despite the strong variation of the ionization parameter as a function of the radial distance from the ionizing stars. These results show the reliability in using the R23 index to characterize the metallicity of HII regions even when only a fraction of the total area is covered by the observations.

We present a multiwavelength (ultraviolet, infrared, optical and CO) study of a set of luminous HII regions in M33: NGC 604, NGC 595, NGC 592, NGC 588 and IC131. We study the emission distribution in the interiors of the HII regions to investigate the relation between the dust emission at 8 micron and 24 micron and the location of the massive stars and gas. We find that the 24 micron emission is closely related to the location of the ionized gas, while the 8 micron emission is more related to the boundaries of the molecular clouds consistently with its expected association with photodissociation regions (PDRs). Ultraviolet emission is generally surrounded by the H-alpha emission. For NGC 604 and NGC 595, where CO data are available, we see a radial gradient of the emission distribution at the wavelengths studied here: from the center to the boundary of the HII regions we observe ultraviolet, H-alpha, 24 micron, 8 micron and CO emission distributions. We quantify the star formation for our HII regions using the integrated fluxes at the set of available wavelengths, assuming an instantaneous burst of star formation. We show that a linear combination of 24 micron and H-alpha emission better describes the star formation for these objects than the dust luminosities by themselves. For NGC 604, we obtain and compare extinction maps derived from the Balmer decrement and from the 24 micron and H-alpha emission line ratio. Although the maps show locally different values in extinction, we find similar integrated extinctions derived from the two methods. We also investigate here the possible existence of embedded star formation within NGC 604.

We present a study of the extinction, traced by the Balmer decrement, in HII regions in the dwarf galaxies NGC 1569 and NGC 4214. We find that the large-scale extinction around the most prominent HII regions in both galaxies forms a shell in which locally the intrinsic extinction can adopt relatively high values ($A_V = 0.8 - 0.9$ mag) despite the low metallicity and thus the low overall dust content. The small-scale extinction (spatial resolution $\sim$0.3'') shows fluctuations that are most likely due to variations in the dust distribution. We compare the distribution of the extinction to that of the dust emission, traced by Spitzer emission at 8 and 24 \mi, and to the emission of cold dust at 850 \mi. We find in general a good agreement between all tracers, expect for the 850 \mi emission in NGC 4214 which is more extended than the extinction and the other emissions. Whereas in NGC 1569 the dust emission at all wavelengths is very similar, NGC 4214 shows spatial variations in the 24-to-850 \mi ratio. We furthermore compared the 24 \mi and the extinction-corrected Halpha emission from HII regions in a sample of galaxies with a wide range of metallicities and found a good correlation between both emissions, independent of metallicity. We suggest that this lack of dependence on metallicity might be due to the formation of dust shells with a relatively constant opacity, like the ones observed here, around ionizing stars.

We investigate the use of the rest-frame 24microns luminosity as an indicator of the star formation rate (SFR) in galaxies with different metallicities by comparing it to the (extinction corrected) Halpha luminosity. We carry out this analysis in 2 steps: First, we compare the emission from HII regions in different galaxies with metallicities between 12+log(O/H) = 8.1 and 8.9. We find that the 24microns and the extinction corrected Halpha luminosities from individual HII regions follow the same correlation for all galaxies, independent of their metallicity. Second, the role of metallicity is explored further for the integrated luminosity in a sample of galaxies with metallicities in the range of 12+log(O/H) = 7.2 - 9.1. For this sample we compare the 24microns and Halpha luminosities integrated over the entire galaxies and find a lack of the 24microns emission for a given Halpha luminosity for low metallicity objects, likely reflecting a low dust content. These results suggest that the 24microns luminosity is a good metallicity independent tracer for the SFR in individual HII regions. On the other hand, metallicity has to be taken into account when using the 24microns luminosity as a tracer for the SFR of entire galaxies.

Large HI shells, with diameters of hundreds of pc and expansion velocities of 10-20kms-1 are well observed features of local gas rich galaxies. These shells could well be predicted as a result of the impact of OB associations on the ISM, but doubt has been cast on this scenario by the apparent absence of OB stars close to the centres of a large fraction of these shells in recent observations of the SMC. Using Fabry-Perot scanned Halpha emission line mapping of nearby galaxy discs we have detected, in all the HII regions where the observations yield sufficient angular resolution and S:N ratio, dominant Halpha shells with radii a few tens of pc, expanding at velocities of 50-100kms-1. We have applied a simple dynamically consistent framework in which we can extrapolate the properties of the observed Halpha shells to a few 10^7yr after the formation of the OB stars. The framework includes the dynamical inputs of both winds and SNe on the surrounding ISM. The results give quantitative statistical support to the hypothesis that the Halpha emitting shells are generic progenitors of the HI shells. During the time taken for an expanding shell to reach the size of a typical HI shell, the OB association may well lose its most luminous stars so the absence of such stars near the centres of many of the HI shells is well explained in this scenario.

We study two dimensional Fabry-Perot interferometric observations of the nearby face-on late-type spiral galaxy, NGC 628, in order to analyse the ionized gas component of the interstellar medium. Covering the galaxy out to a radius larger than 12 kpc, and with a spatial sampling of 1.6 arcsec, we investigate the large-scale dynamics as well as feedback from individual HII regions into their surrounding medium. We study the role of gravitational perturbations along with that of external triggers which can disturb the kinematics and morphology of NGC 628. We verify the presence of an inner rapidly rotating disc-like component in NGC 628, which we interpret as caused by slow secular evolution of the large-scale spiral arms and oval structure. In combination with auxiliary data, we find indication for that gas is falling in from the outer parts towards the central regions, where a nuclear ring has formed at the location of the inner Lindblad resonance radius of an m=2 perturbation which could help build a pseudo-bulge in NGC 628. Moreover, we calculate radial profiles of the emission-line velocity dispersion which we use to study the role of feedback from individual HII regions. The mean velocity dispersion for the ionized gas (even when excluding pixels belonging to individual HII regions) is almost constant out to 12 kpc, although it varies from 14 to 20 km/s, with a steady decline in the outer parts. The current paper demonstrates a number of tools that we have developed for building a solid frame work for studying the evolution of structure in spiral galaxies using two dimensional kinematic observations.

We study two dimensional Fabry-Perot interferometric observations of the nearby face-on late-type spiral galaxy, NGC 628. We investigate the role of the individual HII regions together with the large-scale gravitational mechanisms which govern star formation and overall evolution in spiral galaxies. Our kinematical analysis (reinforced by literature maps in HI and CO at lower angular resolution) enables us to verify the presence of an inner rapidly rotating inner disk-like component which we attribute to long term secular evolution of the large-scale spiral arms and oval structure. We find that gas is falling in from the outer parts towards the bluer central regions. This could be an early phase in the formation of a pseudo-bulge. We find signatures of radial motions caused by an m=2 perturbation, which are likely to be responsible for the inflow of material forming the circumnuclear ring and the rapidly rotating inner structure.

We investigate the distribution of the intrinsic extinction in NGC1569 using an extinction map derived from the H_alpha/H_beta emission line ratio. We compare the extinction distribution to that of the dust emission traced by SPITZER IRAC(8microns) and MIPS (24microns) maps. The intrinsic extinction shows spatial variations, from zones with negligible extinction to zones with values up to A(V)=0.8mag. We find an extinction shell and establish a relation between this shell and the interstellar expanding structure produced by the stellar winds comig from the Super Star Cluster (SSCs) A and B in the center of the galaxy. We suggest that the extinction shell has been produced by the accumulation of dust at the border of the shell. Although we find a good spatial correlation between the Balmer extinction and infrared emission, there is a spatial displacement between the 8microns and 24microns maxima and the maximum in Balmer extinction which needs further investigation.

We investigate spatial the distribution of the intrinsic extinction in the starburst dwarf galaxy NGC1569 creating an extinction map of the whole galaxy derived from the Halpha/Hbeta emission line ratio. We differentiate the extinction in the HII regions from the extinction of the diffuse gas. The intrinsic extinction shows considerable variations over the plane of the galaxy, from negligible extinction up to highest values of A(V)=0.8mag. The extinction map shows small scale clumpy structures possibly due to a clumpy dust distribution. We also identify in this map a shell structure, for which we establish a causal relation with the expanding gas structure produced by the stellar winds coming from the Super Star Clusters (SSC) in the center of the galaxy. The comparison of the spatial profiles of the extinction, dust and gaseous emissions crossing the border of the shell shows a layered structure; the peak of this Halpha distribution lying closest to the SSC A, followed outwards by the peak of the extinction and at a still larger distances by the bulk of the atomic gas. We suggest that the extinction shell has been produced by the SSCs and that it can be explained by the accumulation of dust at the border of this ionized gas structure.

We present an analysis of the integrated Halpha emission line profiles for the HII region population of the spiral galaxies NGC 1530, NGC 6951 and NGC 3359. We show that 70% of the line profiles show two or three Gaussian components. The relations between the Halpha luminosity and non-thermal line width for the HII regions of the three galaxies are studied and compared with the relation found taken all the HII regions of the three galaxies as a single distribution. A clearer envelope in non-thermal line width is found when only those HII regions with non-thermal line width bigger than 13kms are considered. The linear fit for the envelope is logL=36.8+2.0*log(sigma). The masses of the HII regions on the envelope using the virial theorem and the mass estimates from the Halpha luminosity are comparable, which offers evidence that the HII regions on the envelope are virialized systems, while the remaining regions, the majority, are not in virial equilibrium.

We study the presence of low intensity high velocity components, which we have termed wing features in the integrated Halpha emission line profiles of the HII region populations of the spiral barred galaxies NGC 1530, NGC 3359 and NGC 6951. We find that more than a third of the HII region line profiles in each galaxy show these components. The highest fraction is obtained in the galaxy whose line profiles show the best S:N, which suggests that wing features of this type may well exist in most, if not all, HII region line profiles. Applying selection criteria to the wing features, we obtain a sample of HII regions with clearly defined high velocity components in their profiles. Deconvolution of a representative sample of the line profiles eliminates any doubt that the wing features could possibly be due to instrumental effects. We present an analysis of the high velocity low intensity features fitting them with Gaussian functions; the emission measures, central velocities and velocity dispersions for the red and blue features take similar values. We interpret the features as signatures of expanding shells inside the HII regions. Up to a shell radius of R(shell)~0.2R(reg), the stellar winds from the central ionizing stars appear to satisfy the energy and momentum requirements for the formation and driving the shell. Several examples of the most luminous HII regions show that the shells appear to have larger radii; in these cases additional mechanisms may well be needed to explain the kinetic energies and momenta of the shells.

We present emission line mapping of the strongly barred galaxy NGC 1530 obtained from Fabry-Perot interferometry in Halpha, at significantly enhanced angular resolution compared with published studies. The main point of the work is to examine in detail the non-circular components of the velocity field of the gas, presumably induced by the strongly non-axisymmetric gravitational potential of the bar. These reveal quasi-elliptical motions with speeds of order 100 km/s aligned with the bar. It also shows how these flows swing in towards and around the nucleus as they cross a circumnuclear resonance, from the dominant "x1 orbits" outside the resonance to "x2 orbits" within it. Comparing cross-sections of the residual velocity map along and across the bar with the Halpha intensity map indicates a systematic offset between regions of high non-circular velocity and massive star formation. To investigate further we produce velocity gradient maps along and across the bar. These illustrate very nicely the shear compression of the gas, revealed by the location of the dust lanes along loci of maximum velocity gradient perpendicular to the bar. They also show clearly how shear acts to inhibit massive star formation, whereas shocks act to enhance it. Although the inhibiting effect of gas shear flow on star formation has long been predicted, this is the clearest observational illustration so far of the effect. It is also the clearest evidence that dust picks out shock-induced inflow along bars. These observations should be of considerable interest to those modelling massive star formation in general.

We present spherically symmetric and plane parallel photoionization models of NGC 346, an HII region in the Small Magellanic Cloud. The models are based on CLOUDY and on the observations of Peimbert, Peimbert, & Ruiz (2000). We find that approximately 45% of the H ionization photons escape from the HII region providing an important ionizing source for the low density interstellar medium of the SMC. The predicted I(4363)/I(5007) value is smaller than that observed, probably implying that there is an additional source of energy not taken into account by the models. From the ionization structure of the best model and the observed line intensities we determine the abundances of N, Ne, S, Ar, and Fe relative to O.