Search

The Planck mission provided all-sky dust emission maps in the submm to mm range at an angular resolution of 5'. In addition, some specific sources can be observed at long wavelengths and higher resolution using ground-based telescopes. These observations are limited to small scales and require extensive data processing before they become available for scientific analysis. They also suffer from extended emission filtering. At present, we are still unable to fully understand the emissivity variations observed in different astrophysical environments at long wavelengths. It is therefore challenging to estimate any dust emission in the submm-mm at a better resolution than the 5' from Planck. In this analysis, based on supervised deep learning algorithms, we produced dust emission predictions in the two Planck bands centered at 850 mic and 1.38 mm, at the Herschel resolution (37''). Herschel data of Galactic environments, ranging from 160 to 500 mic and smoothed to 5', were used to train the neural network. Then, using Herschel data only, the model was applied to predict dust emission maps at 37''. The neural network is capable of reproducing dust emission maps of various Galactic environments. Remarkably, it also performs well for nearby extragalactic environments. This could indicate that large dust grains have similar properties in both our Galaxy and nearby galaxies, or at least that their spectral behaviors are comparable in Galactic and extragalactic environments. We provide dust emission prediction maps at 850 mic and 1.38 mm at the 37'' of several surveys: Hi-GAL, Gould Belt, Cold Cores, HERITAGE, Helga, HerM33es, KINGFISH, and VNGS. The ratio of these two wavelength brightness bands reveals a derived emissivity spectral index statistically close to 1 for all the surveys, which favors the hypothesis of a flattened dust emission spectrum for wavelengths larger than 850 mic., Comment: Accepted for publication in A&A (18 pages)

The millimeter continuum emission from galaxies provides important information about cold dust, its distribution, heating, and role in their InterStellar Medium (ISM). This emission also carries an unknown portion of the free-free and synchrotron radiation. The IRAM 30m Guaranteed Time Large Project, Interpreting Millimeter Emission of Galaxies with IRAM and NIKA2 (IMEGIN) provides a unique opportunity to study the origin of the millimeter emission on angular resolutions of <18" in a sample of nearby galaxies. As a pilot study, we present millimeter observations of two IMEGIN galaxies, NGC 2146 (starburst) and NGC 2976 (peculiar dwarf) at 1.15 mm and 2 mm. Combined with the data taken with Spitzer, Herschel, Plank, WSRT, and the 100m Effelsberg telescopes, we model the infrared-to-radio Spectral Energy Distribution (SED) of these galaxies, both globally and at resolved scales, using a Bayesian approach to 1) dissect different components of the millimeter emission, 2) investigate the physical properties of dust, and 3) explore correlations between millimeter emission, gas, and Star Formation Rate (SFR). We find that cold dust is responsible for most of the 1.15 mm emission in both galaxies and at 2 mm in NGC 2976. The free-free emission emits more importantly in NGC 2146 at 2 mm. The cold dust emissivity index is flatter in the dwarf galaxy ($\beta = 1.3\pm 0.1$) compared to the starburst galaxy ($\beta = 1.7\pm 0.1$). Mapping the dust-to-gas ratio, we find that it changes between 0.004 and 0.01 with a mean of $0.006\pm0.001$ in the dwarf galaxy. In addition, no global balance holds between the formation and dissociation of H$_2$ in this galaxy. We find tight correlations between the millimeter emission and both the SFR and molecular gas mass in both galaxies., Comment: Accepted for publication in Astronomy & Astrophysics

The JWST has captured the sharpest IR images ever taken of the Horsehead nebula, a prototypical moderately irradiated PDR that is fully representative of most of the UV-illuminated molecular gas in the Milky Way and star-forming galaxies. We investigate the impact of FUV radiation of a molecular cloud and constrain the structure of the edge of the PDR and its illumination conditions. We used NIRCam and MIRI to obtain 17 broadband and 6 narrowband maps from 0.7 to 28 $\mu$m. We mapped the dust emission, scattered light, and several gas phase lines. We also used HST-WFC3 maps at 1.1 and 1. 6 $\mu$m, along with HST-STIS spectroscopic observations of the H$\alpha$ line. We probed the structure of the edge of the Horsehead and resolved its spatial complexity. We detected a network of faint striated features extending perpendicularly to the PDR front into the H\,II region in filters sensitive to nano-grain emission and light scattered by larger grains. This may indeed figure as the first detection of the entrainment of dust particles in the evaporative flow. The map of the 1-0 S(1) line of H$_2$ presents sharp sub-structures on scales as small as 1.5 arcsec. The ionization and dissociation fronts appear at distances 1-2 arcsec behind the edge of the PDR and seem to spatially coincide, indicating a thickness of the neutral atomic layer below 100 au. All broadband maps present strong color variations which can be explained by dust attenuation. Deviations of the emissions in the H$\alpha$, Pa$\alpha,$ and Br$\alpha$ lines also indicate dust attenuation. With a very simple model, we derive the main features of the extinction curve. A small excess of extinction at 3 $\mu$m may be attributed to icy H$_2$O mantles onto grains. In all lines of sight crossing the inner regions of the Horsehead, it appears that dust attenuation is non-negligible over the entire spectral range., Comment: 25 pages, 17 figures

We study the emission of dust grains within the Orion Bar - a well-known, highly far-UV (FUV)-irradiated PDR. The Orion Bar because of its edge-on geometry provides an exceptional benchmark for characterizing dust evolution and the associated driving processes under varying physical conditions. Our goal is to constrain the local properties of dust by comparing its emission to models. Taking advantage of the recent JWST PDRs4All data, we follow the dust emission as traced by JWST NIRCam (at 3.35 and 4.8 micron) and MIRI (at 7.7, 11.3, 15.0, and 25.5 micron), along with NIRSpec and MRS spectroscopic observations. First, we constrain the minimum size and hydrogen content of carbon nano-grains from a comparison between the observed dust emission spectra and the predictions of the THEMIS dust model coupled to the numerical code DustEM. Using this dust model, we then perform 3D radiative transfer simulations of dust emission with the SOC code and compare to data obtained along well chosen profiles across the Orion Bar. The JWST data allows us, for the first time, to spatially resolve the steep variation of dust emission at the illuminated edge of the Orion Bar PDR. By considering a dust model with carbonaceous nano-grains and submicronic coated silicate grains, we derive unprecedented constraints on the properties of across the Orion Bar. To explain the observed emission profiles with our simulations, we find that the nano-grains must be strongly depleted with an abundance (relative to the gas) 15 times less than in the diffuse ISM. The NIRSpec and MRS spectroscopic observations reveal variations in the hydrogenation of the carbon nano-grains. The lowest hydrogenation levels are found in the vicinity of the illuminating stars suggesting photo-processing while more hydrogenated nano-grains are found in the cold and dense molecular region, potentially indicative of larger grains., Comment: 18 pages, 11 figures, accepted for publication in A&A

In the framework of the IMEGIN Large Program, we used the NIKA2 camera on the IRAM 30-m telescope to observe the edge-on galaxy NGC 891 at 1.15 mm and 2 mm and at a FWHM of 11.1" and 17.6", respectively. Multiwavelength data enriched with the new NIKA2 observations fitted by the HerBIE SED code (coupled with the THEMIS dust model) were used to constrain the physical properties of the ISM. Emission originating from the diffuse dust disk is detected at all wavelengths from mid-IR to mm, while mid-IR observations reveal warm dust emission from compact HII regions. Indications of mm excess emission have also been found in the outer parts of the galactic disk. Furthermore, our SED fitting analysis constrained the mass fraction of the small (< 15 Angstrom) dust grains. We found that small grains constitute 9.5% of the total dust mass in the galactic plane, but this fraction increases up to ~ 20% at large distances (|z| > 3 kpc) from the galactic plane., Comment: To appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

As the millimeter wavelength range remains a largely unexplored spectral region for galaxies, the IMEGIN large program aims to map the millimeter continuum emission of 22 nearby galaxies at 1.15 and 2 mm. Using the high-resolution maps produced by the NIKA2 camera, we explore the existence of very cold dust and take possible contamination by free-free and synchrotron emission into account. We study the IR-to-radio emission coming from different regions along the galactic plane and at large vertical distances. New observations of NGC 891, using the NIKA2 camera on the IRAM 30m telescope, along with a suite of observations at other wavelengths were used to perform a multiwavelength study of the spectral energy distribution in the interstellar medium in this galaxy. This analysis was performed globally and locally, using the advanced hierarchical Bayesian fitting code, HerBIE, coupled with the THEMIS dust model. Our dust modeling is able to reproduce the near-IR to millimeter emission of NGC 891, with the exception of an excess at a level of 25% obtained by the NIKA2 observations in the outermost parts of the disk. The radio continuum and thermal dust emission are distributed differently in the disk and galaxy halo. Different dusty environments are also revealed by a multiwavelength investigation of the emission features. Our detailed decomposition at millimeter and centimeter wavelengths shows that emission at 1 mm is purely originated by dust. Radio components become progressively important with increasing wavelengths. Finally, we find that emission arising from small dust grains accounts for ~ 9.5% of the total dust mass, reaching up to 20% at large galactic latitudes. Shock waves in the outflows that shatter the dust grains might explain this higher fraction of small grains in the halo., Comment: 15 pages, 10 figures, Accepted for publication in A&A

Cosmic dust models are key ingredients in advancing our understanding of astronomical environments as diverse as interstellar clouds in galaxies, circumstellar envelopes around evolved and young stars, and protoplanetary disks. Such models consist of several dust populations, each with different compositions and size distributions. They may also consider different grain shapes, although most models assume spherical grains. All include a component of silicate dust. The absorption and emission properties of these dust components are calculated from the optical constants of each dust material which have various experimental, phenomenological, and theoretical origins depending on the model. We aim to provide the community with new sets of optical constants for amorphous silicate dust analogues at low temperatures. The analogues consist of four Mg-rich silicate samples of stoichiometry ranging from enstatite to olivine, and of eight samples of Mg and Fe rich silicates with a pyroxene stoichiometry and differing magnesium and iron content. We calculated the optical constants from transmission measurements using the Kramers-Kronig relations, assuming that the grains are small compared to the wavelength and prolate in shape with axis ratios of 1.5 and 2 for the Mg and Fe rich samples, respectively. New optical constants for silicate dust analogues were calculated over the wavelength range from 5 to 800-1000 microns, depending on the sample, and at temperatures of 10, 30, 100, 200, and 300 K.We determined the uncertainties on the derived optical constants based on the assumptions used to calculate them. To facilitate the use of these data in cosmic dust models, we provide optical constants extrapolated outside the measured spectral range into the UV-NIR and mm-cm wavelength ranges, as well as formulae that can be used to interpolate them at any temperature in the range 10 - 300 K.

Context. Carbonaceous nano-grains play a fundamental role in the physico-chemistry of the interstellar medium (ISM) and especially of photon-dominated regions (PDRs). Their properties vary with the local physical conditions and affect the local chemistry and dynamics. Aims. We aim to highlight the evolution of carbonaceous nano-grains in three different PDRs and propose a scenario of dust evolution as a response to the physical conditions. Methods. We used Spitzer/IRAC (3.6, 4.5, 5.8, and 8 $\mu$m) and Spitzer/MIPS (24 $\mu$m) together with Herschel/PACS (70 $\mu$m) to map dust emission in IC63 and the Orion Bar. To assess the dust properties, we modelled the dust emission in these regions using the radiative transfer code SOC together with the THEMIS dust model. Results. Regardless of the PDR, we find that nano-grains are depleted and that their minimum size is larger than in the diffuse ISM (DISM), which suggests that the mechanisms that lead nano-grains to be photo-destroyed are very efficient below a given critical size limit. The evolution of the nano-grain dust-to-gas mass ratio with both G0 and the effective temperature of the illuminating star indicates a competition between the nano-grain formation through the fragmentation of larger grains and nano-grain photo-destruction. We modelled dust collisions driven by radiative pressure with a classical 1D approach to show that this is a viable scenario for explaining nano-grain formation through fragmentation and, thus, the variations observed in nano-grain dust-to-gas mass ratios from one PDR to another. Conclusions. We find a broad variation in the nano-grain dust properties from one PDR to another, along with a general trend of nano-grain depletion in these regions. We propose a viable scenario of nano-grain formation through fragmentation of large grains due to radiative pressure-induced collisions.

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

The Polarized Instrument for Long-wavelength Observation of the Tenuous interstellar medium (PILOT) is a balloon-borne experiment that aims to measure the polarized emission of thermal dust at a wavelength of 240 um (1.2 THz). The PILOT experiment flew from Timmins, Ontario, Canada in 2015 and 2019 and from Alice Springs, Australia in April 2017. The in-flight performance of the instrument during the second flight was described in Mangilli et al. 2019. In this paper, we present data processing steps that were not presented in Mangilli et al. 2019 and that we have recently implemented to correct for several remaining instrumental effects. The additional data processing concerns corrections related to detector cross-talk and readout circuit memory effects, and leakage from total intensity to polarization. We illustrate the above effects and the performance of our corrections using data obtained during the third flight of PILOT, but the methods used to assess the impact of these effects on the final science-ready data, and our strategies for correcting them will be applied to all PILOT data. We show that the above corrections, and in particular that for the intensity to polarization leakage, which is most critical for accurate polarization measurements with PILOT, are accurate to better than 0.4 % as measured on Jupiter during flight#3.

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

The mm-to-cm range of the Spectral Energy Distribution of spiral galaxies remains largely unexplored. Its coverage is required to disentangle the contribution of dust emission, free-free and synchrotron radiation and can provide constraints on dust models, star-formation rates and ISM properties. We present the case for a synergy between NIKA2 observations of nearby spirals and those from planned and current instrumentation at the Sardinia Radio Telescope, and report on a pilot K-band program to search for Anomalous Microwave Emission, an elusive emission component which is presumably related to dust., Comment: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

In this era of spatially resolved observations of planet forming disks with ALMA and large ground-based telescopes such as the VLT, Keck and Subaru, we still lack statistically relevant information on the quantity and composition of the material that is building the planets, such as the total disk gas mass, the ice content of dust, and the state of water in planetesimals. SPICA is an infrared space mission concept developed jointly by JAXA and ESA to address these questions. The key unique capabilities of SPICA that enable this research are (1) the wide spectral coverage 10-220 micron, (2) the high line detection sensitivity of (1-2) 10-19 W m-2 with R~2000-5000 in the far-IR (SAFARI) and 10-20 W m-2 with R~29000 in the mid-IR (SMI, spectrally resolving line profiles), (3) the high far-IR continuum sensitivity of 0.45 mJy (SAFARI), and (4) the observing efficiency for point source surveys. This paper details how mid- to far-IR infrared spectra will be unique in measuring the gas masses and water/ice content of disks and how these quantities evolve during the planet forming period. These observations will clarify the crucial transition when disks exhaust their primordial gas and further planet formation requires secondary gas produced from planetesimals. The high spectral resolution mid-IR is also unique for determining the location of the snowline dividing the rocky and icy mass reservoirs within the disk and how the divide evolves during the build-up of planetary systems. Infrared spectroscopy (mid- to far-IR) of key solid state bands is crucial for assessing whether extensive radial mixing, which is part of our Solar System history, is a general process occurring in most planetary systems and whether extrasolar planetesimals are similar to our Solar System comets/asteroids. ... (abbreviated), Comment: accepted for publication in PASA

Carbon is one of the most abundant components in the Universe. While silicates have been the main focus of solid phase studies in protoplanetary discs (PPDs), little is known about the solid carbon content especially in the planet-forming regions ($\sim $0.1 to 10 au). Fortunately, several refractory carbonaceous species present C-H bonds (such as hydrogenated nano-diamond and amorphous carbon as well as polycyclic aromatic hydrocarbons (PAHs)), which generate infrared (IR) features that can be used to trace the solid carbon reservoirs. The new mid-IR instrument MATISSE, installed at the Very Large Telescope Interferometer (VLTI), can spatially resolve the inner regions ($\sim$ 1 to 10 au) of PPDs and locate, down to the au-scale, the emission coming from carbon grains. Our aim is to provide a consistent view on the radial structure, down to the au-scale, as well as basic physical properties and the nature of the material responsible for the IR continuum emission in the inner disk region around HD 179218. We implemented a temperature-gradient model to interpret the disk IR continuum emission, based on a multiwavelength dataset comprising a broadband spectral energy distribution (SED) and VLTI H-, L-, and N-bands interferometric data obtained in low spectral resolution. Then, we added a ring-like component, representing the carbonaceous L-band features-emitting region, to assess its detectability in future higher spectral resolution observations employing mid-IR interferometry.

The large disparity in physical conditions from the diffuse interstellar medium (ISM) to denser clouds such as photon-dominated regions (PDRs) triggers an evolution of the dust properties (i.e. composition, size, and shape). The gas physics and chemistry are tightly connected to these dust properties and are therefore affected by dust evolution and especially the nano-grain depletion in the outer irradiated part of PDRs. We highlight the influence of nano-grain depletion on the gas physics and chemistry in the Horsehead nebula, a prototypical PDR. We used a model for atomic and molecular gas in PDRs, the Meudon PDR code, using diffuse ISM-like dust and Horsehead-like dust to study the influence of nano-grain depletion on the gas physics and chemistry, focusing on the impact on photoelectric heating and H2 formation and, therefore, on the H2 gas lines. We find that nano-grain depletion in the Horsehead strongly affects gas heating through the photoelectric effect and thus the gas temperature and the H2 formation, hence the H -> H2 position. Consequently, the first four pure rotational lines of H2 (e.g. 0-0 S(0), S(1), S(2), and S(3)) vary by a factor of 2 to 14. The 0-0 S(3) line that is often underestimated in models is underestimated even more when taking nano-grain depletion into account due to the decrease in gas heating through the photoelectric effect. This strongly suggests that our understanding of the excitation of H2 and/or of heating processes in the Horsehead, and more generally in PDRs, is still incomplete. Nano-grain depletion in the outer part of the Horsehead has a strong influence on several gas tracers that will be prominent in JWST observations of irradiated clouds. We therefore need to take this depletion into account in order to improve our understanding of the Horsehead, and more generally PDRs, and to contribute to the optimal scientific return of the mission., Comment: Accepted in A&A, comments welcome

Light scattering at near-infrared wavelengths has been used to study the optical properties of the interstellar dust grains, but these studies are limited by the assumptions on the strength of the radiation field. On the other hand, thermal dust emission can be used to constrain the properties of the radiation field, although this is hampered by uncertainty about the dust emissivity. We test if current dust models allow us to model a molecular cloud simultaneously in the near infrared (NIR) and far infrared (FIR) wavelengths and compare the results with observations. Our aim is to place constraints on the properties of the dust grains and the strength of the radiation field. We present computations of dust emission and scattered light of a quiescent molecular cloud LDN1512. We construct radiative transfer models for LDN1512 that include an anisotropic radiation field and a three-dimensional cloud model. We are able to reproduce the observed FIR observations, with a radiation field derived from the DIRBE observations, with all of the tested dust models. However, with the same density distribution and the assumed radiation field, the models fail to reproduce the observed NIR scattering in all cases except for models that take into account dust evolution via coagulation and mantle formation. We find that the column densities derived from our radiative transfer modelling can differ by a factor of up to two, compared to the column densities derived from the observations with modified blackbody fits. The discrepancy in the column densities is likely caused because of temperature difference between a modified blackbody fit and the real spectra. We show that the observed dust emission can be reproduced with several different assumptions about the properties of the dust grains. However, in order to reproduce the observed scattered surface brightness dust evolution must be taken into account., Comment: 16 pages, 16 figures. Recommended for publication to A&A

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., Comment: 18 pages, 11 figures, accepted for publication in A&A

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.

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., Comment: 13 pages, 6 figures, 2 tables, A&A accepted

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., Comment: accepted by A&A, abstract modified to match arXiv standards

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., Comment: Accepted for publication in PASA; 37 pages, 16 figures

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., Comment: Accepted for publication in Astronomy and Astrophysics

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., Comment: 24 pages, 15 figures, published in A&A, bibref: 2019A&A...623A...5D, This Version 2 has increased accuracy in Table B1 in appendix B

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., Comment: accepted for publication by A&A

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., Comment: accepted for publication by A&A

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., Comment: Editor in chief: C. Ferrari; Editors: M. Alves, S. Bosse, S. Corbel, A. Ferrari, K. Ferri\`ere, S. Gauffre, E. Josselin, G. Lagache, S. Lambert, G. Marquette, J.-M. Martin, M.-A. Miville-Desch\^enes, L. Montier, B. Semelin, G. Theureau, S. Vergani, N. Vilmer, P. Zarka; Original file with high resolution figures at SKA-France link: https://ska-france.oca.eu/images/SKA-France-Media/FWB_051017.pdf

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., Comment: 33 pages, 9 figures (at low resolution here), 8 tables, Accepted for publication in A&A

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., Comment: 9 pages, 1 figure

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., Comment: Accepted in A&A

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., Comment: 66 pages, accepted for publication in PASP

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., Comment: 25 pages, 19 figures, A&A accepted

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., Comment: 24 pages, 30 figures, accepted for publication in A&A

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., Comment: Accepted for publication in A&A

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., Comment: 12 pages, 15 figures, accepted for publication in A&A along with the companion paper entitled, Mantle formation, coagulation and the origin of cloud/core shine: II Comparison with observations, by Ysard et al. (also accepted for publication in A&A)

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., Comment: Paper accepted for publication in A&A with companion paper "Mantle formation, coagulation and the origin of cloud/core-shine: I. Dust scattering and absorption in the IR", A.P Jones, M. Koehler, N. Ysard, E. Dartois, M. Godard, L. Gavilan

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., Comment: Accepted for publication in MNRAS

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., Comment: Accepted for publication in MNRAS

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

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., Comment: 28 pages, 20 figures. A&A accepted. Corresponding author: J.D. Soler

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., Comment: 28 pages, A&A, in press

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., Comment: Accepted for publication in the Astrophysical Journal

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., Comment: 22 pages, 20 figures

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., Comment: Final version that has appeared in A&A

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., Comment: 5 pages, 2 figures, accepted for publication in MNRAS Letters

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., Comment: 22 pages, 18 figures, accepted to Astronomy & Astrophysics; abstract shortened and figures reduced for astroph