Age Of Stellar Population Research Articles

A systematic study of the inner rotation curves of galaxies observed as part of the GASS and COLD GASS surveys

We present a systematic analysis of the rotation curves of 187 galaxies with masses greater than 10^10 M_sol, with atomic gas masses from the GALEX Arecibo Sloan Survey (GASS), and with follow-up long-slit spectroscopy from the MMT. Our analysis focuses on stellar rotation curves derived by fitting stellar template spectra to the galaxy spectra binned along the slit. In this way, we are able to obtain accurate rotation velocity measurements for a factor of 2 more galaxies than possible with the Halpha line. Galaxies with high atomic gas mass fractions are the most dark-matter dominated galaxies in our sample and have dark matter halo density profiles that are well fit by Navarro, Frenk & White profiles with an average concentration parameter of 10. The inner slopes and of the rotation curves correlate more strongly with stellar population age than with galaxy mass or structural parameters. At fixed stellar mass, the rotation curves of more actively star-forming galaxies have steeper inner slopes than less actively star-forming galaxies. The ratio between the galaxy specific angular momentum and the total specific angular momentum of its dark matter halo, R_j, correlates strongly with galaxy mass, structure and gas content. Low mass, disk-dominated galaxies with atomic gas mass fractions greater than 20% have median values of R_j of around 1, but massive, bulge-dominated galaxies have R_j=0.2-0.3. We argue that these trends can be understood in a picture where gas inflows triggered by disk instabilities lead to the formation of passive, bulge-dominated galaxies with low specific angular momentum.

A MEGACAM SURVEY OF OUTER HALO SATELLITES. IV. TWO FOREGROUND POPULATIONS POSSIBLY ASSOCIATED WITH THE MONOCEROS SUBSTRUCTURE IN THE DIRECTION OF NGC 2419 AND KOPOSOV 2

The origin of the Galactic halo stellar structure known as the Monoceros ring is still under debate. In this work, we study that halo substructure using deep CFHT wide-field photometry obtained for the globular clusters NGC2419 and Koposov2, where the presence of Monoceros becomes significant because of their coincident projected position. Using Sloan Digital Sky Survey photometry and spectroscopy in the area surrounding these globulars and beyond, where the same Monoceros population is detected, we conclude that a second feature, not likely to be associated with Milky Way disk stars along the line-of-sight, is present as foreground population. Our analysis suggests that the Monoceros ring might be composed of an old stellar population of age t ~ 9Gyr and a new component ~ 4Gyr younger at the same heliocentric distance. Alternatively, this detection might be associated with a second wrap of Monoceros in that direction of the sky and also indicate a metallicity spread in the ring. The detection of such a low-density feature in other sections of this halo substructure will shed light on its nature.

EVOLUTION OF MASS FUNCTIONS OF COEVAL STARS THROUGH WIND-MASS LOSS AND BINARY INTERACTIONS

Accurate determinations of stellar mass functions and ages of stellar populations are crucial to much of astrophysics. We analyze the evolution of stellar mass functions of coeval main-sequence stars, including all relevant aspects of single and binary star evolution. We show that the slope of the upper part of the mass function in a stellar cluster can be quite different from the slope of the initial mass function. Wind-mass loss from massive stars leads to an accumulation of stars which is visible as a peak at the high-mass end of mass functions, thereby flattening the mass function slope. Mass accretion and mergers in close binary systems create a tail of rejuvenated binary products. These blue straggler stars extend the single star mass function by up to a factor of 2 in mass and can appear up to 10 times younger than their parent stellar cluster. Cluster ages derived from their most massive stars that are close to the turn-off may thus be significantly biased. To overcome such difficulties, we propose the use of the binary tail of stellar mass functions as an unambiguous clock to derive the cluster age because the location of the onset of the binary tail identifies the cluster turn-off mass. It is indicated by a pronounced jump in the mass function of old stellar populations and by the wind-mass loss peak in young stellar populations. We further characterize the binary induced blue straggler population in star clusters in terms of their frequency, binary fraction, and apparent age.

Open clusters

Fil: Aidelman, Yael Judith. Consejo Nacional de Investigaciones Cientificas y Tecnicas. Centro Cientifico Tecnologico la Plata. Instituto de Astrofisica de la Plata; Argentina. Universidad Nacional de la Plata. Facultad de Ciencias Astronomicas y Geofisicas; Argentina

EVIDENCE FOR (AND AGAINST) PROGENITOR BIAS IN THE SIZE GROWTH OF COMPACT RED GALAXIES

Most massive passive galaxies are compact at high redshifts, but similarly compact massive galaxies are rare in the local universe. The most common interpretation of this phenomenon is that massive galaxies have grown in size by a factor of about five since redshift z=2. An alternative explanation is that recently quenched massive galaxies are larger (a "progenitor bias"). In this paper we explore the importance of progenitor bias by looking for systematic differences in the stellar populations of compact early-type galaxies in the DEEP2 survey as a function of size. Our analysis is based on applying the statistical technique of bootstrap resampling to constrain differences in the median ages of our samples and to begin to characterize the distribution of stellar populations in our co-added spectra. The light-weighted ages of compact early-type galaxies at redshifts 0.5 < z < 1.4 are compared to those of a control sample of larger galaxies at similar redshifts. We find that massive compact early-type galaxies selected on the basis of red color and high bulge-to-total ratio are younger than similarly selected larger galaxies, suggesting that size growth in these objects is not driven mainly by progenitor bias, and that individual galaxies grow as their stellar populations age. However, compact early-type galaxies selected on the basis of image smoothness and high bulge-to-total ratio are older than a control sample of larger galaxies. Progenitor bias will play a significant role in defining the apparent size changes of early-type galaxies if they are selected on the basis of the smoothness of their light distributions.

Role of feedback in AGN-host coevolution: A study from partially obscured active galactic nuclei

Partially obscured AGNs within a redshift range z=0.011∼0.256 are used to re-study the role of feedback in the AGN-host coevolution issue in terms of their [OIII]λ5007 emission line profile. The spectra of these objects enable us to determine the AGN’s accretion properties directly from their broad Hα emission. This is essential for getting rid of the “circular reasoning” in our previous study of narrow emission-line galaxies, in which the [OIII] emission line was used not only as a proxy of AGN’s bolometric luminosity, but also as a diagnostic of outflow. In addition, the measurement of Dn(4000) index is improved by removing an underlying AGN’s continuum according to the corresponding broad Hα emission. With these improvements, we confirm and reinforce the correlation between L/LEdd and stellar population age. More important is that this correlation is found to be related to both [OIII] line blue asymmetry and bulk blueshift velocity, which suggests a linkage between SMBH growth and host star formation through the feedback process. The current sample of partially obscured AGNs shows that the composite galaxies have younger host stellar population, higher Eddington ratio, less significant [OIII] blue wing and smaller bulk [OIII] line shift than do the Seyfert galaxies.

ALMA OBSERVATIONS OF THE HOST GALAXY OF GRB 090423 ATz= 8.23: DEEP LIMITS ON OBSCURED STAR FORMATION 630 MILLION YEARS AFTER THE BIG BANG

We present rest-frame far-infrared (FIR) and optical observations of the host galaxy of GRB090423 at z=8.23 from the Atacama Large Millimeter Array (ALMA) and the Spitzer Space Telescope, respectively. The host remains undetected to 3-sigma limits of Fnu(222 GHz)<33 microJy and Fnu(3.6 micron)<81 nJy. The FIR limit is about 20 times fainter than the luminosity of the local ULIRG Arp220, and comparable to the local starburst M82. Comparing to model spectral energy distributions we place a limit on the IR luminosity of L_IR(8-1000 micron)<3e10 Lsun, corresponding to a limit on the obscured star formation rate of SFR_IR<5 Msun/yr; for comparison, the limit on the unobscured star formation rate from Hubble Space Telescope rest-frame UV observations is SFR_UV<1 Msun/yr. We also place a limit on the host galaxy stellar mass of <5e7 Msun (for a stellar population age of 100 Myr and constant star formation rate). Finally, we compare our millimeter observations to those of field galaxies at z>4 (Lyman break galaxies, Ly-alpha emitters, and submillimeter galaxies), and find that our limit on the FIR luminosity is the most constraining to date, although the field galaxies have much larger rest-frame UV/optical luminosities than the host of GRB090423 by virtue of their selection techniques. We conclude that GRB host galaxies at z>4, especially those with measured interstellar medium metallicities from afterglow spectroscopy, are an attractive sample for future ALMA studies of high redshift obscured star formation.

Modelling the ages and metallicities of early-type galaxies in Fundamental Plane space

Recent observations have probed the formation histories of nearby elliptical galaxies by tracking correlations between the stellar population parameters, age and metallicity, and the structural parameters that enter the Fundamental Plane, radius and velocity dispersion sigma. These studies have found intriguing correlations between these four parameters. In this work, we make use of a semi-analytic model, based on halo merger trees extracted from the Bolshoi cosmological simulation, that predicts the structural properties of spheroid-dominated galaxies based on an analytic model that has been tested and calibrated against an extensive suite of hydrodynamic+N-body binary merger simulations. We predict the radius, sigma, luminosity, age, and metallicity of spheroid-dominated galaxies, enabling us to compare directly to observations. Our model predicts a strong correlation between age and sigma for early-type galaxies, and no significant correlation between age and radius, in agreement with observations. In addition we predict a strong correlation between metallicity and sigma, and a weak correlation between metallicity and radius, in qualitative agreement with observations. We find that the correlations with sigma arise as a result of the strong link between sigma and the galaxy assembly time. Minor mergers produce a large change in radius while leaving sigma nearly the same, which explains the weaker trends with radius.

The dark nature of GRB 130528A and its host galaxy

We study the dark nature of GRB 130528A through multi-wavelength observations and conclude that the main reason for the optical darkness is local extinction inside of the host galaxy. Automatic observations were performed at BOOTES-4/MET robotic telescope. We also triggered target of opportunity (ToO) observation at the OSN, IRAM PdBI and the GTC+OSIRIS. The host galaxy photometric observations in optical to near-infrared (nIR) wavelengths were achieved through large ground-based aperture telescopes, such as the 10.4m GTC, the 4.2m WHT, 6m BTA, and the 2m LT. Based on these observations, spectral energy distributions (SED) for the host galaxy and afterglow were constructed. Thanks to mm observations at PdBI, we confirm the presence of a mm source within the XRT error circle that faded over the course of our observations and identify the host galaxy. However, we do not find any credible optical source within early observations with BOOTES-4/MET and 1.5m OSN telescopes. Spectroscopic observation of this galaxy by GTC showed a single faint emission line that likely corresponds to [OII] 3727\{AA} at a redshift of 1.250+/-0.001 implying a SFR(M_sun/yr) > 6.18 M_sun/yr without correcting for dust extinction. The probable extinction was revealed through analysis of the afterglow SED, resulting in a value of AV >= ~ 0.9 at the rest frame, this is comparable to extinction levels found among other dark GRBs. The SED of the host galaxy is explained well (chi2/d.o.f.=0.564) by a luminous (MB=-21.16), low-extinction (AV =0, rest frame), and aged (2.6 Gyr) stellar population. We can explain this apparent contradiction in global and line-of-sight extinction if the GRB birth place happened to lie in a local dense environment. In light of having relatively small specific SFR (SSFR) ~ 5.3 M_sun/yr (L/L_star)-1, this also could explain the age of the old stellar population of host galaxy.

THE ASSEMBLY HISTORIES OF QUIESCENT GALAXIES SINCEz= 0.7 FROM ABSORPTION LINE SPECTROSCOPY

We present results from modeling the optical spectra of a large sample of quiescent galaxies between 0.1 < z < 0.7 from the Sloan Digital Sky Survey (SDSS) and the AGN and Galaxy Evolution Survey (AGES). We examine how the stellar ages and abundance patterns of galaxies evolve over time as a function of stellar mass from 10^{9.6}-10^{11.8} Msun. Galaxy spectra are stacked in bins of mass and redshift, and modeled over a wavelength range from 4000 A to 5500 A. Full spectrum stellar population synthesis modeling provides estimates of the age and the abundances of the elements Fe, Mg, C, N, and Ca. We find negligible evolution in elemental abundances at fixed stellar mass over roughly 7 Gyr of cosmic time. In addition, the increase in stellar ages with time for massive galaxies is consistent with passive evolution since z = 0.7. Taken together, these results favor a scenario in which the inner ~ 0.3-3 R_e of massive quiescent galaxies have been passively evolving over the last half of cosmic time. Interestingly, the derived stellar ages are considerably younger than the age of the universe at all epochs, consistent with an equivalent single-burst star formation epoch of z < 1.5. These young stellar population ages coupled with the existence of massive quiescent galaxies at z > 1 indicate the inhomogeneous nature of the z < 0.7 quiescent population. The data also permit the addition of newly-quenched galaxies at masses below ~10^{10.5} Msun at z < 0.7. Additionally, we analyze very deep Keck DEIMOS spectra of the two brightest quiescent galaxies in a cluster at z = 0.83. There is tentative evidence that these galaxies are older than their counterparts in low-density environments. In the Appendix, we demonstrate that our full spectrum modeling technique allows for accurate and reliable modeling of galaxy spectra to low S/N (~20 A^{-1}) and/or low spectral resolution (R ~ 500).

The supergiant shell with triggered star formation in the dwarf irregular galaxy IC 2574: neutral and ionized gas kinematics

We analyse the ionized gas kinematics in the star formation regions of the supergiant shell (SGS) of the IC 2574 galaxy using observations with the Fabry-Perot interferometer at the 6-m telescope of SAO RAS; the data of the THINGS survey are used to analyze the neutral gas kinematics in the area. We perform the 'derotation' of the H-alpha and HI data cubes and show its efficiency in kinematics analysis. We confirm the SGS expansion velocity 25 km/s obtained by Walter & Brinks (1999) and conclude that the SGS is located at the far side of the galactic disc plane. We determine the expansion velocities, kinematic ages, and the required mechanical energy input rates for four star formation complexes in the walls of the SGS; for the remaining ones we give the limiting values of the above parameters. A comparison with the age and energy input of the complexes' stellar population shows that sufficient energy is fed to all HII regions except one. We discuss in detail the possible nature of this region and that of another one, which was believed to be an SNR according to radio observations. We measured the expansion velocity of the latter and confirm its identification as an old SNR. Our observations allowed us to identify a faint diffuse H-alpha emission inside the SGS which was never observed before.

Probing the nuclear star cluster of galaxies with extremely large telescopes

The unprecedented sensitivity and spatial resolution of next generation ground-based extremely large telescopes (ELTs) will open a completely new window on the study of resolved stellar populations. In this paper we study the feasibility of the analysis of nuclear star cluster (NSC) stellar populations with ELTs. To date, NSC stellar population studies are based on the properties of their integrated light. NSC are in fact observed as unresolved sources even with the HST. We explore the possibility to obtain direct estimates of the age of NSC stellar populations from photometry of main-sequence turn-off stars. We simulated ELT observations of NSCs at different distances and with different stellar populations. Photometric measurements on each simulated image were analysed in detail and results about photometric accuracy and completeness are reported here. We found that the main-sequence turn-off is detectable -and therefore the age of stellar populations can be directly estimated- up to 2Mpc for old, up to 3 Mpc for intermediate-age and up to 4-5 Mpc for young stellar populations. We found that for this particular science case the performances of TMT and E-ELT are of comparable quality.

Probing the circumnuclear stellar populations of starburst galaxies in the near-infrared

We employ the NASA Infrared Telescope Facility's near-infrared spectrograph SpeX at 0.8-2.4$\mu$m to investigate the spatial distribution of the stellar populations (SPs) in four well known Starburst galaxies: NGC34, NGC1614, NGC3310 and NGC7714. We use the STARLIGHT code updated with the synthetic simple stellar populations models computed by Maraston (2005, M05). Our main results are that the NIR light in the nuclear surroundings of the galaxies is dominated by young/intermediate age SPs ($t \leq 2\times10^9$yr), summing from $\sim$40\% up to 100\% of the light contribution. In the nuclear aperture of two sources (NGC1614 and NGC3310) we detected a predominant old SP component ($t > 2\times10^9$yr), while for NGC34 and NGC7714 the younger component prevails. Furthermore, we found evidence of a circumnuclear star formation ring-like structure and a secondary nucleus in NGC1614, in agreement with previous studies. We also suggest that the merger/interaction experienced by three of the galaxies studied, NGC1614, NGC3310 and NGC7714 can explain the lower metallicity values derived for the young SP component of these sources. In this scenario the fresh unprocessed metal poorer gas from the destroyed/interacting companion galaxy is driven to the centre of the galaxies and mixed with the central region gas, before star formation takes place. In order to deepen our analysis, we performed the same procedure of SP synthesis using Maraston (2011, M11) EPS models. Our results show that the newer and higher resolution M11 models tend to enhance the old/intermediate age SP contribution over the younger ages.

2D stellar population and gas kinematics of the inner kiloparsec of the post-starburst quasar SDSS J0330−0532

We have used optical Integral Field Spectroscopy in order to map the star formation history of the inner kiloparsec of the Post-Starburst Quasar (PSQ) J0330--0532 and to map its gas and stellar kinematics as well as the gas excitation. PSQs are hypothesized to represent a stage in the evolution of galaxies in which the star formation has been recently quenched due to the feedback of the nuclear activity, as suggested by the presence of the post-starburst population at the nucleus. We have found that the old stellar population (age $\ge$ 2.5 Gyr) dominates the flux at 5100 \AA\ in the inner 0.26 kpc, while both the post-starburst (100 Myr $\le$ age $<$ 2.5 Gyr) and starburst (age $<$ 100 Myr) components dominate the flux in a circumnuclear ring at $\approx$0.5 kpc from the nucleus. With our spatially resolved study we do not have found any post-starburst stellar population in the inner 0.26\,kpc. On the other hand, we do see the signature of AGN feedback in this region, which does not reach the circumnuclear ring where the post-starburst population is observed. We thus do not support the quenching scenario for the J0330-0532. In addition, we have concluded that the strong signature of the post-starburst population in larger aperture spectra (e.g. from Sloan Digital Sky Survey) is partially due to the combination of the young and old age components. Based on the M$_{\rm BH}-\sigma_{\rm star}$ relationship and the stellar kinematics we have estimated a mass for the supermassive black hole of 1.48 $\pm$ 0.66 $\times$10$^7$ M$_\odot$.

The beauty of resolution: The SN Ib factory NGC 2770 spatially resolved

AbstractThe late-type spiral NGC 2770 hosted 3 Type Ib supernovae (SNe) in or next to star-forming regions in its outer spiral arms. We study the properties of the SN sites and the galaxy at different spatial resolutions to infer propeties of the SN progenitors and the SF history of the galaxy. Several 3D techniques are used and, for the first time, we present images of metallicity, shocks and stellar population ages from OSIRIS/GTC imaging with tunable narrowband filters.

DYNAMO – II. Coupled stellar and ionized-gas kinematics in two low-redshift clumpy discs

We study the spatially resolved stellar kinematics of two star-forming galaxies at z = 0.1 from the larger DYnamics of Newly Assembled Massive Objects (DYNAMO) sample. These galaxies, which have been characterized by high levels of star formation and large ionized gas velocity dispersions, are considered possible analogs to high-redshift clumpy disks. They were observed using the GMOS instrument in integral field spectroscopy (IFS) mode at the Gemini Observatory with high spectral resolution (R=5400, equivalent to 24 km/s at the observed wavelengths) and 6 hour exposure times in order to measure the resolved stellar kinematics via absorption lines. We also obtain higher-quality emission line kinematics than previous observations. The spatial resolution (1.2 kpc) is sufficient to show that the ionized gas in these galaxies (as traced by H-beta emission) is morphologically irregular, forming multiple giant clumps while stellar continuum light is smooth and well described by an exponential profile. Clumpy gas morphologies observed in IFS data are confirmed by complementary narrow band H-alpha imaging from the Hubble Space Telescope. Morphological differences between the stars and ionized gas are not reflected dynamically as stellar kinematics are found the be closely coupled to the kinematics of the ionized gas: both components are smoothly rotating with large velocity dispersions (~40 km/s) suggesting that the high gas dispersions are not primarily driven by star-formation feedback. In addition, the stellar population ages of these galaxies are estimated to be quite young (60-500 Myr). The large velocity dispersions measured for these young stars suggest that we are seeing the formation of thick disks and/or stellar bulges in support of recent models which produce these from clumpy galaxies at high redshift.

CHARACTERIZING ULTRAVIOLET AND INFRARED OBSERVATIONAL PROPERTIES FOR GALAXIES. II. FEATURES OF ATTENUATION LAW

Variations in the attenuation law have a significant impact on observed spectral energy distributions for galaxies. As one important observational property for galaxies at ultraviolet and infrared wavelength bands, the correlation between infrared-to-ultraviolet luminosity ratio and ultraviolet color index (or ultraviolet spectral slope), i.e., the IRX-UV relation (or IRX-beta relation), offered a widely used recipe for correcting dust attenuation in galaxies, but the usability appears to be in doubt now because of considerable dispersion in this relation found by many studies. In this paper, on the basis of spectral synthesis modeling and spatially resolved measurements of four nearby spiral galaxies, we provide an interpretation of the deviation in the IRX-UV relation with variations in the attenuation law. From both theoretical and observational viewpoints, two components in the attenuation curve, the linear background and the 2175 Angstrom bump, are suggested to be the parameters in addition to the stellar population age (addressed in the first paper of this series) in the IRX-UV function; different features in the attenuation curve are diagnosed for the galaxies in our sample. Nevertheless, it is often difficult to ascertain the attenuation law for galaxies in actual observations. Possible reasons for preventing the successful detection of the parameters in the attenuation curve are also discussed in this paper, including the degeneracy of the linear background and the 2175 Angstrom bump in observational channels, the requirement for young and dust-rich systems to study, and the difficulty in accurate estimates of dust attenuations at different wavelength bands.

Diffuse gas in galaxies sheds new light on the origin of Type Ia supernovae

We measure the strength of HeII$\lambda$4686 nebular emission in passively evolving ("retired") galaxies, aiming to constrain their populations of hot accreting white dwarfs (WDs) in the context of the single degenerate (SD) scenario of Type Ia supernovae (SNe Ia). In the SD scenario, as a WD burns hydrogen-rich material accreted from a companion star, it becomes a powerful source of ionizing UV emission. If significant populations of such sources exist in galaxies, strong emission in the recombination lines of HeII should be expected from the interstellar medium. To explore this conjecture, we select from the Sloan Digital Sky Survey ~11 500 emission line galaxies with stellar ages >1 Gyr showing no signs of AGN activity and co-add their spectra in bins of stellar population age. For the first time, we detect HeII$\lambda$4686 nebular emission in retired galaxies and find it to be significantly weaker than that expected in the SD scenario, especially in the youngest age bin (1-4 Gyr) where the SN Ia rate is the highest. Instead, the strength of the observed HeII$\lambda$4686 nebular emission is consistent with post-asymptotic giant branch stars being the sole ionizing source in all age bins. These results limit populations of accreting WDs with photospheric temperatures ($T_{\rm eff}$) in the range ~(1.5-6)$\cdot 10^5$ K to the level at which they can account for no more than ~5-10% of the observed SN Ia rate. Conversely, should all WD progenitors of SN Ia go through the phase of steady nuclear burning with $T_{\rm eff}$~(1.5-6)$\cdot 10^5$ K, they do not increase their mass by more than ~0.03 $M_\odot$ in this regime.

The effects of metallicity on the Galactic disk population of white dwarfs

It has been known for a long time that stellar metallicity plays a significant role in the determination of the ages of the different Galactic stellar populations, when main sequence evolutionary tracks are employed. Here we analyze the role that metallicity plays on the white dwarf luminosity function of the Galactic disk, which is often used to determine its age. We employ a Monte Carlo population synthesis code that accounts for the properties of the population of Galactic disk white dwarfs. Our code incorporates the most up-to-date evolutionary cooling sequences for white dwarfs with hydrogen-rich and hydrogen-deficient atmospheres for both carbon-oxygen and oxygen-neon cores. We use two different models to assess the evolution of the metallicity, one in which the adopted metallicity is constant with time, but with a moderate dispersion, and a second one in which the metallicity increases with time. We found that our theoretical results are in a very satisfactory agreement with the observational luminosity functions obtained from the Sloan Digital Sky Survey (SDSS) and from the SuperCOSMOS Sky Survey (SSS), independently of the adopted age-metallicity law. In particular, we found that the age-metallicity law has no noticeable impact in shaping the bright branch of the white dwarf luminosity function, and that the position of its cut-off is almost insensitive to the adopoted age-metallicity relationship. Because the shape of the bright branch of the white dwarf luminosity function is insensitive to the age-metallicity law, it can be safely employed to test the theoretical evolutionary sequences, while due to the limited sensitivity of the position of the drop-off to the distribution of metallicities, its location provides a robust indicator of the age of the Galactic disk.

A multiwavelength analysis of the clumpy FIR-bright sources in M33

We present a multiwavelength study of a sample of far-infrared (FIR) sources detected on the Herschel broad--band maps of the nearby galaxy M33. We perform source photometry on the FIR maps as well as mid-infrared (MIR), H$\alpha$, far-ultraviolet and integrated HI and CO line emission maps. By fitting MIR/FIR dust emission spectra, the source dust masses, temperatures and luminosities are inferred. The sources are classified based on their H$\alpha$ morphology (substructured versus not-substructured) and on whether they have a significant CO detection ($S/N>$3$\sigma$). We find that the sources have dust masses in the range 10$^2$-10$^4$~M$_\odot$ and that they present significant differences in their inferred dust/star formation/gas parameters depending on their H$\alpha$ morphology and CO detection classification. The results suggests differences in the evolutionary states or in the number of embedded HII regions between the subsamples. The source background--subtracted dust emission seems to be predominantly powered by local star formation, as indicated by a strong correlation between the dust luminosity and the dust-corrected H$\alpha$ luminosity and the fact that the extrapolated young stellar luminosity is high enough to account for the observed dust emission. Finally, we do not find a strong correlation between the dust-corrected H$\alpha$ luminosity and the dust mass of the sources, consistent with previous results on the breakdown of simple scaling relations at sub-kpc scales. However, the scatter in the relation is significantly reduced by correcting the H$\alpha$ luminosity for the age of the young stellar populations in the star--forming regions.