Source author record

A. Abergel

A. Abergel 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

45works
9topics
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

45 published item(s)

preprint2025arXiv

JWST observations of photodissociation regions III. Dust modelling at the illuminated edge of the Horsehead PDR

Carbonaceous nano-grains are a significant component of interstellar dust and dominate the mid-infrared emission of photodissociation regions (PDRs). We study the evolution of nano-grains across the illuminated edge of the Horsehead PDR, especially their abundance and size properties. This work is part of the Physics and Chemistry of PDR Fronts program studying dust and gas in PDRs with JWST. We use NIRCam+MIRI photometric bands and NIRSpec+MRS spectroscopy to map the illuminated edge. We model dust emission using the THEMIS dust model with the SOC radiative transfer code. Detailed modeling of high angular resolution JWST data allows us to obtain constraints on nano-grain properties. We find that diffuse ISM dust cannot account for the observed data, requiring evolved grains. A sharp density increase is observed at the illuminated edge, consistent with ALMA observations revealing a sharp transition between molecular and ionized gas. Although the PDR length could not be directly determined, we estimate an upper limit of approximately 0.015 pc. This implies a lower limit on small grain abundance (greater than 0.003), showing small grains are not depleted at the Horsehead edge, unlike in the Orion Bar. Our findings indicate a high-density environment and less steep size distribution for nano-grains at the illuminated edge versus the diffuse ISM. This implies nano-grain destruction mechanisms might be less efficient in the Horsehead's moderate-UV field than in more intense PDRs. These results support a model where nano-grain population recovery is slower in moderate-UV environments, leading to a unique dust size distribution at the edge of the Horsehead Nebula.

preprint2025arXiv

JWST observations of photodissociation regions: II. Warm molecular Hydrogen spectroscopy in the Horsehead nebula

H2 is the most abundant molecule in the interstellar medium and is a useful tool to study photodissociation regions, where radiative feedback from massive stars on molecular clouds is dominant. The James Webb Space Telescope, with its high spatial resolution, sensitivity, and wavelength coverage provides unique access to the detection of most of H2 lines and the analysis of its spatial morphology. Our goal is to use H2 line emission detected with the JWST in the Horsehead nebula to constrain the physical parameters (e.g., extinction, gas temperature, thermal pressure) throughout the PDR and its geometry. The study of H2 morphology reveals that FUV-pumped lines peak closer to the edge of the PDR than thermalized lines. From H2 lines, we estimate the value of extinction throughout the PDR. We find that AV is increasing from the edge of the PDR to the second and third H2 filaments. Then, we study the H2 excitation in different regions across the PDR. The temperature profile shows that the observed gas temperature is quite constant throughout the PDR, with a slight decline in each of the dissociation fronts. This study also reveals that the OPR is far from equilibrium. We observe a spatial separation of para and ortho rovibrational levels, indicating that efficient ortho-para conversion and preferential ortho self-shielding are driving the spatial variations of the OPR. Finally, we derive a thermal pressure in the first filament around P > 6x10$^6$ K cm$^{-3}$, about ten times higher than that of the ionized gas. We highlight that template stationary 1D PDR models cannot account for the intrinsic 2D structure and the very high temperature observed in the Horsehead nebula. We argue the highly excited, over-pressurized H2 gas at the edge of the PDR interface could originate from the mixing between the cold and hot phase induced by the photo-evaporation of the cloud.

preprint2020arXiv

Dust evolution across the Horsehead Nebula

Micro-physical processes on interstellar dust surfaces are tightly connected to dust properties (i.e. dust composition, size and shape) and play a key role in numerous phenomena in the interstellar medium (ISM). The large disparity in physical conditions (i.e. density, gas temperature) in the ISM triggers an evolution of dust properties. The analysis of how dust evolves with the physical conditions is a stepping-stone towards a more thorough understanding of interstellar dust. The aim of this paper is to highlight dust evolution in the Horsehead Nebula PDR region. We use Spitzer/IRAC (3.6, 4.5, 5.8 and 8 μm), Spitzer/MIPS (24 μm) together with Herschel/PACS (70 and 160 μm) and Herschel/SPIRE (250, 350 and 500 μm) to map the spatial distribution of dust in the Horsehead over the entire emission spectral range. We model dust emission and scattering using the THEMIS interstellar dust model together with the 3D radiative transfer code SOC. We find that the nano-grains dust-to-gas ratio in the irradiated outer part of the Horsehead is 6 to 10 times lower than in the diffuse ISM. Their minimum size is 2 to 2.25 times larger than in the diffuse ISM and the power-law exponent of their size distribution, 1.1 to 1.4 times lower than in the diffuse ISM. Regarding the denser part of the Horsehead, it is necessary to use evolved grains (i.e. aggregates, with or without an ice mantle). It is not possible to explain the observations using grains from the diffuse medium. We therefore propose the following scenario to explain our results. In the outer part of the Horsehead, all the nano-grains have not yet had time to re-form completely through photo-fragmentation of aggregates and the smallest of the nano-grains that are sensitive to the radiation field are photo-destroyed. In the inner part of the Horsehead, grains most likely consist of multi-compositional, mantled aggregates.

preprint2016arXiv

Comparison of absolute gain photometric calibration between Planck/HFI and Herschel/SPIRE at 545 and 857 GHz

We compare the absolute gain photometric calibration of the Planck/HFI and Herschel/SPIRE instruments on diffuse emission. The absolute calibration of HFI and SPIRE each relies on planet flux measurements and comparison with theoretical far-infrared emission models of planetary atmospheres. We measure the photometric cross calibration between the instruments at two overlapping bands, 545 GHz / 500 $μ$m and 857 GHz / 350 $μ$m. The SPIRE maps used have been processed in the Herschel Interactive Processing Environment (Version 12) and the HFI data are from the 2015 Public Data Release 2. For our study we used 15 large fields observed with SPIRE, which cover a total of about 120 deg^2. We have selected these fields carefully to provide high signal-to-noise ratio, avoid residual systematics in the SPIRE maps, and span a wide range of surface brightness. The HFI maps are bandpass-corrected to match the emission observed by the SPIRE bandpasses. The SPIRE maps are convolved to match the HFI beam and put on a common pixel grid. We measure the cross-calibration relative gain between the instruments using two methods in each field, pixel-to-pixel correlation and angular power spectrum measurements. The SPIRE / HFI relative gains are 1.047 ($\pm$ 0.0069) and 1.003 ($\pm$ 0.0080) at 545 and 857 GHz, respectively, indicating very good agreement between the instruments. These relative gains deviate from unity by much less than the uncertainty of the absolute extended emission calibration, which is about 6.4% and 9.5% for HFI and SPIRE, respectively, but the deviations are comparable to the values 1.4% and 5.5% for HFI and SPIRE if the uncertainty from models of the common calibrator can be discounted. Of the 5.5% uncertainty for SPIRE, 4% arises from the uncertainty of the effective beam solid angle, which impacts the adopted SPIRE point source to extended source unit conversion factor (Abridged)

preprint2016arXiv

Deep Herschel PACS point spread functions

The knowledge of the point spread function (PSF) of imaging instruments represents a fundamental requirement for astronomical observations. The Herschel PACS PSFs delivered by the instrument control centre are obtained from observations of the Vesta asteroid, providing a characterisation of the central part therefore excluding fainter features. However, in many cases information on both the core and the wings of the PSFs is needed. With this aim, we combine Vesta and Mars dedicated observations and obtain PACS PSFs with an unprecedented dynamic range $(\sim 10^6)$, at slow and fast scan speeds and for the three photometric bands.

preprint2016arXiv

Three-dimensional distribution of hydrogen fluoride gas toward NGC6334 I and I(N)

Aims. We investigate the spatial distribution of a collection of absorbing gas clouds, some associated with the dense, massive star-forming core NGC6334 I, and others with diffuse foreground clouds. For the former category, we aim to study the dynamical properties of the clouds in order to assess their potential to feed the accreting protostellar cores. Methods. We use spectral imaging from the Herschel SPIRE iFTS to construct a map of HF absorption at 243 micron in a 6x3.5 arcmin region surrounding NGC6334 I and I(N). Results. The combination of new, spatially fully sampled, but spectrally unresolved mapping with a previous, single-pointing, spectrally resolved HF signature yields a 3D picture of absorbing gas clouds in the direction of NGC6334. Toward core I, the HF equivalent width matches that of the spectrally resolved observation. The distribution of HF absorption is consistent with three of the seven components being associated with this dense star-forming envelope. For two of the remaining four components, our data suggest that these clouds are spatially associated with the larger scale filamentary star-forming complex. Our data also implies a lack of gas phase HF in the envelope of core I(N). Using a simple description of adsorption onto and desorption from dust grain surfaces, we show that the overall lower temperature of the envelope of source I(N) is consistent with freeze-out of HF, while it remains in the gas phase in source I. Conclusions. We use the HF molecule as a tracer of column density in diffuse gas (n(H) ~ 10^2 - 10^3 cm^-3), and find that it may uniquely trace a relatively low density portion of the gas reservoir available for star formation that otherwise escapes detection. At higher densities prevailing in protostellar envelopes (>10^4 cm^-3), we find evidence of HF depletion from the gas phase under sufficiently cold conditions.

preprint2015arXiv

$Herschel$ SPIRE-FTS observations of RCW 120

The expansion of Galactic HII regions can trigger the formation of a new generation of stars. However, little is know about the physical conditions that prevail in these regions. We study the physical conditions that prevail in specific zones towards expanding HII regions that trace representative media such as the photodissociation region, the ionized region, and condensations with and without ongoing star formation. We use the SPIRE Fourier Transform Spectrometer (FTS) on board $Herschel$ to observe the HII region RCW 120. Continuum and lines are observed in the $190-670\,μ$m range. Line intensities and line ratios are obtained and used as physical diagnostics of the gas. We used the Meudon PDR code and the RADEX code to derive the gas density and the radiation field at nine distinct positions including the PDR surface and regions with and without star-formation activity. For the different regions we detect the atomic lines [NII] at $205\,μ$m and [CI] at $370$ and $609\,μ$m, the $^{12}{\rm CO}$ ladder between the $J=4$ and $J=13$ levels and the $^{13}{\rm CO}$ ladder between the $J=5$ and $J=14$ levels, as well as CH$ ^{+} $ in absorption. We find gas temperatures in the range $45-250\,$K for densities of $10^4-10^6\,{\rm cm}^{-3}$, and a high column density on the order of $N_{\rm H}\sim10^{22}\,{\rm cm}^{-2}$ that is in agreement with dust analysis. The ubiquitousness of the atomic and CH$ ^{+} $ emission suggests the presence of a low-density PDR throughout RCW 120. High-excitation lines of CO indicate the presence of irradiated dense structures or small dense clumps containing young stellar objects, while we also find a less dense medium ($N_{\rm H}\sim10^{20}\,{\rm cm}^{-2}$) with high temperatures ($80-200\,$K).

preprint2015arXiv

Detection of two power-law tails in the probability distribution functions of massive GMCs

We report the novel detection of complex high-column density tails in the probability distribution functions (PDFs) for three high-mass star-forming regions (CepOB3, MonR2, NGC6334), obtained from dust emission observed with Herschel. The low column density range can be fit with a lognormal distribution. A first power-law tail starts above an extinction (Av) of ~6-14. It has a slope of alpha=1.3-2 for the rho~r^-alpha profile for an equivalent density distribution (spherical or cylindrical geometry), and is thus consistent with free-fall gravitational collapse. Above Av~40, 60, and 140, we detect an excess that can be fitted by a flatter power law tail with alpha>2. It correlates with the central regions of the cloud (ridges/hubs) of size ~1 pc and densities above 10^4 cm^-3. This excess may be caused by physical processes that slow down collapse and reduce the flow of mass towards higher densities. Possible are: 1. rotation, which introduces an angular momentum barrier, 2. increasing optical depth and weaker cooling, 3. magnetic fields, 4. geometrical effects, and 5. protostellar feedback. The excess/second power-law tail is closely linked to high-mass star-formation though it does not imply a universal column density threshold for the formation of (high-mass) stars.

preprint2015arXiv

Dust variations in the diffuse interstellar medium: constraints on Milky Way dust from Planck-HFI observations

The Planck-HFI all-sky survey from 353 to 857GHz combined with the 100 microns IRAS show that the dust properties vary in the diffuse ISM at high Galactic latitude (1e19<NH<2.5e20 H/cm2). Our aim is to explain these variations with changes in the ISM properties and grain evolution. Our starting point is the latest core-mantle dust model. It consists of small aromatic-rich carbon grains, larger amorphous carbon grains with aliphatic-rich cores and aromatic-rich mantles, and amorphous silicates with Fe/FeS nano-inclusions covered by aromatic-rich carbon mantles. We explore whether variations in the radiation field or in the gas density distribution in the diffuse ISM could explain the observations. The dust properties are also varied in terms of mantle thickness, Fe/FeS inclusions, carbon abundance, and size distribution. Variations in the radiation field intensity and gas density distribution cannot explain the observed variations but radiation fields harder than the standard ISRF may participate in creating part of them. We further show that variations in the grain mantle thickness coupled with changes in the grain size distribution can reproduce most of the observations. We put a limit on the mantle thickness of the silicates (~10-15nm), and find that aromatic-rich mantles are needed for the carbon grains (at least 5-7.5nm thick). We also find that changes in the carbon abundance in the grains could explain part of the observed variations. Finally, we show that varying the composition of Fe/FeS inclusions in the silicates cannot account for the variations. With small variations in the dust properties, we are able to explain most of the variations in the dust emission observed by Planck-HFI in the diffuse ISM. We also find that the small realistic changes in the dust properties that we consider almost perfectly match the anti-correlation and scatter in the observed beta-T relation.

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

Physical structure of the photodissociation regions in NGC 7023. Observations of gas and dust emission with Herschel

The determination of the physical conditions in molecular clouds is a key step towards our understanding of their formation and evolution of associated star formation. We investigate the density, temperature, and column density of both dust and gas in the photodissociation regions (PDRs) located at the interface between the atomic and cold molecular gas of the NGC 7023 reflection nebula. We study how young stars affect the gas and dust in their environment. Our approach combining both dust and gas delivers strong constraints on the physical conditions of the PDRs. We find dense and warm molecular gas of high column density in the PDRs.

preprint2014arXiv

Planck 2013 results. XI. All-sky model of thermal dust emission

This paper presents an all-sky model of dust emission from the Planck 857, 545 and 353 GHz, and IRAS 100 micron data. Using a modified black-body fit to the data we present all-sky maps of the dust optical depth, temperature, and spectral index over the 353-3000 GHz range. This model is a tight representation of the data at 5 arc min. It shows variations of the order of 30 % compared with the widely-used model of Finkbeiner, Davis, and Schlegel. The Planck data allow us to estimate the dust temperature uniformly over the whole sky, providing an improved estimate of the dust optical depth compared to previous all-sky dust model, especially in high-contrast molecular regions. An increase of the dust opacity at 353 GHz, tau_353/N_H, from the diffuse to the denser interstellar medium (ISM) is reported. It is associated with a decrease in the observed dust temperature, T_obs, that could be due at least in part to the increased dust opacity. We also report an excess of dust emission at HI column densities lower than 10^20 cm^-2 that could be the signature of dust in the warm ionized medium. In the diffuse ISM at high Galactic latitude, we report an anti-correlation between tau_353/N_H and T_obs while the dust specific luminosity, i.e., the total dust emission integrated over frequency (the radiance) per hydrogen atom, stays about constant. The implication is that in the diffuse high-latitude ISM tau_353 is not as reliable a tracer of dust column density as we conclude it is in molecular clouds where the correlation of tau_353 with dust extinction estimated using colour excess measurements on stars is strong. To estimate Galactic E(B-V) in extragalactic fields at high latitude we develop a new method based on the thermal dust radiance, instead of the dust optical depth, calibrated to E(B-V) using reddening measurements of quasars deduced from Sloan Digital Sky Survey data.

preprint2014arXiv

Planck intermediate results. XVII. Emission of dust in the diffuse interstellar medium from the far-infrared to microwave frequencies

The dust-HI correlation is used to characterize the emission properties of dust in the diffuse interstellar medium. We cross-correlate sky maps from Planck, WMAP, and DIRBE, at 17 frequencies from 23 to 3000 GHz, with the Parkes survey of the 21-cm line emission of neutral atomic hydrogen, over a contiguous area of 7500 deg$^2$ centred on the southern Galactic pole. Our analysis yields four specific results. (1) The dust temperature is observed to be anti-correlated with the dust emissivity and opacity. We interpret this result as evidence for dust evolution within the diffuse ISM. The mean dust opacity is measured to be $(7.1 \pm 0.6) 10^{-27} cm^2/H \times (ν/353\, GHz)^{1.53\pm0.03}$ for $100 < ν<353$GHz. (2) We map the spectral index of dust emission at millimetre wavelengths, which is remarkably constant at $β_{mm} = 1.51\pm 0.13$. We compare it with the far infrared spectral index beta_FIR derived from greybody fits at higher frequencies, and find a systematic difference, $β_{mm}-β_{FIR} = -0.15$, which suggests that the dust SED flattens at $ν< 353\,$GHz. (3) We present spectral fits of the microwave emission correlated with HI from 23 to 353 GHz, which separate dust and anomalous microwave emission. The flattening of the dust SED can be accounted for with an additional component with a blackbody spectrum, which accounts for $(26 \pm 6)$% of the dust emission at 100 GHz and could represent magnetic dipole emission. Alternatively, it could account for an increasing contribution of carbon dust, or a flattening of the emissivity of amorphous silicates, at millimetre wavelengths. These interpretations make different predictions for the dust polarization SED. (4) We identify a Galactic contribution to the residuals of the dust-HI correlation, which we model with variations of the dust emissivity on angular scales smaller than that of our correlation analysis.

preprint2013arXiv

First results from the Herschel Gould Belt Survey in Taurus

The whole of the Taurus region (a total area of 52 sq. deg.) has been observed by the Herschel SPIRE and PACS instruments at wavelengths of 70, 160, 250, 350 and 500 μm as part of the Herschel Gould Belt Survey. In this paper we present the first results from the part of the Taurus region that includes the Barnard 18 and L1536 clouds. A new source-finding routine, the Cardiff Source-finding AlgoRithm (CSAR), is introduced, which is loosely based on CLUMPFIND, but that also generates a structure tree, or dendrogram, which can be used to interpret hierarchical clump structure in a complex region. Sources were extracted from the data using the hierarchical version of CSAR and plotted on a mass-size diagram. We found a hierarchy of objects with sizes in the range 0.024-2.7 pc. Previous studies showed that gravitationally bound prestellar cores and unbound starless clumps appeared in different places on the mass-size diagram. However, it was unclear whether this was due to a lack of instrumental dynamic range or whether they were actually two distinct populations. The excellent sensitivity of Herschel shows that our sources fill the gap in the mass-size plane between starless and pre-stellar cores, and gives the first clear supporting observational evidence for the theory that unbound clumps and (gravitationally bound) prestellar cores are all part of the same population, and hence presumably part of the same evolutionary sequence (c.f. Simpson et al. 2011).

preprint2013arXiv

Reconstructing the density and temperature structure of prestellar cores from $Herschel$ data: A case study for B68 and L1689B

Utilizing multi-wavelength dust emission maps acquired with $Herschel$, we reconstruct local volume density and dust temperature profiles for the prestellar cores B68 and L1689B using inverse-Abel transform based technique. We present intrinsic radial dust temperature profiles of starless cores directly from dust continuum emission maps disentangling the effect of temperature variations along the line of sight which was previously limited to the radiative transfer calculations. The reconstructed dust temperature profiles show a significant drop in core center, a flat inner part, and a rising outward trend until the background cloud temperature is reached. The central beam-averaged dust temperatures obtained for B68 and L1689B are 9.3 $\pm$ 0.5 K and 9.8 $\pm$0.5 K, respectively, which are lower than the temperatures of 11.3 K and 11.6 K obtained from direct SED fitting. The best mass estimates derived by integrating the volume density profiles of B68 and L1689B are 1.6 M_sol and 11 M_sol, respectively. Comparing our results for B68 with the near-infrared extinction studies, we find that the dust opacity law adopted by the HGBS project, $κ_λ =0.1(λ/300 μm)^{-2}$, agrees to within 50% with the dust extinction constraints

preprint2012arXiv

Estimating hyperparameters and instrument parameters in regularized inversion. Illustration for SPIRE/Herschel map making

We describe regularized methods for image reconstruction and focus on the question of hyperparameter and instrument parameter estimation, i.e. unsupervised and myopic problems. We developed a Bayesian framework that is based on the \post density for all unknown quantities, given the observations. This density is explored by a Markov Chain Monte-Carlo sampling technique based on a Gibbs loop and including a Metropolis-Hastings step. The numerical evaluation relies on the SPIRE instrument of the Herschel observatory. Using simulated and real observations, we show that the hyperparameters and instrument parameters are correctly estimated, which opens up many perspectives for imaging in astrophysics.

preprint2012arXiv

Interstellar matter and star formation in W5-E - A Herschel view

W5-E has been observed with the Herschel-PACS and -SPIRE photometers, at 100, 160, 250, 350, and 500 microns. The dust temperature map shows a rather uniform temperature, in the range 17.5-20 K in the dense condensations or filaments, 21-22 K in the photodissociation regions, and 24-31 K in the direction of the ionized regions. The column densities are rather low, everywhere lower than 10^23 cm-2, and of the order of a few 10^21 cm-2 in the PDRs. About 8000 solar masses of neutral material surrounds the ionized region, which is low with respect to the volume of this HII region; we suggest that the exciting stars of the W5-E, W5-W, Sh~201, A and B HII regions formed along a dense filament or sheet rather than inside a more spherical cloud. Fifty point sources have been detected at 100 microns. Most of them are Class 0/I YSOs. The SEDs of their envelopes have been fitted using a modified blackbody model. These envelopes are cold, with a mean temperature of 15.7+-1.8K. Their masses are in the range 1.3-47 solar masses. Eleven of these point sources are candidate Class 0 YSOs. Twelve of these point sources are possibly at the origin of bipolar outflows detected in this region. None of the YSOs contain a massive central object, but a few may form a massive star as they have both a massive envelope and also a high envelope accretion rate. Most of the Class 0/I YSOs are observed in the direction of high column density material, for example in the direction of the massive condensations present at the waist of the bipolar Sh 201 HII region or enclosed by the bright-rimmed cloud BRC14. The overdensity of Class 0/I YSOs on the borders of the HII regions strongly suggests that triggered star formation is at work in this region but, due to insufficient resolution, the exact processes at the origin of the triggering are difficult to determine.

preprint2012arXiv

Modelling the dust emission from dense interstellar clouds: disentangling the effects of radiative transfer and dust properties

With Planck and Herschel, we now have the spectral coverage and angular resolution required to observe dense and cold molecular clouds. As these clouds are optically thick at short wavelength but optically thin at long wavelength, it is tricky to conclude anything about dust properties without a proper treatment of the radiative transfer (RT). Our aim is to disentangle the effects of RT and of dust properties on the variations in the dust emission to provide observers with keys to analyse the emission arising from dense clouds. We model cylindrical clouds, illuminated by the ISRF, and carry out full RT calculations. Dust temperatures are solved using DustEM for amorphous carbons and silicates, grains coated with carbon mantles, and mixed aggregates of carbon and silicate. We allow variations of the grain optical properties with wavelength and temperature. We determine observed colour temperatures, T, and emissivity spectral indices, beta, by fitting the dust emission with modified blackbodies, to compare our models with observations. RT effects can neither explain the low T nor the increased submm emissivity measured at the centre of dense clouds, nor the observed beta-T anti-correlation. Adding noise to the modelled data, we show that it is not likely to be the unique explanation for the beta-T anti-correlation observed in starless clouds. It may be explained by intrinsic variations in the grain optical properties with temperature. As for the increased submm emissivity and the low T, they have to originate in variations in the grain optical properties, probably caused by their growth to form porous aggregates. We find it difficult to track back the nature of grains from the spectral variations in their emission. Finally, the column density is underestimated when determined with blackbody fitting because of the discrepancy between T and the true dust temperature at the cloud centre.

preprint2012arXiv

The Dust Properties of Bubble HII Regions as seen by Herschel

Because of their relatively simple morphology, "bubble" HII regions have been instrumental to our understanding of star formation triggered by HII regions. With the far-infrared (FIR) spectral coverage of the Herschel satellite, we can access the wavelengths where these regions emit the majority of their energy through their dust emission. At Herschel wavelengths 70 micron to 500 micron, the emission associated with HII regions is dominated by the cool dust in their photodissociation regions (PDRs). We find average dust temperatures of 26K along the PDRs, with little variation between the HII regions in the sample, while local filaments and infrared dark clouds average 19K and 15K respectively. Higher temperatures lead to higher values of the Jeans mass, which may affect future star formation. The mass of the material in the PDR, collected through the expansion of the HII region, is between ~300 and ~10,000 Solar masses for the HII regions studied here. These masses are in rough agreement with the expected masses swept up during the expansion of the \hii\ regions. Approximately 20% of the total FIR emission is from the direction of the bubble central regions. This suggests that we are detecting emission from the "near-side" and "far-side" PDRs along the line of sight and that bubbles are three-dimensional structures. We find only weak support for a relationship between dust temperature and beta, of a form similar to that caused by noise and calibration uncertainties alone.

preprint2011arXiv

Galactic cold cores II. Herschel study of the extended dust emission around the first Planck detections

Within the project Galactic Cold Cores we are carrying out Herschel photometric observations of cold interstellar clouds detected with the Planck satellite. The three fields observed as part of the Herschel science demonstration phase (SDP) provided the first glimpse into the nature of these sources. We examine the properties of the dust emission within the fields. We determine the dust sub-millimetre opacity, look for signs of spatial variations in the dust spectral index, and estimate how the apparent variations of the parameters could be affected by different sources of uncertainty. We use the Herschel observations where the zero point of the surface brightness scale is set with the help of the Planck satellite data. We derive the colour temperature and column density maps of the regions and determine the dust opacity by a comparison with extinction measurements. By simultaneously fitting the colour temperature and the dust spectral index values we look for spatial variations in the apparent dust properties. With a simple radiative transfer model we estimate to what extent these can be explained by line-of-sight temperature variations, without changes in the dust grain properties. The analysis of the dust emission reveals cold and dense clouds that coincide with the Planck sources and confirm those detections. The derived dust opacity varies in the range kappa(250um) ~ 0.05-0.2 cm^2/g, higher values being observed preferentially in regions of high column density. The average dust spectral index beta is ~ 1.9-2.2. There are indications that beta increases towards the coldest regions. The spectral index decreases strongly near internal heating sources but, according to radiative transfer models, this can be explained by the line-of-sight temperature variations without a change in the dust properties.

preprint2011arXiv

Planck Early Results XXIV: Dust in the diffuse interstellar medium and the Galactic halo

This paper presents the first results of comparison of Planck along with IRAS data with Green Bank Telescope 21-cm observations in 14 fields covering more than 800 square degrees at high Galactic latitude. Galactic dust emission for fields with average HI column density lower than 2 x 10^20 cm^-2 is well correlated with 21-cm emission. The residual emission in these fields, once the HI-correlated emission is removed, is consistent with the expected statistical properties of the cosmic infrared background fluctuations. Fields with larger column densities show significant excess dust emission compared to the HI column density. Regions of excess lie in organized structures that suggest the presence of hydrogen in molecular form, though they are not always correlated with CO emission. Dust emission from intermediate-velocity clouds is detected with high significance. Its spectral properties are consistent with, compared to the local ISM values, significantly hotter dust (T~20 K), lower sub-millimeter dust opacity, and a relative abundance of very small grains to large grains about four times higher. These results are compatible with expectations for clouds that are part of the Galactic fountain in which there is dust shattering and fragmentation. Correlated dust emission in HVCs is not detected; the average of the 99.9% confidence upper limits to the emissivity is 0.15 times the local ISM value at 857 and 3000 GHz, in accordance with gas phase evidence for lower metallicity and depletion in these clouds. Unexpected anti-correlated variations of the dust temperature and emission cross-section per H atom are identified in the local ISM and IVCs, a trend that continues into molecular environments. This suggests that dust growth through aggregation, seen in molecular clouds, is active much earlier in the cloud condensation and star formation processes.

preprint2011arXiv

Planck Early Results: Thermal dust in Nearby Molecular Clouds

Planck allows unbiased mapping of Galactic sub-millimetre and millimetre emission from the most diffuse regions to the densest parts of molecular clouds. We present an early analysis of the Taurus molecular complex, on line-of-sight-averaged data and without component separation. The emission spectrum measured by Planck and IRAS can be fitted pixel by pixel using a single modified blackbody. Some systematic residuals are detected at 353 GHz and 143 GHz, with amplitudes around -7 % and +13 %, respectively, indicating that the measured spectra are likely more complex than a simple modified blackbody. Significant positive residuals are also detected in the molecular regions and in the 217 GHz and 100 GHz bands, mainly caused by to the contribution of the J=2-1 and J=1-0 12CO and 13CO emission lines. We derive maps of the dust temperature T, the dust spectral emissivity index beta, and the dust optical depth at 250 microns tau. The temperature map illustrates the cooling of the dust particles in thermal equilibrium with the incident radiation field, from 16-17 K in the diffuse regions to 13-14 K in the dense parts. The distribution of spectral indices is centred at 1.78, with a standard deviation of 0.08 and a systematic error of 0.07. We detect a significant T-beta anti-correlation. The dust optical depth map reveals the spatial distribution of the column density of the molecular complex from the densest molecular regions to the faint diffuse regions. We use near-infrared extinction and HI data at 21-cm to perform a quantitative analysis of the spatial variations of the measured dust optical depth at 250 microns per hydrogen atom tau/NH. We report an increase of tau/NH by a factor of about 2 between the atomic phase and the molecular phase, which has a strong impact on the equilibrium temperature of the dust particles.

preprint2011arXiv

Planck early results. XXI. Properties of the interstellar medium in the Galactic plane

(abridged) Planck has observed the entire sky from 30 GHz to 857GHz. The observed foreground emission contains contributions from different phases of the interstellar medium (ISM). We have separated the observed Galactic emission into the different gaseous components (atomic, molecular and ionised) in each of a number of Galactocentric rings. Templates are created for various Galactocentric radii using velocity information from atomic (neutral hydrogen) and molecular (12CO) observations. The ionised template is assumed to be traced by free-free emission as observed by WMAP, while 408 MHz emission is used to trace the synchrotron component. Gas emission not traced by the above templates, namely "ark gas", as evidenced using Planck data, is included as an additional template, the first time such a component has been used in this way. These templates are then correlated with each of the Planck frequency bands, as well as other ancillary data. The emission per column density of the gas templates allows us to create distinct spectral energy distributions (SEDs) per Galactocentric ring and in each of the gaseous tracers from 1.4 GHz to 25 THz (12μm). Apart from the thermal dust and free-free emission, we have probed the Galaxy for anomalous (e.g., spinning) dust as well as synchrotron emission. We show that anomalous dust emission is present in the atomic, molecular and dark gas phases throughout the Galactic disk. The derived dust propeties associated with the dark gas phase are derived but do not allow us to reveal the nature of this phase. For all environments, the anomalous emission is consistent with rotation from polycyclic aromatic hydrocarbons (PAHs) and, according to our simple model, accounts for $(25\pm5)%$ (statistical) of the total emission at 30 GHz.

preprint2011arXiv

Spatial variation of the cooling lines in the Orion Bar from Herschel/PACS

We present spatially resolved Herschel/PACS observations of the Orion Bar. We have characterise the emission of the far-infrared fine-structure lines of [CII] (158um), [OI] (63 and 145um), and [NII] (122um) that trace the gas local conditions. The observed distribution and variation of the lines are discussed in relation to the underlying geometry and linked to the energetics associated with the Trapezium stars. These observations enable us to map the spatial distribution of these fine-structure lines with a spatial resolution between 4" and 11" and covering a total square area of about 120"x105". The spatial profile of the emission lines are modelled using the radiative transfer code Cloudy. We find that the spatial distribution of the [CII] line coincides with that of the [OI] lines. The [NII] line peaks closer to the ionising star than the other three lines, but with a small region of overlap. We can distinguish several knots of enhanced emission within the Bar indicating the presence of an inhomogenous and structured medium. The emission profiles cannot be reproduced by a single photo-dissociation region, clearly indicating that, besides the Bar, there is a significant contribution from additional photo-dissociation region(s) over the area studied. The combination of both the [NII] and [OI] 145um lines can be used to estimate the [CII] emission and distinguish between its ionised or neutral origin. We have calculated how much [CII] emission comes from the neutral and ionised region, and find that at least 82% originates from the photo-dissocciation region. Together, the [CII] 158um and [OI] 63 and 145um lines account for 90% of the power emitted by the main cooling lines in the Bar (including CO, H2, etc...), with [OI] 63um alone accounting for 72% of the total.

preprint2011arXiv

Super-resolution in map-making based on a physical instrument model and regularized inversion. Application to SPIRE/Herschel

We investigate super-resolution methods for image reconstruction from data provided by a family of scanning instruments like the Herschel observatory. To do this, we constructed a model of the instrument that faithfully reflects the physical reality, accurately taking the acquisition process into account to explain the data in a reliable manner. The inversion, ie the image reconstruction process, is based on a linear approach resulting from a quadratic regularized criterion and numerical optimization tools. The application concerns the reconstruction of maps for the SPIRE instrument of the Herschel observatory. The numerical evaluation uses simulated and real data to compare the standard tool (coaddition) and the proposed method. The inversion approach is capable to restore spatial frequencies over a bandwidth four times that possible with coaddition and thus to correctly show details invisible on standard maps. The approach is also applied to real data with significant improvement in spatial resolution.

preprint2010arXiv

A Herschel study of the properties of starless cores in the Polaris Flare dark cloud region using PACS and SPIRE

The Polaris Flare cloud region contains a great deal of extended emission. It is at high declination and high Galactic latitude. It was previously seen strongly in IRAS Cirrus emission at 100 microns. We have detected it with both PACS and SPIRE on Herschel. We see filamentary and low-level structure. We identify the five densest cores within this structure. We present the results of a temperature, mass and density analysis of these cores. We compare their observed masses to their virial masses, and see that in all cases the observed masses lie close to the lower end of the range of estimated virial masses. Therefore, we cannot say whether they are gravitationally bound prestellar cores. Nevertheless, these are the best candidates to be potentialprestellar cores in the Polaris cloud region.

preprint2010arXiv

Clouds, filaments and protostars: the Herschel Hi-GAL Milky Way

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2° x 2° tiles approximately centered at l=30° and l=59°. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59° field; a A_V value between 5 and 10 is found for the l=30° field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

preprint2010arXiv

Filamentary structures and compact objects in the Aquila and Polaris clouds observed by Herschel

Our PACS and SPIRE images of the Aquila Rift and part of the Polaris Flare regions, taken during the science demonstration phase of Herschel discovered fascinating, omnipresent filamentary structures that appear to be physically related to compact cores. We briefly describe a new multi-scale, multi-wavelength source extraction method used to detect objects and measure their parameters in our Herschel images. All of the extracted starless cores (541 in Aquila and 302 in Polaris) appear to form in the long and very narrow filaments. With its combination of the far-IR resolution and sensitivity, Herschel directly reveals the filaments in which the dense cores are embedded; the filaments are resolved and have deconvolved widths of 35 arcsec in Aquila and 59 arcsec in Polaris (9000 AU in both regions). Our first results of observations with Herschel enable us to suggest that in general dense cores may originate in a process of fragmentation of complex networks of long, thin filaments, likely formed as a result of an interplay between gravity, interstellar turbulence, and magnetic fields. To unravel the roles of the processes, one has to obtain additional kinematic and polarization information; these follow-up observations are planned.

preprint2010arXiv

First detection of the Methylidyne cation (CH+) fundamental rotational line with the Herschel/SPIRE FTS

Aims. To follow the species chemistry arising in diverse sources of the Galaxy with Herschel. Methods. SPIRE FTS sparse sampled maps of the Orion bar & compact HII regions G29.96-0.02 and G32.80+0.19 have been analyzed. Results. Beyond the wealth of atomic and molecular lines detected in the high-resolution spectra obtained with the FTS of SPIRE in the Orion Bar, one emission line is found to lie at the position of the fundamental rotational transition of CH+ as measured precisely in the laboratory (Pearson & Drouion 2006). This coincidence suggests that it is the first detection of the fundamental rotational transition of CH+. This claim is strengthened by the observation of the lambda doublet transitions arising from its relative, CH, which are also observed in the same spectrum. The broad spectral coverage of the SPIRE FTS allows for the simultaneous measurement of these closely related chemically species, under the same observing conditions. The importance of these lines are discussed and a comparison with results obtained from models of the Photon Dominated Region (PDR) of Orion are presented. The CH+ line also appears in absorption in the spectra of the two galactic compact HII regions G29.96-0.02 and G32.80+0.19, which is likely due to the presence of CH+ in the the Cold Neutral Medium of the galactic plane. These detections will shed light on the formation processes and on the existence of CH+, which are still outstanding questions in astrophysics.

preprint2010arXiv

From filamentary clouds to prestellar cores to the stellar IMF: Initial highlights from the Herschel Gould Belt survey

We summarize the first results from the Gould Belt survey, obtained toward the Aquila Rift and Polaris Flare regions during the 'science demonstration phase' of Herschel. Our 70-500 micron images taken in parallel mode with the SPIRE and PACS cameras reveal a wealth of filamentary structure, as well as numerous dense cores embedded in the filaments. Between ~ 350 and 500 prestellar cores and ~ 45-60 Class 0 protostars can be identified in the Aquila field, while ~ unbound starless cores and no protostars are observed in the Polaris field. The prestellar core mass function (CMF) derived for the Aquila region bears a strong resemblance to the stellar initial mass function (IMF), already confirming the close connection between the CMF and the IMF with much better statistics than earlier studies. Comparing and contrasting our Herschel results in Aquila and Polaris, we propose an observationally-driven scenario for core formation according to which complex networks of long, thin filaments form first within molecular clouds, and then the densest filaments fragment into a number of prestellar cores via gravitational instability.

preprint2010arXiv

Herschel observations of embedded protostellar clusters in the Rosette Molecular Cloud

The Herschel OB young stellar objects survey (HOBYS) has observed the Rosette molecular cloud, providing an unprecedented view of its star formation activity. These new far-infrared data reveal a population of compact young stellar objects whose physical properties we aim to characterise. We compiled a sample of protostars and their spectral energy distributions that covers the near-infrared to submillimetre wavelength range. These were used to constrain key properties in the protostellar evolution, bolometric luminosity, and envelope mass and to build an evolutionary diagram. Several clusters are distinguished including the cloud centre, the embedded clusters in the vicinity of luminous infrared sources, and the interaction region. The analysed protostellar population in Rosette ranges from 0.1 to about 15 Msun with luminosities between 1 and 150 Lsun, which extends the evolutionary diagram from low-mass protostars into the high-mass regime. Some sources lack counterparts at near- to mid-infrared wavelengths, indicating extreme youth. The central cluster and the Phelps & Lada 7 cluster appear less evolved than the remainder of the analysed protostellar population. For the central cluster, we find indications that about 25% of the protostars classified as Class I from near- to mid-infrared data are actually candidate Class 0 objects. As a showcase for protostellar evolution, we analysed four protostars of low- to intermediate-mass in a single dense core, and they represent different evolutionary stages from Class 0 to Class I. Their mid- to far-infrared spectral slopes flatten towards the Class I stage, and the 160 to 70um flux ratio is greatest for the presumed Class 0 source. This shows that the Herschel observations characterise the earliest stages of protostellar evolution in detail.

preprint2010arXiv

Herschel-SPIRE observations of the Polaris flare : structure of the diffuse interstellar medium at the sub-parsec scale

We present a power spectrum analysis of the Herschel-SPIRE observations of the Polaris flare, a high Galactic latitude cirrus cloud midway between the diffuse and molecular phases. The SPIRE images of the Polaris flare reveal for the first time the structure of the diffuse interstellar medium down to 0.01 parsec over a 10 square degrees region. These exceptional observations highlight the highly filamentary and clumpy structure of the interstellar medium even in diffuse regions of the map. The power spectrum analysis shows that the structure of the interstellar medium is well described by a single power law with an exponent of -2.7 +- 0.1 at all scales from 30" to 8 degrees. That the power spectrum slope of the dust emission is constant down to the SPIRE angular resolution is an indication that the inertial range of turbulence extends down to the 0.01 pc scale. The power spectrum analysis also allows the identification of a Poissonian component at sub-arcminute scales in agreement with predictions of the cosmic infrared background level at SPIRE wavelengths. Finally, the comparison of the SPIRE and IRAS 100 micron data of the Polaris flare clearly assesses the capability of SPIRE in maping diffuse emission over large areas.

preprint2010arXiv

Herschel-SPIRE spectroscopy of G29.96-0.02: fitting the full SED

We use the SPIRE Fourier-Transform Spectrometer (FTS) on-board the ESA Herschel Space Telescope to analyse the submillimetre spectrum of the Ultra-compact HII region G29.96-0.02. Spectral lines from species including 13CO, CO, [CI], and [NII] are detected. A sparse map of the [NII] emission shows at least one other HII region neighbouring the clump containing the UCHII. The FTS spectra are combined with ISO SWS and LWS spectra and fluxes from the literature to present a detailed spectrum of the source spanning three orders of magnitude in wavelength. The quality of the spectrum longwards of 100 μm allows us to fit a single temperature greybody with temperature 80.3\pm0.6K and dust emissivity index 1.73\pm0.02, an accuracy rarely obtained with previous instruments. We estimate a mass of 1500 Msol for the clump containing the HII region. The clump's bolometeric luminosity of 4 x 10^6 Lsol is comparable to, or slightly greater than, the known O-star powering the UCHII region.

preprint2010arXiv

Herschel-SPIRE spectroscopy of the DR21 molecular cloud core

We present far-infrared spectra and maps of the DR21 molecular cloud core between 196 and 671 microns, using the Herschel-SPIRE spectrometer. Nineteen molecular lines originating from CO, 13CO, HCO+ and H2O, plus lines of [N II] and [CI] were recorded, including several transitions not previously detected. The CO lines are excited in warm gas with Tkin ~ 125 K and nH2 ~ 7 x 10^4 cm-3, CO column density N(CO) ~ 3.5 x 10^18 cm^-2 and a filling factor of ~ 12%, and appear to trace gas associated with an outflow. The rotational temperature analysis incorporating observations from ground-based telescopes reveals an additional lower excitation CO compoment which has a temperature ~ 78 K and N(CO) ~ 4.5 x 10^21 cm^-2. Astronomy & Astrophysics HERSCHEL special Issue, in press.

preprint2010arXiv

Optical performance of the JWST MIRI flight model: characterization of the point spread function at high-resolution

The Mid Infra Red Instrument (MIRI) is one of the four instruments onboard the James Webb Space Telescope (JWST), providing imaging, coronagraphy and spectroscopy over the 5-28 microns band. To verify the optical performance of the instrument, extensive tests were performed at CEA on the flight model (FM) of the Mid-InfraRed IMager (MIRIM) at cryogenic temperatures and in the infrared. This paper reports on the point spread function (PSF) measurements at 5.6 microns, the shortest operating wavelength for imaging. At 5.6 microns the PSF is not Nyquist-sampled, so we use am original technique that combines a microscanning measurement strategy with a deconvolution algorithm to obtain an over-resolved MIRIM PSF. The microscanning consists in a sub-pixel scan of a point source on the focal plane. A data inversion method is used to reconstruct PSF images that are over-resolved by a factor of 7 compared to the native resolution of MIRI. We show that the FWHM of the high-resolution PSFs were 5-10% wider than that obtained with Zemax simulations. The main cause was identified as an out-of-specification tilt of the M4 mirror. After correction, two additional test campaigns were carried out, and we show that the shape of the PSF is conform to expectations. The FWHM of the PSFs are 0.18-0.20 arcsec, in agreement with simulations. 56.1-59.2% of the total encircled energy (normalized to a 5 arcsec radius) is contained within the first dark Airy ring, over the whole field of view. At longer wavelengths (7.7-25.5 microns), this percentage is 57-68%. MIRIM is thus compliant with the optical quality requirements. This characterization of the MIRIM PSF, as well as the deconvolution method presented here, are of particular importance, not only for the verification of the optical quality and the MIRI calibration, but also for scientific applications.

preprint2010arXiv

Physical properties of the Sh2-104 HII region as seen by Herschel

Context: Sh2-104 is a Galactic H ii region with a bubble morphology, detected at optical and radio wavelengths. It is considered the first observational confirmation of the collect-and-collapse model of triggered star-formation. Aims: We aim to analyze the dust and gas properties of the Sh2-104 region to better constrain its effect on local future generations of stars. In addition, we investigate the relationship between the dust emissivity index β and the dust temperature, T_dust. Methods: Using Herschel PACS and SPIRE images at 100, 160, 250, 350 and 500 μm we determine T_dust and β throughout Sh2-104, fitting the spectral energy distributions (SEDs) obtained from aperture photometry. With the SPIRE Fourier transform spectrometer (FTS) we obtained spectra at different positions in the Sh2-104 region. We detect J-ladders of CO and 13CO, with which we derive the gas temperature and column density. We also detect proxies of ionizing flux as the [NII] 3P1-3P0 and [CI] 3P2-3P1 transitions. Results: We find an average value of β ~ 1.5 throughout Sh2-104, as well as a T dust difference between the photodissociation region (PDR, ~ 25 K) and the interior (~ 40 K) of the bubble. We recover the anti-correlation between β and dust temperature reported numerous times in the literature. The relative isotopologue abundances of CO appear to be enhanced above the standard ISM values, but the obtained value is very preliminary and is still affected by large uncertainties.

preprint2010arXiv

Small-scale structure in the Rosette molecular cloud revealed by Herschel

We present a preliminary analysis of the small-scale structure found in new 70-520 micron continuum maps of the Rosette molecular cloud (RMC), obtained with the SPIRE and PACS instruments of the Herschel Space Observatory. We find 473 clumps within the RMC using a new structure identification algorithm, with sizes up to ~1.0 pc in diameter. A comparison with recent Spitzer maps reveals that 371 clumps are "starless" (without an associated young stellar object), while 102 are "protostellar." Using the respective values of dust temperature, we determine the clumps have masses (M_C) over the range -0.75 <= log (M_C/M_sun) <= 2.50. Linear fits to the high-mass tails of the resulting clump mass spectra (CMS) have slopes that are consistent with those found for high-mass clumps identified in CO emission by other groups.

preprint2010arXiv

Star formation triggered by the Galactic HII region RCW 120: First results from the Herschel Space Observatory

By means of different physical mechanisms, the expansion of HII regions can promote the formation of new stars of all masses. RCW 120 is a nearby Galactic HII region where triggered star formation occurs. This region is well-studied - there being a wealth of existing data - and is nearby. However, it is surrounded by dense regions for which far infrared data is essential to obtain an unbiased view of the star formation process and in particular to establish whether very young protostars are present. We attempt to identify all Young Stellar Objects (YSOs), especially those previously undetected at shorter wavelengths, to derive their physical properties and obtain insight into the star formation history in this region. We use Herschel-PACS and -SPIRE images to determine the distribution of YSOs observed in the field. We use a spectral energy distribution fitting tool to derive the YSOs physical properties. Herschel-PACS and -SPIRE images confirm the existence of a young source and allow us to determine its nature as a high-mass (8-10 MSun) Class 0 object (whose emission is dominated by a massive envelope) towards the massive condensation 1 observed at (sub)-millimeter wavelengths. This source was not detected at 24 micron and only barely seen in the MISPGAL 70 micron data. Several other red sources are detected at Herschel wavelengths and coincide with the peaks of the millimeter condensations. SED fitting results for the brightest Herschel sources indicate that, apart from the massive Class 0 that forms in condensation 1, young low mass stars are forming around RCW 120. The YSOs observed on the borders of RCW 120 are younger than its ionizing star, which has an age of about 2.5 Myr.

preprint2010arXiv

The Aquila prestellar core population revealed by Herschel

The origin and possible universality of the stellar initial mass function (IMF) is a major issue in astrophysics. One of the main objectives of the Herschel Gould Belt Survey is to clarify the link between the prestellar core mass function (CMF) and the IMF. We present and discuss the core mass function derived from Herschel data for the large population of prestellar cores discovered with SPIRE and PACS in the Aquila Rift cloud complex at d ~ 260 pc. We detect a total of 541 starless cores in the entire ~11 deg^2 area of the field imaged at 70-500 micron with SPIRE/PACS. Most of these cores appear to be gravitationally bound, and thus prestellar in nature. Our Herschel results confirm that the shape of the prestellar CMF resembles the stellar IMF, with much higher quality statistics than earlier submillimeter continuum ground-based surveys.

preprint2010arXiv

The Herschel first look at protostars in the Aquila Rift

As part of the science demonstration phase of the Herschel mission of the Gould Belt Key Program, the Aquila Rift molecular complex has been observed. The complete ~ 3.3deg x 3.3deg imaging with SPIRE 250/350/500 micron and PACS 70/160 micron allows a deep investigation of embedded protostellar phases, probing of the dust emission from warm inner regions at 70 and 160 micron to the bulk of the cold envelopes between 250 and 500 micron. We used a systematic detection technique operating simultaneously on all Herschel bands to build a sample of protostars. Spectral energy distributions are derived to measure luminosities and envelope masses, and to place the protostars in an M_env - L_bol evolutionary diagram. The spatial distribution of protostars indicates three star-forming sites in Aquila, with W40/Sh2-64 HII region by far the richest. Most of the detected protostars are newly discovered. For a reduced area around the Serpens South cluster, we could compare the Herschel census of protostars with Spitzer results. The Herschel protostars are younger than in Spitzer with 7 Class 0 YSOs newly revealed by Herschel. For the entire Aquila field, we find a total of ~ 45-60 Class 0 YSOs discovered by Herschel. This confirms the global statistics of several hundred Class~0 YSOs that should be found in the whole Gould Belt survey.

preprint2010arXiv

The Herschel view of star formation in the Rosette molecular cloud under the influence of NGC 2244

The Rosette molecular cloud is promoted as the archetype of a triggered star-formation site. This is mainly due to its morphology, because the central OB cluster NGC 2244 has blown a circular-shaped cavity into the cloud and the expanding HII-region now interacts with the cloud. Studying the spatial distribution of the different evolutionary states of all star-forming sites in Rosette and investigating possible gradients of the dust temperature will help to test the 'triggered star-formation' scenario in Rosette. We use continuum data obtained with the PACS (70 and 160 micron) and SPIRE instruments (250, 350, 500 micron) of the Herschel telescope during the Science Demonstration Phase of HOBYS. Three-color images of Rosette impressively show how the molecular gas is heated by the radiative impact of the NGC 2244 cluster. A clear negative temperature gradient and a positive density gradient (running from the HII-region/molecular cloud interface into the cloud) are detected. Studying the spatial distribution of the most massive dense cores (size scale 0.05 to 0.3 pc), we find an age-sequence (from more evolved to younger) with increasing distance to the cluster NGC 2244. No clear gradient is found for the clump (size-scale up to 1 pc) distribution. The existence of temperature and density gradients and the observed age-sequence imply that star formation in Rosette may indeed be influenced by the radiative impact of the central NGC 2244 cluster. A more complete overview of the prestellar and protostellar population in Rosette is required to obtain a firmer result.

preprint2010arXiv

The Herschel-SPIRE instrument and its in-flight performance

The Spectral and Photometric Imaging Receiver (SPIRE), is the Herschel Space Observatory`s submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 microns, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 microns (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4' x 8', observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6'. The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the Herschel cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5-2.

preprint2010arXiv

The physical properties of the dust in the RCW 120 HII region as seen by Herschel

Context. RCW 120 is a well-studied, nearby Galactic HII region with ongoing star formation in its surroundings. Previous work has shown that it displays a bubble morphology at mid-infrared wavelengths and has a massive layer of collected neutral material seen at sub-mm wavelengths. Given the well-defined photo-dissociation region (PDR) boundary and collected layer, it is an excellent laboratory to study the "collect and collapse" process of triggered star formation. Using Herschel Space Observatory data at 100, 160, 250, 350, and 500 micron, in combination with Spitzer and APEX-LABOCA data, we can for the first time map the entire spectral energy distribution of an HII region at high angular resolution. Aims. We seek a better understanding of RCW120 and its local environment by analysing its dust temperature distribution. Additionally, we wish to understand how the dust emissivity index, beta, is related to the dust temperature. Methods. We determine dust temperatures in selected regions of the RCW 120 field by fitting their spectral energy distribution (SED), derived using aperture photometry. Additionally, we fit the SED extracted from a grid of positions to create a temperature map. Results. We find a gradient in dust temperature, ranging from >30 K in the interior of RCW 120, to ~20K for the material collected in the PDR, to ~10K toward local infrared dark clouds and cold filaments. Our results suggest that RCW 120 is in the process of destroying the PDR delineating its bubble morphology. The leaked radiation from its interior may influence the creation of the next generation of stars. We find support for an anti-correlation between the fitted temperature and beta, in rough agreement with what has been found previously. The extended wavelength coverage of the Herschel data greatly increases the reliability of this result.

preprint2009arXiv

Data inversion for over-resolved spectral imaging in astronomy

We present an original method for reconstructing a three-dimensional object having two spatial dimensions and one spectral dimension from data provided by the infrared slit spectrograph on board the Spitzer Space Telescope. During acquisition, the light flux is deformed by a complex process comprising four main elements (the telescope aperture, the slit, the diffraction grating and optical distortion) before it reaches the two-dimensional sensor. The originality of this work lies in the physical modelling, in integral form, of this process of data formation in continuous variables. The inversion is lso approached with continuous variables in a semi-parametric format decomposing the object into a family of Gaussian functions. The estimate is built in a deterministic regularization framework as the minimizer of a quadratic criterion. These specificities give our method the power to over-resolve. Its performance is illustrated using real and simulated data. We also present a study of the resolution showing a 1.5-fold improvement relative to conventional methods.

preprint2007arXiv

Analysis of the emission of very small dust particles from Spitzer spectro-imagery data using blind signal separation methods

This work was conducted as part of the SPECPDR program, dedicated to the study of very small particles and astrochemistry, in Photo-Dissociation Regions (PDRs). We present the analysis of the mid-IR spectro-imagery observations of Ced 201, NCG 7023 East and North-West and rho-Ophiuchi West filament. Using the data from all four modules of the InfraRed Spectrograph onboard the Spitzer Space Telescope, we produced a spectral cube ranging from 5 to 35 um for each one of the observed PDRs. The resulting cubes were analysed using Blind Signal Separation methods (NMF and FastICA). For Ced 201, rho-Ophiuchi West filament and NGC 7023 East, we find that two signals can be extracted from the original data cubes, which are 5 to 35 um spectra. The main features of the first spectrum are a strong continuum emission at long wavelengths, and a broad 7.8 um band. On the contrary, the second spectrum exhibits the classical Aromatic Infrared Bands (AIBs) and no continuum. The reconstructed spatial distribution maps show that the latter spectrum is mainly present at the cloud surface, close to the star whereas the first one is located slightly deeper inside the PDR. The study of the spectral energy distribution of Ced 201 up to 100 um suggests that, in cool PDRs, the 5-25 um continuum is carried by Very Small Grains (VSGs). The AIB spectra in the observed objects can be interpreted as the contribution of neutral and positively-charged Polycyclic Aromatic Hydrocarbons (PAHs). We extracted the 5 to 25 um emission spectrum of VSGs in cool PDRs, these grains being most likely carbonaceous. We show that the variations of the mid-IR (5-35 um) spectra of PDRs can be explained by the photo-chemical processing of these VSGs and PAHs, VSGs being the progenitors of free PAH.