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.

Context: Dust in late-type galaxies in the local Universe is responsible for absorbing approximately one third of the energy emitted by stars. It is often assumed that dust heating is mainly attributable to the absorption of UV and optical photons emitted by the youngest (<= 100 Myr) stars. Consequently, thermal re-emission by dust at FIR wavelengths is often linked to the star-formation activity of a galaxy. However, several studies argue that the contribution to dust heating by much older stars might be more significant. Advances in radiation transfer (RT) simulations finally allow us to actually quantify the heating mechanisms of diffuse dust by the stellar radiation field. Aims: As one of the main goals in the DustPedia project, we have constructed detailed 3D stellar and dust RT models for nearby galaxies. We analyse the contribution of the different stellar populations to the dust heating in four face-on barred galaxies: NGC1365, M83, M95, and M100. We aim to quantify the fraction directly related to young stars, both globally and on local scales, and to assess the influence of the bar on the heating fraction. Results: We derive global attenuation laws for each galaxy and confirm that galaxies of high sSFR have shallower attenuation curves and weaker UV bumps. On average, 36.5% of the bolometric luminosity is absorbed by dust. We report a clear effect of the bar structure on the radial profiles of the dust-heating fraction by the young stars, and the dust temperature. We find that the young stars are the main contributors to the dust heating, donating, on average ~59% of their luminosity to this purpose throughout the galaxy. This dust-heating fraction drops to ~53% in the bar region and ~38% in the bulge region where the old stars are the dominant contributors to the dust heating. We also find a strong link between the heating fraction by the young stars and the sSFR.

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 {\mu}m), Spitzer/MIPS (24 {\mu}m) together with Herschel/PACS (70 and 160 {\mu}m) and Herschel/SPIRE (250, 350 and 500 {\mu}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.

The size and composition of dust grains are critical in setting the dynamical, physical and chemical evolution of the media in which they are present. Thanks to facilities such as ALMA and in the future the SKA, their thermal emission in the (sub)mm to cm has become a convenient way to trace grain properties. Our aim is to understand the influence of the grain composition and size distribution on the shape of their SED in dense ISM regions such as molecular clouds, prestellar cores, YSOs and protoplanetary discs. Starting from the optical constants defined in the THEMIS model for amorphous hydrogenated carbon and silicate grains in addition to water ice, we define 6 material mixtures representative of the expected dust composition in dense ISM regions. The optical properties of 0.01 micron to 10 cm grains are then calculated with effective medium and Mie theories. The corresponding SEDs are calculated for isolated clouds either externally heated by the ISRF alone or in addition to an internal source. The 3 main outcomes of this study are as follows. First, the dust mass absorption coefficient strongly depends on both its composition and size distribution potentially leading to errors in dust mass estimates by factors up to 3 and 20, respectively. Second, it appears almost impossible to retrieve the grain composition from the (sub)mm to cm SED shape alone as its spectral index for lambda > 3 mm does not depend on composition. Third, using the true dust opacity spectral index to estimate grain sizes may lead to erroneous findings as the observed spectral index can be highly modified by the dust temperature distribution along the line-of-sight, which depends on the specific heating source and on the geometry of the studied region. Based on the interpretation of only the spectral shape of (sub)mm to cm SEDs, the determination of the dust masses, compositions and sizes are highly uncertain.

Observations of evolution in the dust-to-metal ratio allow us to constrain the dominant dust processing mechanisms. In this work, we present a study of the dust-to-metal and dust-to-gas ratios in a sub-sample of ~500 DustPedia galaxies. Using literature and MUSE emission line fluxes, we derived gas-phase metallicities (oxygen abundances) for over 10000 individual regions and determine characteristic metallicities for each galaxy. We study how the relative dust, gas, and metal contents of galaxies evolve by using metallicity and gas fraction as proxies for evolutionary state. The global oxygen abundance and nitrogen-to-oxygen ratio are found to increase monotonically as galaxies evolve. Additionally, unevolved galaxies (gas fraction > 60%, metallicity 12 + log(O/H) < 8.2) have dust-to-metal ratios that are about a factor of 2.1 lower (a factor of six lower for galaxies with gas fraction > 80%) than the typical dust-to-metal ratio (Md/MZ ~ 0.214) for more evolved sources. However, for high gas fractions, the scatter is larger due to larger observational uncertainties as well as a potential dependence of the dust grain growth timescale and supernova dust yield on local conditions and star formation histories. We find chemical evolution models with a strong contribution from dust grain growth describe these observations reasonably well. The dust-to-metal ratio is also found to be lower for low stellar masses and high specific star formation rates (with the exception of some sources undergoing a starburst). Finally, the metallicity gradient correlates weakly with the HI-to-stellar mass ratio, the effective radius and the dust-to-stellar mass ratio, but not with stellar mass.

The dust mass absorption coefficient, $\kappa_{d}$, is the conversion function used to infer physical dust masses from observations of dust emission. However, it is notoriously poorly constrained, and it is highly uncertain how it varies, either between or within galaxies. Here we present the results of a proof-of concept study, using the DustPedia data for two nearby face-on spiral galaxies M74 (NGC 628) and M83 (NGC 5236), to create the first ever maps of $\kappa_{d}$ in galaxies. We determine $\kappa_{d}$ using an empirical method that exploits the fact that the dust-to-metals ratio of the interstellar medium is constrained by direct measurements of the depletion of gas-phase metals. We apply this method pixel-by-pixel within M74 and M83, to create maps of $\kappa_{d}$. We also demonstrate a novel method of producing metallicity maps for galaxies with irregularly-sampled measurements, using the machine learning technique of Gaussian process regression. We find strong evidence for significant variation in $\kappa_{d}$. We find values of $\kappa_{d}$ at 500 $\mu$m spanning the range 0.11-0.25 ${\rm m^{2}\,kg^{-1}}$ in M74, and 0.15-0.80 ${\rm m^{2}\,kg^{-1}}$ in M83. Surprisingly, we find that $\kappa_{d}$ shows a distinct inverse correlation with the local density of the interstellar medium. This inverse correlation is the opposite of what is predicted by standard dust models. However, we find this relationship to be robust against a large range of changes to our method - only the adoption of unphysical or highly unusual assumptions would be able to suppress it.

Aims: We compare the far-infrared to sub-millimetre dust emission properties measured in high Galactic latitude cirrus with those determined in a sample of 204 late-type DustPedia galaxies. The aim is to verify if it is appropriate to use Milky Way dust properties to derive dust masses in external galaxies. Methods: We used Herschel observations and atomic and molecular gas masses to estimate the disc-averaged dust emissivity at 250 micrometres, and from this, the absorption cross section per H atom and per dust mass. The emissivity requires one assumption, which is the CO-to-H_2 conversion factor, and the dust temperature is additionally required for the absorption cross section per H atom; yet another constraint on the dust-to-hydrogen ratio D/H, depending on metallicity, is required for the absorption cross section dust mass. Results: We find epsilon(250) = 0.82 +/- 0.07 MJy sr^-1 (1E20 H cm^-2)^-1 for galaxies with 4 < F(250)/F(500) < 5. This depends only weakly on the adopted CO-to-H_2 conversion factor. The value is almost the same as that for the Milky Way at the same colour ratio. Instead, for F(250)/F(500) > 6, epsilon(250) is lower than predicted by its dependence on the heating conditions. The reduction suggests a variation in dust emission properties for spirals of earlier type, higher metallicity, and with a higher fraction of molecular gas. When the standard emission properties of Galactic cirrus are used for these galaxies, their dust masses might be underestimated by up to a factor of two. Values for the absorption cross sections at the Milky Way metallicity are also close to those of the cirrus. Mild trends of the absorption cross sections with metallicity are found, although the results depend on the assumptions made.

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.

We use a sub-set of the DustPedia galaxy sample (461 galaxies) to investigate the effect the environment has had on galaxies. We consider Virgo cluster and field samples and also assign a density contrast parameter to each galaxy, as defined by the local density of SDSS galaxies. We consider their chemical evolution (using M_{Dust}/M_{Baryon} and M_{Gas}/M_{Baryon}), their specific star formation rate (SFR/M_{Stars}), star formation efficiency (SFR/M_{Gas}), stars-to-dust mass ratio (M_{Stars}/M_{Dust}), gas-to-dust mass ratio (M_{Gas}/M_{Dust}) and the relationship between star formation rate per unit mass of dust and dust temperature (SFR/M_{Dust} and T_{Dust}). Late type galaxies (later than Sc) in all of the environments can be modelled using simple closed box chemical evolution and a simple star formation history (SFR(t) \propto t\exp{-t/\tau}). For earlier type galaxies the physical mechanisms that give rise to their properties are clearly much more varied and require a more complicated model (mergers, gas in or outflow). However, we find little or no difference in the properties of galaxies of the same morphological type within the cluster, field or with different density contrasts. It appears that it is morphology, how and whenever this is laid down, and consistent internal physical processes that primarily determine the derived properties of galaxies in the DustPedia sample and not processes related to differences in the local environment.

Dust offers a unique probe of the interstellar medium (ISM) across multiple size, density, and temperature scales. Dust is detected in outflows of evolved stars, star-forming molecular clouds, planet-forming disks, and even in galaxies at the dawn of the Universe. These grains also have a profound effect on various astrophysical phenomena from thermal balance and extinction in galaxies to the building blocks for planets, and changes in dust grain properties will affect all of these phenomena. A full understanding of dust in all of its forms and stages requires a multi-disciplinary investigation of the dust life cycle. Such an investigation can be achieved with a statistical study of dust properties across stellar evolution, star and planet formation, and redshift. Current and future instrumentation will enable this investigation through fast and sensitive observations in dust continuum, polarization, and spectroscopy from near-infrared to millimeter wavelengths.

In the interstellar medium, carbon (nano-)grains are a major component of interstellar dust. This solid phase is more vulnerable to processing and destruction than its silicate counterpart. It exhibits a complex, size-dependent evolution due to interactions within different radiative and dynamical environments. Infrared signatures of these nanocarbon grains are seen in a large number of disks around Herbig HAeBe stars. We probe the composition and evolution of carbon nano-grains at the surface of (pre-)transitional protoplanetary disks around Herbig stars. We present spatially resolved infrared emission spectra obtained with NAOS CONICA at the VLT in the 3-4 $\mu$m range with a spatial resolution of 0.1", which allow us to trace aromatic, olefinic and aliphatic bands which are attributed to sub-nanometer hydrocarbon grains. We apply a gaussian fitting to analyse the observed spectral signatures. Finally, we propose an interpretation in the framework of the The Heterogeneous dust Evolution Model of Interstellar Solids (THEMIS). We show the presence of several spatially extended spectral features, related to aromatic and aliphatic hydrocarbon material in disks around Herbig stars, from ~ 10 to 50-100 au, and even in inner gaps devoided of large grains. The correlation and constant intensity ratios between aliphatic and aromatic CH stretching bands suggest a common nature of the carriers. Given their expected high destruction rates due to UV photons, our observations suggest that they are continuously replenished at the disk surfaces.

The dust sub-millimetre polarisation of star-forming clouds carries information on dust and the role of magnetic fields in cloud evolution. With observations of a dense filamentary cloud G035.39-00.33, we aim to characterise the dust emission properties and the variations of the polarisation fraction. JCMT SCUBA-2/POL-2 data at 850um are combined with Planck 850um (353GHz) data to map polarisation fractions. With previous SCUBA-2 observations (450um and 850um) and Herschel data, the column densities are determined via modified blackbody fits and via radiative transfer modelling. Models are constructed to examine how the polarisation angles and fractions depend on potential magnetic field geometries and grain alignment. POL-2 data show clear changes in the magnetic field orientation. The filament has a peak column density of N(H2)~7 10^22 cm-2, a minimum dust temperature of T~12 K, and a mass of some 4300Msun for the area N(H2)> 5 10^21 cm-2. The estimated average value of the dust opacity spectral index is beta ~ 1.9. The ratio of sub-millimetre and J band optical depths is tau(250 um)/tau(J) ~ 2.5 10^-3, more than four times the typical values for diffuse medium. The polarisation fraction decreases as a function of column density to p ~ 1% in the central filament. Because of noise, the observed decrease of p(N) is significant only at N(H2)>2 10^22 cm-2. The observations suggest that the grain alignment is not constant. Although the data can be explained with a complete loss of alignment at densities above ~ 10^4 cm-3 or using the predictions of radiative torques alignment, the uncertainty of the field geometry and the spatial filtering of the SCUBA-2 data prevent strong conclusions. G035.39-00.33 shows strong signs of dust evolution and the low polarisation fraction is suggestive of a loss of polarised emission from its densest parts.

Grains are determinant for setting the chemical, physical, dynamical, and radiative properties of all the media in which they are present. Their influence depends on the grain composition, size, and geometrical structure which vary throughout the dust lifecycle. Grain growth arises in dense media as traced by an enhancement of the FIR emissivity and by cloudshine and coreshine. We investigate the imprint of the grain characteristics on the dust unpolarised optical properties from the visible to the FIR for isolated grains and aggregates. Using optical constants for both carbonaceous and silicate materials, we derive the absorption and scattering efficiencies, the asymmetry factor of the phase function, the albedo, and the mass opacity, using either Mie theory or DDA. We study the effects of the size, porosity and shape of the grains, and of the monomers constituting the aggregates, on the optical properties. For aggregates we study the influence of the number of monomers and of mixing monomer sizes. Grain structure changes result in variations at all wavelengths. Porosity, grain elongation, and aggregation all produce an increase in the FIR opacity. The spectral dependence of this increase depends on the nature of the material composing the grain: it is independent of the wavelength for insulators but not for conductors. For aggregates, the FIR increase does not depend on the monomer size and saturates for aggregates containing 6 or more monomers. For silicates, the MIR spectral feature at 18 microns is more sensitive to the details of the grain structure than the 10 microns feature. This study provides a basis for understanding the range of variations achievable as a result of varying the grain characteristics. It emphasises the importance of considering the detailed grain structure in determining the dust optical properties and of using exact methods.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

The dust properties inferred from the analysis of Planck observations in total and polarized emission challenge current dust models. We propose new dust models compatible with polarized and unpolarized data in extinction and emission for translucent lines of sight ($0.5 < A_V < 2.5$). We amended the DustEM tool to model polarized extinction and emission. We fit the spectral dependence of the mean extinction, polarized extinction, SED and polarized SED with PAHs, astrosilicates and amorphous carbon (a-C). The astrosilicate population is aligned along the magnetic field lines, while the a-C population may be aligned or not. With their current optical properties, oblate astrosilicate grains are not emissive enough to reproduce the emission to extinction polarization ratio $P_{353}/p_V$ derived with Planck data. Models using prolate astrosilicate grains with an elongation $a/b=3$ and an inclusion of 20% of porosity succeed. The spectral dependence of the polarized SED is steeper in our models than in the data. Models perform slightly better when a-C grains are aligned. A small (6%) volume inclusion of a-C in the astrosilicate matrix removes the need for porosity and perfect grain alignment, and improves the fit to the polarized SED. Dust models based on astrosilicates can be reconciled with Planck data by adapting the shape of grains and adding inclusions of porosity or a-C in the astrosilicate matrix.

The purpose of this work is the characterization of the radial distribution of dust, stars, gas, and star-formation rate (SFR) in a sub-sample of 18 face-on spiral galaxies extracted from the DustPedia sample. This study is performed by exploiting the multi-wavelength, from UV to sub-mm bands, DustPedia database, in addition to molecular (12CO) and atomic (HI) gas maps and metallicity abundance information available in the literature. We fitted the surface brightness profiles of the tracers of dust and stars, the mass surface density profiles of dust, stars, molecular gas, and total gas, and the SFR surface density profiles with an exponential curve and derived their scale-lengths. We also developed a method to solve for the CO-to-H2 conversion factor (alpha_CO) per galaxy by using dust and gas mass profiles. Although each galaxy has its own peculiar behaviour, we identified a common trend of the exponential scale-lengths vs. wavelength. On average, the scale-lengths normalized to the B-band 25 mag/arcsec^2 radius decrease from UV to 70 micron, from 0.4 to 0.2, and then increase back up to 0.3 at 500 microns. The main result is that, on average, the dust mass surface density scale-length is about 1.8 times the stellar one derived from IRAC data and the 3.6 micron surface brightness, and close to that in the UV. We found a mild dependence of the scale-lengths on the Hubble stage T: the scale-lengths of the Herschel bands and the 3.6 micron scale-length tend to increase from earlier to later types, the scale-length at 70 micron tends to be smaller than that at longer sub-mm wavelength with ratios between longer sub-mm wavelengths and 70 micron that decrease with increasing T. The scale-length ratio of SFR and stars shows a weak increasing trend towards later types.

Here we introduce the interstellar dust modelling framework THEMIS (The Heterogeneous dust Evolution Model for Interstellar Solids), which takes a global view of dust and its evolution in response to the local conditions in interstellar media. This approach is built upon a core model that was developed to explain the dust extinction and emission in the diffuse interstellar medium. The model was then further developed to self-consistently include the effects of dust evolution in the transition to denser regions. The THEMIS approach is under continuous development and currently we are extending the framework to explore the implications of dust evolution in HII regions and the photon-dominated regions associated with star formation. We provide links to the THEMIS, DustEM and DustPedia websites where more information about the model, its input data and applications can be found.

Filaments are key for star formation models. As part of the study carried out by the Herschel GCC Programme, here we study the filament properties presented in GCC.VII in context with theoretical models of filament formation and evolution. A conservative sample of filaments at a distance D<500pc was extracted with the Getfilaments algorithm. Their physical structure was quantified according to two main components: the central (Gaussian) region (core component), and the power-law like region dominating the filament column density profile at larger radii (wing component). The properties and behaviour of these components relative to the total linear mass density of the filament and its environmental column density were compared with theoretical models describing the evolution of filaments under gravity-dominated conditions. The feasibility of a transition to supercritical state by accretion is dependent on the combined effect of filament intrinsic properties and environmental conditions. Reasonably self-gravitating (high Mline-core) filaments in dense environments (av\sim3mag) can become supercritical in timescales of t\sim1Myr by accreting mass at constant or decreasing width. The trend of increasing Mline-tot (Mline-core and Mline-wing), and ridge Av with background also indicates that the precursors of star-forming filaments evolve coevally with their environment. The simultaneous increase of environment and filament Av explains the association between dense environments and high Mline-core values, and argues against filaments remaining in constant single-pressure equilibrium states. The simultaneous growth of filament and background in locations with efficient mass assembly, predicted in numerical models of collapsing clouds, presents a suitable scenario for the fulfillment of the combined filament mass-environment criterium that is in quantitative agreement with Herschel observations.

Dust modeling is crucial to infer dust properties and budget for galaxy studies. However, there are systematic disparities between dust grain models that result in corresponding systematic differences in the inferred dust properties of galaxies. Quantifying these systematics requires a consistent fitting analysis. We compare the output dust parameters and assess the differences between two dust grain models, Compiegne et al (2011), and THEMIS (Jones et al 2013, Kohler et al 2015). In this study, we use a single fitting method applied to all the models to extract a coherent and unique statistical analysis. We fit the models to the dust emission seen by Spitzer and Herschel in the Small and Large Magellanic Clouds (SMC and LMC). The observations cover the infrared (IR) spectrum from a few microns to the sub-millimeter range. For each fitted pixel, we calculate the full n-D likelihood, based on the method described in Gordon et al (2014). The free parameters are both environmental ($U$, the interstellar radiation field strength; $\alpha_\mathrm{ISRF}$, power-law coefficient for a multi-U environment; $\Omega^*$, the starlight strength) and intrinsic to the model ($Y_\mathrm{i}$: abundances of the grain species $i$; $\alpha_\mathrm{sCM20}$, coefficient in the small carbon grain size distribution). Fractional residuals of 5 different sets of parameters show that fitting THEMIS brings a more accurate reproduction of the observations than the Compiegne model. However, independent variations of the dust species show strong model-dependencies. We find that the abundance of silicates can only be constrained to an upper-limit and the silicate/carbon ratio is different than that seen in our Galaxy. In the LMC, our fits result in dust masses slightly lower than those found in literature , by a factor lower than 2. In the SMC, we find dust masses in agreement with previous studies.

The European Space Agency has invested heavily in two cornerstones missions; Herschel and Planck. The legacy data from these missions provides us with an unprecedented opportunity to study cosmic dust in galaxies so that we can answer fundamental questions about, for example: the origin of the chemical elements, physical processes in the interstellar medium (ISM), its effect on stellar radiation, its relation to star formation and how this relates to the cosmic far infrared background. In this paper we describe the DustPedia project, which is enabling us to develop tools and computer models that will help us relate observed cosmic dust emission to its physical properties (chemical composition, size distribution, temperature), to its origins (evolved stars, super novae, growth in the ISM) and the processes that destroy it (high energy collisions and shock heated gas). To carry out this research we will combine the Herschel/Planck data with that from other sources of data, providing observations at numerous wavelengths (< 41) across the spectral energy distribution, thus creating the DustPedia database. To maximise our spatial resolution and sensitivity to cosmic dust we limit our analysis to 4231 local galaxies (v < 3000 km/s) selected via their near infrared luminosity (stellar mass). To help us interpret the data we have developed a new physical model for dust (THEMIS), a new Bayesian method of fitting and interpreting spectral energy distributions (HerBIE) and a state-of-the-art Monte Carlo photon tracing radiative transfer model (SKIRT). In this the first of the DustPedia papers we describe the project objectives, data sets used and provide an insight into the new scientific methods we plan to implement.

We study the correlation between far-infared/submm dust emission and atomic gas column density in order to derive the properties of the high Galactic latitude, low density, Milky Way cirrus in the foreground of the Virgo cluster of galaxies. Dust emission maps from 60 to 850 um are obtained from SPIRE observations carried out within the Herschel Virgo Cluster Survey, complemented by IRAS-IRIS and Planck-HFI maps. Data from the Arecibo legacy Fast ALFA Survey is used to derive atomic gas column densities for two broad velocity components, low and intermediate velocity clouds. Dust emissivities are derived for each gas component and each far-infared/submm band. For the low velocity clouds, we measure an average emissivity 0.79 +/- 0.08 times 1E-20 MJy sr^-1 cm^2 at 250 um. After fitting a modified blackbody to the available bands, we estimated a dust absorption cross-section 0.49 +/- 0.13 times 1E-25 cm^2 H^-1 at 250 um (with dust temperature T = 20.4 +/- 1.5 K and spectral index beta = 1.53 +/- 0.17). The results are in excellent agreement with those obtained by Planck over a much larger coverage of the high Galactic latitude cirrus (50% of the sky vs 0.2% in our work). For dust associated with intermediate velocity gas, we confirm earlier Planck results and find a higher temperature and lower emissivity and cross-section. After subtracting the modelled components, we find regions at scales smaller than 20' where the residuals deviate significantly from the average, cosmic-infrared-background dominated, scatter. These large residuals are most likely due to local variations in the cirrus dust properties (and/or the dust/atomic-gas correlation) or to high-latitude molecular clouds with average N_H2 <~ 1E20 cm^-2. We find no conclusive evidence for intracluster dust emission in Virgo.

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.

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.

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.

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.

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 {\mu}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.

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.

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\'enic or sub-Alfv\'enic 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.

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.

We compare the performance of several dust models in reproducing the dust spectral energy distribution (SED) per unit extinction in the diffuse interstellar medium (ISM). We use our results to constrain the variability of the optical properties of big grains in the diffuse ISM, as published by the Planck collaboration. We use two different techniques to compare the predictions of dust models to data from the Planck HFI, IRAS and SDSS surveys. First, we fit the far-infrared emission spectrum to recover the dust extinction and the intensity of the interstellar radiation field (ISRF). Second, we infer the ISRF intensity from the total power emitted by dust per unit extinction, and then predict the emission spectrum. In both cases, we test the ability of the models to reproduce dust emission and extinction at the same time. We identify two issues. Not all models can reproduce the average dust emission per unit extinction: there are differences of up to a factor $\sim2$ between models, and the best accord between model and observation is obtained with the more emissive grains derived from recent laboratory data on silicates and amorphous carbons. All models fail to reproduce the variations in the emission per unit extinction if the only variable parameter is the ISRF intensity: this confirms that the optical properties of dust are indeed variable in the diffuse ISM. Diffuse ISM observations are consistent with a scenario where both ISRF intensity and dust optical properties vary. The ratio of the far-infrared opacity to the $V$ band extinction cross-section presents variations of the order of $\sim20\%$ ($40-50\%$ in extreme cases), while ISRF intensity varies by $\sim30\%$ ($\sim60\%$ in extreme cases). This must be accounted for in future modelling.

Variations in the observed dust emission and extinction indicate a systematic evolution of grain properties in the transition from the diffuse interstellar medium (ISM) to denser molecular clouds. The differences in the dust spectral energy distribution (SED) observed from the diffuse ISM to denser regions, namely an increase in the spectral index at long wavelengths, an increase in the FIR opacity, and a decrease in temperature, are usually assumed to be the result of changes in dust properties. We investigate if evolutionary processes, such as coagulation and accretion, are able to change the dust properties of grains in a way that is consistent with observations. We use a core-mantle grain model to describe diffuse ISM-type grains, and using DDA we calculate how the accretion of mantles and coagulation into aggregates vary the grain optical properties. We calculate the dust SED and extinction using DustEM and the radiative transfer code CRT. We show that the accretion of an aliphatic carbon mantle on diffuse ISM-type dust leads to an increase in the FIR opacity by a factor of about 2 and in the FIR/submm spectral index from 1.5 to 1.8, and to a decrease in the temperature by about 2 K. We also show that the coagulation of these grains into aggregates further decreases the temperature by 3 K and increases the spectral index up to a value of $\sim$2. The FIR opacity is increased by a factor of 3 (7) for these aggregates (with an additional ice-mantle) compared to the diffuse ISM-dust. Dust evolution in the ISM resulting from coagulation and accretion, leads to significant changes in the optical properties of the grains that can explain the observed variations in the dust SED in the transition from the diffuse ISM to denser regions.

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.

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.

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.

In this chapter, we will outline the scientific motivation for studying Anomalous Microwave Emission (AME) with the SKA. AME is thought to be due to electric dipole radiation from small spinning dust grains, although thermal fluctuations of magnetic dust grains may also contribute. Studies of this mysterious component would shed light on the emission mechanism, which then opens up a new window onto the interstellar medium (ISM). AME is emitted mostly in the frequency range $\sim 10$--100\,GHz, and thus the SKA has the potential of measuring the low frequency side of the AME spectrum, particularly in band 5. Science targets include dense molecular clouds in the Milky Way, as well as extragalactic sources. We also discuss the possibility of detecting rotational line emission from Poly-cyclic Aromatic Hydrocarbons (PAHs), which could be the main carriers of AME. Detecting PAH lines of a given spacing would allow for a definitive identification of specific PAH species.

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.

Observational evidence seems to indicate that the depletion of interstellar carbon into dust shows rather wide variations and that carbon undergoes rather rapid recycling in the interstellar medium (ISM). Small hydrocarbon grains are processed in photo-dissociation regions by UV photons, by ion and electron collisions in interstellar shock waves and by cosmic rays. A significant fraction of hydrocarbon dust must therefore be re-formed by accretion in the dense, molecular ISM. A new dust model (Jones et al., Astron. Astrophys., 2013, 558, A62) shows that variations in the dust observables in the diffuse interstellar medium (nH = 1000 cm^3), can be explained by systematic and environmentally-driven changes in the small hydrocarbon grain population. Here we explore the consequences of gas-phase carbon accretion onto the surfaces of grains in the transition regions between the diffuse ISM and molecular clouds (e.g., Jones, Astron. Astrophys., 2013, 555, A39). We find that significant carbonaceous dust re-processing and/or mantle accretion can occur in the outer regions of molecular clouds and that this dust will have significantly different optical properties from the dust in the adjacent diffuse ISM. We conclude that the (re-)processing and cycling of carbon into and out of dust is perhaps the key to advancing our understanding of dust evolution in the ISM.

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.

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.

Anomalous microwave emission (AME) is believed to be due to electric dipole radiation from small spinning dust grains. The aim of this paper is a statistical study of the basic properties of AME regions and the environment in which they emit. We used WMAP and Planck maps, combined with ancillary radio and IR data, to construct a sample of 98 candidate AME sources, assembling SEDs for each source using aperture photometry on 1deg-smoothed maps from 0.408 GHz up to 3000 GHz. Each spectrum is fitted with a simple model of free-free, synchrotron (where necessary), cosmic microwave background (CMB), thermal dust, and spinning dust components. We find that 42 of the 98 sources have significant (>5sigma) excess emission at frequencies between 20 and 60 GHz. An analysis of the potential contribution of optically thick free-free emission from ultra-compact HII regions, using IR colour criteria, reduces the significant AME sample to 27 regions. The spectrum of the AME is consistent with model spectra of spinning dust. The AME regions tend to be more spatially extended than regions with little or no AME. The AME intensity is strongly correlated with the submillimetre/IR flux densities and comparable to previous AME detections in the literature. AME emissivity, defined as the ratio of AME to dust optical depth, varies by an order of magnitude for the AME regions. The AME regions tend to be associated with cooler dust in the range 14-20 K and an average emissivity index of +1.8, while the non-AME regions are typically warmer, at 20-27 K. In agreement with previous studies, the AME emissivity appears to decrease with increasing column density. The emerging picture is that the bulk of the AME is coming from the polycyclic aromatic hydrocarbons and small dust grains from the colder neutral interstellar medium phase (Abridged).

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 $\mu$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.

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.

We observed the L1506 filament, which is located in the Taurus molecular complex, with the Herschel PACS and SPIRE instruments. Our aim is to prove the variation in grain properties along the entire length of the filament. In particular, we want to determine above which gas density this variation arises and what changes in the grain optical properties/size distribution are required. We use the 3D radiative transfer code CRT, coupled to the dust emission and extinction code DustEM, to model the emission and extinction of the dense filament. We test a range of optical properties and size distributions for the grains: dust of the diffuse interstellar medium (interstellar PAHs and amorphous carbons and silicates) and both compact and fluffy aggregates. We find that the grain opacity has to increase across the filament to fit simultaneously the near-IR extinction and Herschel emission profiles of L1506. We interpret this change to be a consequence of the coagulation of dust grains to form fluffy aggregates. Grains similar to those in the diffuse medium have to be present in the outer layers of the cloud, whereas aggregates must prevail above gas densities of a few 1000 H/cm3. This corresponds to line-of-sights with visual extinction in the V band of the order of 2 to 3. The dust opacity at 250 microns is raised by a factor of 1.8 to 2.2, while the grain average size is increased by a factor of 5. These exact numbers depend naturally on the dust model chosen to fit the data. Our findings agree with the constraints given by the study of the gas molecular lines. Using a simple approach, we show that the aggregates may have time to form inside the filament within the cloud lifetime. Our model also characterises the density structure of the filament, showing that the filament width is not constant along L1506 but instead varies by a factor of the order of 4.

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 (\nu < 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$\alpha$ 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.

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.

We use the most up-to-date cosmological evolution models of star-forming (SF) galaxies and radio sources to compute the extragalactic number counts and the cosmic background from 408MHz to 12THz. The model of SF galaxies reproduces the constraints obtained by Spitzer, Herschel, and ground-based submm/mm experiments: number counts, redshift distribution of galaxies, cosmic background intensity and anisotropies. The template SEDs of this model are extrapolated to the radio adding synchrotron, free-free, and spinning dust emissions. To fix the synchrotron intensity, we use the IR/radio flux ratio, q70, and a spectral index beta=-3. For a constant q70, our model added to the AGN contribution provides a good fit to the number counts from 12THz to 408MHz and to the CIB. Spinning dust accounts for up to 20% of the cosmic microwave background produced by SF galaxies, but for less than 10% of the total background when AGN are included. The SF galaxies account for 77.5% of the number counts at 1.4GHz for a flux of 1e-4Jy. However, the model does not explain the CRB measured with the ARCADE2 experiment. Considering the case when q70 decreases strongly with redshift, this still does not explain the ARCADE2 results. It also yields to an overestimate of the low-flux number counts in the radio. Thus, we rule out a steep variation of q70 with redshift at least for z<3.5. Adding a population of faint SF galaxies at high redshift (Lir<1e11Lsun and 4<z<6), which would reproduce the ARCADE2 results, leads to predictions of the CIB much higher than what is observed, ruling out this as the explanation for the ARCADE2 results. Considering our findings and previous studies, we conclude that if the radio emission measured by ARCADE2 is astrophysical in origin, it has to originate in the Galaxy or in a new kind of radio sources (with no mid- to far-IR counterparts) or emission mechanism still to be discovered.

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.

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.

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.

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.

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.

(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\mu 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.

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.

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.

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.

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.

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.

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 {\beta} 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 {\beta} over the all temperature range. {\beta} ranges from 1.4 to 2.8 with a mean value around 2.1, and several possible scenarios are possible, including {\beta}(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.

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

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.

Dust emission at submm to cm wavelengths is often simply the Rayleigh-Jeans tail of dust particles at thermal equilibrium and is used as a cold mass tracer in various environments including nearby galaxies. However, well-sampled spectral energy distributions of the nearby, star-forming Magellanic Clouds have a pronounced (sub-)millimeter excess (Israel et al., 2010). This study attempts to confirm the existence of such a millimeter excess above expected dust, free-free and synchrotron emission and to explore different possibilities for its origin. We model NIR to radio spectral energy distributions of the Magellanic Clouds with dust, free-free and synchrotron emission. A millimeter excess emission is confirmed above these components and its spectral shape and intensity are analysed in light of different scenarios: very cold dust, Cosmic Microwave Background (CMB) fluctuations, a change of the dust spectral index and spinning dust emission. We show that very cold dust or CMB fluctuations are very unlikely explanations for the observed excess in these two galaxies. The excess in the LMC can be satisfactorily explained either by a change of the spectral index due to intrinsic properties of amorphous grains, or by spinning dust emission. In the SMC however, due to the importance of the excess, the dust grain model including TLS/DCD effects cannot reproduce the observed emission in a simple way. A possible solution was achieved with spinning dust emission, but many assumptions on the physical state of the interstellar medium had to be made. Further studies, using higher resolution data from Planck and Herschel, are needed to probe the origin of this observed submm-cm excess more definitely. Our study shows that the different possible origins will be best distinguished where the excess is the highest, as is the case in the SMC.

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.

The emission of cold dust grains at long wavelengths will soon be observed by the Planck and Herschel satellites and provide new constraints on the nature of interstellar dust. The microwave anomalous emission, proposed to be due to spinning PAHs, should help to better define these species. Moreover, understanding the fluctuations of the anomalous emission over the sky is crucial for CMB studies. We focus on the long wavelength emission of interstellar PAHs in their rovibrational and rotational transitions. The PAH emission spectrum from the IR to the microwave range is presented and compared to anomalous emission observations. To model their long wavelength emission, we treat PAHs as isolated systems and follow consistently their IR and rotational emissions. We consider several interstellar phases and discuss how the anomalous emission may constrain their size distribution. Our model of PAH emission accounts for the mid-IR spectra of the diffuse interstellar medium and of the Orion Bar. For lambda<3mm the PAH IR emission does not scale any more with the radiation field intensity (Go) unlike the mid-IR part of the spectrum. This emission represents less than 10% of the total dust emission at 100 GHz. We find the broadband emissivity of spinning PAHs per carbon atom to be rather constant for Go<100 and nH<100cm-3. Observations of anomalous emission in the Perseus molecular cloud are explained by our model of spinning PAH emission with a standard abundance, supporting PAHs as the origin of the anomalous emission. The behaviour of the spinning dust emissivity with Go provides a clear test that can be tested against observations of anomalous and dust mid-IR emissions. Comparison of these emissions will provide constraints on the size and/or electric dipole moment of interstellar PAHs.

The galactic anomalous microwave emission detected between 10 and 90 GHz is a major foreground to CMB fluctuations. Well correlated to dust emission at 100 $\mu$m, the anomalous emission is interstellar but its origin is still debated. Some possible explanations relate it to dust: emission of spinning, small (nanometric) grains carrying a permanent electric dipole or magnetic fluctuations in larger (submicronic) grains. To probe the origin of the anomalous emission, we compare microwave data to dust IR emission and search for specific signatures predicted by models of spinning dust. For the anomalous emission, we use the 23 GHz all-sky map deduced from WMAP data by Miville-Deschenes et al. (2008). The dust emission is traced by IRAS data. Models show that spinning dust emission is little sensitive to the intensity of the radiation field (Go) for 10<nu<30 GHz while the corresponding mid-IR emission is proportional to Go. To test this behaviour in our comparison, we derive Go from the dust temperature maps of Schlegel et al. (1998). From all-sky maps, we show that the anomalous emission is better correlated to the emission of small grains (at 12 $\mu$m) than to that of big grains (at 100 $\mu$m). In addition we show that the former correlation is significantly improved when the 12 $\mu$m flux is divided by Go, as predicted by current models of spinning dust.

The main goals of this study is to use the information from both WMAP intensity and polarization data to do a separation of the Galactic components, with a focus on the synchrotron and anomalous emissions. Our analysis is made at 23 GHz where the signal-to-noise ratio is the highest and the estimate of the CMB map is less critical. Our estimate of the synchrotron intensity is based on an extrapolation of the Haslam 408 MHz data with a spatially varying spectral index constrained by the WMAP 23 GHz polarization data and a bi-symmetrical spiral model of the galactic magnetic field with a turbulent part following a -5/3 power law spectrum. The 23 GHz polarization data are found to be compatible with a magnetic field with a pitch angle p=-8.5 degrees and an amplitude of the turbulent part of the magnetic field 0.57 times the local value of the field, in agreement with what is found using rotation measures of pulsars and polarized extinction by dust. The synchrotron spectral index between 408 MHz and 23 GHz obtained from polarization data and our model of the magnetic field has a mean value of Beta=-3.00 with a limited spatial variation with a standard deviation of 0.06. When thermal dust, free-free and synchrotron are removed from the WMAP intensity data, the residual anomalous emission is highly correlated with thermal dust emission with a spectrum in agreement with spinning dust models. Considering a classical model of the large scale Galactic magnetic field, we show that the polarization data of WMAP are in favor of a soft synchrotron intensity highly correlated with the 408 MHz data. Furthermore the combination of the WMAP polarization and intensity data brings strong evidence for the presence of unpolarized spinning dust emission in the 20-60 GHz range.