Source author record

N. Ysard

N. Ysard 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

43works
4topics
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

43 published item(s)

preprint2026arXiv

Resolved ISM properties and scaling relations in the barred galaxy NGC 3627: constraints from NIKA2 observations

We investigate the interplay between star formation, interstellar medium (ISM) components, and dust properties in NGC 3627 using new NIKA2 1.15 and 2 mm observations from the IMEGIN Large Program. Our goal is to analyze dust and radio emission, decompose contributions in the millimeter-centimeter regime, and explore ISM properties within the galaxy. We perform spectral energy distribution fitting, at both global and spatial scales, using the THEMIS dust model within the HerBIE code, applied to data from 3.4 $μ$m to 6 cm. We decompose emission into dust, free-free, and synchrotron components, and examine correlations with gas surface density and star formation activity. Additionally, we analyze the small dust grain fraction and its variation across the galaxy. We find $\sim$10% radio emission at 2 mm, peaking at 18% in the southern bar-end, which hosts the highest star formation activity. However, an isolated star-forming region beyond this bar-end is the most efficient, as indicated by its elevated dust production efficiency and effective yield, predicted by our simplistic dust evolution model. The 160 $μ$m emission shows the strongest correlation with molecular gas, while 1.15 mm better traces the dust mass surface density. Small grains, which make up $\sim$13% of dust mass (2 $\times$ 10$^{7}$ M$_{\odot}$), are depleted in intense radiation fields, with a notable deficit in the southern tidal tail. ISM properties and chemical evolution indicate that dynamical processes, such as bar-driven gas flows and tidal interactions, are crucial in shaping the galactic structure, influencing star formation efficiency, and dust distribution.

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.

preprint2025arXiv

The essential elements of dust evolution: a-C(:H) nanoparticle sub-structures and photo-fragmentation

Hydrogenated amorphous carbon materials, a-C(:H), are heterogeneous structures consisting of carbon atoms in different hybridisation states and bonding configurations and are thought to constitute a significant and observationally important fraction of the interstellar dust material. This work aims to characterise semi-conducting a-C(:H) nanoparticle structures and, in particular, their property-characterising aromatic domain size distribution and so predict how they will behave in intense UV radiation fields that can fragment them through dissociative and charge effects as a result of carbon-carbon bond-breaking. Using a statistical approach we determine the typical sizes of the aromatic domains, their size distribution, how they are network-bonded, and where they are to be found within the structure. We consider the effects of thermal excitation, photo-dissociation and charging of a-C(:H) nanoparticles, and the products of their fragmentation. The derived UV photon-induced fragmentation lifetimes for nanometre-sized a-C(:H) nanoparticles, with radii ~0.4-0.5nm radius and containing ~40-60 carbon atoms, are of the order of 10^6-10^7yr in the diffuse interstellar medium and likely 10^2-10^4 times shorter in photodissociation regions, depending on the local radiation field intensity. Grains larger than this are stable against photodissociation. In H{\footnotesize II} regions only a-C(:H) nanoparticles with radii greater than 0.7nm (> 150 carbon atoms) are likely to survive. The photon-driven fragmentation of sub-nanometre a-C(:H) particles was determined to be important in the diffuse interstellar medium and also in high excitation regions, such as photodissociation and HII regions. However, in these same regions Coulomb fragmentation is unlikely to be an important dust destruction process.

preprint2022arXiv

Dust emissivity in resolved spiral galaxies

Context: The far-infrared (FIR) and sub-millimeter (submm) emissivity of the Milky Way (MW) cirrus is an important benchmark for dust grain models. Dust masses in other galaxies are generally derived from the FIR/submm using the emission properties of these MW-calibrated models. Aims: We seek to derive the FIR/submm emissivity in nine nearby spiral galaxies to check its compatibility with MW cirrus measurements. Methods: We obtained values of the emissivity at 70 to 500 um, using maps of dust emission from the Herschel satellite and of gas surface density from the THINGS and HERACLES surveys on a scale generally corresponding to 440 pc. We studied the variation of the emissivity with the surface brightness ratio I(250um)/I(500um), a proxy for the intensity of the interstellar radiation field heating the dust. Results: We find that the average value of the emissivity agrees with MW estimates for pixels sharing the same color as the cirrus, namely, for I(250um)/I(500um) = 4.5. For I(250um)/I(500um) > 5, the measured emissivity is instead up to a factor ~2 lower than predicted from MW dust models heated by stronger radiation fields. Regions with higher I(250um)/I(500um) are preferentially closer to the galactic center and have a higher overall (stellar+gas) surface density and molecular fraction. The results do not depend strongly on the adopted CO-to-molecular conversion factor and do not appear to be affected by the mixing of heating conditions. Conclusions: Our results confirm the validity of MW dust models at low density, but are at odds with predictions for grain evolution in higher density environments. If the lower-than-expected emissivity at high I(250um)/I(500um) is the result of intrinsic variations in the dust properties, it would imply an underestimation of the dust mass surface density of up to a factor ~2 when using current dust models.

preprint2022arXiv

Searching for Anomalous Microwave Emission in nearby galaxies. K-band observations with the Sardinia Radio Telescope

We observed four nearby spiral galaxies (NGC 3627, NGC 4254, NGC 4736 and NGC 5055) in the K band with the 64-m Sardinia Radio Telescope, with the aim of detecting the Anomalous Microwave Emission (AME), a radiation component presumably due to spinning dust grains, observed so far in the Milky Way and in a handful of other galaxies only (most notably, M 31). We mapped the galaxies at 18.6 and 24.6 GHz and studied their global photometry together with other radio-continuum data from the literature, in order to find AME as emission in excess of the synchrotron and thermal components. We only find upper limits for AME. These non-detections, and other upper limits in the literature, are nevertheless consistent with the average AME emissivity from the few detections: it is $ε^\mathrm{AME}_{\mathrm{30~GHz}} = 2.4\pm0.4 \times 10^{-2}$ MJy sr$^{-1}$ (M$_\odot$ pc$^{-2}$)$^{-1}$ in units of dust surface density (equivalently, $1.4\pm0.2 \times 10^{-18}$ Jy sr$^{-1}$ (H cm$^{-2}$)$^{-1}$ in units of H column density). We finally suggest to search for AME in quiescent spirals with relatively low radio luminosity, such as M~31.

preprint2021arXiv

The mm-to-cm SED of spiral galaxies. Synergies between NIKA2 and SRT instruments

The mm-to-cm range of the Spectral Energy Distribution of spiral galaxies remains largely unexplored. Its coverage is required to disentangle the contribution of dust emission, free-free and synchrotron radiation and can provide constraints on dust models, star-formation rates and ISM properties. We present the case for a synergy between NIKA2 observations of nearby spirals and those from planned and current instrumentation at the Sardinia Radio Telescope, and report on a pilot K-band program to search for Anomalous Microwave Emission, an elusive emission component which is presumably related to dust.

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.

preprint2020arXiv

High-resolution, 3D radiative transfer modelling IV. AGN-powered dust heating in NGC 1068

Dust emission, an important diagnostic of star formation and ISM mass throughout the Universe, can be powered by sources unrelated to ongoing star formation. In the framework of the DustPedia project we have set out to disentangle the radiation of the ongoing star formation from that of the older stellar populations. This is done through detailed, 3D radiative transfer simulations of face-on spiral galaxies. In this particular study, we focus on NGC 1068, which contains an active galactic nucleus (AGN). The effect of diffuse dust heating by AGN (beyond the torus) was so far only investigated for quasars. This additional dust heating source further contaminates the broadband fluxes on which classic galaxy modelling tools rely to derive physical properties. We aim to fit a realistic model to the observations of NGC 1068 and quantify the contribution of the several dust heating sources. Our model is able to reproduce the global spectral energy distribution of the galaxy. It matches the resolved optical and infrared images fairly well, but deviates in the UV and the submm. We find a strong wavelength dependency of AGN contamination to the broadband fluxes. It peaks in the MIR, drops in the FIR, but rises again at submm wavelengths. We quantify the contribution of the dust heating sources in each 3D dust cell and find a median value of 83% for the star formation component. The AGN contribution is measurable at the percentage level in the disc, but quickly increases in the inner few 100 pc, peaking above 90%. This is the first time the phenomenon of an AGN heating the diffuse dust beyond its torus is quantified in a nearby star-forming galaxy. NGC 1068 only contains a weak AGN, meaning this effect can be stronger in galaxies with a more luminous AGN. This could significantly impact the derived star formation rates and ISM masses for such systems.

preprint2019arXiv

Probing the cold magnetized Universe with SPICA-POL (B-BOP)

SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars.

preprint2016arXiv

Galactic Cold Cores VII: Filament Formation and Evolution - Methods & Observational Constraints

The association of filaments with protostellar objects has made these structures a priority target in star formation studies. The datasets of the Herschel Galactic Cold Cores Key Programme allow for a statistical study of filaments with a wide range of intrinsic and environmental characteristics. Characterisation of this sample can be used to identify key physical parameters and quantify the role of environment in the formation of supercritical filaments. Filaments were extracted from fields at D<500pc with the getfilaments algorithm and characterised according to their column density profiles and intrinsic properties. Each profile was fitted with a beam-convolved Plummer-like function and quantified based on the relative contributions from the filament 'core', represented by a Gaussian, and 'wing' component, dominated by the power-law of the Plummer-like function. These parameters were examined for populations associated with different background levels. We find that filaments increase their core (Mcore) and wing (Mwing) contributions while increasing their total linear mass density (Mtot). Both components appear to be linked to the local environment, with filaments in higher backgrounds having systematically more massive Mcore and Mwing. This dependence on the environment supports an accretion-based model for filament evolution in the local neighbourhood (D<500pc). Structures located in the highest backgrounds develop the highest central Av, Mcore, and Mwing as Mtot increases with time, favoured by the local availability of material and the enhanced gravitational potential. Our results indicate that filaments acquiring a significantly massive central region with Mcore>Mcrit/2 may become supercritical and form stars. This translates into a need for filaments to become at least moderately self-gravitating in order to undergo localised star formation or become star-forming filaments.

preprint2016arXiv

Mantle formation, coagulation and the origin of cloud/core shine: I. Modelling dust scattering and absorption in the infra-red

Context. The observed cloudshine and coreshine (C-shine) have been explained in terms of grain growth leading to enhanced scatter- ing from clouds in the J, H and K photometric bands and the Spitzer IRAC 3.6 and 4.5 μm bands. Aims. Using our global dust modelling approach THEMIS (The Heterogeneous dust Evolution Model at the IaS) we explore the effects of dust evolution in dense clouds, through aliphatic-rich carbonaceous mantle formation and grain-grain coagulation. Methods. We model the effects of wide band gap a-C:H mantle formation and the low-level aggregation of diffuse interstellar medium dust in the moderately-extinguished outer regions of molecular clouds. Results. The formation of wide band gap a-C:H mantles on amorphous silicate and amorphous carbon (a-C) grains leads to a decrease in their absorption cross-sections but no change in their scattering cross-sections at near-IR wavelengths, resulting in higher albedos. Conclusions. The evolution of dust, with increasing density and extinction in the diffuse to dense molecular cloud transition, through mantle formation and grain aggregation, appears to be a likely explanation for the observed C-shine.

preprint2016arXiv

Mantle formation, coagulation and the origin of cloud/core-shine: II. Comparison with observations

Many dense interstellar clouds are observable in emission in the near-IR, commonly referred to as "Cloudshine", and in the mid-IR, the so-called "Coreshine". These C-shine observations have usually been explained with grain growth but no model has yet been able to self-consistently explain the dust spectral energy distribution from the near-IR to the submm. We want to demonstrate the ability of our new core/mantle evolutionary dust model THEMIS (The Heterogeneous dust Evolution Model at the IaS), which has been shown to be valid in the far-IR and submm, to reproduce the C-shine observations. Our starting point is a physically motivated core/mantle dust model. It consists of 3 dust populations: small aromatic-rich carbon grains; bigger core/mantle grains with mantles of aromatic-rich carbon and cores either made of amorphous aliphatic-rich carbon or amorphous silicate. We assume an evolutionary path where these grains, when entering denser regions, may first form a second aliphatic-rich carbon mantle (coagulation of small grains, accretion of carbon from the gas phase), second coagulate together to form large aggregates, and third accrete gas phase molecules coating them with an ice mantle. To compute the corresponding dust emission and scattering, we use a 3D Monte-Carlo radiative transfer code. We show that our global evolutionary dust modelling approach THEMIS allows us to reproduce C-shine observations towards dense starless clouds. Dust scattering and emission is most sensitive to the cloud central density and to the steepness of the cloud density profile. Varying these two parameters leads to changes, which are stronger in the near-IR, in both the C-shine intensity and profile. With a combination of aliphatic-rich mantle formation and low-level coagulation into aggregates, we can self-consistently explain the observed C-shine and far-IR/submm emission towards dense starless clouds.

preprint2016arXiv

Planck intermediate results. XXIX. All-sky dust modelling with Planck, IRAS, and WISE observations

We present all-sky modelling of the high resolution Planck, IRAS, and WISE infrared (IR) observations using the physical dust model presented by Draine and Li in 2007 (DL). We study the performance and results of this model, and discuss implications for future dust modelling. The present work extends the DL dust modelling carried out on nearby galaxies using Herschel and Spitzer data to Galactic dust emission. We employ the DL dust model to generate maps of the dust mass surface density, the optical extinction Av, and the starlight intensity parametrized by Umin. The DL model reproduces the observed spectral energy distribution (SED) satisfactorily over most of the sky, with small deviations in the inner Galactic disk and in low ecliptic latitude areas. We compare the DL optical extinction Av for the diffuse interstellar medium with optical estimates for 2 10^5 quasi-stellar objects (QSOs) observed in the Sloan digital sky survey. The DL Av estimates are larger than those determined towards QSOs by a factor of about 2, which depends on Umin. The DL fitting parameter Umin, effectively determined by the wavelength where the SED peaks, appears to trace variations in the far-IR opacity of the dust grains per unit Av, and not only in the starlight intensity. To circumvent the model deficiency, we propose an empirical renormalization of the DL Av estimate, dependent of Umin, which compensates for the systematic differences found with QSO observations. This renormalization also brings into agreement the DL Av estimates with those derived for molecular clouds from the near-IR colours of stars in the 2 micron all sky survey. The DL model and the QSOs data are used to compress the spectral information in the Planck and IRAS observations for the diffuse ISM to a family of 20 SEDs normalized per Av, parameterized by Umin, which may be used to test and empirically calibrate dust models.

preprint2016arXiv

Stellar classification of CoRoT targets

The CoRoT faint stars channel observed about 163 600 targets to detect transiting planetary companions. Because CoRoT targets are faint (11< r <16) and close to the galactic plane, only a small subsample has been observed spectroscopically. We describe the latest classification scheme used to derive the spectral type of CoRoT targets, which is based on broadband multi-colour photometry. We assess the accuracy of this spectral classification for the first time. We find that the classification method performs better for stars that were observed during the mission-dedicated photometric ground-based campaigns.The luminosity class is wrong for less than 7% of the targets. Generally, the effective temperature of stars classified as early type (O, B, and A) is overestimated. Conversely, the temperature of stars classified as later type tends to be underestimated. This is mainly due to the adverse effect of interstellar reddening. We find that the median error on the effective temperature is less than 5% for dwarf stars classified with a spectral later than F0, but it is worse for earlier type stars, with up to 20\% error for A and late-B dwarfs, and up to 70% for early-B and O-type dwarfs. Similar results are found for giants, with a median error that is lower than 7% for G- and later type giants, but greater than 25% for earlier types. Overall, we find an average median absolute temperature difference |Delta Teff| = 533+\-6 K for the whole sample of stars classified as dwarfs and |Delta Teff| = 280+\-3 K for the whole sample of giant stars. The corresponding standard deviation is of about 92+\-5 K for dwarfs and 304+\-4 K for giants. Typically for late-type stars, this means that the classification is accurate to about half a class.

preprint2015arXiv

CARMA Observations of Galactic Cold Cores: Searching for Spinning Dust Emission

We present the first search for spinning dust emission from a sample of 34 Galactic cold cores, performed using the CARMA interferometer. For each of our cores we use photometric data from the Herschel Space Observatory to constrain N_{H}, T_{d}, n_{H}, and G_{0}. By computing the mass of the cores and comparing it to the Bonnor-Ebert mass, we determined that 29 of the 34 cores are gravitationally unstable and undergoing collapse. In fact, we found that 6 cores are associated with at least one young stellar object, suggestive of their proto-stellar nature. By investigating the physical conditions within each core, we can shed light on the cm emission revealed (or not) by our CARMA observations. Indeed, we find that only 3 of our cores have any significant detectable cm emission. Using a spinning dust model, we predict the expected level of spinning dust emission in each core and find that for all 34 cores, the predicted level of emission is larger than the observed cm emission constrained by the CARMA observations. Moreover, even in the cores for which we do detect cm emission, we cannot, at this stage, discriminate between free-free emission from young stellar objects and spinning dust emission. We emphasise that, although the CARMA observations described in this analysis place important constraints on the presence of spinning dust in cold, dense environments, the source sample targeted by these observations is not statistically representative of the entire population of Galactic cores.

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

Galactic cold cores V. Dust opacity

The project Galactic Cold Cores has made Herschel observations of interstellar clouds where the Planck satellite survey has located cold and compact clumps. The sources range from starless clumps to protostellar cores. We examine 116 Herschel fields to estimate the submillimetre dust opacity and its variations. The submillimetre dust opacity was derived from Herschel data, and near-infrared observations of the reddening of background stars are converted into near-infrared optical depth. We studied the systematic errors affecting these parameters and used modelling to correct for the expected biases. The ratio of 250um and J band opacities is correlated with the cloud location and star formation activity. We find a median ratio of tau(250um)/tau(J)= (1.6+-0.2)*10^-3, which is more than three times the mean value in diffuse medium. Assuming a spectral index beta=1.8 instead of beta=2.0, the value would be lower by ~30%. No significant systematic variation is detected with Galactocentric distance or with Galactic height. The tau(250um)/tau(J) maps reveal six fields with clear increase of submillimetre opacity of up to tau(250um)/tau(J) ~ 4*10^-3. These are all nearby fields with spatially resolved clumps of high column density. We interpret the increase in the far-infrared opacity as a sign of grain growth in the densest and coldest regions of interstellar clouds.

preprint2015arXiv

New radio observations of anomalous microwave emission in the HII region RCW175

We have observed the HII region RCW175 with the 64m Parkes telescope at 8.4GHz and 13.5GHz in total intensity, and at 21.5GHz in both total intensity and polarization. High angular resolution, high sensitivity, and polarization capability enable us to perform a detailed study of the different constituents of the HII region. For the first time, we resolve three distinct regions at microwave frequencies, two of which are part of the same annular diffuse structure. Our observations enable us to confirm the presence of anomalous microwave emission (AME) from RCW175. Fitting the integrated flux density across the entire region with the currently available spinning dust models, using physically motivated assumptions, indicates the presence of at least two spinning dust components: a warm component with a relatively large hydrogen number density n_H=26.3/cm^3 and a cold component with a hydrogen number density of n_H=150/cm^3. The present study is an example highlighting the potential of using high angular-resolution microwave data to break model parameter degeneracies. Thanks to our spectral coverage and angular resolution, we have been able to derive one of the first AME maps, at 13.5GHz, showing clear evidence that the bulk of the AME arises in particular from one of the source components, with some additional contribution from the diffuse structure. A cross-correlation analysis with thermal dust emission has shown a high degree of correlation with one of the regions within RCW175. In the center of RCW175, we find an average polarized emission at 21.5GHz of 2.2\pm0.2(rand.)\pm0.3(sys.)% of the total emission, where we have included both systematic and statistical uncertainties at 68% CL. This polarized emission could be due to sub-dominant synchrotron emission from the region and is thus consistent with very faint or non-polarized emission associated with AME.

preprint2015arXiv

Planck intermediate results. XXXV. Probing the role of the magnetic field in the formation of structure in molecular clouds

Within ten nearby (d < 450 pc) Gould Belt molecular clouds we evaluate statistically the relative orientation between the magnetic field projected on the plane of sky, inferred from the polarized thermal emission of Galactic dust observed by Planck at 353 GHz, and the gas column density structures, quantified by the gradient of the column density, $N_H$. The selected regions, covering several degrees in size, are analyzed at an effective angular resolution of 10' FWHM, thus sampling physical scales from 0.4 to 40 pc in the nearest cloud. The column densities in the selected regions range from $N_H \approx 10^{21}$ to $10^{23}$ cm$^{-2}$, and hence they correspond to the bulk of the molecular clouds. The relative orientation is evaluated pixel by pixel and analyzed in bins of column density using the novel statistical tool called "Histogram of Relative Orientations". Throughout this study, we assume that the polarized emission observed by Planck at 353 GHz is representative of the projected morphology of the magnetic field in each region, i.e., we assume a constant dust grain alignment efficiency, independent of the local environment. Within most clouds we find that the relative orientation changes progressively with increasing $N_H$, from preferentially parallel or having no preferred orientation to preferentially perpendicular. In simulations of magnetohydrodynamic turbulence in molecular clouds this trend in relative orientation is a signature of Alfvénic or sub-Alfvénic turbulence, implying that the magnetic field is significant for the gas dynamics at the scales probed by Planck. We compare the deduced magnetic field strength with estimates we obtain from other methods and discuss the implications of the Planck observations for the general picture of molecular cloud formation and evolution.

preprint2015arXiv

Using cm Observations to Constrain the Abundance of Very Small Dust Grains in Galactic Cold Cores

In this analysis we illustrate how the relatively new emission mechanism known as spinning dust can be used to characterize dust grains in the interstellar medium. We demonstrate this by using spinning dust emission observations to constrain the abundance of very small dust grains (a $\lesssim$ 10nm) in a sample of Galactic cold cores. Using the physical properties of the cores in our sample as inputs to a spinning dust model, we predict the expected level of emission at a wavelength of 1cm for four different very small dust grain abundances, which we constrain by comparing to 1cm CARMA observations. For all of our cores we find a depletion of very small grains, which we suggest is due to the process of grain growth. This work represents the first time that spinning dust emission has been used to constrain the physical properties of interstellar dust grains.

preprint2014arXiv

A hidden reservoir of Fe/FeS in interstellar silicates?

The depletion of iron and sulphur into dust in the interstellar medium and the exact nature of interstellar amorphous silicate grains is still an open question. We study the incorporation of iron and sulphur into amorphous silicates of olivine- and pyroxene-type and their effects on the dust spectroscopy and thermal emission. We used the Maxwell-Garnett effective-medium theory to construct the optical constants for a mixture of silicates, metallic iron, and iron sulphide. We also studied the effects of iron and iron sulphide in aggregate grains. Iron sulphide inclusions within amorphous silicates that contain iron metal inclusions shows no strong differences in the optical properties of the grains. A mix of amorphous olivine- and pyroxene-type silicate broadens the silicate features. An amorphous carbon mantle with a thickness of 10 nm on the silicate grains leads to an increase in absorption on the short-wavelength side of the 10 $μ$m silicate band. The assumption of amorphous olivine-type and pyroxene-type silicates and a 10 nm thick amorphous carbon mantle better matches the interstellar silicate band profiles. Including iron nano-particles leads to an increase in the mid-IR extinction, while up to 5 ppm of sulphur can be incorporated as Fe/FeS nano inclusions into silicate grains without leaving a significant trace of its presence.

preprint2014arXiv

Multiwavelength study of the high-latitude cloud L1642: chain of star formation

L1642 is one of the two high galactic latitude (|b| > 30deg) clouds confirmed to have active star formation. We examine the properties of this cloud, especially the large-scale structure, dust properties, and compact sources in different stages of star formation. We present high-resolution far-infrared and submm observations with the Herschel and AKARI satellites and mm observations with the AzTEC/ASTE telescope, which we combined with archive data from near- and mid-infrared (2MASS, WISE) to mm observations (Planck). The Herschel observations, combined with other data, show a sequence of objects from a cold clump to young stellar objects at different evolutionary stages. Source B-3 (2MASS J04351455-1414468) appears to be a YSO forming inside the L1642 cloud, instead of a foreground brown dwarf, as previously classified. Herschel data reveal striation in the diffuse dust emission around L1642. The western region shows striation towards NE and has a steeper column density gradient on its southern side. The densest central region has a bow-shock like structure showing compression from the west and a filamentary tail extending towards east. The differences suggest that these may be spatially distinct structures, aligned only in projection. We derive values of the dust emission cross-section per H nucleon for different regions of the cloud. Modified black-body fits to the spectral energy distribution of Herschel and Planck data give emissivity spectral index beta values 1.8-2.0 for the different regions. The compact sources have lower beta values and show an anticorrelation between T and beta. Markov chain Monte Carlo calculations demonstrate the strong anticorrelation between beta and T errors and the importance of mm Planck data in constraining the estimates. L1642 reveals a more complex structure and sequence of star formation than previously known.

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

The dust masses of powerful radio galaxies: clues to the triggering of their activity

We use deep Herschel Space Observatory observations of a 90% complete sample of 32 intermediate-redshift 2Jy radio galaxies (0.05 < z < 0.7) to estimate the dust masses of their host galaxies and thereby investigate the triggering mechanisms for their quasar-like AGN. The dust masses derived for the radio galaxies (7.2x10^5 < M_d < 2.6x10^8 M_sun) are intermediate between those of quiescent elliptical galaxies on the one hand, and ultra luminous infrared galaxies (ULIRGs) on the other. Consistent with simple models for the co-evolution of supermassive black holes and their host galaxies, these results suggest that most of the radio galaxies represent the late time re-triggering of AGN activity via mergers between the host giant elliptical galaxies and companion galaxies with relatively low gas masses. However, a minority of the radio galaxies in our sample (~20%) have high, ULIRG-like dust masses, along with evidence for prodigious star formation activity. The latter objects are more likely to have been triggered in major, gas-rich mergers that represent a rapid growth phase for both their host galaxies and their supermassive black holes.

preprint2014arXiv

The evolution of amorphous hydrocarbons in the ISM: dust modelling from a new vantage point

Context. The evolution of amorphous hydrocarbon materials, a-C(:H), principally resulting from ultraviolet (UV) photon absorption- induced processing, are likely at the heart of the variations in the observed properties of dust in the interstellar medium. Aims. The consequences of the size-dependent and compositional variations in a-C(:H), from aliphatic-rich a-C:H to aromatic-rich a-C, are studied within the context of the interstellar dust extinction and emission. Methods. Newly-derived optical property data for a-C(:H) materials, combined with that for an amorphous forsterite-type silicate with iron nano-particle inclusions, a-SilFe, are used to explore dust evolution in the interstellar medium. Results. We present a new dust model that consists of a power-law distribution of small a-C grains and log-normal distributions of large a-SilFe and a-C(:H) grains. The model, which is firmly anchored by laboratory-data, is shown to quite naturally explain the variations in the infrared (IR) to far-ultraviolet (FUV) extinction, the 217 nm UV bump, the IR absorption and emission bands and the IR-mm dust emission. Conclusions. The major strengths of the new model are its inherent simplicity and built-in capacity to follow dust evolution in interstellar media. We show that mantle accretion in molecular clouds and UV photo-processing in photo-dominated regions are likely the major drivers of dust evolution.

preprint2013arXiv

Planck intermediate results. XII: Diffuse Galactic components in the Gould Belt System

We perform an analysis of the diffuse low-frequency Galactic components in the Southern part of the Gould Belt system (130^\circ\leq l\leq 230^\circ and -50^\circ\leq b\leq -10^\circ). Strong ultra-violet (UV) flux coming from the Gould Belt super-association is responsible for bright diffuse foregrounds that we observe from our position inside the system and that can help us improve our knowledge of the Galactic emission. Free-free emission and anomalous microwave emission (AME) are the dominant components at low frequencies (ν< 40 GHz), while synchrotron emission is very smooth and faint. We separate diffuse free-free emission and AME from synchrotron emission and thermal dust emission by using Planck data, complemented by ancillary data, using the "Correlated Component Analysis" (CCA) component separation method and we compare with the results of cross-correlation of foreground templates with the frequency maps. We estimate the electron temperature T_e from H$α$ and free-free emission using two methods (temperature-temperature plot and cross-correlation) and we obtain T_e ranging from 3100 to 5200 K, for an effective fraction of absorbing dust along the line of sight of 30% (f_d=0.3). We estimate the frequency spectrum of the diffuse AME and we recover a peak frequency (in flux density units) of 25.5 \pm 1.5 GHz. We verify the reliability of this result with realistic simulations that include the presence of biases in the spectral model for the AME and in the free-free template. By combining physical models for vibrational and rotational dust emission and adding the constraints from the thermal dust spectrum from Planck and IRAS we are able to get a good description of the frequency spectrum of the AME for plausible values of the local density and radiation field.

preprint2013arXiv

The degeneracy between dust colour temperature and spectral index. Comparison of methods for estimating the beta(T) relation

Sub-millimetre dust emission provides information on the physics of interstellar clouds and dust. Noise can produce anticorrelation between the colour temperature T_C and the spectral index beta. This must be separated from the intrinsic beta(T) relation of dust. We compare methods for the analysis of the beta(T) relation. We examine sub-millimetre observations simulated as simple modified black body emission or using 3D radiative transfer modelling. In addition to chi^2 fitting, we examine the results of the SIMEX method, basic Bayesian model, hierarchical models, and one method that explicitly assumes a functional form for beta(T). All methods exhibit some bias. Bayesian method shows significantly lower bias than direct chi^2 fits. The same is true for hierarchical models that also result in a smaller scatter in the temperature and spectral index values. However, significant bias was observed in cases with high noise levels. Beta and T estimates of the hierarchical model are biased towards the relation determined by the data with the highest S/N ratio. This can alter the recovered beta(T) function. With the method where we explicitly assume a functional form for the beta(T) relation, the bias is similar to the Bayesian method. In the case of an actual Herschel field, all methods agree, showing some degree of anticorrelation between T and beta. The Bayesian method and the hierarchical models can both reduce the noise-induced parameter correlations. However, all methods can exhibit non-negligible bias. This is particularly true for hierarchical models and observations of varying signal-to-noise ratios and must be taken into account when interpreting the results.

preprint2012arXiv

Galactic cold cores III. General cloud properties

We examine the cloud structure around the Planck detections in 71 fields observed with the Herschel SPIRE instrument. We wish to determine the general physical characteristics of the fields and to examine the morphology of the clouds where the cold high column density clumps are found. We derive colour temperature and column density maps of the fields. We examine the infrared spectral energy distributions of the main clumps. The clouds are categorised according to their large scale morphology. With the help of recently released WISE satellite data, we look for signs of enhanced mid-infrared scattering (coreshine), an indication of growth of the dust grains, and examine the star formation activity associated with the cold clumps. The mapped clouds have distances ranging from ~100pc to several kiloparsecs and cover a range of sizes and masses from cores of less than 10 solar masses to clouds with masses in excess of 10000 solar mass. Most fields contain some filamentary structures and in about half of the cases a filament or a few filaments dominate the morphology. In one case out of ten, the clouds show a cometary shape or have sharp boundaries indicative of compression by an external force. The width of the filaments is typically ~0.2-0.3pc. However, there is significant variation from 0.1pc to 1pc and the estimates are sensitive to the methods used and the very definition of a filament. Enhanced mid-infrared scattering, coreshine, was detected in four clouds with six additional tentative detections. The cloud LDN183 is included in our sample and remains the best example of this phenomenon. About half of the fields are associated with active star formation as indicated by the presence of mid-infrared point sources. The mid-infrared sources often coincide with structures whose sub-millimetre spectra are still dominated by the cold dust.

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

Reliability of NH3 as the temperature probe of cold cloud cores

The temperature is a central parameter affecting the chemical and physical properties of dense cores of interstellar clouds and their evolution to star formation. The chemistry and the dust properties are temperature dependent and the interpretation of observation requires the knowledge of the temperature and its variations. Measurement of the gas kinetic temperature is possible with molecular line spectroscopy, the ammonia molecule, NH3, being the most commonly used tracer. We want to determine the accuracy of the temperature estimates derived from ammonia spectra. The normal interpretation of NH3 observations assumes that all the hyperfine line components are tracing the same gas volume. In the case of temperature gradients they may be sensitive to different layers and cause errors in the optical depth and gas temperature estimates. We examine a series of spherical cloud models, 1.0 and 0.5 M_Sun Bonnor-Ebert spheres, with different radial temperature profiles. We calculate synthetic NH3 spectra and compare the derived column densities and temperatures to the true values. For high signal-to-noise observations, the estimated gas kinetic temperatures are within ~0.3 K of the real mass averaged temperature and the column densities are correct to within ~10%. When the S/N ratio of the (2,2) spectrum decreases below 10, the temperature errors are of the order of 1K but without a significant bias. When the density of the models is increased by a factor of a few, the results begin to show significant bias because of the saturation of the (1,1) main group. The ammonia spectra are found to be a reliable tracer of the mass averaged gas temperature. Because the radial temperature profiles of the cores are not well constrained, the central temperature could still differ from this value. If the cores are optically very thick, there are no guarantees of the accuracy.

preprint2012arXiv

The degeneracy between the dust colour temperature and the spectral index. The problem of multiple chi^2 minima

Because of the Herschel and Planck satellite missions, there is strong interest in the interpretation the sub-millimetre dust spectra from interstellar clouds. Much work has been done to understand the dependence between the spectral index beta_Obs and the colour temperature T_C that is partly caused by the noise. The (T_C, beta_Obs) confidence regions are elongated, banana-shaped structures. We studied under which conditions these exhibit anomalous, strongly non-Gaussian behaviour that could affect the interpretation of the observed (T_C, beta_Obs) relations. We examined modified black body spectra and spectra calculated from radiative transfer models of filamentary clouds at wavelengths 100um-850um. We performed modified black body fits and examined the structure of the chi^2(T_, beta_Obs) function. We show cases where, when the signal-to-noise ratio is low, the chi^2 has multiple local minima in the (T_C, beta_Obs) plane. A small change in the weighting of the data points can cause the solution to jump to completely different values. In particular, noise can lead to the appearance of a separate population of solutions with low colour temperatures and high spectral indices. The anomalies are caused by the noise but the tendency to show multiple chi^2 minima depends on the model and the wavelengths analysed. Deviations from the assumed single modified black body spectrum are not important. The presence of local minima implies that the results obtained from the chi^2 minimisation depend on the starting point of the optimisation and may correspond to non-global minima. The (T_C,beta_Obs) distributions may be contaminated by a few solutions with unrealistically low colour temperatures and high spectral indices. Proper weighting must be applied to avoid the determination of the beta_Obs(T_C) relation to be unduly affected by these measurements.

preprint2012arXiv

The effect of temperature mixing on the observable (T,beta)-relation of interstellar dust clouds

Detailed studies of the shape of dust emission spectra are possible thanks to the current instruments capable of observations in several sub-millimetre bands (e.g., Herschel and Planck). However, some controversy remains even on the basic effects resulting from the mixing of temperatures along the line-of-sight. Studies have suggested either a positive or a negative correlation between the colour temperature T_C and the observed spectral index beta_Obs. Our aim is to show that both cases are possible and to determine the factors leading to either behaviour. We start by studying the sum of two or three modified black bodies of different temperature. With radiative transfer modelling, we examine the probability distributions of the dust mass as a function of the physical dust temperature. With these results as a guideline, we examine the (T_C, beta_Obs) relations for different sets of clouds. Even in the case of modified blackbodies at temperatures T_0 and T_0+ Delta T_0, the correlation between T_C and beta_Obs can be either positive or negative. If one compares models where Delta T_0 is varied, the correlation is negative. If the models differ in their mean temperature T_0 rather than in Delta T_0, the correlation remains positive. Radiative transfer models show that externally heated clouds have different mean temperatures but the widths of their temperature distributions are rather similar. Thus, the correlation between T_C and beta_Obs is expected to be positive. The same result applies to clouds illuminated by external radiation fields of different intensity. For internally heated clouds a negative correlation is the more likely alternative. If the signal-to-noise ratio is high, the observed negative correlation could be explained by the temperature dependence of the dust optical properties but that intrinsic dependence could be even steeper than the observed one.

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

Modelling the spinning dust emission from dense interstellar clouds

Electric dipole emission arising from PAHs is often invoked to explain the anomalous microwave emission (AME). This assignation is based on an observed tight correlation between the mid-IR emission of PAHs and the AME; and a good agreement between models of spinning dust and the broadband AME spectrum. So far often detected at large scale in the diffuse interstellar medium, the AME has recently been studied in detail in well-known dense molecular clouds with the help of Planck data. While much attention has been given to the physics of spinning dust emission, the impact of varying local physical conditions has not yet been considered in detail. Our aim is to study the emerging spinning dust emission from interstellar clouds with realistic physical conditions and radiative transfer. We use the DustEM code from Compiegne et al. to describe the extinction and IR emission of all dust populations. The spinning dust emission is obtained with SpDust, as described by Silsbee et al., that we have coupled to DustEM. We carry out full radiative transfer simulations and carefully estimate the local gas state as a function of position within interstellar clouds. We show that the spinning dust emission is sensitive to the abundances of the major ions and we propose a simple scheme to estimate these abundances. We also investigate the effect of changing the cosmic-ray rate. In dense media, where radiative transfer is mandatory, we show that the relationship between the spinning and mid-IR emissivities of PAHs is no longer linear and that the spinning dust emission may actually be strong at the centre of clouds where the mid-IR PAH emission is weak. These results provide new ways to trace grain growth from diffuse to dense medium and will be useful for the analysis of AME at the scale of interstellar clouds.

preprint2011arXiv

On The Gas Temperature of Molecular Cloud Cores

We investigate the uncertainties affecting the temperature profiles of dense cores of interstellar clouds. In regions shielded from external ultraviolet radiation, the problem is reduced to the balance between cosmic ray heating, line cooling, and the coupling between gas and dust. We show that variations in the gas phase abundances, the grain size distribution, and the velocity field can each change the predicted core temperatures by one or two degrees. We emphasize the role of non-local radiative transfer effects that often are not taken into account, for example, when modelling the core chemistry. These include the radiative coupling between regions of different temperature and the enhanced line cooling near the cloud surface. The uncertainty of the temperature profiles does not necessarily translate to a significant error in the column density derived from observations. However, depletion processes are very temperature sensitive and a two degree difference can mean that a given molecule no longer traces the physical conditions in the core centre.

preprint2011arXiv

Planck Early Results: Origin of the submm excess dust emission in the Magellanic Clouds

The integrated Spectral Energy Distributions of the Large and Small Magellanic Cloud appear significantly flatter than expected from dust models based on their FIR and radio emission. The origin of this millimetre excess is still unexplained, and is here investigated using the Planck data. The background CMB contribution is subtracted using an ILC method performed locally around the galaxies. The foreground emission from the Milky Way is subtracted. After subtraction, the emission of both galaxies correlates closely with the gas emission of the LMC and SMC. The millimetre excess in the LMC can be explained by CMB fluctuations, but a significant excess is still present in the SMC SED. The Planck and IRIS data at 100 micron are combined to produce thermal dust temperature and optical depth maps of the two galaxies. The LMC temperature map shows the presence of a warm inner arm already found with the Spitzer data, but also shows the existence of a previously unidentified cold outer arm. Several cold regions are found along this arm, some of which are associated with known molecular clouds. The average emissivity spectral index is found to be consistent with beta=1.5 and beta=1.2 below 500 microns for the LMC and SMC respectively, significantly flatter than the values observed in the Milky Way. Furthermore, there is evidence in the SMC for a further flattening of the SED in the sub-mm. The spatial distribution of the millimetre dust excess in the SMC follows the gas and thermal dust distribution. Different models are explored in order to fit the dust emission in the SMC. It is concluded that the millimetre excess is unlikely to be caused by very cold dust emission and that it could be due to a combination of spinning dust. emission and thermal dust emission by more amorphous dust grains than those present in our Galaxy.

preprint2011arXiv

Planck Early Results: The Galactic Cold Core Population revealed by the first all-sky survey

We present the statistical properties of the first version of the Cold Core Catalogue of Planck Objects (C3PO), in terms of their spatial distribution, temperature, distance, mass, and morphology. We also describe the statistics of the Early Cold Core Catalogue (ECC, delivered with the Early Release Compact Source Catalogue, ERCSC) that is the subset of the 915 most reliable detections of the complete catalogue. We have used the CoCoCoDeT algorithm to extract 10783 cold sources. Temperature and dust emission spectral index β values are derived using the fluxes in the IRAS 100 \mum band and the three highest frequency Planck bands. Temperature spans from 7K to 17K, and peaks around 13K. Data are not consistent with a constant value of β over the all temperature range. β ranges from 1.4 to 2.8 with a mean value around 2.1, and several possible scenarios are possible, including β(T) and the effect of multiple T components folded into the measurements. For one third of the objects the distances are obtained. Most of the detections are within 2 kpc in the Solar neighbourhood, but a few are at distances greater than 4 kpc. The cores are distributed from the deep Galactic plane, despite the confusion, to high latitudes (>30$^{\circle}$). The associated mass estimates range from 1 to $10^5$ solar masses. Using their physical properties these cold sources are shown to be cold clumps, defined as the intermediate cold sub-structures between clouds and cores. These cold clumps are not isolated but mostly organized in filaments associated with molecular clouds. The Cold Core Catalogue of Planck Objects (C3PO) is the first unbiased all-sky catalogue of cold objects. It gives an unprecedented statistical view to the properties of these potential pre-stellar clumps and offers a unique possibility for their classification in terms of their intrinsic properties and environment.

preprint2011arXiv

Planck Early Results: The submillimetre properties of a sample of Galactic cold clumps

(abridged) We perform a detailed investigation of sources from the Cold Cores Catalogue of Planck Objects (C3PO). Our goal is to probe the reliability of the detections, validate the separation between warm and cold dust emission components, provide the first glimpse at the nature, internal morphology and physical characterictics of the Planck-detected sources. We focus on a sub-sample of ten sources from the C3PO list, selected to sample different environments, from high latitude cirrus to nearby (150pc) and remote (2kpc) molecular complexes. We present Planck surface brightness maps and derive the dust temperature, emissivity spectral index, and column densities of the fields. With the help of higher resolution Herschel and AKARI continuum observations and molecular line data, we investigate the morphology of the sources and the properties of the substructures at scales below the Planck beam size.

preprint2011arXiv

Planck early results. XX. New light on anomalous microwave emission from spinning dust grains

Anomalous microwave emission (AME) has been observed by numerous experiments in the frequency range ~10-60 GHz. Using Planck maps and multi-frequency ancillary data, we have constructed spectra for two known AME regions: the Perseus and Rho Ophiuchi molecular clouds. The spectra are well fitted by a combination of free-free radiation, cosmic microwave background, thermal dust, and electric dipole radiation from small spinning dust grains. The spinning dust spectra are the most precisely measured to date, and show the high frequency side clearly for the first time. The spectra have a peak in the range 20-40 GHz and are detected at high significances of 17.1 sigma for Perseus and 8.4 sigma for Rho Ophiuchi. In Perseus, spinning dust in the dense molecular gas can account for most of the AME; the low density atomic gas appears to play a minor role. In Rho Ophiuchi, the ~30 GHz peak is dominated by dense molecular gas, but there is an indication of an extended tail at frequencies 50-100 GHz, which can be accounted for by irradiated low density atomic gas. The dust parameters are consistent with those derived from other measurements. We have also searched the Planck map at 28.5 GHz for candidate AME regions, by subtracting a simple model of the synchrotron, free-free, and thermal dust. We present spectra for two of the candidates; S140 and S235 are bright HII regions that show evidence for AME, and are well fitted by spinning dust models.

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.

preprint2010arXiv

Sub-millimeter to centimeter excess emission from the Magellanic Clouds. I. Global spectral energy distribution

In order to reconstruct the global SEDs of the Magellanic Clouds over eight decades in spectral range, we combined literature flux densities representing the entire LMC and SMC respectively, and complemented these with maps extracted from the WMAP and COBE databases covering the missing the 23--90 GHz (13--3.2 mm) and the poorly sampled 1.25--250 THz (240--1.25 micron). We have discovered a pronounced excess of emission from both Magellanic Clouds, but especially the SMC, at millimeter and sub-millimeter wavelengths. We also determined accurate thermal radio fluxes and very low global extinctions for both LMC and SMC. Possible explanations are briefly considered but as long as the nature of the excess emission is unknown, the total dust masses and gas-to-dust ratios of the Magellanic Clouds cannot reliably be determined.

preprint2010arXiv

The global dust SED: Tracing the nature and evolution of dust with DustEM

The Planck and Herschel missions are currently measuring the farIR-mm emission of dust, which combined with existing IR data, will for the first time provide the full SED of the galactic ISM dust emission with an unprecedented sensitivity and angular resolution. It will allow a systematic study of the dust evolution processes that affect the SED. Here we present a versatile numerical tool, DustEM, that predicts the emission and extinction of dust given their size distribution and their optical and thermal properties. In order to model dust evolution, DustEM has been designed to deal with a variety of grain types, structures and size distributions and to be able to easily include new dust physics. We use DustEM to model the dust SED and extinction in the diffuse interstellar medium at high-galactic latitude (DHGL), a natural reference SED. We present a coherent set of observations for the DHGL SED. The dust components in our DHGL model are (i) PAHs, (ii) amorphous carbon and (iii) amorphous silicates. We use amorphous carbon dust, rather than graphite, because it better explains the observed high abundances of gas-phase carbon in shocked regions of the interstellar medium. Using the DustEM model, we illustrate how, in the optically thin limit, the IRAS/Planck HFI (and likewise Spitzer/Herschel for smaller spatial scales) photometric band ratios of the dust SED can disentangle the influence of the exciting radiation field intensity and constrain the abundance of small grains relative to the larger grains. We also discuss the contributions of the different grain populations to the IRAS, Planck and Herschel channels. Such information is required to enable a study of the evolution of dust as well as to systematically extract the dust thermal emission from CMB data and to analyze the emission in the Planck polarized channels. The DustEM code described in this paper is publically available.