Source author record

M. Boquien

M. Boquien appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

91works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

91 published item(s)

preprint2025arXiv

Stochastic star formation activity of galaxies within the first billion years probed by JWST

In this work, we aim at characterizing the burstiness level of high-redshift galaxy SFHs and its evolution. We implement a stochastic SFH in CIGALE using PSD, to estimate the burstiness level of star formation in galaxies at 6<z<12. We find that SFHs with a high level of stochasticity better reproduce the SEDs of z>6 galaxies, while smoother assumptions introduce biases when applied to galaxies with bursty star-formation activity. The assumed stochasticity level of the SFH also affects the constraints on galaxies' physical properties, including the main sequence. Successively assuming different levels of burstiness, we determined the best-suited SFH for each 6<z<12 galaxy in the JADES sample from a Bayes Factor analysis. Galaxies are classified according to their level of burstiness, and the corresponding physical properties are associated to them. For massive galaxies (8.8< logM*/Msun<9.5), the fraction of bursty galaxies increases from 0.38+/-0.08 to 0.77+/-0.2 at z~6 and z~12, respectively. At all redshifts, only <20% of low-mass galaxies are classified as bursty; although, this estimate is uncertain because their faintness leads to a low S/N. For bursty galaxies, the log10(SFR10/SFR100) ratio, another indicator of bursty star formation, does not evolve with redshift, but the fraction of galaxies with a high log10(SFR10/SFR100) slightly increases from 0.28+/-0.06 to 0.38+/-0.11 between z~6 and z~9. We include additional constraints from observations on sigmaUV, the dispersion of the UV magnitude distribution, and SFE, finding a maximum of 0.72+/-0.02 mag and 0.06+/-0.01 for sigmaUV and SFE, respectively. This confirms that neither alone is responsible for the weak evolution of the UVLF at z>10. Our results add further evidence that a combination with other mechanisms is likely responsible for the high-z UVLF. The stochastic SFH module is public as part of CIGALE version 2025.1.

preprint2022arXiv

A Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE).XIV. The main sequence relation in a rich environment down to M_star ~ 10^6 Mo

Using a compilation of Halpha fluxes for 384 star forming galaxies detected during the VESTIGE survey, we study several important scaling relations for a complete sample of galaxies in a rich environment. The extraordinary sensitivity of the data allows us to sample the whole dynamic range of the Halpha luminosity function, from massive (M*~10^11 Mo) to dwarf systems (M*~10^6 Mo). This extends previous works to a dynamic range in stellar mass and star formation rate (10^-4<SFR<10 Mo yr^-1) never explored so far. The main sequence (MS) relation derived for all star forming galaxies within one virial radius of the Virgo cluster has a slope comparable to that observed in other nearby samples of isolated objects, but has a dispersion ~3 times larger. The dispersion is tightly connected to the available amount of HI gas, with gas-poor systems located far below objects of similar stellar mass but with a normal HI content. When measured on unperturbed galaxies with a normal HI gas content, the relation has a slope a=0.92, an intercept b=-1.57, and a scatter ~0.40. We compare these observational results to the prediction of models. The observed scatter in the MS relation can be reproduced only after a violent and active stripping process such as ram-pressure that removes gas from the disc and quenches star formation on short (<1 Gyr) timescales. This rules out milder processes such as starvation. This interpretation is also consistent with the position of galaxies of different star formation activity and gas content within the phase-space diagram. We also show that the star forming regions formed in the stripped material outside perturbed galaxies are located well above the MS relation drawn by unperturbed systems. These HII regions, which might be at the origin of compact sources typical in rich environments, are living a starburst phase lasting only <50 Myr, later becoming quiescent systems.

preprint2022arXiv

Comparison of the star formation of X-ray selected active galactic nuclei in eFEDS with star-forming galaxies

We use $\sim 1800$ X-ray Active Galactic Nuclei (AGN) in the eROSITA Final Equatorial-Depth Survey (eFEDS), that span over two orders of magnitude in X-ray luminosity, $\rm L_{X,2-10keV} \approx 10^{43-45}\,ergs^{-1}$, and compare their star-formation rate (SFR) relative to that of non-AGN star-forming systems, at $\rm 0.5<z<1.5$. For that purpose, we compile a reference galaxy catalogue with $\sim 17000$ sources. Both samples have the same photometric coverage, from optical to far-infrared. We construct the spectral energy distributions (SEDs) of all sources and fit them using the CIGALE code, utilizing the same templates and parametric grid for both samples. We account for the mass incompleteness and exclude quiescent systems from both datasets. These allow us to compare the SFR of the two populations in a uniform manner, minimising systematic effects. Based on our analysis, AGN at low and moderate L$_X$ ($\rm L_{X,2-10keV} < 10^{44}\,ergs^{-1}$), have SFR that is lower, or at most, equal to that of star-forming galaxies, in agreement with previous studies. The large number of luminous X-ray AGN available in the eFEDS dataset, enable us to expand our investigations at higher L$_X$ to test previous, tentative results. At $\rm L_{X,2-10keV} > 10^{44.2}\,ergs^{-1}$, the SFR of AGN appears enhanced, by $\sim 30\%$, compared to that of star-forming sources, for systems with stellar mass, $\rm 10.5 < log\,[M_*(M_\odot)] < 11.5$, confirming indications found in previous studies. The most massive sources $\rm log\,[M_*(M_\odot)] > 11.5$, present a flat SFR$_{norm}$-L$_X$ relation up to $\rm L_{X,2-10keV} \sim 10^{44.5}\,ergs^{-1}$, with SFR similar to that of star-forming galaxies. However, at higher L$_X$ ($\rm L_{X,2-10keV} \sim 10^{45}\,ergs^{-1}$), we find indications that the SFR of these massive AGN hosts may be enhanced compared to that of non-AGN systems.

preprint2022arXiv

IR characteristic emission and dust properties of star-forming galaxies at 4.5 $<$ z $<$ 6.2

The luminosity functions at z < 4 - 5 suggest that most galaxies have a relatively low stellar mass (logM_star = 10) and a low dust attenuation (A_FUV = 1.0). The physical properties of these objects are quite homogeneous. We used an approach where we combined their rest-frame far-infrared and submillimeter emissions and utilized the universe and the redshift as a spectrograph to increase the amount of information in a collective way. From a subsample of 27 ALMA-detected galaxies at z > 4.5, we built an infrared spectral energy distribution composite template. It was used to fit, with CIGALE, the 105 galaxies (detections and upper limits) in the sample from the FUV to the FIR. The derived physical parameters provide information to decipher the nature of the dust cycle and of the stellar populations in these galaxies. The derived IR composite template is consistent with the galaxies in the studied sample. A delayed star formation history with tau_main = 500 Myrs is slightly favored by the statistical analysis as compared to a delayed with a final burst or a continuous star formation history. The position of the sample in the star formation rate (SFR)- M_star diagram is consistent with previous papers. The redshift evolution of the log M_star versus A_FUV relation is in agreement with evolution in the redshift of this relation. This evolution is necessary to explain the cosmic evolution of the average dust attenuation of galaxies. Evolution is also observed in the L_dust/ L_FUV (IRX) versus UV slope beta_FUV diagram: younger galaxies have bluer beta_FUV. We modeled the shift of galaxies in the IRX versus the beta_FUV diagram with the mass-weighted age as a free parameter, and we provide an equation to make predictions.

preprint2022arXiv

Linking stellar populations to HII regions across nearby galaxies: I. Constraining pre-supernova feedback from young clusters in NGC1672

One of the fundamental factors regulating the evolution of galaxies is stellar feedback. However, we still do not have strong observational constraints on the relative importance of the different feedback mechanisms (e.g. radiation, ionised gas pressure, stellar winds) in driving HII region evolution and molecular cloud disruption. In this letter, we constrain the relative importance of the various feedback mechanisms from young massive star populations by resolving HII regions across the disk of the nearby star-forming galaxy NGC 1672. We combine measurements of ionised gas nebular lines obtained by PHANGS-MUSE, with high-resolution imaging from the HST in both the narrow-band Hα and broad-band filters. We identify a sample of 40 isolated, compact HII regions in the HST Hα image, for which we measure the sizes that were previously unresolved in seeing-limited ground-based observations. Additionally, we identify the ionisation source(s) for each HII region from catalogues produced as part of the PHANGS-HST survey. We find that the HII regions investigated are mildly dominated by thermal or wind pressure, yet their elevation above the radiation pressure is within the expected uncertainty range. We see that radiation pressure provides a substantially higher contribution to the total pressure than previously found in the literature over similar size scales. In general, we find higher pressures within more compact HII regions, which is driven by the inherent size scaling relations of each pressure term, albeit with significant scatter introduced by the variation in the stellar population properties (e.g. luminosity, mass, age, metallicity). For nearby galaxies, here we provide a promising approach that could yield the statistics required to map out how the importance of different stellar feedback mechanisms evolve over the lifetime of an HII region.

preprint2022arXiv

The ALPINE-ALMA [CII] Survey: Investigation of 10 Galaxies at $z\sim4.5$ with [OII] and [CII] Line Emission $-$ ISM Properties and [OII]$-$SFR Relation

We present $10$ main-sequence ALPINE galaxies (log($M/M_{\odot}$) = 9.2-11.1 and ${\rm SFR}=23-190\,{\rm M_{\odot}\,yr^{-1}}$) at $z\sim4.5$ with optical [OII] measurements from Keck/MOSFIRE spectroscopy and Subaru/MOIRCS narrow-band imaging. This is the largest such multi-wavelength sample at these redshifts, combining various measurements in the ultra-violet, optical, and far-infrared including [CII]$_{158{\rm μm}}$ line emission and dust continuum from ALMA and H$α$ emission from Spitzer photometry. For the first time, this unique sample allows us to analyze the relation between [OII] and total star-formation rate (SFR) and the interstellar medium (ISM) properties via [OII]/[CII] and [OII]/\halpha luminosity ratios at $z\sim4.5$. The [OII]$-$SFR relation at $z\sim4.5$ cannot be described using standard local descriptions, but is consistent with a metal-dependent relation assuming metallicities around $50\%$ solar. To explain the measured dust-corrected luminosity ratios of $L[OII]/L[CII] \sim 0.98^{+0.21}_{-0.22}$ and $L[OII]/LHa \sim -0.22^{+0.13}_{-0.15}$ for our sample, ionization parameters $\log(U)< -2$ and electron densities $\log(\rm n_e / {\rm [cm^{-3}]}) \sim 2.5-3$ are required. The former is consistent with galaxies at $z\sim2-3$, however lower than at $z>6$. The latter may be slightly higher than expected given the galaxies' specific SFR. The analysis of this pilot sample suggests that typical log($ M/M_{\odot})$ > 9 galaxies at $z\sim4.5$ to have broadly similar ISM properties as their descendants at $z\sim2$ and suggest a strong evolution of ISM properties since the Epoch of Reionization at $z>6$.

preprint2022arXiv

The ALPINE-ALMA [CII] Survey: the population of [CII]-undetected galaxies and their role in the $\mathrm{L_{[CII]}}$-SFR relation

The [CII] 158$~μ$m emission line represents so far one of the most profitable tools for the investigation of the high-redshift galaxies in the early Universe. Being one of the brightest cooling lines in the rest-frame far-infrared regime of star-forming galaxies, it has been successfully exploited as a tracer of star-formation rate (SFR) in local sources. The picture is more complex at higher redshifts, where its usability in this context is still under investigation. Recent results from the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey suggest that there is no (or weak) evolution of the L$\mathrm{_{[CII]}}$-SFR relation up to $z\sim6$ but their reliability is hampered by the presence of a large population of [CII] non-detected galaxies. In this work, we characterize the population of [CII] non-detections in ALPINE. By stacking their ALMA spectra, we obtain a signal detected at $\sim5.1σ$, resulting in a [CII] luminosity of $\mathrm{log(L_\mathrm{[CII]}}/\mathrm{L_{\odot}})$ $\sim7.8$. When combining this value with those from the [CII] detections, we find a $\mathrm{L_{[CII]}}$-SFR relation with a slope $b=1.14\pm0.11$, in agreement within the uncertainties both with the linear relation found in the local Universe, and with the previous findings from ALPINE at $z\sim5$. This suggests that the [CII] line can be considered a good tracer of star formation up to the distant Universe. Finally, we show that the galaxies of our sample that most deviate from the observed L$_\mathrm{[CII]}$-SFR relation could suffer from a less precise redshift estimation, perhaps artificially reducing their [CII] luminosity. In this respect, we claim that there is no evidence in favour of a deficit of [CII] content in high-z galaxies, in contrast with earlier studies.

preprint2021arXiv

A Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE).XI. Two dimensional H$α$ kinematics of the edge-on ram pressure stripped galaxy NGC 4330

Using the VESTIGE survey, a deep narrow-band H$α$ imaging survey of the Virgo cluster carried on at the CFHT with MegaCam, we discovered a long diffuse tail of ionised gas in the edge-on late-type galaxy NGC 4330. This peculiar feature witnesses an ongoing ram pressure stripping (RPS) event able to remove the gas in the outer disc region. Tuned hydrodynamic simulations suggest that the RPS event is occurring almost face-on, making NGC 4330 the ideal candidate to study the effects of the perturbation in the direction perpendicular to the disc plane. We present here two new independent sets of Fabry-Perot observations (R$\simeq$10000) in order to understand the effects of the RPS process on the ionised gas kinematics. Despite their limited sensitivity to the diffuse gas emission, the data allowed us to measure the velocity and the velocity dispersion fields over the galaxy disc and in several features at the edges or outside the stellar disc formed after the RPS event. We have constructed the position-velocity diagrams and the rotation curves of the galaxy using three different techniques. The data show, consistent with the hydrodynamic simulations, that the galaxy has an inner solid-body rotation up to $\sim$2.4 kpc, with non-circular streaming motions outwards the disc and in the several external features formed during the interaction of the galaxy with the surrounding intracluster medium. The data also indicate a decrease of the rotational velocity of the gas with increasing distance from the galaxy disc along the tails, suggesting a gradual but not linear loss of angular momentum in the stripped gas. Consistent with a RPS scenario, the $i$-band image shows a boxy shape at the southwest edge of the disc, where the stellar orbits might have been perturbed by the modification of the gravitational potential well of the galaxy due to the displacement of the gas in the $z$-direction.

preprint2021arXiv

A Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE).XII. Ionised gas emission in the inner regions of lenticular galaxies

As part of the VESTIGE survey, a blind narrow-band Ha+[NII] imaging survey of the Virgo cluster carried out with MegaCam at the CFHT, we discovered 8 massive lenticular galaxies with prominent ionised gas emission features in their inner (few kpc) regions. These features are either ionised gas filaments similar to those observed in cooling flows (2 gal), or thin discs with sizes 0.7<R(Ha)<2.0 kpc (6 gal), thus significantly smaller than those of the stellar disc. These discs have morphological properties similar to those of the dust seen in absorption in high-resolution HST images. Using a unique set of multifrequency data we show that while the gas located within these inner discs is photoionised by young stars, signaling ongoing star formation, the gas in the filamentary structures is shock-ionised. These discs have a star formation surface brightness similar to those observed in late-type galaxies. Because of their reduced size, however, these lenticular galaxies are located below the main sequence of unperturbed or cluster star-forming systems. By comparing the dust masses measured from absorption maps in optical images, from the Balmer decrement, or estimated by fitting the UV-to-far-IR spectral energy distribution of the target galaxies, we confirm that those derived from optical attenuation maps are heavily underestimated because of geometrical effects due to the relative distribution of the absorbing dust and the emitting stars. We have also shown that these galaxies have gas-to-dust ratios of G/D~80, and that the star formation within these discs follows the Schmidt relation, albeit with an efficiency reduced by a factor of ~ 2.5. Using our unique set of multifrequency data, we discuss the possible origin of the ionised gas in these objects, which suggests multiple and complex formation scenarios for massive lenticular galaxies in clusters.

preprint2020arXiv

A Cautionary Tale of Attenuation in Star Forming Regions

The attenuation of light in star forming galaxies is correlated with a multitude of physical parameters including star formation rate, metallicity and total dust content. This variation in attenuation is even more prevalent on the kiloparsec scale, which is relevant to many current spectroscopic integral field unit surveys. To understand the cause of this variation, we present and analyse \textit{Swift}/UVOT near-UV (NUV) images and SDSS/MaNGA emission-line maps of 29 nearby ($z<0.084$) star forming galaxies. We resolve kiloparsec-sized star forming regions within the galaxies and compare their optical nebular attenuation (i.e., the Balmer emission line optical depth, $τ^l_B\equivτ_{\textrm{H}β}-τ_{\textrm{H}α}$) and NUV stellar continuum attenuation (via the NUV power-law index, $β$) to the attenuation law described by Battisti et al. The data agree with that model, albeit with significant scatter. We explore the dependence of the scatter of the $β$-$τ^l_B$ measurements from the star forming regions on different physical parameters, including distance from the nucleus, star formation rate and total dust content. Finally, we compare the measured $τ^l_B$ and $β$ between the individual star forming regions and the integrated galaxy light. We find a strong variation in $β$ between the kiloparsec scale and the larger galaxy scale not seen in $τ^l_B$. We conclude that the sight-line dependence of UV attenuation and the reddening of $β$ due to the light from older stellar populations could contribute to the $β$-$τ^l_B$ discrepancy.

preprint2020arXiv

A hyper luminous starburst at z=4.72 magnified by a lensing galaxy pair at z=1.48

[Abridged] We discovered in the Herschel Reference Survey an extremely bright IR source with $S_{500}$~120mJy (Red Virgo 4 - RV4). Based on IRAM/EMIR and IRAM/NOEMA detections of the CO(5-4), CO(4-3), and [CI] lines, RV4 is located at z=4.724, yielding a total observed L$_{IR}$ of 1.1+/-0.6x0$^{14}$L$_{\odot}$. At the position of the Herschel emission, three blobs are detected with the VLA at 10cm. The CO(5-4) line detection of each blob confirms that they are at the same redshift with the same line width, indicating that they are multiple images of the same source. In Spitzer and deep optical observations, two sources, High-z Lens 1 (HL1) West and HL1 East, are detected at the center of the three VLA/NOEMA blobs. These two sources are placed at z=1.48 with XSHOOTER spectra, suggesting that they could be merging and gravitationally lensing the emission of RV4. HL1 is the second most distant lens known to date in strong lensing systems. The Einstein radius of the lensing system is 2.2"+/-0.2 (20kpc). The high redshift of HL1 and the large Einstein radius are highly unusual for a strong lensing system. We present the ISM properties of the background source RV4. Different estimates of the gas depletion time yield low values suggesting that RV4 is a SB galaxy. Among all high-z SMGs, this source exhibits one of the lowest L$_{[CI]}$ to L$_{IR}$ ratios, 3.2+/-0.9x10$^{-6}$, suggesting an extremely short gas tdepl of only 14+/-5Myr. It also shows a relatively high L$_{[CI]}$ to L$_{CO(4-3)}$ ratio (0.7+/-0.2) and low L$_{CO(5-4)}$ to L$_{IR}$ ratio (only ~50% of the value expected for normal galaxies) hinting a low density of gas. Finally, we discuss that the short tdepl of RV4 can be explained by either a very high SFE, which is difficult to reconcile with major mergers simulations of high-z galaxies, or a rapid decrease of SF, which would bias the estimate of tdepl toward low value.

preprint2020arXiv

A Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE). VIII. Modeling ram pressure stripping of diffuse gas in the Virgo cluster spiral galaxy NGC 4330

NGC 4330 is one of the Virgo galaxies whose UV distribution shows a tail structure. An associated tail structure is also observed in the HI and H$α$ emission distributions. Previous dynamical modeling showed that the galaxy is approaching the cluster center and is therefore undergoing increasing ram pressure stripping. Recent stellar population fitting of deep optical spectra together with multiband photometry lead to the determination of the time when star formation was quenched in the galactic disk. We introduce a new version of the dynamical model that includes the diffuse ionized gas and aim to reproduce the HI, H$α$, UV distributions together with the star formation histories of the outer gas-free parts of the galactic disk. The results of 50 simulations with five different Lorentzian temporal ram-pressure profiles and five different delays between the simulation onset and peak ram pressure are presented. The inclusion of diffuse gas stripping changes significantly the HI, UV, and H$α$ emission distributions. The simulations with diffuse gas stripping naturally lead to vertical low surface density filaments in the downwind region of the galactic disk. These filaments occur less frequently in the simulations without diffuse gas stripping. The simulations with diffuse gas stripping lead to better joint fits to the SEDs and optical spectra. The HI, NUV, and H$α$ morphologies of the model snapshots which best reproduce the SEDs and optical spectra are sufficiently different to permit a selection of best-fit models. We conclude that the inclusion of diffuse gas stripping significantly improves the resemblance between the model and observations. Our preferred model yields a time to peak ram pressure of 140 Myr in the future. The spatial coincidence of the radio continuum and diffuse H$α$ tails suggests that both gas phases are stripped together.

preprint2020arXiv

A Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE).VI. Environmental quenching on HII region scales

The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band Halpha+[NII] imaging survey of the Virgo cluster carried out with MegaCam at the Canada-French-Hawaii telescope (CFHT). We use a new set of data extracted from VESTIGE to study the impact of the hostile cluster environment on the star formation process down to the scale of HII regions (~ 50 pc). HII regions are identified and their parameters measured using the HIIphot code on a sample of 114 late-type galaxies spanning a wide range in morphological type (Sa-Sd, Im, BCD), stellar mass (10^6.5 <= M_star <= 10^11 Mo), and star formation activity (10^-3 <= SFR <= 10 Mo yr^-1). Owing to the exquisite average resolution of the VESTIGE data (0.65 arcsec), we detect 11302 HII regions with an Halpha luminosity L(Halpha) >= 10^37 erg s^-1. We show that the typical number of HII regions in gas-stripped objects is significantly lower than in healthy late-types of similar stellar mass. We also show that in these gas-stripped galaxies the number of HII regions significantly drops outside the effective radius, suggesting that the quenching process occurs outside-in, in agreement with other multifrequency observations. These new results consistently confirm that the main mechanism responsible for the decrease of the star formation activity observed in cluster galaxies is ram pressure, allowing us to discard other milder processes such as starvation or strangulation unable to reproduce the observed radially truncated profiles.

preprint2020arXiv

A Virgo Environmental Survey Tracing Ionised Gas Emission. VESTIGE VIII. Bridging the cluster-ICM-galaxy evolution at small scales

We measure FIR emission from tails of stripped dust following the ionised and atomic gas components in galaxies undergoing ram pressure stripping. We study the dust-to-gas relative distribution and mass ratio in the stripped interstellar medium and relate them to those of the intra-cluster medium, thus linking the cluster-ICM-galaxy evolution at small-scales. The galaxy sample consists of three Scd Virgo galaxies with stellar masses in the range $10^9\lesssim \mathrm{M_{*}} \lesssim 10^{10}\, \mathrm{M_{\odot}}$, and within 1 Mpc from the cluster centre, namely NGC 4330, NGC 4522, and NGC 4654. Through the analysis of VESTIGE H$α$, $Herschel$ SPIRE far-infrared, and VIVA HI data, we trace the spatial distribution of the tails and infer the dust and gas masses from the measured far-infrared 250 $μ$m and HI flux densities. Dust-to-gas mass ratios (DGRs) in the tails are analysed as a function of the galaxy mass, metallicity, and dust temperature. Along the stripped component, the dust distribution closely follows the HI and H$α$ emitting gas, all extending beyond the optical disc. In these regions, the DGRs are $2.0\pm0.6\times10^{-3}$, $0.7\pm0.1\times10^{-3}$, and $0.4\pm0.03\times10^{-3}$, for NGC 4330, NGC 4522, and NGC 4654, respectively, i.e. up to a factor of 15 less than the values measured in the main body of nearby galaxies. We also find a negative trend in the DGR as a function of the metallicity that can be explained in terms of a dust component more centrally concentrated in more metal-rich systems. Together with the finding that the stripped dust is cold, $T_{d} \lesssim 25\, K$, our results support an outside-in stripping scenario of the galaxy interstellar medium. This study shows that ram pressure stripping is a key mechanism in the building up of the Virgo intra-cluster component injecting dust grains into the ICM, thus contributing to its metal enrichment.

preprint2020arXiv

Differences and similarities of stellar populations in LAEs and LBGs at $z\sim$ 3.4 - 6.8

The differences between the inherent stellar populations (SPs) of LAEs and LBGs are a key factor in understanding early galaxy formation and evolution. We have run a set of SP burst-like models for a sample of 1,558 sources at $3.4<z<6.8$ from the Survey for High-$z$ Absorption Red and Dead Sources (SHARDS) over the GOODS-N field. This work focuses on the differences between the three different observational subfamilies of our sample: LAE-LBGs, no-Ly$α$ LBGs and pure LAEs. Single and double SP synthetic spectra were used to model the SEDs, adopting a Bayesian information criterion to analyse under which situations a second SP is required. We find that the sources are well modelled using a single SP in $\sim79\%$ of the cases. The best models suggest that pure LAEs are typically young low mass galaxies ($t\sim26^{+41}_{-25}$ Myr; $M_{\mathrm{star}}\sim5.6^{+12.0}_{-5.5}\times10^{8}\ M_{\odot}$), undergoing one of their first bursts of star formation. On the other hand, no-Ly$α$ LBGs require older SPs ($t\sim71\pm12$ Myr), and they are substantially more massive ($M_{\mathrm{star}}\sim3.5\pm1.1\times10^{9}\ M_{\odot}$). LAE-LBGs appear as the subgroup that more frequently needs the addition of a second SP, representing an old and massive galaxy caught in a strong recent star-forming episode. The relative number of sources found from each subfamily at each $z$ supports an evolutionary scenario from pure LAEs and single SP LAE-LBGs to more massive LBGs. Stellar Mass Functions are also derived, finding an increase of $M^{*}$ with cosmic time and a possible steepening of the low mass slope from $z\sim6$ to $z\sim5$ with no significant change to $z\sim4$. Additionally, we have derived the SFR-$M_{\mathrm{star}}$ relation, finding a $\mathrm{SFR}\propto M_{\mathrm{star}}^β$ behaviour with negligible evolution from $z\sim4$ to $z\sim6$.

preprint2020arXiv

SDSS IV MaNGA -- Metallicity and ionisation parameter in local star-forming galaxies from Bayesian fitting to photoionisation models

We measured gas-phase metallicity, ionisation parameter and dust extinction for 1795 representative local star-forming galaxies using integral field spectroscopy from the SDSS-IV MaNGA survey. We self-consistently derive these quantities by comparing observed line fluxes with photoionisation models using a Bayesian framework. We also present the first comprehensive study of the [SIII]$λλ$9069,9532 nebular lines, which have long been predicted to be ideal tracers of the ionisation parameter. Unfortunately, we find that current photoionisation models substantially over-predict [SIII] lines intensity, while broadly reproducing other observed optical line ratios. We discuss how to nonetheless make use of the information provided by [SIII] lines by setting a prior on the ionisation parameter. Following this approach, we derive spatially-resolved maps and radial profiles of metallicity and ionisation parameter. The metallicity radial profiles are comparable with previous works, with metallicity declining toward the outer parts and a flattening in the central regions, in agreement with infall models of galaxy formation, that predict that spiral discs build up through accretion of material, which leads to an inside-out growth. On the other hand, ionisation parameter radial profiles are flat for low-mass galaxies, while their slope becomes positive as galaxy mass increases. However, the ionisation parameter maps we obtain are clumpy, especially for low-mass galaxies. Ionisation parameter is tightly correlated with the H$α$ equivalent width [EW(H$α$)], following a nearly universal relation, which we attribute to the change of the spectral shape of ionising sources due to ageing of HII regions. We derive a positive correlation between ionisation parameter and metallicity at fixed EW(H$α$), in disagreement with previous theoretical works expecting an anti-correlation.

preprint2020arXiv

The ALPINE-ALMA [CII] Survey: Multi-Wavelength Ancillary Data and Basic Physical Measurements

We present the ancillary data and basic physical measurements for the galaxies in the ALMA Large Program to Investigate C+ at Early Times (ALPINE) survey - the first large multi-wavelength survey which aims at characterizing the gas and dust properties of 118 main-sequence galaxies at redshifts 4.4 < z < 5.9 via the measurement of [CII]-emission at 158 micro-meter (64% at >3.5$σ$) and the surrounding far-infrared (FIR) continuum in conjunction with a wealth of optical and near-infrared data. We outline in detail the spectroscopic data and selection of the galaxies as well as the ground- and space-based imaging products. In addition, we provide several basic measurements including stellar masses, star formation rates (SFR), rest-frame ultra-violet (UV) luminosities, UV continuum slopes ($β$), and absorption line redshifts, as well as H-alpha emission derived from Spitzer colors. We find that the ALPINE sample is representative of the 4 < z < 6 galaxy population selected by photometric methods and only slightly biased towards bluer colors ($Δβ$ ~ 0.2). Using [CII] as tracer of the systemic redshift (confirmed for one galaxy at z=4.5 out of 118 for which we obtained optical [OII]$λ$3727A emission), we confirm red shifted Ly-alpha emission and blue shifted absorption lines similar to findings at lower redshifts. By stacking the rest-frame UV spectra in the [CII] rest-frame we find that the absorption lines in galaxies with high specific SFR are more blue shifted, which could be indicative of stronger winds and outflows.

preprint2020arXiv

The ALPINE-ALMA [CII] Survey: On the nature of an extremely obscured serendipitous galaxy

We report the serendipitous discovery of a bright galaxy (Gal-A) observed as part of the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). While this galaxy is detected both in line and continuum emission in ALMA Band 7, it is completely dark in UV/optical filters and only presents a marginal detection in the UltraVISTA Ks band. We discuss the nature of the observed ALMA line, i.e. whether the emission comes from [CII] at z~4.6, or from high-J CO transitions at z~2.2. In the first case we find a [CII]-to-FIR luminosity ratio of log(L_[CII]/L_FIR)=-2.5, consistent with the average value for local star-forming galaxies (SFGs); in the second case, instead, the source would lie outside of the empirical relations between L_CO and L_FIR found in the literature. At both redshifts, we derive the star-formation rate (SFR) from the ALMA continuum, and the stellar mass (M*) by using stellar population synthesis models as input for LePHARE spectral energy distribution (SED) fitting. Exploiting our results, we believe that Gal-A is a Main-Sequence (MS), dusty SFG at z=4.6 (i.e. [CII] emitter) with log(SFR/[M/yr])~1.4 and log(M*/M)~9.7. This work underlines the crucial role of the ALPINE survey in making a census of this class of objects, in order to unveil their contribution to the global star-formation rate density (SFRD) of the Universe at the end of the Reionisation epoch.

preprint2020arXiv

X-CIGALE: fitting AGN/galaxy SEDs from X-ray to infrared

CIGALE is a powerful multiwavelength spectral energy distribution (SED) fitting code for extragalactic studies. However, the current version of CIGALE is not able to fit X-ray data, which often provide unique insights into AGN intrinsic power. We develop a new X-ray module for CIGALE, allowing it to fit SEDs from the X-ray to infrared (IR). We also improve the AGN fitting of CIGALE from UV-to-IR wavelengths. We implement a modern clumpy two-phase torus model, SKIRTOR. To account for moderately extincted type 1 AGNs, we implement polar-dust extinction. We publicly release the source code (named X-CIGALE). We test X-CIGALE with X-ray detected AGNs in SDSS, COSMOS, and AKARI-NEP. The fitting quality (as indicated by reduced $χ^2$) is good in general, indicating that X-CIGALE is capable of modelling the observed SED from X-ray to IR. We discuss constrainability and degeneracy of model parameters in the fitting of AKARI-NEP, for which excellent mid-IR photometric coverage is available. We also test fitting a sample of AKARI-NEP galaxies for which only X-ray upper limits are available from Chandra observations, and find that the upper limit can effectively constrain the AGN SED contribution for some systems. Finally, using X-CIGALE, we assess the ability of Athena to constrain the AGN activity in future extragalactic studies.

preprint2020arXiv

X-ray flux in the SED modelling: An application of X-CIGALE in the XMM-XXL field

X-CIGALE, built upon the spectral energy distribution (SED) code of CIGALE, implements important new features: The code accounts for obscuring material in the polars of the AGN and has the ability to fit X-ray fluxes. In this work, we use ~2500 spectroscopic, X-ray AGN from the XMM-XXL-North field and examine the improvements the new features bring in the SED modelling analysis. Based on our results, X-CIGALE successfully connects the X-ray with the UV luminosity in the whole range spanned by our sample (log Lx(2-10 keV) = (42 - 46) erg/s). The addition of the new features globally improves the efficiency of X-CIGALE in the estimation and characterization of the AGN component. The classification into type 1 and type 2 based on their inclination angle is improved, especially at redshifts lower than 1. The statistical significance of the AGN fraction, fracAGN, measurements is increased, in particular for luminous X-ray sources (log LX > 45 erg/s). These conclusions hold under the condition that (mid-) IR photometry is available in the SED fitting process. The addition of polar dust increases the AGN fraction and the efficiency of the SED decomposition to detect AGN among X-ray selected sources. X-CIGALE estimates a strong AGN (fracAGN > 0.3) in more than 90% of the infrared selected AGN and 75% of X-ray detected AGN not selected by IR colour criteria. The latter drops to ~50% when polar dust is not included. The ability of X-CIGALE to include X-ray information in the SED fitting process can be instrumental in the optimal exploitation of the wealth of data that current (eROSITA) and future (ATHENA) missions will provide us.

preprint2019arXiv

Modeling Dust and Starlight in Galaxies Observed by Spitzer and Herschel: The KINGFISH Sample

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.

preprint2019arXiv

Spatial segregation of massive clusters in dwarf galaxies

The relative average minimum projected separations of star clusters in the Legacy ExtraGalactic UV Survey (LEGUS) and in tidal dwarfs around the interacting galaxy NGC 5291 are determined as a function of cluster mass to look for cluster-cluster mass segregation. Class 2 and 3 LEGUS clusters, which have a more irregular internal structure than the compact and symmetric class 1 clusters, are found to be mass segregated in low mass galaxies, which means that the more massive clusters are systematically bunched together compared to the lower mass clusters. This mass segregation is not present in high-mass galaxies nor for class 1 clusters. We consider possible causes for this segregation including differences in cluster formation and scattering in the shallow gravitational potentials of low mass galaxies.

preprint2019arXiv

The ALPINE-ALMA [CII] survey: Star-formation-driven outflows and circumgalactic enrichment in the early Universe

We study the efficiency of galactic feedback in the early Universe by stacking the [C II] 158 um emission in a large sample of normal star-forming galaxies at 4 < z < 6 from the ALMA Large Program to INvestigate [C II] at Early times (ALPINE) survey. Searching for typical signatures of outflows in the high-velocity tails of the stacked [C II] profile, we observe (i) deviations from a single-component Gaussian model in the combined residuals and (ii) broad emission in the stacked [C II] spectrum, with velocities of |v|<~ 500 km/s. The significance of these features increases when stacking the subset of galaxies with star formation rates (SFRs) higher than the median (SFRmed = 25 Msun/yr), thus confirming their star-formation-driven nature. The estimated mass outflow rates are comparable to the SFRs, yielding mass-loading factors of the order of unity (similarly to local star-forming galaxies), suggesting that star-formation-driven feedback may play a lesser role in quenching galaxies at z > 4. From the stacking analysis of the datacubes, we find that the combined [C II] core emission (|v|< 200 km/s) of the higher-SFR galaxies is extended on physical sizes of ~ 30 kpc (diameter scale), well beyond the analogous [C II] core emission of lower-SFR galaxies and the stacked far-infrared continuum. The detection of such extended metal-enriched gas, likely tracing circumgalactic gas enriched by past outflows, corroborates previous similar studies, confirming that baryon cycle and gas exchanges with the circumgalactic medium are at work in normal star-forming galaxies already at early epochs.

preprint2016arXiv

Comparing [CII], HI, and CO dynamics of nearby galaxies

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 $μ$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.

preprint2016arXiv

Dust Properties in HII Regions in M33

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.

preprint2016arXiv

Millimeter and Submillimeter Excess Emission in M33 revealed by Planck and LABOCA

Previous studies have shown the existence of an excess of emission at submillimeter (submm) and millimeter (mm) wavelengths in the spectral energy distribution (SED) of many low-metallicity galaxies. The goal of the present study is to model separately the emission from the star forming (SF) component and the emission from the diffuse interstellar medium (ISM) in the nearby spiral galaxy M33. We decomposed the observed SED of M33 into its SF and diffuse components. Mid-infrared (MIR) and far-infrared (FIR) fluxes were extracted from Spitzer and Herschel data. At submm and mm wavelengths, we used ground-based observations from APEX to measure the emission from the SF component and data from the Planck space telescope to estimate the diffuse emission. Both components were separately fitted using radiation transfer models based on standard dust properties and a realistic geometry. Both modeled SEDs were combined to build the global SED of M33. In addition, the radiation field necessary to power the dust emission in our modeling was compared with observations from GALEX, Sloan, and Spitzer. Our modeling is able to reproduce the observations at MIR and FIR wavelengths, but we found a strong excess of emission at submm and mm wavelengths, where the model expectations severely underestimate the LABOCA and Planck fluxes. We also found that the ultraviolet (UV) radiation escaping the galaxy is 70% higher than the model predictions. We determined a gas-to-dust mass ratio Gdust~100, significantly lower than the value expected from the sub-solar metallicity of M33. We discussed different hypotheses to explain the discrepancies found in our study (i.e., excess of emission at submm and mm wavelengths, deficit of UV attenuation, and abnormally low value of Gdust), concluding that different dust properties in M33 is the most plausible explanation.

preprint2016arXiv

Molecular gas and star formation in the Tidal Dwarf Galaxy VCC 2062

The physical mechanisms driving star formation (SF) in galaxies are still not fully understood. Tidal dwarf galaxies (TDGs), made of gas ejected during galaxy interactions, seem to be devoid of dark matter and have a near-solar metallicity. The latter makes it possible to study molecular gas and its link to SF using standard tracers (CO, dust) in a peculiar environment. We present a detailed study of a nearby TDG in the Virgo Cluster, VCC 2062, using new high-resolution CO(1--0) data from the Plateau de Bure, deep optical imaging from the Next Generation Virgo Cluster Survey (NGVS), and complementary multiwavelength data. Until now, there was some doubt whether VCC 2062 was a true TDG, but the new deep optical images from the NGVS reveal a stellar bridge between VCC 2062 and its parent galaxy, NGC 4694, which is clear proof of its tidal origin. Several high-resolution tracers (\halpha, UV, 8~\mi, and 24~\mi) of the star formation rate (SFR) are compared to the molecular gas distribution as traced by the CO(1-0). Coupled with the SFR tracers, the NGVS data are used with the CIGALE code to model the stellar populations throughout VCC 2062, yielding a declining SFR in the recent past, consistent with the low \halpha/UV ratio, and a high burst strength. HI emission covers VCC 2062, whereas the CO is concentrated near the HI maxima. The CO peaks correspond to two very distinct regions: one with moderate SF to the NE and one with only slightly weaker CO emission but with nearly no SF. Even where SF is clearly present, the SFR is below the value expected from the surface density of the molecular and the total gas as compared to spiral galaxies and other TDGs. After discussing different possible explanations, we conclude that the low surface brightness is a crucial parameter to understand the low SFR.

preprint2016arXiv

The quenching of the star formation activity in cluster galaxies

We study the star formation quenching mechanism in cluster galaxies by fitting the SED of the Herschel Reference Survey, a complete volume-limited K-band-selected sample of nearby galaxies including objects in different density regions, from the core of the Virgo cluster to the general field. The SED are fitted using the CIGALE SED modelling code. The truncated activity of cluster galaxies is parametrised using a specific SFH with 2 free parameters, the quenching age QA and the quenching factor QF. These 2 parameters are crucial for the identification of the quenching mechanism which acts on long timescales if starvation while rapid and efficient if ram pressure. To be sensitive to an abrupt and recent variation of the star formation activity, we combine in a new way 20 UV to FIR photometric bands with 3 age-sensitive Balmer line absorption indices extracted from available medium-resolution integrated spectroscopy and with Halpha narrow band imaging data. The use of a truncated SFH significantly increases the quality of the fit in those objects whose atomic gas content has been removed during the interaction with the hostile cluster environment. The typical QA of the perturbed late-type galaxies is QA < 300 Myr whenever the activity of star formation is reduced by 50% < QF <= 80% and QA < 500 Myr for QF > 80%, while that of the quiescent early-types is QA ~ 1-3 Gyr. The fraction of late-types with a star formation activity reduced by QF > 80% and with an HI-deficiency parameter HI-def > 0.4 drops by a factor of ~ 5 from the inner half virial radius of the Virgo cluster, where the hot diffuse X-ray emitting gas of the cluster is located, to the outer regions. The efficient quenching of the star formation activity observed in Virgo suggests that the dominant stripping process is ram pressure. We discuss the implication of this result in the cosmological context of galaxy evolution.

preprint2016arXiv

The selective effect of environment on the atomic and molecular gas-to-dust ratio of nearby galaxies in the Herschel Reference Survey

We combine dust, atomic (HI) and molecular (H$_{2}$) hydrogen mass measurements for 176 galaxies in the Herschel Reference Survey to investigate the effect of environment on the gas-to-dust mass ($M_{\rm gas}/M_{\rm dust}$) ratio of nearby galaxies. We find that, at fixed stellar mass, the average $M_{\rm gas}/M_{\rm dust}$ ratio varies by no more than a factor of $\sim$2 when moving from field to cluster galaxies, with Virgo galaxies being slightly more dust rich (per unit of gas) than isolated systems. Remarkably, once the molecular and atomic hydrogen phases are investigated separately, we find that \hi-deficient galaxies have at the same time lower $M_{\rm HI}/M_{\rm dust}$ ratio but higher $M_{\rm H_{2}}/M_{\rm dust}$ ratio than \hi-normal systems. In other words, they are poorer in atomic but richer in molecular hydrogen if normalized to their dust content. By comparing our findings with the predictions of theoretical models, we show that the opposite behavior observed in the $M_{\rm HI}/M_{\rm dust}$ and $M_{\rm H_{2}}/M_{\rm dust}$ ratios is fully consistent with outside-in stripping of the interstellar medium (ISM), and is simply a consequence of the different distribution of dust, \hi\ and H$_{2}$ across the disk. Our results demonstrate that the small environmental variations in the total $M_{\rm gas}/M_{\rm dust}$ ratio, as well as in the gas-phase metallicity, do not automatically imply that environmental mechanisms are not able to affect the dust and metal content of the ISM in galaxies.

preprint2016arXiv

Towards universal hybrid star formation rate estimators

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.

preprint2015arXiv

Cold gas properties of the Herschel Reference Survey. I. 12CO(1-0) and HI data

We present new 12CO(1-0) observations of 59 late-type galaxies belonging to the Herschel Reference Survey (HRS), a complete K-band-selected, volume-limited (15 <D< 25 Mpc) sample of nearby galaxies spanning a wide range in morphological type and luminosity. We studied different recipes to correct single-beam observations of nearby galaxies of different sizes and inclinations for aperture effects. This was done by comparing single-beam and multiple-beam observations along the major axis, which were corrected for aperture effects using different empirical or analytical prescriptions, to integrated maps of several nearby galaxies, including edge-on systems observed by different surveys. The resulting recipe is an analytical function determined by assuming that late-type galaxies are 3D exponentially declining discs with a characteristic scale length r_CO} = 0.2 r_24.5, where r_24.5 is the optical, g- (or B-) band isophotal radius at the 24.5 mag arcsec^-2 (25 mag arcsec^-2), as well as a scale height z_CO = 1/100 r_24.5. Our new CO data are then combined with those available in the literature to produce the most updated catalogue of CO observations for the HRS, now including 225 out of the 322 galaxies of the complete sample. The 3D exponential disc integration is applied to all the galaxies of the sample to measure their total CO fluxes, which are later transformed into molecular gas masses using a constant and a luminosity-dependent X_CO conversion factor. We also collect HI data for 315 HRS galaxies from the literature and present it in a homogenised form.

preprint2015arXiv

Gas dynamics in tidal dwarf galaxies: disc formation at z=0

Tidal dwarf galaxies (TDGs) are recycled objects that form within the collisional debris of interacting/merging galaxies. They are expected to be devoid of non-baryonic dark matter, since they can form only from dissipative material ejected from the discs of the progenitor galaxies. We investigate the gas dynamics in a sample of six bona-fide TDGs around three interacting and post-interacting systems: NGC 4694, NGC 5291, and NGC 7252 ("Atoms for Peace"). For NGC 4694 and NGC 5291 we analyse existing HI data from the Very Large Array (VLA), while for NGC 7252 we present new HI observations from the Jansky VLA together with long-slit and integral-field optical spectroscopy. For all six TDGs, the HI emission can be described by rotating disc models. These HI discs, however, have undergone less than a full rotation since the time of the interaction/merger event, raising the question of whether they are in dynamical equilibrium. Assuming that these discs are in equilibrium, the inferred dynamical masses are consistent with the observed baryonic masses, implying that TDGs are devoid of dark matter. This puts constraints on putative "dark discs" (either baryonic or non-baryonic) in the progenitor galaxies. Moreover, TDGs seem to systematically deviate from the baryonic Tully-Fisher relation. These results provide a challenging test for alternative theories like MOND.

preprint2015arXiv

Measuring star formation with resolved observations: the test case of M33

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.

preprint2015arXiv

NGC4370: a case study for testing our ability to infer dust distribution and mass in nearby galaxies

A fraction of the early-type galaxy population hosts a prominent dust lane. Methods to quantify the dust content of these systems based on optical imaging data usually yield dust masses which are an order of magnitude lower than dust masses derived from the observed FIR emission. High-quality optical data from the Next Generation Virgo cluster Survey (NGVS) and FIR/submm observations from the Herschel Virgo Cluster Survey (HeViCS) allow us to revisit previous methods to determine the dust content in galaxies and explore new ones. We aim to derive the dust mass in NGC 4370 from both optical and FIR data, and investigate the need to invoke a putative diffuse dust component. We create color and attenuation maps, which are converted to approximate dust mass maps based on simple dust geometries. Dust masses are also derived from SED fits to FIR/submm observations. Finally, inverse radiative transfer fitting is performed to investigate more complex dust geometries. The empirical methods applied to the optical data yield lower limits of 3.4e5 solar masses, an order of magnitude below the total dust masses derived from SED fitting. In contrast, radiative transfer models yield dust masses which are slightly lower, but fully consistent with the FIR-derived mass. Dust is more likely to be distributed in a ring around the centre of NGC 4370 as opposed to an exponential disc or a simple foreground screen. Moreover, using inverse radiative transfer fitting, we are able to constrain most of the parameters describing these geometries. The resulting dust masses are high enough to account for the dust observed at FIR/submm wavelengths, so that no diffuse dust component needs to be invoked. We furthermore caution for the interpretation of dust masses and optical depths based on optical data alone, using overly simplistic star-dust geometries and the neglect of scattering effects. [ABRIDGED]

preprint2015arXiv

The bolometric and UV attenuation in normal spiral galaxies of the Herschel Reference Survey

The dust in nearby galaxies absorbs a fraction of the UV-optical-near-infrared radiation produced by stars. This energy is consequently re-emitted in the infrared. We investigate the portion of the stellar radiation absorbed by spiral galaxies from the HRS by modelling their UV-to-submillimetre spectral energy distributions. Our models provide an attenuated and intrinsic SED from which we find that on average 32 % of all starlight is absorbed by dust. We define the UV heating fraction as the percentage of dust luminosity that comes from absorbed UV photons and find that this is 56 %, on average. This percentage varies with morphological type, with later types having significantly higher UV heating fractions. We find a strong correlation between the UV heating fraction and specific star formation rate and provide a power-law fit. Our models allow us to revisit the IRX-AFUV relations, and derive these quantities directly within a self-consistent framework. We calibrate this relation for different bins of NUV-r colour and provide simple relations to relate these parameters. We investigated the robustness of our method and we conclude that the derived parameters are reliable within the uncertainties which are inherent to the adopted SED model. This calls for a deeper investigation on how well extinction and attenuation can be determined through panchromatic SED modelling.

preprint2015arXiv

The imprint of rapid star formation quenching on the spectral energy distributions of galaxies

[Abridged] In high density environment, the gas content of galaxies is stripped, leading to a rapid quenching of their star formation activity. This dramatic environmental effect is generally not taken into account in the SFHs usually assumed to perform spectral energy distribution (SED) fitting of these galaxies, yielding to a poor fit of their stellar emission and, consequently, a biased estimate of the SFR. We aim at reproducing the SFH of galaxies that underwent a rapid star formation quenching using a truncated delayed SFH that we implemented in the SED fitting code CIGALE. We show that the ratio between the instantaneous SFR and the SFR just before the quenching ($r_{SFR}$) is well constrained as long as rest frame UV data are available. This SED modelling is applied to the Herschel Reference Survey (HRS) containing isolated galaxies and sources falling in the dense environment of the Virgo cluster. The latter are HI-deficient due to ram pressure stripping. We show that the truncated delayed SFH successfully reproduces their SED while typical SFH assumptions fail. A good correlation is found between $r_{SFR}$ and HI-def, the parameter quantifying the gas deficiency of cluster galaxies, meaning that SED fitting results can be used to provide a tentative estimate of the gas deficiency of galaxies for which HI observations are not available. The HRS galaxies are placed on the SFR-$M_*$ diagram showing that the HI-deficient sources lie in the quiescent region confirming previous studies. Using the $r_{SFR}$ parameter, we derive the SFR of these sources before quenching and show that they were previously on the main sequence relation. We show that the $r_{SFR}$ parameter is also well recovered for deeply obscured high redshift sources, as well as in absence of IR data. SED fitting is thus a powerful tool to identify galaxies that underwent a rapid star formation quenching.

preprint2014arXiv

[CII] 158 $μ$m Emission as a Star Formation Tracer

The [CII] 157.74 $μ$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 $Σ_{\rm [CII]}-Σ_{\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 $Σ_{\rm[CII]}-Σ_{\rm SFR}$ correlation is valid over almost 5 orders of magnitude in $Σ_{\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.

preprint2014arXiv

A resolved analysis of cold dust and gas in the nearby edge-on spiral NGC 891

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

preprint2014arXiv

Cold gas properties of the Herschel Reference Survey. II. Molecular and total gas scaling relations

We study the properties of the cold gas component of the interstellar medium of the Herschel Reference Survey, a complete volume-limited (15<D<25 Mpc), K-band-selected sample of galaxies spanning a wide range in morphological type (from E to Im) and stellar mass (10^9<M*<10^11 Mo). The multifrequency data in our hands are used to trace the molecular gas mass distribution and the main scaling relations of the sample, which put strong constraints on galaxy formation simulations. We extend the main scaling relations concerning the total and the molecular gas component determined for massive galaxies (M* > 10^10 Mo) from the COLD GASS survey down to stellar masses M* ~ 10^9 Mo. As scaling variables we use M*, the stellar surface density mu*, the specific star formation rate SSFR, and the metallicity of the target galaxies. By comparing molecular gas masses determined using a constant or a luminosity dependent conversion factor, we estimate the robustness of these scaling relations on the very uncertain assumptions used to transform CO line intensities into molecular gas masses. The molecular gas distribution of a K-band-selected sample is different from that of a far-infrared-selected sample since it includes a significantly smaller number of objects with M(H2) < 6 10^9 Mo. In spiral galaxies the molecular gas phase is only 25-30% of the atomic gas. The analysis also indicates that the slope of the main scaling relations depends on the adopted conversion factor. Among the sampled relations, all those concerning M(gas)/M* are statistically significant and show little variation with X_CO. We observe a significant correlation between M(H2)/M* and SSFR, M(H2)/M(HI) and mu*, M(H2)/M(HI), and 12+log(O/H) regardless of the adopted X_CO. The total and molecular gas consumption timescales are anticorrelated with the SSFR.

preprint2014arXiv

Cold gas properties of the Herschel Reference Survey. III. Molecular gas stripping in cluster galaxies

The HRS is a complete volume-limited sample of nearby objects including Virgo cluster and isolated objects. Using a recent compilation of HI and CO data we study the effects of the cluster on the molecular gas content of spiral galaxies. We first identify M* as the scaling variable that traces the total H2 mass of galaxies better. We show that, on average, HI-deficient galaxies are significantly offset from the M(H2) vs. M* relation for HI-normal galaxies. We use the M(H2) vs. M* scaling relation to define the H2-deficiency parameter. This parameter shows a weak and scattered relation with the HI-def, here taken as a proxy for galaxy interactions with the cluster environment. We also show that, as for the HI, the extent of the H2 disc decreases with increasing HI-deficiency. These results show that cluster galaxies have, on average, a lower H2 content than similar objects in the field. The slope of the H2-def vs. HI-def relation is less than 1, while the D(HI)/D(i) vs. HI-def relation is steeper than the D(CO)/D(i) vs. HI-def relation, thereby indicating that the H2 gas is removed less efficiently than the HI. This result can be understood if the HI is distributed on a flat disc more extended than the stellar disc, thus less anchored to the gravitational potential well of the galaxy than the H2. There is a clear trend between the NUV-i colour and H2-def, which suggests that H2 removal quenches the activity of star formation. This causes galaxies migrate from the blue cloud to the green valley and, eventually, to the red sequence. The total gas-consumption timescale of gas deficient cluster galaxies is comparable to that of isolated systems, and is significantly larger than the typical timescale for total gas removal in a ram pressure stripping process, thus suggesting that ram pressure, rather than starvation, is the dominant process driving the evolution of these cluster galaxies.

preprint2014arXiv

Cool dust heating and temperature mixing in nearby star-forming galaxies

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 β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.

preprint2014arXiv

Dust Spectral Energy Distributions of Nearby Galaxies: an Insight from the Herschel Reference Survey

We gather infrared (IR) photometric data from 8 to 500 microns (Spitzer, WISE, IRAS and Herschel) for all of the HRS galaxies. Draine & Li (2007) models are fit to the data from which the stellar contribution has been carefully removed. We find that our photometric coverage is sufficient to constrain all of the models parameters and that a strong constraint on the 20-60 microns range is mandatory to estimate the relative contribution of the photo-dissociation regions to the IR SED. The SED models tend to systematically under-estimate the observed 500 microns flux densities, especially for low mass systems. We provide the output parameters for all of the galaxies: the minimum intensity of the interstellar radiation field (ISRF), the fraction of PAH, the relative contribution of PDR and evolved stellar population to the dust heating, the $M_{dust}$ and the $L_{IR}$. For a subsample of gas-rich galaxies, we analyze the relations between these parameters and the integrated properties of galaxies, such as $M_*$, SFR, metallicity, H$α$ and H-band surface brightness, and the FUV attenuation. A good correlation between the fraction of PAH and the metallicity is found implying a weakening of the PAH emission in galaxies with low metallicities. The intensity of the IRSF and the H-band and H$α$ surface brightnesses are correlated, suggesting that the diffuse dust component is heated by both the young stars in star forming regions and the diffuse evolved population. We use these results to provide a new set of IR templates calibrated with Herschel observations on nearby galaxies and a mean SED template to provide the z=0 reference for cosmological studies. For the same purpose, we put our sample on the SFR-$M_*$ diagram. The templates are compared to the most popular IR SED libraries, enlightening a large discrepancy between all of them in the 20-100 microns range.

preprint2014arXiv

PACS photometry of the Herschel Reference Survey - Far-infrared/sub-millimeter colours as tracers of dust properties in nearby galaxies

We present Herschel/PACS 100 and 160 micron integrated photometry for the 323 galaxies in the Herschel Reference Survey (HRS), a K-band-, volume-limited sample of galaxies in the local Universe. Once combined with the Herschel/SPIRE observations already available, these data make the HRS the largest representative sample of nearby galaxies with homogeneous coverage across the 100-500 micron wavelength range. In this paper, we take advantage of this unique dataset to investigate the properties and shape of the far-infrared/sub-millimeter spectral energy distribution in nearby galaxies. We show that, in the stellar mass range covered by the HRS (8<log(M*/Msun)<12), the far-infrared/sub-millimeter colours are inconsistent with a single modified black-body having the same dust emissivity index beta for all galaxies. In particular, either beta decreases, or multiple temperature components are needed, when moving from metal-rich/gas-poor to metal-poor/gas-rich galaxies. We thus investigate how the dust temperature and mass obtained from a single modified black-body depend on the assumptions made on beta. We show that, while the correlations between dust temperature, galaxy structure and star formation rate are strongly model dependent, the dust mass scaling relations are much more reliable, and variations of beta only change the strength of the observed trends.

preprint2014arXiv

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

With the newly available SPIRE images at 250 and 500 micron from Herschel Space Observatory, we study quantitative correlations over a sub-kpc scale among three distinct emission components in the interstellar medium of the nearby spiral galaxy M 81 (NGC 3031): (a) $I_{8}$ or $I_{24}$, the surface brightness of the mid-infrared emission observed in the Spitzer IRAC 8 or MIPS 24 micron band, with $I_8$ and $I_{24}$ being dominated by the emissions from Polycyclic Aromatic Hydrocarbons (PAHs) and very small grains (VSGs) of dust, respectively; (b) $I_{500}$, that of the cold dust continuum emission in the Herschel SPIRE 500 micron band, dominated by the emission from large dust grains heated by evolved stars, and (c) $I_{{\rm H}α}$, a nominal surface brightness of the H$α$ line emission, from gas ionized by newly formed massive stars. The results from our correlation study, free from any assumption on or modeling of dust emissivity law or dust temperatures, present solid evidence for significant heating of PAHs and VSGs by evolved stars. In the case of M 81, about 67% (48%) of the 8 micron (24 micron) emission derives its heating from evolved stars, with the remainder attributed to radiation heating associated with ionizing stars.

preprint2014arXiv

The Heating of Mid-Infrared Dust in the Nearby Galaxy M33: A Testbed for Tracing Galaxy Evolution

Because the 8 μm polycyclic aromatic hydrocarbon (PAH) emission has been found to correlate with other well-known star formation tracers, it has widely been used as a star formation rate (SFR) tracer. There are, however, studies that challenge the accuracy and reliability of the 8 μm emission as a SFR tracer. Our study, part of the Herschel M33 Extended Survey (HERM33ES) open time key program, aims at addressing this issue by analyzing the infrared emission from the nearby spiral galaxy M33 at the high spatial scale of 75 pc. Combining data from the Herschel Space Observatory and the Spitzer Space Telescope we find that the 8 μm emission is better correlated with the 250 μm emission, which traces cold interstellar gas, than with the 24 μm emission. The L(8)/L(24) ratio is highly depressed in 24 μm luminous regions, which correlate with known HII regions. We also compare our results with the dust emission models by Draine & Li (2007). We confirm that the depression of 8 μm PAH emission near star-forming regions is higher than what is predicted by models; this is possibly an effect of increased stellar radiation from young stars destroying the dust grains responsible for the 8 μm emission as already suggested by other authors. We find that the majority of the 8 μm emission is fully consistent with heating by the diffuse interstellar medium, similar to what recently determined for the dust emission in M31 by Draine at al. (2013). We also find that the fraction of 8 μm emission associated with the diffuse interstellar radiation field ranges between 60% and 80% and is 40% larger than the diffuse fraction at 24 μm.

preprint2014arXiv

The Herschel Exploitation of Local Galaxy Andromeda (HELGA): IV. Dust scaling relations at sub-kpc resolution

The imprints of dust-starlight interactions are visible in scaling relations between stellar characteristics, star formation parameters and dust properties. We aim to examine dust scaling relations on a sub-kpc resolution in the Andromeda galaxy (M31) by comparing the properties on a local and global scale to other galaxies of the local universe. New Herschel observations are combined with available data from GALEX, SDSS, WISE and Spitzer to construct a dataset covering UV to submm wavelengths. We work at the resolution of the SPIRE $500\; μ$m beam, with pixels corresponding to physical regions of 137 x 608 pc in the galaxy's disk. A panchromatic spectral energy distribution was modelled for each pixel and several dust scaling relations are investigated. We find, on a sub-kpc scale, strong correlations between $M_d/M_\star$ and NUV-r, and between $M_d/M_\star$ and $μ_\star$ (the stellar mass surface density). Striking similarities with corresponding relations based on integrated galaxies are found. We decompose M31 in four macro-regions based on their FIR morphology; the bulge, inner disk, star forming ring and the outer disk. All regions closely follow the galaxy-scale average trends. The specific star formation characteristics we derive for these macro-regions give strong hints of an inside-out formation of the bulge-disk geometry, as well as an internal downsizing process. However, within each macro-region, a great diversity in individual micro-regions is found. Furthermore, we confirm that dust in the bulge of M31 is heated only by the old stellar populations. In general, the local dust scaling relations indicate that the dust content in M31 is maintained by a subtle interplay of past and present star formation. The similarity with galaxy-based relations strongly suggests that they are in situ correlations, with underlying processes that must be local in nature. (Abriged)

preprint2014arXiv

The Herschel Virgo Cluster Survey. XV. Planck submillimetre sources in the Virgo Cluster

We cross-correlate the Planck Catalogue of Compact Sources (PCCS) with the fully sampled 84 deg2 Herschel Virgo Cluster Survey (HeViCS) fields. We search for and identify the 857 and 545 GHz PCCS sources in the HeViCS fields by studying their FIR/submm and optical counterparts. We find 84 and 48 compact Planck sources in the HeViCS fields at 857 and 545 GHz, respectively. Almost all sources correspond to individual bright Virgo Cluster galaxies. The vast majority of the Planck detected galaxies are late-type spirals, with the Sc class dominating the numbers, while early-type galaxies are virtually absent from the sample, especially at 545 GHz. We compare the HeViCS SPIRE flux densities for the detected galaxies with the four different PCCS flux density estimators and find an excellent correlation with the aperture photometry flux densities, even at the highest flux density levels. We find only seven PCCS sources in the HeViCS fields without a nearby galaxy as obvious counterpart, and conclude that all of these are dominated by Galactic cirrus features or are spurious detections. No Planck sources in the HeViCS fields seem to be associated to high-redshift proto-clusters of dusty galaxies or strongly lensed submm sources. Finally, our study is the first empirical confirmation of the simulation-based estimated completeness of the PCCS, and provides a strong support of the internal PCCS validation procedure.

preprint2014arXiv

The Herschel Virgo Cluster Survey. XVIII. Star-forming dwarf galaxies in a cluster environment

To assess the effects of the cluster environment on the different components of the interstellar medium, we analyse the FIR-submm properties of a sample of star-forming dwarf (SFD) galaxies detected by the Herschel Virgo Cluster Survey (HeViCS). We determine dust masses and dust temperatures by fitting a modified black body (MBB) function to the spectral energy distributions (SEDs). Stellar and gas masses, star formation rates (SFRs), and metallicities are obtained from the analysis of a set of ancillary data. Dust is detected in 49 out of 140 optically identified dwarfs covered by the HeViCS field; considering only dwarfs brighter than $m_B$ = 18 mag, this gives a detection rate of 43%. After evaluating different emissivity indices, we find that the FIR-submm SEDs are best-fit by $β$=1.5, with a median dust temperature $T_d$ = 22.4 K. Assuming $β$=1.5, 67% of the 23 galaxies detected in all five Herschel bands show emission at 500 $μ$m in excess of the MBB model. The excess is inversely correlated with SFR and stellar masses. To study the variations in the global properties of our sample due to environmental effects, we compare the Virgo SFDs to other Herschel surveys, such as KINGFISH, the Dwarf Galaxy Survey (DGS), and the HeViCS bright galaxy catalogue (BGC). We explore the relations between stellar mass and HI fraction, specific SFR, dust fraction, gas-to-dust ratio over a wide range of stellar masses. Highly HI-deficient Virgo dwarf galaxies are mostly characterised by quenched star formation activity and lower dust fractions giving hints for dust stripping in cluster dwarfs. However, we find that the fraction of dust removed has to be less than that of the HI component. Since the Virgo SFDs are likely to be crossing the cluster for the first time, a longer timescale might be necessary to strip the more centrally concentrated dust distribution.

preprint2014arXiv

The identification of dust heating mechanisms in nearby galaxies using Herschel 160/250 and 250/350 micron surface brightness ratios

We examined variations in the 160/250 and 250/350 micron surface brightness ratios within 24 nearby (<30 Mpc) face-on spiral galaxies observed with the Herschel Space Observatory to identify the heating mechanisms for dust emitting at these wavelengths. The analysis consisted of both qualitative and quantitative comparisons of the 160/250 and 250/350 micron ratios to H alpha and 24 micron surface brightnesses, which trace the light from star forming regions, and 3.6 micron emission, which traces the light from the older stellar populations of the galaxies. We find broad variations in the heating mechanisms for the dust. In one subset of galaxies, we found evidence that emission at <=160 microns (and in rare cases potentially at <=350 microns) originates from dust heated by star forming regions. In another subset, we found that the emission at >=250 microns (and sometimes at >=160 microns) originates from dust heated by the older stellar population. In the rest of the sample, either the results are indeterminate or both of these stellar populations may contribute equally to the global dust heating. The observed variations in dust heating mechanisms does not necessarily match what has been predicted by dust emission and radiative transfer models, which could lead to overestimated dust temperatures, underestimated dust masses, false detections of variability in dust emissivity, and inaccurate star formation rate measurements.

preprint2014arXiv

The IRAM M33 CO(2-1) Survey - A complete census of the molecular gas out to 7 kpc

In order to study the ISM and the interplay between the atomic and molecular components in a low-metallicity environment, we present a complete high angular and spectral resolution map and data cube of the 12CO(2-1) emission from the Local Group galaxy M33. Its metallicity is roughly half-solar, such that we can compare its ISM with that of the Milky Way with the main changes being the metallicity and the gas mass fraction. The data have a 12" angular resolution (50pc) with a spectral resolution of 2.6 km/s and a mean noise level of 20 mK per channel in antenna temperature. A radial cut along the major axis was also observed in the 12CO(1-0) line. The CO data cube and integrated intensity map are optimal when using HI data to define the baseline window and the velocities over which the CO emission is integrated. Great care was taken when building these maps, testing different windowing and baseline options and investigating the effect of error beam pickup. The total CO(2-1) luminosity is 2.8e7 K km/s pc2, following the spiral arms in the inner disk. There is no clear variation in the CO(2-1/1-0) intensity ratio with radius and the average value is roughly 0.8. The total molecular gas mass is estimated, using a N(H2)/Ico(1-0)=4e20 cm-2/(K km/s) conversion factor, to be 3.1e8 Msol. The CO spectra in the cube were shifted to zero velocity by subtracting the velocity of the HI peak from the CO spectra. Hence, the velocity dispersion between the atomic and molecular components is extremely low, independently justifying the use of the HI line in building our maps. Stacking the spectra in concentric rings shows that the CO linewidth and possibly the CO-HI velocity dispersion decrease in the outer disk. Using the CO(2-1) emission to trace the molecular gas, the probability distribution function of the H2 column density shows an excess at high column density above a log normal distribution.

preprint2014arXiv

The physical characteristics of the gas in the disk of Centaurus A using the Herschel Space Observatory

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

preprint2014arXiv

The relationship between polycyclic aromatic hydrocarbon emission and far-infrared dust emission from NGC 2403 and M83

We examine the relation between polycyclic aromatic hydrocarbon (PAH) emission at 8 microns and far-infrared emission from hot dust grains at 24 microns and from large dust grains at 160 and 250 microns in the nearby spiral galaxies NGC 2403 and M83 using data from the Spitzer Space Telescope and Herschel Space Observatory. We find that the PAH emission in NGC 2403 is better correlated with emission at 250 microns from dust heated by the diffuse interstellar radiation field (ISRF) and that the 8/250 micron surface brightness ratio is well-correlated with the stellar surface brightness as measured at 3.6 microns. This implies that the PAHs in NGC 2403 are intermixed with cold large dust grains in the diffuse interstellar medium (ISM) and that the PAHs are excited by the diffuse ISRF. In M83, the PAH emission appears more strongly correlated with 160 micron emission originating from large dust grains heated by star forming regions. However, the PAH emission in M83 is low where the 24 micron emission peaks within star forming regions, and enhancements in the 8/160 micron surface brightness ratios appear offset relative to the dust and the star forming regions within the spiral arms. This suggests that the PAHs observed in the 8 micron band are not excited locally within star forming regions but either by light escaping non-axisymmetrically from star forming regions or locally by young, non-photoionising stars that have migrated downstream from the spiral density waves. The results from just these two galaxies show that PAHs may be excited by different stellar populations in different spiral galaxies.

preprint2013arXiv

An Overview of the Dwarf Galaxy Survey

The Dwarf Galaxy Survey (DGS) program is studying low-metallicity galaxies using 230h of far-infrared (FIR) and submillimetre (submm) photometric and spectroscopic observations of the Herschel Space Observatory and draws to this a rich database of a wide range of wavelengths tracing the dust, gas and stars. This sample of 50 galaxies includes the largest metallicity range achievable in the local Universe including the lowest metallicity (Z) galaxies, 1/50 Zsun, and spans 4 orders of magnitude in star formation rates. The survey is designed to get a handle on the physics of the interstellar medium (ISM) of low metallicity dwarf galaxies, especially on their dust and gas properties and the ISM heating and cooling processes. The DGS produces PACS and SPIRE maps of low-metallicity galaxies observed at 70, 100, 160, 250, 350, and 500 mic with the highest sensitivity achievable to date in the FIR and submm. The FIR fine-structure lines, [CII] 158 mic, [OI] 63 mic, [OI] 145 mic, [OIII] 88 mic, [NIII] 57 mic and [NII] 122 and 205 mic have also been observed with the aim of studying the gas cooling in the neutral and ionized phases. The SPIRE FTS observations include many CO lines (J=4-3 to J=13-12), [NII] 205 mic and [CI] lines at 370 and 609 mic. This paper describes the sample selection and global properties of the galaxies, the observing strategy as well as the vast ancillary database available to complement the Herschel observations. The scientific potential of the full DGS survey is described with some example results included.

preprint2013arXiv

Calibration of the total infrared luminosity of nearby galaxies from Spitzer and Herschel bands

We present new empirical calibrations to estimate resolved and integrated total infrared luminosities from Spitzer and Herschel bands used as monochromatic or combined tracers. We base our calibrations on resolved elements of nearby galaxies (3 to 30 Mpc) observed with Herschel. We perform a resolved SED modelling of these objects using the Draine and Li (2007) dust models and investigate the influence of the addition of SPIRE measurements in the estimation of LTIR. We find that using data up to 250 um leads to local LTIR values consistent with those obtained with a complete coverage (up to 500 um) within 10 per cent for most of our resolved elements. We then study the distribution of energy in the resolved SEDs of our galaxies. The bulk of energy (30-50 per cent) is contained in the (70-160 um) band. The (24-70 um) fraction decreases with increasing metallicity. The (160-1100 um) submillimeter band can account for up to 25 per cent of the LTIR in metal-rich galaxies. We investigate the correlation between TIR surface brightnesses/luminosities and monochromatic Spitzer and Herschel surface brightnesses/luminosities. The three PACS bands can be used as reliable monochromatic estimators of the LTIR, the 100 um band being the most reliable monochromatic tracer. There is also a strong correlation between the SPIRE 250um and LTIR, although with more scatter than for the PACS relations. We also study the ability of our monochromatic relations to reproduce integrated LTIR of nearby galaxies as well as LTIR of z=1-3 sources. Finally, we provide calibration coefficients that can be used to derive TIR surface brightnesses/luminosities from a combination of Spitzer and Herschel surface brightnesses/fluxes and analyse the associated uncertainties.

preprint2013arXiv

Cold dust but warm gas in the unusual elliptical galaxy NGC 4125

Data from the Herschel Space Observatory have revealed an unusual elliptical galaxy, NGC 4125, which has strong and extended submillimeter emission from cold dust but only very strict upper limits to its CO and HI emission. Depending on the dust emissivity, the total dust mass is 2-5x10^6 Msun. While the neutral gas-to-dust mass ratio is extremely low (< 12-30), including the ionized gas traced by [CII] emission raises this limit to < 39-100. The dust emission follows a similar r^{1/4} profile to the stellar light and the dust to stellar mass ratio is towards the high end of what is found in nearby elliptical galaxies. We suggest that NGC 4125 is currently in an unusual phase where evolved stars produced in a merger-triggered burst of star formation are pumping large amounts of gas and dust into the interstellar medium. In this scenario, the low neutral gas-to-dust mass ratio is explained by the gas being heated to temperatures >= 10^4 K faster than the dust is evaporated. If galaxies like NGC 4125, where the far-infrared emission does not trace neutral gas in the usual manner, are common at higher redshift, this could have significant implications for our understanding of high redshift galaxies and galaxy evolution.

preprint2013arXiv

Dense Gas in M33 (HerM33es)

We aim to better understand the emission of molecular tracers of the diffuse and dense gas in giant molecular clouds and the influence that metallicity, optical extinction, density, far-UV field, and star formation rate have on these tracers. Using the IRAM 30m telescope, we detected HCN, HCO+, 12CO, and 13CO in six GMCs along the major axis of M33 at a resolution of ~ 114pc and out to a radial distance of 3.4kpc. Optical, far-infrared, and submillimeter data from Herschel and other observatories complement these observations. To interpret the observed molecular line emission, we created two grids of models of photon-dominated regions, one for solar and one for M33-type subsolar metallicity. The observed HCO+/HCN line ratios range between 1.1 and 2.5. Similarly high ratios have been observed in the Large Magellanic Cloud. The HCN/CO ratio varies between 0.4% and 2.9% in the disk of M33. The 12CO/13CO line ratio varies between 9 and 15 similar to variations found in the diffuse gas and the centers of GMCs of the Milky Way. Stacking of all spectra allowed HNC and C2H to be detected. The resulting HCO+/HNC and HCN/HNC ratios of ~ 8 and 6, respectively, lie at the high end of ratios observed in a large set of (ultra-)luminous infrared galaxies. HCN abundances are lower in the subsolar metallicity PDR models, while HCO+ abundances are enhanced. For HCN this effect is more pronounced at low optical extinctions. The observed HCO+/HCN and HCN/CO line ratios are naturally explained by subsolar PDR models of low optical extinctions between 4 and 10 mag and of moderate densities of n = 3x10^3 - 3x10^4 cm^-3, while the FUV field strength only has a small effect on the modeled line ratios. The line ratios are almost equally well reproduced by the solar-metallicity models, indicating that variations in metallicity only play a minor role in influencing these line ratios.

preprint2013arXiv

Gas and dust cooling along the major axis of M33 (HerM33es): ISO/LWS CII observations

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%.

preprint2013arXiv

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

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

preprint2013arXiv

Regional variations in the dense gas heating and cooling in M51 from Herschel far-infrared spectroscopy

We present Herschel PACS and SPIRE spectroscopy of the most important far-infrared cooling lines in M51, [CII](158 μm), [NII](122 & 205 μm), [OI](63 and 145 μm) and [OIII](88 μm). We compare the observed flux of these lines with the predicted flux from a photon dominated region model to determine characteristics of the cold gas such as density, temperature and the far-ultraviolet radiation field, G_0, resolving details on physical scales of roughly 600 pc. We find an average [CII]/F_TIR of 4 x 10^{-3}, in agreement with previous studies of other galaxies. A pixel-by-pixel analysis of four distinct regions of M51 shows a radially decreasing trend in both the far-ultraviolet (FUV) radiation field, G_0 and the hydrogen density, n, peaking in the nucleus of the galaxy, then falling off out to the arm and interarm regions. We see for the first time that the FUV flux and gas density are similar in the differing environments of the arm and interarm regions, suggesting that the inherent physical properties of the molecular clouds in both regions are essentially the same.

preprint2013arXiv

Revealing the cold dust in low-metallicity environments: I - Photometry analysis of the Dwarf Galaxy Survey with Herschel

We present new photometric data from our Herschel Key Programme, the Dwarf Galaxy Survey (DGS), dedicated to the observation of the gas and dust in 48 low-metallicity environments. They were observed with PACS and SPIRE onboard Herschel at 70,100,160,250,350, and 500 microns. We focus on a systematic comparison of the derived FIR properties (FIR luminosity, dust mass, dust temperature and emissivity index) with more metal-rich galaxies and investigate the detection of a potential submm excess. The data reduction method is adapted for each galaxy to derive the most reliable photometry from the final maps. PACS flux densities are compared with the MIPS 70 and 160 microns bands. We use colour-colour diagrams and modified blackbody fitting procedures to determine the dust properties of the DGS galaxies. We also include galaxies from the Herschel KINGFISH sample, containing more metal-rich environments, totalling 109 galaxies. The location of the DGS galaxies on Herschel colour-colour diagrams highlights the differences in global environments of low-metallicity galaxies. The dust in DGS galaxies is generally warmer than in KINGFISH galaxies (T_DGS~32 K, T_KINGFISH~23 K). The emissivity index, beta, is ~1.7 in the DGS, but metallicity does not make a strong effect on beta. The dust-to-stellar mass ratio is lower in low-metallicity galaxies: M_dust/M_star~0.02% for the DGS vs 0.1% for KINGFISH. Per unit dust mass, dwarf galaxies emit ~6 times more in the FIR than higher metallicity galaxies. Out of the 22 DGS galaxies detected at 500 micron, 41% present an excess in the submm not explained by our dust SED model. The excess mainly appears in lower metallicity galaxies (12+log(O/H) < 8.3), and the strongest excesses are detected in the most metal-poor galaxies. We stress the need for observations longwards of the Herschel wavelengths to detect any submm excess appearing beyond 500 micron.

preprint2013arXiv

Spectral Energy Distributions of HII regions in M33 (HerM33es)

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.

preprint2013arXiv

Star Formation & Dust Heating in the FIR Bright Sources of M83

We investigate star formation and dust heating in the compact FIR bright sources detected in the Herschel maps of M83. We use the source extraction code GETSOURCES to detect and extract sources in the FIR, as well as their photometry in the MIR and Hα. By performing infrared SED fitting and applying an Hα based star formation rate (SFR) calibration, we derive the dust masses and temperatures, SFRs, gas masses and star formation efficiencies (SFEs). The detected sources lie exclusively on the spiral arms and represent giant molecular associations (GMAs), with gas masses and sizes of 10^6-10^8 solar masses and 200-300 pc, respectively. The inferred parameters show little to no radial dependence and there is only a weak correlation between the SFRs and gas masses, which suggests that more massive clouds are less efficient at forming stars. Dust heating is mainly due to local star formation. However, although the sources are not optically thick, the total intrinsic young stellar population luminosity can almost completely account for the dust luminosity. This suggests that other radiation sources contribute to the dust heating as well and approximately compensate for the unabsorbed fraction of UV light.

preprint2013arXiv

Towards understanding the relation between the gas and the attenuation in galaxies at kpc scales

[abridged] Aims. The aim of the present paper is to provide new and more detailed relations at the kpc scale between the gas surface density and the face-on optical depth directly calibrated on galaxies, in order to compute the attenuation not only for semi-analytic models but also observationally as new and upcoming radio observatories are able to trace gas ever farther in the Universe. Methods. We have selected a sample of 4 nearby resolved galaxies and a sample of 27 unresolved galaxies from the Herschel Reference Survey and the Very Nearby Galaxies Survey, for which we have a large set of multi-wavelength data from the FUV to the FIR including metallicity gradients for resolved galaxies, along with radio HI and CO observations. For each pixel in resolved galaxies and for each galaxy in the unresolved sample, we compute the face-on optical depth from the attenuation determined with the CIGALE SED fitting code and an assumed geometry. We determine the gas surface density from HI and CO observations with a metallicity-dependent XCO factor. Results. We provide new, simple to use, relations to determine the face-on optical depth from the gas surface density, taking the metallicity into account, which proves to be crucial for a proper estimate. The method used to determine the gas surface density or the face-on optical depth has little impact on the relations except for galaxies that have an inclination over 50d. Finally, we provide detailed instructions on how to compute the attenuation practically from the gas surface density taking into account possible information on the metallicity. Conclusions. Examination of the influence of these new relations on simulated FUV and IR luminosity functions shows a clear impact compared to older oft-used relations, which in turn could affect the conclusions drawn from studies based on large scale cosmological simulations.

preprint2013arXiv

Ultraviolet to infrared emission of z>1 galaxies: Can we derive reliable star formation rates and stellar masses?

We seek to derive star formation rates (SFR) and stellar masses (M_star) in distant galaxies and to quantify the main uncertainties affecting their measurement. We explore the impact of the assumptions made in their derivation with standard calibrations or through a fitting process, as well as the impact of the available data, focusing on the role of IR emission originating from dust. We build a sample of galaxies with z>1, all observed from the UV to the IR (rest frame). The data are fitted with the code CIGALE, which is also used to build and analyse a catalogue of mock galaxies. Models with different SFHs are introduced. We define different set of data, with or without a good sampling of the UV range, NIR, and thermal IR data. The impact of these different cases on the determination of M_star and SFR are analysed. Exponentially decreasing models with a redshift formation of the stellar population z ~8 cannot fit the data correctly. The other models fit the data correctly at the price of unrealistically young ages when the age of the single stellar population is taken to be a free parameter. The best fits are obtained with two stellar populations. As long as one measurement of the dust emission continuum is available, SFR are robustly estimated whatever the chosen model is, including standard recipes. M_star measurement is more subject to uncertainty, depending on the chosen model and the presence of NIR data, with an impact on the SFR-M_star scatter plot. Conversely, when thermal IR data from dust emission are missing, the uncertainty on SFR measurements largely exceeds that of stellar mass. Among all physical properties investigated here, the stellar ages are found to be the most difficult to constrain and this uncertainty acts as a second parameter in SFR measurements and as the most important parameter for M_star measurements.

preprint2013arXiv

Variation in the dust emissivity index across M33 with Herschel and Spitzer (HerM33es)

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.

preprint2012arXiv

Cool and warm dust emission from M33 (HerM33es)

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).

preprint2012arXiv

Dust and gas power-spectrum in M33 (HERM33ES)

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.

preprint2012arXiv

FIR colours of nearby late-type galaxies in the Herschel Reference Survey

We study the far infrared (60-500 $μ$m) colours of late-type galaxies in the $Herschel$ Reference Survey, a K-band selected, volume limited sample of nearby galaxies. The far infrared colours are correlated with each other, with tighter correlations for the indices that are closer in wavelength. We also compare the different colour indices to various tracers of the physical properties of the target galaxies, such as the surface brightness of the ionising and non-ionising stellar radiation, the dust attenuation and the metallicity. The emission properties of the cold dust dominating the far infrared spectral domain are regulated by the properties of the interstellar radiation field. Consistent with that observed in nearby, resolved galaxies, our analysis shows that the ionising and the non-ionising stellar radiation, including that emitted by the most evolved, cold stars, both contribute to the heating of the cold dust component. This work also shows that metallicity is another key parameter characterising the cold dust emission of normal, late-type galaxies. A single modified black body with a grain emissivity index $β$=1.5 better fits the observed SPIRE flux density ratios $S250/S350$ vs. $S350/S500$ than $β$=2, although values of $β$ $\simeq$ 2 are possible in metal rich, high surface brightness galaxies. Values of $β$ $\lesssim$ 1.5 better represent metal poor, low surface brightness objects. This observational evidence provides strong constraints for dust emission models of normal, late type galaxies.

preprint2012arXiv

Herschel and JCMT observations of the early-type dwarf galaxy NGC 205

We present Herschel dust continuum, James Clerk Maxwell Telescope CO(3-2) observations and a search for [CII] 158 micron and [OI] 63 micron spectral line emission for the brightest early-type dwarf satellite of Andromeda, NGC 205. While direct gas measurements (Mgas ~ 1.5e+6 Msun, HI + CO(1-0)) have proven to be inconsistent with theoretical predictions of the current gas reservoir in NGC 205 (> 1e+7 Msun), we revise the missing interstellar medium mass problem based on new gas mass estimates (CO(3-2), [CII], [OI]) and indirect measurements of the interstellar medium content through dust continuum emission. Based on Herschel observations, covering a wide wavelength range from 70 to 500 micron, we are able to probe the entire dust content in NGC 205 (Mdust ~ 1.1-1.8e+4 Msun at Tdust ~ 18-22 K) and rule out the presence of a massive cold dust component (Mdust ~ 5e+5 Msun, Tdust ~ 12 K), which was suggested based on millimeter observations from the inner 18.4 arcsec. Assuming a reasonable gas-to-dust ratio of ~ 400, the dust mass in NGC 205 translates into a gas mass Mgas ~ 4-7e+6 Msun. The non-detection of [OI] and the low L_[CII]-to-L_CO(1-0) line intensity ratio (~ 1850) imply that the molecular gas phase is well traced by CO molecules in NGC 205. We estimate an atomic gas mass of 1.5e+4 Msun associated with the [CII] emitting PDR regions in NGC 205. From the partial CO(3-2) map of the northern region in NGC 205, we derive a molecular gas mass of M_H2 ~ 1.3e+5 Msun. [abridged]

preprint2012arXiv

Spectrally resolved CII emission in M~33 ({\tt HerM33es}): Physical conditions and kinematics around BCLMP 691

This work presents high spectral resolution observations of the \CII\ line at 158 \micron, one of the major cooling lines of the interstellar medium, taken with the HIFI heterodyne spectrometer on the Herschel satellite. In BCLMP 691, an \HII\ region far north (3.3 kpc) in the disk of M 33, the \CII\ and CO line profiles show similar velocities within $0.5 \kms$, while the \HI\ line velocities are systematically shifted towards lower rotation velocities by $\sim 5\kms$. Observed at the same $12"$ angular resolution, the \CII\ lines are broader than those of CO by about 50% but narrower than the \HI\ lines. The \CII\ line to far-infrared continuum ratio suggests a photoelectric heating efficiency of 1.1%. The data, together with published models indicate a UV field $G_0 \sim 100$ in units of the solar neighborhood value, a gas density $n_H \sim 1000 \cc$, and a gas temperature $T\sim 200$ K. Adopting these values, we estimate the C$^+$ column density to be $N_{C^+} \approx 1.3 \times 10^{17} \cmt$. The \CII\ emission comes predominantly from the warm neutral region between the \HII\ region and the cool molecular cloud behind it. From published abundances, the inferred C$^+$ column corresponds to a hydrogen column density of $N_H \sim 2 \times 10^{21} \cmt$. The CO observations suggest that $N_H = 2 N_{H_2} \sim 3.2 \times 10^{21} \cmt$ and 21cm measurements, also at $12"$ resolution, yield $N_\HI \approx 1.2 \times 10^{21} \cmt$ within the \CII\ velocity range. Thus, some H$_2$ not detected in CO must be present, in agreement with earlier findings based on the SPIRE 250 -- 500 $μ$m emission.

preprint2012arXiv

Submillimetre Photometry of 323 Nearby Galaxies from the Herschel Reference Survey

The Herschel Reference Survey (HRS) is a guaranteed time Herschel key project aimed at studying the physical properties of the interstellar medium in galaxies of the nearby universe. This volume limited, K-band selected sample is composed of galaxies spanning the whole range of morphological types (from ellipticals to late-type spirals) and environments (from the field to the centre of the Virgo Cluster). We present flux density measurements of the whole sample of 323 galaxies of the HRS in the three bands of the Spectral and Photometric Imaging Receiver (SPIRE), at 250, 350 and 500 microns. Aperture photometry is performed on extended galaxies and point spread function (PSF) fitting on timeline data for unresolved objects; we carefully estimate errors and upper limits. The flux densities are found to be in good agreement with those of the HeViCS and KINGFISH key projects in all SPIRE bands, and of the Planck consortium at 350 and 550 microns, for the galaxies in common. This submillimetre catalogue of nearby galaxies is a benchmark for the study of the dust properties in the local universe, giving the zero redshift reference for any cosmological survey.

preprint2012arXiv

The Dust & Gas Properties of M83

We examine the dust and gas properties of the nearby, barred galaxy M83, which is part of the Very Nearby Galaxy Survey. Using images from the PACS and SPIRE instruments of Herschel, we examine the dust temperature and dust mass surface density distribution. We find that the nuclear, bar and spiral arm regions exhibit higher dust temperatures and masses compared to interarm regions. However, the distribution of dust temperature and mass are not spatially coincident. Assuming a trailing spiral structure, the dust temperature peaks in the spiral arms lie ahead of the dust surface density peaks. The dust mass surface density correlates well with the distribution of molecular gas as traced by CO (J=3-2) images (JCMT) and the star formation rate as traced by H?2 with a correction for obscured star formation using 24 micron emission. Using HI images from THINGS to trace the atomic gas component, we make total gas mass surface density maps and calculate the gas-to-dust ratio. We find a mean gas-to-dust ratio of 84 \pm 4 with higher values in the inner region assuming a constant CO-to-H2 conversion factor. We also examine the gas-to-dust ratio using CO-to-H2 conversion factor that varies with metallicity.

preprint2012arXiv

The dust scaling relations of the Herschel Reference Survey

We combine new Herschel/SPIRE sub-millimeter observations with existing multiwavelength data to investigate the dust scaling relations of the Herschel Reference Survey, a magnitude-, volume-limited sample of ~300 nearby galaxies in different environments. We show that the dust-to-stellar mass ratio anti-correlates with stellar mass, stellar mass surface density and NUV-r colour across the whole range of parameters covered by our sample. Moreover, the dust-to-stellar mass ratio decreases significantly when moving from late- to early-type galaxies. These scaling relations are similar to those observed for the HI gas-fraction, supporting the idea that the cold dust is tightly coupled to the cold atomic gas component in the interstellar medium. We also find a weak increase of the dust-to-HI mass ratio with stellar mass and colour but no trend is seen with stellar mass surface density. By comparing galaxies in different environments we show that, although these scaling relations are followed by both cluster and field galaxies, HI-deficient systems have, at fixed stellar mass, stellar mass surface density and morphological type systematically lower dust-to-stellar mass and higher dust-to-HI mass ratios than HI-normal/field galaxies. This provides clear evidence that dust is removed from the star-forming disk of cluster galaxies but the effect of the environment is less strong than what is observed in the case of the HI disk. Such effects naturally arise if the dust disk is less extended than the HI and follows more closely the distribution of the molecular gas phase, i.e., if the dust-to-atomic gas ratio monotonically decreases with distance from the galactic center.

preprint2012arXiv

The gas-to-dust mass ratio of Centaurus A as seen by Herschel

We present photometry of the nearby galaxy NGC 5128 (Centaurus A) observed with the PACS and SPIRE instruments on board the Herschel Space Observatory, at 70, 160, 250, 350 and 500 μm, as well as new CO J = 3-2 observations taken with the HARP-B instrument on the JCMT. Using a single component modified blackbody, we model the dust spectral energy distribution within the disk of the galaxy using all five Herschel wavebands, and find dust temperatures of ~30 K towards the centre of the disk and a smoothly decreasing trend to ~20 K with increasing radius. We find a total dust mass of (1.59 \pm 0.05) \times 10^7 M\odot, and a total gas mass of (2.7 \pm 0.2) \times 10^9 M\odot. The average gas-to-dust mass ratio is 103 \pm 8 but we find an interesting increase in this ratio to approximately 275 toward the centre of Cen A. We discuss several possible physical processes that may be causing this effect, including dust sputtering, jet entrainment and systematic variables such as the XCO factor. Dust sputtering by X-rays originating in the AGN or the removal of dust by the jets are our most favoured explanations.

preprint2012arXiv

The Herschel Exploitation of Local Galaxy Andromeda (HELGA) II: Dust and Gas in Andromeda

We present an analysis of the dust and gas in Andromeda, using Herschel images sampling the entire far-infrared peak. We fit a modified-blackbody model to ~4000 quasi-independent pixels with spatial resolution of ~140pc and find that a variable dust-emissivity index (beta) is required to fit the data. We find no significant long-wavelength excess above this model suggesting there is no cold dust component. We show that the gas-to-dust ratio varies radially, increasing from ~20 in the center to ~70 in the star-forming ring at 10kpc, consistent with the metallicity gradient. In the 10kpc ring the average beta is ~1.9, in good agreement with values determined for the Milky Way (MW). However, in contrast to the MW, we find significant radial variations in beta, which increases from 1.9 at 10kpc to ~2.5 at a radius of 3.1kpc and then decreases to 1.7 in the center. The dust temperature is fairly constant in the 10kpc ring (ranging from 17-20K), but increases strongly in the bulge to ~30K. Within 3.1kpc we find the dust temperature is highly correlated with the 3.6 micron flux, suggesting the general stellar population in the bulge is the dominant source of dust heating there. At larger radii, there is a weak correlation between the star formation rate and dust temperature. We find no evidence for 'dark gas' in M31 in contrast to recent results for the MW. Finally, we obtained an estimate of the CO X-factor by minimising the dispersion in the gas-to-dust ratio, obtaining a value of (1.9+/-0.4)x10^20 cm^-2 [K kms^-1]^-1.

preprint2012arXiv

The Herschel Exploitation of Local Galaxy Andromeda (HELGA). I: Global far-infrared and sub-mm morphology

We have obtained Herschel images at five wavelengths from 100 to 500 micron of a ~5.5x2.5 degree area centred on the local galaxy M31 (Andromeda), our nearest neighbour spiral galaxy, as part of the Herschel guaranteed time project "HELGA". The main goals of HELGA are to study the characteristics of the extended dust emission, focusing on larger scales than studied in previous observations of Andromeda at an increased spatial resolution, and the obscured star formation. In this paper we present data reduction and Herschel maps, and provide a description of the far-infrared morphology, comparing it with features seen at other wavelengths. We use high--resolution maps of the atomic hydrogen, fully covering our fields, to identify dust emission features that can be associated to M31 with confidence, distinguishing them from emission coming from the foreground Galactic cirrus. Thanks to the very large extension of our maps we detect, for the first time at far-infrared wavelengths, three arc-like structures extending out to ~21, ~26 and ~31 kpc respectively, in the south-western part of M31. The presence of these features, hosting ~2.2e6 Msol of dust, is safely confirmed by their detection in HI maps. Overall, we estimate a total dust mass of ~5.8e7 Msol, about 78% of which is contained in the two main ring-like structures at 10 and 15 kpc, at an average temperature of 16.5 K. We find that the gas-to-dust ratio declines exponentially as a function of the galacto-centric distance, in agreement with the known metallicity gradient, with values ranging from 66 in the nucleus to ~275 in the outermost region. [Abridged]

preprint2012arXiv

The Herschel Virgo Cluster Survey XII: FIR properties of optically-selected Virgo Cluster galaxies

The Herschel Virgo Cluster Survey (HeViCS) is the deepest, confusion-limited survey of the Virgo Cluster at far-infrared (FIR) wavelengths. The entire survey at full depth covers $\sim$55 sq. deg. in 5 bands (100-500 \micron), encompassing the areas around the central dominant elliptical galaxies (M87, M86 & M49) and extends as far as the NW cloud, the W cloud and the Southern extension. The survey extends beyond this region with lower sensitivity so that the total area covered is 84 sq. deg. In this paper we describe the data, the data acquisition techniques and present the detection rates of the optically selected Virgo Cluster Catalogue (VCC). We detect 254 (34%) of 750 VCC galaxies found within the survey boundary in at least one band and 171 galaxies are detected in all five bands. For the remainder of the galaxies we have measured strict upper limits for their FIR emission. The population of detected galaxies contains early- as well as late-types although the latter dominate the detection statistics. We have modelled 168 galaxies, showing no evidence of a strong synchrotron component in their FIR spectra, using a single-temperature modified blackbody spectrum with a fixed emissivity index ($β= 2$). A study of the $χ^2$ distribution indicates that this model is not appropriate in all cases, and this is supported by the FIR colours which indicate a spread in $β$=1--2. Statistical comparison of the dust mass and temperature distributions from 140 galaxies with $χ^2_{dof=3} < 7.8$ (95% confidence level) shows that late-types have typically colder, more massive dust reservoirs; the early-type dust masses have a mean of ${\rm log}(<M> / M_{\sun}) = 6.3 \pm 0.3 $, while for late-types ${\rm log}(<M> / M_{\sun}) =7.1 \pm 0.1$... (abridged)

preprint2012arXiv

The IRX-beta relation on sub-galactic scales in star-forming galaxies of the Herschel Reference Survey

UV and optical surveys are essential to gain insight into the processes driving galaxy formation and evolution. The rest-frame UV emission is key to measure the cosmic SFR. However, UV light is strongly reddened by dust. In starburst galaxies, the UV colour and the attenuation are linked, allowing to correct for dust extinction. Unfortunately, evidence has been accumulating that the relation between UV colour and attenuation is different for normal star-forming galaxies when compared to starburst galaxies. It is still not understood why star-forming galaxies deviate from the UV colour-attenuation relation of starburst galaxies. Previous work and models hint that the role of the shape of the attenuation curve and the age of stellar populations have an important role. In this paper we aim at understanding the fundamental reasons to explain this deviation. We have used the CIGALE SED fitting code to model the far UV to the far IR emission of a set of 7 reasonably face-on spiral galaxies from the HRS. We have explored the influence of a wide range of physical parameters to quantify their influence and impact on the accurate determination of the attenuation from the UV colour, and why normal galaxies do not follow the same relation as starburst galaxies. We have found that the deviation can be best explained by intrinsic UV colour differences between different regions in galaxies. Variations in the shape of the attenuation curve can also play a secondary role. Standard age estimators of the stellar populations prove to be poor predictors of the intrinsic UV colour. These results are also retrieved on a sample of 58 galaxies when considering their integrated fluxes. When correcting the emission of normal star-forming galaxies for the attenuation, it is crucial to take into account possible variations in the intrinsic UV colour as well as variations of the shape of the attenuation curve.

preprint2011arXiv

Dust heating sources in galaxies: the case of M33 (HERM33ES)

Dust emission is one of the main windows to the physics of galaxies and to star formation as the radiation from young, hot stars is absorbed by the dust and reemitted at longer wavelengths. The recently launched Herschel satellite now provides a view of dust emission in the far-infrared at an unequaled resolution and quality up to 500 μm. In the context of the Herschel HERM33ES open time key project, we are studying the moderately inclined Scd local group galaxy M33 which is located only 840 kpc away. In this article, using Spitzer and Herschel data ranging from 3.6 μm to 500 μm, along with HI, Hα maps, and GALEX ultraviolet data we have studied the emission of the dust at the high spatial resolution of 150 pc. Combining Spitzer and Herschel bands, we have provided new, inclination corrected, resolved estimators of the total infrared brightness and of the star formation rate from any combination of these bands. The study of the colors of the warm and cold dust populations shows that the temperature of the former is, at high brightness, dictated by young massive stars but, at lower brightness, heating is taken over by the evolved populations. Conversely, the temperature of the cold dust is tightly driven by the evolved stellar populations.

preprint2011arXiv

Spectral Energy Distributions of a set of HII regions in M33 (HerM33es)

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.

preprint2011arXiv

Star formation in M33 (HerM33es)

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.

preprint2011arXiv

Status of GDL - GNU Data Language

GNU Data Language (GDL) is an open-source interpreted language aimed at numerical data analysis and visualisation. It is a free implementation of the Interactive Data Language (IDL) widely used in Astronomy. GDL has a full syntax compatibility with IDL, and includes a large set of library routines targeting advanced matrix manipulation, plotting, time-series and image analysis, mapping, and data input/output including numerous scientific data formats. We will present the current status of the project, the key accomplishments, and the weaknesses - areas where contributions are welcome !

preprint2011arXiv

The Herschel M33 extended survey (HerM33es): PACS spectroscopy of the star forming region BCLMP 302

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.

preprint2011arXiv

Variation in the dust spectral index across M33

Using the Herschel PACS and SPIRE FIR/submm data, we investigate variations in the dust spectral index $β$ in the nearby spiral galaxy M33 at a linear resolution of 160 pc. We use an iteration method in two different approaches, single and two-component modified black body models. In both approaches, $β$ is higher in the central disk than in the outer disk similar to the dust temperature. There is a positive correlation between $β$ and Halpha as well as with the molecular gas traced by CO(2-1). A Monte-Carlo simulation shows that the physical parameters are better constrained when using the two-component model.

preprint2010arXiv

100 mum and 160 mum emission as resolved star-formation rate estimators in M33 (HERM33ES)

Over the past few years several studies have provided estimates of the SFR (star-formation rate) or the total infrared luminosity from just one infrared band. However these relations are generally derived for entire galaxies, which are known to contain a large scale diffuse emission that is not necessarily related to the latest star-formation episode. We provide new relations to estimate the SFR from resolved star-forming regions at 100 mum and 160 mum. We select individual star-forming regions in the nearby (840 kpc) galaxy M33. We estimate the SFR combining the emission in Halpha and at 24 mum to calibrate the emission at 100 mum and 160 mum as SFR estimators, as mapped with PACS/Herschel. The data are obtained in the framework of the HERM33ES open time key project. There is less emission in the HII regions at 160 mum than at 100 mum. Over a dynamic range of almost 2 dex in Sigma(SFR) we find that the 100 mum emission is a nearly linear estimator of the SFR, whereas that at 160 mum is slightly superlinear. The behaviour of individual star-forming regions is surprisingly similar to that of entire galaxies. At high Sigma(SFR), star formation drives the dust temperature, whereas uncertainties and variations in radiation-transfer and dust-heated processes dominate at low Sigma(SFR). Detailed modelling of both galaxies and individual star forming regions will be needed to interpret similarities and differences between the two and assess the fraction of diffuse emission in galaxies.

preprint2010arXiv

Cool gas and dust in M33: Results from the Herschel M33 extended survey (HERM33ES)

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 $σ$ using SPIRE and recent CO and \HI\ observations; a variation in $σ$ 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.

preprint2010arXiv

PACS and SPIRE photometer maps of M33: First results of the Herschel M33 extended survey (HERM33ES)

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.

preprint2010arXiv

Properties of compact 250 μm emission and HII regions in M33 (HERM33ES)

Within the framework of the HERM33ES Key Project, using the high resolution and sensitivity of the Herschel photometric data, we study the compact emission in the Local Group spiral galaxy M33 to investigate the nature of the compact SPIRE emission sources. We extracted a catalogue of sources at 250um in order to investigate the nature of this compact emission. Taking advantage of the unprecedented Herschel resolution at these wavelengths, we also focus on a more precise study of some striking Halpha shells in the northern part of the galaxy. We present a catalogue of 159 compact emission sources in M33 identified by SExtractor in the 250um SPIRE band that is the one that provides the best spatial resolution. We also measured fluxes at 24um and Halpha for those 159 extracted sources. The morphological study of the shells also benefits from a multiwavelength approach including Halpha, far-UV from GALEX, and infrared from both Spitzer IRAC 8um and MIPS 24um in order to make comparisons. For the 159 compact sources selected at 250um, we find a very strong Pearson correlation coefficient with the MIPS 24um emission (r24 = 0.94) and a rather strong correlation with the Halpha emission, although with more scatter (rHa = 0.83). The morphological study of the Halpha shells shows a displacement between far-ultraviolet, Halpha, and the SPIRE bands. The cool dust emission from SPIRE clearly delineates the Halpha shell structures. The very strong link between the 250um compact emission and the 24um and Halpha emissions, by recovering the star formation rate from standard recipes for HII regions, allows us to provide star formation rate calibrations based on the 250um compact emission alone. The different locations of the Halpha and far-ultraviolet emissions with respect to the SPIRE cool dust emission leads to a dynamical age of a few Myr for the Halpha shells and the associated cool dust.

preprint2010arXiv

Star Formation in Collision Debris: Insights from the modeling of their Spectral Energy Distribution

During galaxy-galaxy interactions, massive gas clouds can be injected into the intergalactic medium which in turn become gravitationally bound, collapse and form stars, star clusters or even dwarf galaxies. The objects resulting from this process are both "pristine", as they are forming their first generation of stars, and chemically evolved because the metallicity inherited from their parent galaxies is high. Such characteristics make them particularly interesting laboratories to study star formation. After having investigated their star-forming properties, we use photospheric, nebular and dust modeling to analyze here their spectral energy distribution (SED) from the far-ultraviolet to the mid-infrared regime for a sample of 7 star-forming regions. Our analysis confirms that the intergalactic star forming regions in Stephan's Quintet, around Arp 105, and NGC 5291, appear devoid of stellar populations older than 10^9 years. We also find an excess of light in the near-infrared regime (from 2 to 4.5 microns) which cannot be attributed to stellar photospheric or nebular contributions. This excess is correlated with the star formation rate intensity suggesting that it is probably due to emission by very small grains fluctuating in temperature as well as the polycyclic aromatic hydrocarbons (PAH) line at 3.3 micron. Comparing the attenuation via the Balmer decrement to the mid-infrared emission allows us to check the reliability of the attenuation estimate. It suggests the presence of embedded star forming regions in NGC 5291 and NGC 7252. Overall the SED of star-forming regions in collision debris (and Tidal Dwarf Galaxies) resemble more that of dusty star-forming regions in galactic disks than to that of typical star-forming dwarf galaxies.

preprint2010arXiv

Total Infrared Luminosity Estimation of Resolved and Unresolved Galaxies

The total infrared (TIR) luminosity from galaxies can be used to examine both star formation and dust physics. We provide here new relations to estimate the TIR luminosity from various Spitzer bands, in particular from the 8 micron and 24 micron bands. To do so, we use 45" subregions within a subsample of nearby face-on spiral galaxies from the Spitzer Infrared Nearby Galaxies Survey (SINGS) that have known oxygen abundances as well as integrated galaxy data from the SINGS, the Local Volume Legacy Survey (LVL) and Engelbracht et al. (2008) samples. Taking into account the oxygen abundances of the subregions, the star formation rate intensity, and the relative emission of the polycyclic aromatic hydrocarbons at 8 micron, the warm dust at 24 micron and the cold dust at 70 micron and 160 micron we derive new relations to estimate the TIR luminosity from just one or two of the Spitzer bands. We also show that the metallicity and the star formation intensity must be taken into account when estimating the TIR luminosity from two wave bands, especially when data longward of 24 micron are not available.

preprint2009arXiv

Star-Forming or Starbursting? The Ultraviolet Conundrum

Compared to starburst galaxies, normal star forming galaxies have been shown to display a much larger dispersion of the dust attenuation at fixed reddening through studies of the IRX-beta diagram (the IR/UV ratio "IRX" versus the UV color "beta"). To investigate the causes of this larger dispersion and attempt to isolate second parameters, we have used GALEX UV, ground-based optical, and Spitzer infrared imaging of 8 nearby galaxies, and examined the properties of individual UV and 24 micron selected star forming regions. We concentrated on star-forming regions, in order to isolate simpler star formation histories than those that characterize whole galaxies. We find that 1) the dispersion is not correlated with the mean age of the stellar populations, 2) a range of dust geometries and dust extinction curves are the most likely causes for the observed dispersion in the IRX-beta diagram 3) together with some potential dilution of the most recent star-forming population by older unrelated bursts, at least in the case of star-forming regions within galaxies, 4) we also recover some general characteristics of the regions, including a tight positive correlation between the amount of dust attenuation and the metal content. Although generalizing our results to whole galaxies may not be immediate, the possibility of a range of dust extinction laws and geometries should be accounted for in the latter systems as well.