Source author record

V. Lebouteiller

V. Lebouteiller 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

49works
7topics
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

49 published item(s)

preprint2026arXiv

Calibrating Mid-Infrared Emission Features As Diagnostics of Star Formation in Infrared-Luminous Galaxies via Radiative Transfer Modeling

Luminous infrared galaxies are key sites of obscured stellar mass assembly at z > 0.5. Their star formation rates (SFRs) are often estimated using the luminosities of the 6.2 micron and 11.2 micron polycyclic aromatic hydrocarbon (PAH) features, or those of the [Ne II] and [Ne III] fine-structure lines, as they are minimally affected by obscuration. It is uncertain whether the calibration of these features as SFR tracers depends on the starburst bolometric luminosity or the level of Active Galactic Nucleus (AGN) activity. We here investigate the relationship between the luminosities of PAH and Neon lines with star formation rate for highly luminous objects using radiative transfer modeling and archival observations of 42 local Ultraluminous (>= 10^12 L_sun) Infrared Galaxies (ULIRGs). We find that PAH and [Ne II] features arise mainly in star-forming regions, with small contributions from the AGN or host, but that the [Ne III] line has a mixed contribution from both star formation and AGN activity. We present relations between L_PAH and L_NeII, and both starburst luminosity and SFR. We find relations for lower luminosity (L_IR ~= 10^10-10^12 L_sun) systems underestimate the SFRs in local ULIRGs by up to ~1 dex. The 6.2 micron and 11.2 micron PAH features, and the [Ne II] line, are thus good tracers of SFR in ULIRGs. We do not find that a more luminous AGN affects the relationship between SFR and PAH or Neon luminosity, but that it can make PAH emission harder to discern. Our results and derived relations are relevant to studies of star-forming and composite galaxies at z < 3 with the James Webb Space Telescope.

preprint2022arXiv

METAL: The Metal Evolution, Transport, and Abundance in the Large Magellanic Cloud Hubble program. III. Interstellar Depletions, Dust-to-Metal, and Dust-to-Gas Ratios Versus Metallicity

The metallicity and gas density dependence of interstellar depletions, the dust-to-gas (D/G), and dust-to-metal (D/M) ratios have important implications for how accurately we can trace the chemical enrichment of the universe; either by using FIR dust emission as a tracer of the ISM; or by using spectroscopy of damped Lyman-alpha systems (DLAs) to measure chemical abundances over a wide range of redshifts. We collect and compare large samples of depletion measurements in the Milky Way (MW), LMC (Z=0.5 Zsun), and SMC (Z=0.2 Zsun). The relation between the depletions of different elements do not strongly vary between the three galaxies, implying that abundance ratios should trace depletions accurately down to 20% solar metallicity. From the depletions, we derive D/G and D/M. The D/G increases with density, consistent with the more efficient accretion of gas-phase metals onto dust grains in the denser ISM. For log N(H) > 21 cm^-2, the depletion of metallicity tracers (S, Zn) exceeds -0.5 dex, even at 20% solar metallicity. The gas fraction of metals increases from the MW to the LMC (factor 3) and SMC (factor 6), compensating the reduction in total heavy element abundances and resulting in those three galaxies having the same neutral gas-phase metallicities. The D/G derived from depletions are a factor of 2 (LMC) and 5 (SMC) higher than the D/G derived from FIR, 21 cm, and CO emission, likely due to the combined uncertainties on the dust FIR opacity and on the depletion of carbon and oxygen.

preprint2022arXiv

The Infrared Database of Extragalactic Observables from Spitzer. -- II. The Database & The Diagnostic Power of Crystalline Silicate Features in Galaxy Spectra

We present the Infrared Database of Extragalactic Observables from Spitzer (IDEOS), a homogeneous, publicly available, database of 77 fitted mid-infrared observables in the 5.4-36um range, comprising measurements for 3335 galaxies observed in the low-resolution staring mode of the Infrared Spectrometer onboard the Spitzer Space Telescope. Among the included observables are PAH fluxes and their equivalent widths, the strength of the 9.8um silicate feature, emission line fluxes, solid-state features, rest frame continuum fluxes, synthetic photometry, and a mid-infrared spectral classification. The IDEOS spectra were selected from the Cornell Atlas of Spitzer-IRS Sources. To our surprise we have detected at a >95% confidence level crystalline silicates in the spectra of 786 IDEOS galaxies. The detections range from single band detections to detections of all fitted crystalline bands (16, 19, 23, 28 and 33um). We find the strength of the crystalline silicate bands to correlate with the amorphous silicate strength, and the change from an emission to an absorption feature to occur at higher obscuration as the wavelength of the crystalline silicate band is longer. These observed characteristics are consistent with an origin for the amorphous and crystalline silicate features in a centrally heated dust geometry, either an edge-on disk or a cocoon. We find the 23 and 33um crystalline silicate bands to be well-suited to classify the obscuration level of galactic nuclei, even in the presence of strong circumnuclear star formation. Based on our detection statistics, we conclude that crystalline silicates are a common component of the interstellar medium of galactic nuclei.

preprint2022arXiv

Trends in Silicates in the $β$ Pictoris Disk

While beta Pic is known to host silicates in ring-like structures, whether the properties of these silicate dust vary with stellocentric distance remains an open question. We re-analyze the beta Pictoris debris disk spectrum from the Spitzer Infrared Spectrograph (IRS) and a new IRTF/SpeX spectrum to investigate trends in Fe/Mg ratio, shape, and crystallinity in grains as a function of wavelength, a proxy for stellocentric distance. By analyzing a re-calibrated and re-extracted spectrum, we identify a new 18 micron forsterite emission feature and recover a 23 micron forsterite emission feature with a substantially larger line-to-continuum ratio than previously reported. We find that these prominent spectral features are primarily produced by small submicron-sized grains, which are continuously generated and replenished from planetesimal collisions in the disk and can elucidate their parent bodies' composition. We discover three trends about these small grains: as stellocentric distance increases, (1) small silicate grains become more crystalline (less amorphous), (2) they become more irregular in shape, and (3) for crystalline silicate grains, the Fe/Mg ratio decreases. Applying these trends to beta Pic's planetary architecture, we find that the dust population exterior to the orbits of beta Pic b and c differs substantially in crystallinity and shape. We also find a tentative 3-5 micron dust excess due to spatially unresolved hot dust emission close to the star. From our findings, we infer that the surfaces of large planetesimals are more Fe-rich and collisionally-processed closer to the star but more Fe-poor and primordial farther from the star.

preprint2016arXiv

A milestone toward understanding PDR properties in the extreme environment of LMC-30Dor

More complete knowledge of galaxy evolution requires understanding the process of star formation and interaction between the interstellar radiation field and the interstellar medium in galactic environments traversing a wide range of physical parameter space. Here we focus on the impact of massive star formation on the surrounding low metallicity ISM in 30 Doradus in the Large Magellanic Cloud. A low metal abundance, as is the case of some galaxies of the early universe, results in less ultra-violet shielding for the formation of the molecular gas necessary for star formation to proceed. The half-solar metallicity gas in this region is strongly irradiated by the super star cluster R136, making it an ideal laboratory to study the structure of the ISM in an extreme environment. Our spatially resolved study investigates the gas heating and cooling mechanisms, particularly in the photo-dissociation regions where the chemistry and thermal balance are regulated by far-ultraviolet photons (6 eV< hν<13.6 eV). We present Herschel observations of far-infrared fine-structure lines obtained with PACS and SPIRE/FTS. We have combined atomic fine-structure lines from Herschel and Spitzer observations with ground-based CO data to provide diagnostics on the properties and the structure of the gas by modeling it with the Meudon PDR code. We derive the spatial distribution of the radiation field, the pressure, the size, and the filling factor of the photodissociated gas and molecular clouds. We find a range of pressure of ~ 10^5 - 1.7x10^6 cm^{-3} K and a range of incident radiation field G_UV ~ 10^2 - 2.5x10^4 through PDR modeling. Assuming a plane-parallel geometry and a uniform medium, we find a total extinction of 1-3 mag , which correspond to a PDR cloud size of 0.2 to 3pc, with small CO depth scale of 0.06 to 0.5pc. We also determine the three dimensional structure of the gas. (Abridged)

preprint2015arXiv

CASSIS: The Cornell Atlas of Spitzer/Infrared Spectrograph Sources. II. High-resolution observations

The Infrared Spectrograph (IRS) on board the Spitzer Space Telescope observed about 15,000 objects during the cryogenic mission lifetime. Observations provided low-resolution (R~60-127) spectra over ~5-38um and high-resolution (R~600) spectra over ~10-37um. The Cornell Atlas of Spitzer/IRS Sources (CASSIS) was created to provide publishable quality spectra to the community. Low-resolution spectra have been available in CASSIS since 2011, and we present here the addition of the high-resolution spectra. The high-resolution observations represent approximately one third of all staring observations performed with the IRS instrument. While low-resolution observations are adapted to faint objects and/or broad spectral features (e.g., dust continuum, molecular bands), high-resolution observations allow more accurate measurements of narrow features (e.g., ionic emission lines) as well as a better sampling of the spectral profile of various features. Given the narrow aperture of the two high-resolution modules, cosmic ray hits and spurious features usually plague the spectra. Our pipeline is designed to minimize these effects through various improvements. A super sampled point-spread function was created in order to enable the optimal extraction in addition to the full aperture extraction. The pipeline selects the best extraction method based on the spatial extent of the object. For unresolved sources, the optimal extraction provides a significant improvement in signal-to-noise ratio over a full aperture extraction. We have developed several techniques for optimal extraction, including a differential method that eliminates low-level rogue pixels (even when no dedicated background observation was performed). The updated CASSIS repository now includes all the spectra ever taken by the IRS, with the exception of mapping observations.

preprint2015arXiv

From forced collapse to H ii region expansion in Mon R2: Envelope density structure and age determination with Herschel

The surroundings of HII regions can have a profound influence on their development, morphology, and evolution. This paper explores the effect of the environment on H II regions in the MonR2 molecular cloud. We aim to investigate the density structure of envelopes surrounding HII regions and to determine their collapse and ionisation expansion ages. The Mon R2 molecular cloud is an ideal target since it hosts an H II region association. Column density and temperature images derived from Herschel data were used together to model the structure of HII bubbles and their surrounding envelopes. The resulting observational constraints were used to follow the development of the Mon R2 ionised regions with analytical calculations and numerical simulations. The four hot bubbles associated with H II regions are surrounded by dense, cold, and neutral gas envelopes. The radial density profiles are reminiscent of those of low-mass protostellar envelopes. The inner parts of envelopes of all four HII regions could be free-falling because they display shallow density profiles. As for their outer parts, the two compact HII regions show a density profile, which is typical of the equilibrium structure of an isothermal sphere. In contrast, the central UCHii region shows a steeper outer profile, that could be interpreted as material being forced to collapse. The size of the heated bubbles, the spectral type of the irradiating stars, and the mean initial neutral gas density are used to estimate the ionisation expansion time, texp, 0.1Myr,for the dense UCHII and compact HII regions and 0.35 Myr for the extended HII region. The envelope transition radii between the shallow and steeper density profiles are used to estimate the time elapsed since the formation of the first proto stellar embryo, Tinf : 1Myr, for the ultra-compact, 1.5 / 3Myr for the compact, and greater than 6Myr for the extended HII regions.

preprint2015arXiv

Linking dust emission to fundamental properties in galaxies: The low-metallicity picture

In this work, we aim at providing a consistent analysis of the dust properties from metal-poor to metal-rich environments by linking them to fundamental galactic parameters. We consider two samples of galaxies: the Dwarf Galaxy Survey (DGS) and KINGFISH, totalling 109 galaxies, spanning almost 2 dex in metallicity. We collect infrared (IR) to submillimetre (submm) data for both samples and present the complete data set for the DGS sample. We model the observed spectral energy distributions (SED) with a physically-motivated dust model to access the dust properties. Using a different SED model (modified blackbody), dust composition (amorphous carbon), or wavelength coverage at submm wavelengths results in differences in the dust mass estimate of a factor two to three, showing that this parameter is subject to non-negligible systematic modelling uncertainties. For eight galaxies in our sample, we find a rather small excess at 500 microns (< 1.5 sigma). We find that the dust SED of low-metallicity galaxies is broader and peaks at shorter wavelengths compared to more metal-rich systems, a sign of a clumpier medium in dwarf galaxies. The PAH mass fraction and the dust temperature distribution are found to be driven mostly by the specific star-formation rate, SSFR, with secondary effects from metallicity. The correlations between metallicity and dust mass or total-IR luminosity are direct consequences of the stellar mass-metallicity relation. The dust-to-stellar mass ratios of metal-rich sources follow the well-studied trend of decreasing ratio for decreasing SSFR. The relation is more complex for highly star-forming low-metallicity galaxies and depends on the chemical evolutionary stage of the source (i.e., gas-to-dust mass ratio). Dust growth processes in the ISM play a key role in the dust mass build-up with respect to the stellar content at high SSFR and low metallicity. (abridged)

preprint2015arXiv

Star-formation rates from young-star counts and the structure of the ISM across the NGC346/N66 complex in the SMC

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

preprint2015arXiv

The dust properties and physical conditions of the interstellar medium in the LMC massive star forming complex N11

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

preprint2015arXiv

The Herschel Dwarf Galaxy Survey: I. Properties of the low-metallicity ISM from PACS spectroscopy

The far-infrared (FIR) lines are key tracers of the physical conditions of the interstellar medium (ISM) and are becoming workhorse diagnostics for galaxies throughout the universe. Our goal is to explain the differences and trends observed in the FIR line emission of dwarf galaxies compared to more metal-rich galaxies. We present Herschel PACS spectroscopic observations of the CII157um, OI63 and 145um, OIII88um, NII122 and 205um, and NIII57um fine-structure cooling lines in a sample of 48 low-metallicity star-forming galaxies of the guaranteed time key program Dwarf Galaxy Survey. We correlate PACS line ratios and line-to-LTIR ratios with LTIR, LTIR/LB, metallicity, and FIR color, and interpret the observed trends in terms of ISM conditions and phase filling factors with Cloudy radiative transfer models. We find that the FIR lines together account for up to 3 percent of LTIR and that star-forming regions dominate the overall emission in dwarf galaxies. Compared to metal-rich galaxies, the ratios of OIII/NII122 and NIII/NII122 are high, indicative of hard radiation fields. In the photodissociation region (PDR), the CII/OI63 ratio is slightly higher than in metal-rich galaxies, with a small increase with metallicity, and the OI145/OI63 ratio is generally lower than 0.1, demonstrating that optical depth effects should be small on the scales probed. The OIII/OI63 ratio can be used as an indicator of the ionized gas/PDR filling factor, and is found ~4 times higher in the dwarfs than in metal-rich galaxies. The high CII/LTIR, OI/LTIR, and OIII/LTIR ratios, which decrease with increasing LTIR and LTIR/LB, are interpreted as a combination of moderate FUV fields and low PDR covering factor. Harboring compact phases of low filling factor and a large volume filling factor of diffuse gas, the ISM of low-metallicity dwarf galaxies has a more porous structure than that in metal-rich galaxies.

preprint2015arXiv

The Spatial Variation of the Cooling Lines in the Reflection Nebula NGC7023

Context: The north-west photo-dissociation region (PDR) in the reflection nebula NGC 7023 displays a complex structure. Filament-like condensations at the edge of the cloud can be traced via the emission of the main cooling lines, offering a great opportunity to study the link between the morphology and energetics of these regions. Aims: We study the spatial variation of the far-infrared fine-structure lines of [C II] (158 um) and [O I] (63 and 145 um). These lines trace the local gas conditions across the PDR. Methods: We used observations from the Herschel/PACS instrument to map the spatial distribution of these fine-structure lines. The observed region covers a square area of about 110" x 110" with an angular resolution that varies from 4" to 11". We compared this emission with ground-based and Spitzer observations of H2 lines, Herschel/SPIRE observations of CO lines, and Spitzer/IRAC 3.6 um images that trace the emission of polycyclic aromatic hydrocarbons. Results: The [C II] (158 um) and [O I] (63 and 145 um) lines arise from the warm cloud surface where the PDR is located and the gas is warm, cooling the region. We find that although the relative contribution to the cooling budget over the observed region is dominated by [O I]63 um (>30%), H2 contributes significantly in the PDR (35%), as does [C II]158 um outside the PDR (30%). Other species contribute little to the cooling ([O I]145 um 9%, and CO 4%). The [O I] maps resolve these condensations into two structures and show that the peak of [O I] is slightly displaced from the molecular H2 emission. The size of these structures is about 8" (0.015 pc) and in surface cover about 9% of the PDR emission. Finally, we did not detect emission from [N II]122 um, suggesting that the cavity is mostly filled with non-ionised gas.

preprint2014arXiv

A stubbornly large mass of cold dust in the ejecta of Supernova 1987A

We present new Herschel photometric and spectroscopic observations of Supernova 1987A, carried out in 2012. Our dedicated photometric measurements provide new 70 micron data and improved imaging quality at 100 and 160 micron compared to previous observations in 2010. Our Herschel spectra show only weak CO line emission, and provide an upper limit for the 63 micron [O I] line flux, eliminating the possibility that line contaminations distort the previously estimated dust mass. The far-infrared spectral energy distribution (SED) is well fitted by thermal emission from cold dust. The newly measured 70 micron flux constrains the dust temperature, limiting it to nearly a single temperature. The far-infrared emission can be fitted by 0.5+-0.1 Msun of amorphous carbon, about a factor of two larger than the current nucleosynthetic mass prediction for carbon. The observation of SiO molecules at early and late phases suggests that silicates may also have formed and we could fit the SED with a combination of 0.3 Msun of amorphous carbon and 0.5 Msun of silicates, totalling 0.8 Msun of dust. Our analysis thus supports the presence of a large dust reservoir in the ejecta of SN 1987A. The inferred dust mass suggests that supernovae can be an important source of dust in the interstellar medium, from local to high-redshift galaxies.

preprint2014arXiv

The molecular gas reservoir of 6 low-metallicity galaxies from the Herschel Dwarf Galaxy Survey: A ground-based follow-up survey of CO(1-0), CO(2-1), and CO(3-2)

We aim to quantify the molecular gas reservoir in a subset of 6 low-metallicity galaxies from the Herschel Dwarf Galaxy Survey with newly acquired CO data, and link this reservoir to the observed star formation activity. We present CO(1-0), CO(2-1), and CO(3-2) observations obtained at the ATNF Mopra 22-m, APEX, and IRAM 30-m telescopes, as well as [CII] 157um and [OI] 63um observations obtained with the Herschel/PACS spectrometer in the 6 galaxies: Haro11, Mrk1089, Mrk930, NGC4861, NGC625, and UM311. We derive molecular gas mass from several methods including the use of the CO-to-H2 conversion factor Xco (both Galactic and metallicity-scaled values) and of dust measurements. The molecular and atomic gas reservoirs are compared to the star formation activity. We also constrain the physical conditions of the molecular clouds using the non-LTE code RADEX and the spectral synthesis code Cloudy. We detect CO in 5 of the 6 galaxies, including first detections in Haro11 (Z~0.4 Zsun), Mrk930 (0.2 Zsun), and UM311 (0.5 Zsun), but CO remains undetected in NGC4861 (0.2 Zsun). The CO luminosities are low while [CII] is bright in these galaxies, resulting in [CII]/CO(1-0)>10000. Our dwarf galaxies are in relatively good agreement with the Schmidt-Kennicutt relation for total gas. They show short molecular depletion time scales, even when considering metallicity-scaled Xco factors. Those galaxies are dominated by their HI gas, except Haro11 which has high star formation efficiency and is dominated by ionized and molecular gas. We determine the mass of each ISM phase in Haro11 using Cloudy and estimate an equivalent Xco factor which is 10 times higher than the Galactic value. Overall, our results confirm the emerging picture that CO suffers from significant selective photodissociation in low-metallicity dwarf galaxies.

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

A multiwavelength study of the Magellanic-type galaxy NGC 4449 - I. Modelling the spectral energy distribution, the ionization structure and the star formation history

[Abridged] We present an integrated photometric spectral energy distribution (SED) of the Magellanic-type galaxy NGC 4449 from the far-ultraviolet (UV) to the submillimetre, including new observations acquired by the Herschel Space Observatory. We include integrated UV photometry from the Swift Ultraviolet and Optical Telescope using a measurement technique which is appropriate for extended sources with coincidence loss. In this paper, we examine the available multiwavelength data to infer a range of ages, metallicities and star formation rates for the underlying stellar populations, as well as the composition and the total mass of dust in NGC 4449. We present an iterative scheme, which allows us to build an in-depth and multicomponent representation of NGC 4449 `bottom-up', taking advantage of the broad capabilities of the photoionization and radiative transfer code MOCASSIN (MOnte CArlo SimulationS of Ionized Nebulae). We fit the observed SED, the global ionization structure and the emission line intensities, and infer a recent SFR of 0.4 Msolar/yr and a total stellar mass of approximately 1e9 Msolar emitting with a bolometric luminosity of 5.7e9 Lsolar. Our fits yield a total dust mass of 2.9e6 Msolar including 2 per cent attributed to polycyclic aromatic hydrocarbons. We deduce a dust to gas mass ratio of 1/190 within the modelled region. While we do not consider possible additional contributions from even colder dust, we note that including the extended HI envelope and the molecular gas is likely to bring the ratio down to as low as ~ 1/800.

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

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

Diagnostics of AGN-driven Molecular Outflows in ULIRGs from Herschel-PACS Observations of OH at 119um

We report on our observations of the 79 and 119um doublet transitions of OH for 24 local (z<0.262) ULIRGs observed with Herschel-PACS as part of the Herschel ULIRG Survey (HERUS). Some OH119 profiles display a clear P-Cygni shape and therefore imply outflowing OH gas, other profiles are predominantly in absorption or are completely in emission. We find that the relative strength of the OH emission component decreases as the silicate absorption increases. This locates the OH outflows inside the obscured nuclei. The maximum outflow velocities for our sources range from less than 100 to 2000 km/s, with 15/24 (10/24) sources showing OH absorption at velocities exceeding 700 km/s (1000 km/s). Three sources show maximum OH outflow velocities exceeding that of Mrk231. Since outflow velocities above 500-700 km/s are thought to require an active galactic nucleus (AGN) to drive them, about 2/3 of our ULIRG sample may host AGN-driven molecular outflows. This finding is supported by the correlation we find between the maximum OH outflow velocity and the IR-derived bolometric AGN luminosity. No such correlation is found with the IR-derived star formation rate. The highest outflow velocities are found among sources which are still deeply embedded. We speculate that the molecular outflows in these sources may be in an early phase of disrupting the nuclear dust veil before these sources evolve into less obscured AGN. Four of our sources show high-velocity wings in their [C II] fine-structure line profiles implying neutral gas outflow masses of at least 2-4.5 x 10^8 Msun.

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

Herschel observations and a model for IRAS 08572+3915: a candidate for the most luminous infrared galaxy in the local (z < 0.2) Universe

We present Herschel photometry and spectroscopy, carried out as part of the Herschel ULIRG survey (HERUS), and a model for the infrared to submillimetre emission of the ultraluminous infrared galaxy IRAS 08572+3915. This source shows one of the deepest known silicate absorption features and no polycyclic aromatic hydrocarbon (PAH) emission. The model suggests that this object is powered by an active galactic nucleus (AGN) with a fairly smooth torus viewed almost edge-on and a very young starburst. According to our model the AGN contributes about 90% of the total luminosity of 1.1 x 10^13 Lo, which is about a factor of five higher than previous estimates. The large correction of the luminosity is due to the anisotropy of the emission of the best fit torus. Similar corrections may be necessary for other local and high-z analogs. This correction implies that IRAS 08572+3915 at a redshift of 0.05835 may be the nearest hyperluminous infrared galaxy and probably the most luminous infrared galaxy in the local (z < 0.2) Universe. IRAS 08572+3915 shows a low ratio of [CII] to IR luminosity (log L_[CII]/L_{IR} < -3.8) and a [OI]63um to [CII]158um line ratio of about 1 that supports the model presented in this letter.

preprint2013arXiv

Learning the Fundamental MIR Spectral Components of Galaxies with Non-Negative Matrix Factorisation

The mid-infrared (MIR) spectra observed with the \textit{Spitzer} Infrared Spectrograph (IRS) provide a valuable dataset for untangling the physical processes and conditions within galaxies. This paper presents the first attempt to blindly learn fundamental spectral components of MIR galaxy spectra, using non-negative matrix factorisation (NMF). NMF is a recently developed multivariate technique shown to be successful in blind source separation problems. Unlike the more popular multivariate analysis technique, principal component analysis, NMF imposes the condition that weights and spectral components are non-negative. This more closely resembles the physical process of emission in the mid-infrared, resulting in physically intuitive components. By applying NMF to galaxy spectra in the Cornell Atlas of Spitzer/IRS sources (CASSIS), we find similar components amongst different NMF sets. These similar components include two for AGN emission and one for star formation. [... ABBREVIATED...] We show an NMF set with seven components can reconstruct the general spectral shape of a wide variety of objects, though struggle to fit the varying strength of emission lines. We also show that the seven components can be used to separate out different types of objects. We model this separation with Gaussian Mixtures modelling and use the result to provide a classification tool. We also show the NMF components can be used to separate out the emission from AGN and star formation regions and define a new star formation/AGN diagnostic which is consistent with all mid-infrared diagnostics already in use but has the advantage that it can be applied to mid-infrared spectra with low signal to noise or with limited spectral range. The 7 NMF components and code for classification are made public on arxiv and are available at: \url{https://github.com/pdh21/NMF_software/}

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

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

The thermal dust emission in the N158-N159-N160 (LMC) star forming complex mapped by Spitzer, Herschel and LABOCA

We present a study of the infrared/submm emission of the LMC star forming complex N158-N159-N160. Combining observations from the Spitzer Space Telescope (3.6-70um), the Herschel Space Observatory (100-500um) and LABOCA (870um) allows us to work at the best angular resolution available now for an extragalactic source. We observe a remarkably good correlation between SPIRE and LABOCA emission and resolve the low surface brightnesses emission. We use the Spitzer and Herschel data to perform a resolved Spectral Energy Distribution (SED) modelling of the complex. Using MBB, we derive a global emissivity index beta_c of 1.47. If beta cold is fixed to 1.5, we find an average temperature of 27K. We also apply the Galliano et al. (2011) modelling technique (and amorphous carbon to model carbon dust) to derive maps of the star formation rate, the mean starlight intensity, the fraction of PAHs or the dust mass surface density of the region. We observe that the PAH fraction strongly decreases in the HII regions. This decrease coincides with peaks in the mean radiation field intensity map. The dust surface densities follow the FIR distribution, with a total dust mass of 2.1x10^4 Msolar (2.8 times less than when using graphite grains) in the resolved elements we model. We find a non-negligible amount of dust in the molecular cloud N159 South (showing no massive SF). We also investigate the drivers of the Herschel/PACS and SPIRE submm colours as well as the variations in the gas-to-dust mass ratio (G/D) and the XCO conversion factor in the region N159. We finally model individual regions to analyse variations in the SED shape across the complex and the 870um emission in more details. No measurable submm excess emission at 870um seems to be detected in these regions.

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.

preprint2012arXiv

[CII] 158 micron Luminosities and Star Formation Rate in Dusty Starbursts and AGN

Results are presented for [CII] 158 micron line fluxes observed with the Herschel PACS instrument in 112 sources with both starburst and AGN classifications, of which 102 sources have confident detections. Results are compared with mid-infrared spectra from the Spitzer Infrared Spectrometer and with L(IR) from IRAS fluxes; AGN/starburst classifications are determined from equivalent width of the 6.2 micron PAH feature. It is found that the [CII] line flux correlates closely with the flux of the 11.3 micron PAH feature independent of AGN/starburst classification, log [f([CII] 158 micron)/f(11.3 micron PAH)] = -0.22 +- 0.25. It is concluded that [CII] line flux measures the photodissociation region associated with starbursts in the same fashion as the PAH feature. A calibration of star formation rate for the starburst component in any source having [CII] is derived comparing [CII] luminosity L([CII]) to L(IR) with the result that log SFR = log L([CII)]) - 7.08 +- 0.3, for SFR in solar masses per year and L([CII]) in solar luminosities. The decreasing ratio of L([CII]) to L(IR) in more luminous sources (the "[CII] deficit") is shown to be a consequence of the dominant contribution to L(IR) arising from a luminous AGN component because the sources with largest L(IR) and smallest L([CII])/L(IR) are AGN.

preprint2012arXiv

Globules and pillars seen in the [CII] 158 micron line with SOFIA

Molecular globules and pillars are spectacular features, found only in the interface region between a molecular cloud and an HII-region. Impacting Far-ultraviolet (FUV) radiation creates photon dominated regions (PDRs) on their surfaces that can be traced by typical cooling lines. With the GREAT receiver onboard SOFIA we mapped and spectrally resolved the [CII] 158 micron atomic fine-structure line and the highly excited 12CO J=11-10 molecular line from three objects in Cygnus X (a pillar, a globule, and a strong IRAS source). We focus here on the globule and compare our data with existing Spitzer data and recent Herschel Open-Time PACS data. Extended [CII] emission and more compact CO-emission was found in the globule. We ascribe this emission mainly to an internal PDR, created by a possibly embedded star-cluster with at least one early B-star. However, external PDR emission caused by the excitation by the Cyg OB2 association cannot be fully excluded. The velocity-resolved [CII] emission traces the emission of PDR surfaces, possible rotation of the globule, and high-velocity outflowing gas. The globule shows a velocity shift of ~2 km/s with respect to the expanding HII-region, which can be understood as the residual turbulence of the molecular cloud from which the globule arose. This scenario is compatible with recent numerical simulations that emphazise the effect of turbulence. It is remarkable that an isolated globule shows these strong dynamical features traced by the [CII]-line, but it demands more observational studies to verify if there is indeed an embedded cluster of B-stars.

preprint2012arXiv

Herschel-SPIRE Imaging Spectroscopy of Molecular Gas in M82

We present new Herschel-SPIRE imaging spectroscopy (194-671 microns) of the bright starburst galaxy M82. Covering the CO ladder from J=4-3 to J=13-12, spectra were obtained at multiple positions for a fully sampled ~ 3 x 3 arcminute map, including a longer exposure at the central position. We present measurements of 12CO, 13CO, [CI], [NII], HCN, and HCO+ in emission, along with OH+, H2O+ and HF in absorption and H2O in both emission and absorption, with discussion. We use a radiative transfer code and Bayesian likelihood analysis to model the temperature, density, column density, and filling factor of multiple components of molecular gas traced by 12CO and 13CO, adding further evidence to the high-J lines tracing a much warmer (~ 500 K), less massive component than the low-J lines. The addition of 13CO (and [CI]) is new and indicates that [CI] may be tracing different gas than 12CO. No temperature/density gradients can be inferred from the map, indicating that the single-pointing spectrum is descriptive of the bulk properties of the galaxy. At such a high temperature, cooling is dominated by molecular hydrogen. Photon-dominated region (PDR) models require higher densities than those indicated by our Bayesian likelihood analysis in order to explain the high-J CO line ratios, though cosmic-ray enhanced PDR models can do a better job reproducing the emission at lower densities. Shocks and turbulent heating are likely required to explain the bright high-J emission.

preprint2012arXiv

Oxygen-rich dust production in IC 10

We report the detection of oxygen-rich circumstellar envelopes in stars of the nearby (700 kpc) starburst galaxy IC 10. The star formation history and the chemical environment of this galaxy makes it an ideal target to observe dust production by high-mass stars in a low-metallicity environment. The goal of this study is to identify oxygen-rich stars in IC 10 and to constrain their nature between asymptotic giant branch stars (AGBs), red supergiants (RSGs), and other infrared bright sources. We examine the mass-loss rate of the stars and compare to results obtained for the Magellanic Clouds. Our objectives are to (1) assess whether RSGs can be significant dust producers in IC 10, and (2), solve the discrepancy between the star formation history of IC 10 and the relatively low number of RSGs detected in the optical. We search for silicate dust in emission by using the spectral map observed with the Infrared Spectrograph on board the Spitzer Space Telescope. The optical (UBVRI) and infrared (JHK, Spitzer/IRAC and Spitzer/MIPS) photometry is used to assert the membership of the stars to IC 10 and disentangle between AGBs and RSGs. Radiative models are used to infer mass-loss rates and stellar luminosities. The luminosity and colors of at least 9 silicate emission sources are consistent with stars within IC 10. Furthermore, the photometry of 2 of these sources is consistent with RSGs. We derive dust mass-loss rates similar to the values found in the Magellanic Clouds. Accounting for the sample completeness, RSGs are not important contributors to the dust mass budget in IC 10.

preprint2012arXiv

Physical conditions in the gas phases of the giant HII region LMC-N11 unveiled by Herschel - I. Diffuse [CII] and [OIII] emission in LMC-N11B

(Abridged) The Magellanic Clouds provide a nearby laboratory for metal-poor dwarf galaxies. The low dust abundance enhances the penetration of UV photons into the interstellar medium (ISM), resulting in a relatively larger filling factor of the ionized gas. Furthermore, there is likely a hidden molecular gas reservoir probed by the [CII]157um line. We present Herschel/PACS maps in several tracers, [CII], [OI]63um,145um, [NII]122um, [NIII]57um, and [OIII]88um in the HII region N11B in the Large Magellanic Cloud. Halpha and [OIII]5007A images were used as complementary data to investigate the effect of dust extinction. Observations were interpreted with photoionization models to infer the gas conditions and estimate the ionized gas contribution to the [CII] emission. Photodissociation regions (PDRs) are probed through polycyclic aromatic hydrocarbons (PAHs). We first study the distribution and properties of the ionized gas. We then constrain the origin of [CII]157um by comparing to tracers of the low-excitation ionized gas and of PDRs. [OIII] is dominated by extended emission from the high-excitation diffuse ionized gas; it is the brightest far-infrared line, ~4 times brighter than [CII]. The extent of the [OIII] emission suggests that the medium is rather fragmented, allowing far-UV photons to permeate into the ISM to scales of >30pc. Furthermore, by comparing [CII] with [NII], we find that 95% of [CII] arises in PDRs, except toward the stellar cluster for which as much as 15% could arise in the ionized gas. We find a remarkable correlation between [CII]+[OI] and PAH emission, with [CII] dominating the cooling in diffuse PDRs and [OI] dominating in the densest PDRs. The combination of [CII] and [OI] provides a proxy for the total gas cooling in PDRs. Our results suggest that PAH emission describes better the PDR gas heating as compared to the total infrared 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 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 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.

preprint2012arXiv

The nature of the interstellar medium of the starburst low-metallicity galaxy Haro11: a multi-phase model of the infrared emission

(abridged) Our goal is to describe the multi-phase ISM of the IR bright low-metallicity galaxy Haro 11, dissecting the photoionised and photodissociated gas components. We present observations of the mid- and far-IR fine-structure cooling lines obtained with the Spitzer/IRS and Herschel/PACS spectrometers. We use the spectral synthesis code Cloudy to methodically model the ionised and neutral gas from which these lines originate. We find that the mid- and far-IR lines account for ~1% of the total IR luminosity L_TIR. Haro 11 is undergoing a phase of intense star formation, as traced by the brightest line [OIII] 88um, with L_[OIII]/L_TIR ~0.3%, and high ratios of [NeIII]/[NeII] and [SIV]/[SIII]. Due to their different origins, the observed lines require a multi-phase modeling comprising: a compact HII region, dense fragmented photodissociation regions (PDRs), a diffuse extended low-ionisation/neutral gas which has a volume filling factor of at least 90%, and porous warm dust in proximity to the stellar source. For a more realistic picture of the ISM of Haro 11 we would need to model the clumpy source and gas structures. We combine these 4 model components to explain the emission of 17 spectral lines, investigate the global energy balance of the galaxy through its spectral energy distribution, and establish a phase mass inventory. While the ionic emission lines of Haro 11 essentially originate from the dense HII region component, a diffuse low-ionisation gas is needed to explain the [NeII], [NII], and [CII] line intensities. The [OIII] 88um line intensity is not fully reproduced by our model, hinting towards the possible presence of yet another low-density high-ionisation medium. The [OI] emission is consistent with a dense PDR of low covering factor, and we find no evidence for an X-ray dominated component. The PDR component accounts for only 10% of the [CII] emission.

preprint2011arXiv

CASSIS: The Cornell Atlas of Spitzer/IRS Sources

We present the spectral atlas of sources observed in low resolution with the Infrared Spectrograph on board the Spitzer Space Telescope. More than 11,000 distinct sources were extracted using a dedicated algorithm based on the SMART software with an optimal extraction (AdOpt package). These correspond to all 13,000 low resolution observations of fixed objects (both single source and cluster observations). The pipeline includes image cleaning, individual exposure combination, and background subtraction. A particular attention is given to bad pixel and outlier rejection at the image and spectra levels. Most sources are spatially unresolved so that optimal extraction reaches the highest possible signal-to-noise ratio. For all sources, an alternative extraction is also provided that accounts for all of the source flux within the aperture. CASSIS provides publishable quality spectra through an online database together with several important diagnostics, such as the source spatial extent and a quantitative measure of detection level. Ancillary data such as available spectroscopic redshifts are also provided. The database interface will eventually provide various ways to interact with the spectra, such as on-the-fly measurements of spectral features or comparisons among spectra.

preprint2011arXiv

Investigations of dust heating in M81, M83, and NGC 2403 with the Herschel Space Observatory

We use Spitzer Space Telescope and Herschel Space Observatory far-infrared data along with ground-based optical and near-infrared data to understand how dust heating in the nearby face-on spiral galaxies M81, M83, and NGC 2403 is affected by the starlight from all stars and by the radiation from star forming regions. We find that 70/160 micron surface brightness ratios tend to be more strongly influenced by star forming regions. However, the 250/350 micron and 350/500 micron surface brightness ratios are more strongly affected by the light from the total stellar populations, suggesting that the dust emission at >250 microns originates predominantly from a component that is colder than the dust seen at <160 microns and that is relatively unaffected by star formation activity. We conclude by discussing the implications of this for modelling the spectral energy distributions of both nearby and more distant galaxies and for using far-infrared dust emission to trace star formation.

preprint2011arXiv

Non-Standard Grain Properties, Dark Gas Reservoir, and Extended Submillimeter Excess, Probed by Herschel in the Large Magellanic Cloud

Aims: In this paper, we perform detailed modelling of the Spitzer and Herschel observations of the LMC, in order to: (i) systematically study the uncertainties and biases affecting dust mass estimates; and to (ii) explore the peculiar ISM properties of the LMC. Methods: To achieve these goals, we have modelled the spatially resolved SEDs with two alternate grain compositions, to study the impact of different submillimetre opacities on the dust mass. We have rigorously propagated the observational errors (noise and calibration) through the entire fitting process, in order to derive consistent parameter uncertainties. Results: First, we show that using the integrated SED leads to underestimating the dust mass by ~50 % compared to the value obtained with sufficient spatial resolution, for the region we studied. This might be the case, in general, for unresolved galaxies. Second, we show that Milky Way type grains produce higher gas-to-dust mass ratios than what seems possible according to the element abundances in the LMC. A spatial analysis shows that this dilemma is the result of an exceptional property: the grains of the LMC have on average a larger intrinsic submm opacity (emissivity index beta~1.7 and opacity kappa_abs(160 microns)=1.6 m2/kg) than those of the Galaxy. By studying the spatial distribution of the gas-to-dust mass ratio, we are able to constrain the fraction of unseen gas mass between ~10, and ~100 % and show that it is not sufficient to explain the gas-to-dust mass ratio obtained with Milky Way type grains. Finally, we confirm the detection of a 500 microns extended emission excess with an average relative amplitude of ~15 %, varying up to 40 %. This excess anticorrelates well with the dust mass surface density. Although we do not know the origin of this excess, we show that it is unlikely the result of very cold dust, or CMB fluctuations.

preprint2010arXiv

Influence of the Environment on PAH Emission in Star-Forming Regions

We investigate the emission properties of polycyclic aromatic hydrocarbons (PAHs) in various metallicity environments with the Infrared Spectrograph on board Spitzer. Local giant HII regions are used as references as they enable access to the distinct interstellar medium components that contribute to the mid-infrared spectrum of star-forming galaxies: photodissociation regions (PDRs), photoionized gas, stellar clusters, and embedded regions. Three objects are considered, NGC3603 in the Milky Way, 30Doradus in the Large Magellanic Cloud, and N66 in the Small Magellanic Cloud. From the variations of the PAH/14um ratio, we find that PAHs are destroyed in the ionized gas for a radiation field such that [NeIII]/[NeII]>3. From the variations of the PAH/Hu-alpha ratio, we find that the PAH emission sources in the giant HII regions follow the same photodestruction law regardless of metallicity. We then compare these results with observations of starburst galaxies, HII galaxies, and blue compact dwarf galaxies (BCDs). While the integrated mid-infrared spectra of BCDs are reminiscent of a warm dusty ionized gas, we observe a significant contribution to the PAH emission in starburst galaxies that is not arising from PDRs.

preprint2010arXiv

The extraordinary mid-infrared spectral properties of FeLoBAL Quasars

We present mid-infrared spectra of six FeLoBAL QSOs at 1<z<1.8, taken with the Spitzer space telescope. The spectra span a range of shapes, from hot dust dominated AGN with silicate emission at 9.7 microns, to moderately obscured starbursts with strong Polycyclic Aromatic Hydrocarbon (PAH) emission. The spectrum of one object, SDSS 1214-0001, shows the most prominent PAHs yet seen in any QSO at any redshift, implying that the starburst dominates the mid-IR emission with an associated star formation rate of order 2700 solar masses per year. With the caveats that our sample is small and not robustly selected, we combine our mid-IR spectral diagnostics with previous observations to propose that FeLoBAL QSOs are at least largely comprised of systems in which (a) a merger driven starburst is ending, (b) a luminous AGN is in the last stages of burning through its surrounding dust, and (c) which we may be viewing over a restricted line of sight range.

preprint2009arXiv

A Spitzer high resolution mid-infrared spectral atlas of starburst galaxies

We present an atlas of Spitzer/IRS high resolution (R~600) 10-37um spectra for 24 well known starburst galaxies. The spectra are dominated by fine-structure lines, molecular hydrogen lines, and emission bands of polycyclic aromatic hydrocarbons. Six out of the eight objects with a known AGN component show emission of the high excitation [NeV] line. This line is also seen in one other object (NGC4194) with, a priori, no known AGN component. In addition to strong polycyclic aromatic hydrocarbon emission features in this wavelength range (11.3, 12.7, 16.4um), the spectra reveal other weak hydrocarbon features at 10.6, 13.5, 14.2um, and a previously unreported emission feature at 10.75um. An unidentified absorption feature at 13.7um is detected in many of the starbursts. We use the fine-structure lines to derive the abundance of neon and sulfur for 14 objects where the HI 7-6 line is detected. We further use the molecular hydrogen lines to sample the properties of the warm molecular gas. Several basic diagrams characterizing the properties of the sample are also shown. We have combined the spectra of all the pure starburst objects to create a high S/N template, which is available to the community.

preprint2009arXiv

Advanced optimal extraction for the Spitzer/IRS

We present new advances in the spectral extraction of point-like sources adapted to the Infrared Spectrograph onboard the Spitzer Space Telescope. For the first time, we created a super-sampled point spread function of the low-resolution modules. We describe how to use the point spread function to perform optimal extraction of a single source and of multiple sources within the slit. We also examine the case of the optimal extraction of one or several sources with a complex background. The new algorithms are gathered in a plugin called Adopt which is part of the SMART data analysis software.

preprint2009arXiv

On the nature of the first galaxies selected at 350 microns

[abridged] We present constraints on the nature of the first galaxies selected at 350 microns. The sample includes galaxies discovered in the deepest blank-field survey at 350 microns (in the Bootes Deep Field) and also later serendipitous detections in the Lockman Hole. Spectral energy distribution templates are fit to identified counterparts, and the sample is found to comprise IR-luminous galaxies at 1<z<3 predominantly powered by star formation. The first spectrum of a 350-micron-selected galaxy provides an additional confirmation, showing prominent dust grain features typically associated with star-forming galaxies. Compared to submillimeter galaxies selected at 850 and 1100 microns, galaxies selected at 350 microns have a similar range of far-infrared color temperatures. However, no 350-micron-selected sources are reliably detected at 850 or 1100 microns. Galaxies in our sample with redshifts 1<z<2 show a tight correlation between the far- and mid-infrared flux densities, but galaxies at higher redshifts show a large dispersion in their mid- to far-infrared colors. The 350-micron flux densities (15<S(350)<40 mJy) place these objects near the Herschel/SPIRE 350-micron confusion threshold, with the lower limit on the star formation rate density suggesting the bulk of the 350-micron contribution will come from less luminous infrared sources and normal galaxies.

preprint2009arXiv

Probing the Excitation of Extreme Starbursts: High Resolution Mid-IR Spectroscopy of Blue Compact Dwarfs

We present an analysis of the mid-infrared emission lines for a sample of 12 low metallicity Blue Compact Dwarf (BCD) galaxies based on high resolution observations obtained with Infrared Spectrograph on board the {\rm Spitzer} Space Telescope. We compare our sample with a local sample of typical starburst galaxies and active galactic nuclei (AGNs), to study the ionization field of starbursts over a broad range of physical parameters and examine its difference from the one produced by AGN. The high-ionization line [OIV]25.89$μ$m is detected in most of the BCDs, starbursts, and AGNs in our sample. We propose a diagnostic diagram of the line ratios [OIV]25.89$μ$m/[SIII]33.48$μ$m as a function of [NeIII]15.56$μ$m/[NeII]12.81$μ$m which can be useful in identifying the principal excitation mechanism in a galaxy. Galaxies in this diagram split naturally into two branches. Classic AGNs as well as starburst galaxies with an AGN component populate the upper branch, with stronger AGNs displaying higher [NeIII]/[NeII] ratios. BCDs and pure starbursts are located in the lower branch. We find that overall the placement of galaxies on this diagram correlates well with their corresponding locations in the log([NII]/H$α$) vs. log([OIII]/H$β$) diagnostic diagram, which has been widely used in the optical. The two diagrams provide consistent classifications of the excitation mechanism in a galaxy. On the other hand, the diagram of [NeIII]15.56$μ$m/[NeII]12.81$μ$m vs. [SIV]10.51$μ$m/[SIII]18.71$μ$m is not as efficient in separating AGNs from BCDs and pure starbursts. (abridged)

preprint2008arXiv

The Mid-Infrared Properties of Blue Compact Dwarf Galaxies

The unprecedented sensitivity of the Spitzer Space Telescope has enabled us for the first time to detect a large sample of Blue Compact Dwarf galaxies (BCDs), which are intrinsically faint in the infrared. In the present paper we present a summary of our findings which providing essential information on the presence/absence of the Polycyclic Aromatic Hydrocarbon features in metal-poor environments. In addition, using Spitzer/IRS high-resolution spectroscopy, we study the elemental abundances of neon and sulfur in BCDs and compare with the results from optical studies. Finally, we present an analysis of the mid- and far-infrared to radio correlation in low luminosity low metallicity galaxies.