Showing 1–41 of 41 results for author: Guillet, V

THEMIS 2.0: A self-consistent model for dust extinction, emission, and polarisation

Authors: Nathalie Ysard, Anthony Peter Jones, Vincent Guillet, Karine Demyk, Marjorie Decleir, Laurent Verstraete, Ilyes Choubani, Marc-Antoine Miville-Deschênes, Lapo Fanciullo

Abstract: Recent observations in emission, extinction, and polarisation have at least partially invalidated most of the astronomical standard grain models for the diffuse ISM. Moreover, lab measurements on interstellar silicate analogues have shown differences with the optical properties used in these standard models. To address these issues, our objective is twofold: (i) to update the optical properties of… ▽ More Recent observations in emission, extinction, and polarisation have at least partially invalidated most of the astronomical standard grain models for the diffuse ISM. Moreover, lab measurements on interstellar silicate analogues have shown differences with the optical properties used in these standard models. To address these issues, our objective is twofold: (i) to update the optical properties of silicates and (ii) to develop the THEMIS dust model to allow the calculation of polarised extinction and emission. Based on optical constants measured in the lab for amorphous silicates and on observational constraints in mid-IR extinction and X-ray scattering, we defined new optical constants for the THEMIS silicates. Absorption and scattering efficiencies for spheroidal grains were then derived with the discrete dipole approximation. These new optical properties make it possible to explain the dust emission and extinction, both total and polarised. The model is not yet pushed to its limits since it does not require the perfect alignment of all grains to explain the observations and it therefore has the potential to accommodate the highest polarisation levels inferred from extinction measures. Moreover, the dispersion of the optical properties of the different lab silicates naturally explain the variations in both the total and polarised emission and extinction observed in the diffuse ISM. A single, invariant model calibrated on one single set of observations is obsolete for explaining contemporary observations. We are proposing a completely flexible dust model based entirely on lab measurements that has the potential to make major advances in understanding the nature of ISM grains and how they evolve as a function of their environment. Even if challenging, this is also relevant for future missions that will aim to perform precise measurements of the CMB spectral distortions and polarisation. △ Less

Submitted 30 May, 2024; v1 submitted 15 January, 2024; originally announced January 2024.

Comments: Accepted for publication by A&A ; abstract modified to match the arxiv standard ; typo corrected on 19/01/2024 in Tab. 3 ; units corrected on 07/02/2024 in Figs. 9 and C.2

Fast methods for tracking grain coagulation and ionization. III. Protostellar collapse with non-ideal MHD

Authors: Pierre Marchand, Ugo Lebreuilly, Mordecai-Mark Mac Low, Vincent Guillet

Abstract: Dust grains influence many aspects of star formation, including planet formation, opacities for radiative transfer, chemistry, and the magnetic field via Ohmic, Hall, and ambipolar diffusion. The size distribution of the dust grains is the primary characteristic influencing all these aspects. Grain size increases by coagulation throughout the star formation process. We describe here numerical simu… ▽ More Dust grains influence many aspects of star formation, including planet formation, opacities for radiative transfer, chemistry, and the magnetic field via Ohmic, Hall, and ambipolar diffusion. The size distribution of the dust grains is the primary characteristic influencing all these aspects. Grain size increases by coagulation throughout the star formation process. We describe here numerical simulations of protostellar collapse using methods described in earlier papers of this series. We compute the evolution of the grain size distribution from coagulation and the non-ideal magnetohydrodynamics effects self-consistently and at low numerical cost. We find that the coagulation efficiency is mostly affected by the time spent in high-density regions. Starting from sub-micron radii, grain sizes reach more than 100 μm in an inner protoplanetary disk that is only 1000 years old. We also show that the growth of grains significantly affects the resistivities, and indirectly the dynamics and angular momentum of the disk. △ Less

Submitted 9 January, 2023; v1 submitted 4 January, 2023; originally announced January 2023.

Comments: 11 pages, 9 figures, accepted in A&A

Journal ref: A&A 670, A61 (2023)

Protostellar collapse simulations in spherical geometry with dust coagulation and fragmentation

Authors: Ugo Lebreuilly, Valentin Vallucci-Goy, Vincent Guillet, Maxime Lombart, Pierre Marchand

Abstract: We model the coagulation and fragmentation of dust grains during the protostellar collapse with our newly developed shark code. It solves the gas-dust hydrodynamics in a spherical geometry and the coagulation/fragmentation equation. It also computes the ionization state of the cloud and the Ohmic, ambipolar and Hall resistivities. We find that the dust size distribution evolves significantly durin… ▽ More We model the coagulation and fragmentation of dust grains during the protostellar collapse with our newly developed shark code. It solves the gas-dust hydrodynamics in a spherical geometry and the coagulation/fragmentation equation. It also computes the ionization state of the cloud and the Ohmic, ambipolar and Hall resistivities. We find that the dust size distribution evolves significantly during the collapse, large grain formation being controlled by the turbulent differential velocity. When turbulence is included, only ambipolar diffusion remains efficient at removing the small grains from the distribution, brownian motion is only efficient as a standalone process. The macroscopic gas-dust drift is negligible for grain growth and only dynamically significant near the first Larson core. At high density, we find that the coagulated distribution is unaffected by the initial choice of dust distribution. Strong magnetic fields are found to enhance the small grains depletion, causing an important increase of the ambipolar diffusion. This hints that the magnetic field strength could be regulated by the small grain population during the protostellar collapse. Fragmentation could be effective for bare silicates, but its modeling relies on the choice of ill-constrained parameters. It is also found to be negligible for icy grains. When fragmentation occurs, it strongly affects the magnetic resistivities profiles. Dust coagulation is a critical process that needs to be fully taken into account during the protostellar collapse. The onset and feedback of fragmentation remains uncertain and its modeling should be further investigated. △ Less

Submitted 3 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS

Frequency dependence of the thermal dust $E/B$ ratio and $EB$ correlation: insights from the spin-moment expansion

Authors: Léo Vacher, Jonathan Aumont, François Boulanger, Ludovic Montier, Vincent Guillet, Alessia Ritacco, Jens Chluba

Abstract: The change of physical conditions across the turbulent and magnetized interstellar medium (ISM) induces a 3D spatial variation of the properties of Galactic polarized emission. The observed signal results from the averaging of different spectral energy distributions (SED) and polarization angles, along and between lines of sight. As a consequence, the total Stokes parameters $Q$ and $U$ will have… ▽ More The change of physical conditions across the turbulent and magnetized interstellar medium (ISM) induces a 3D spatial variation of the properties of Galactic polarized emission. The observed signal results from the averaging of different spectral energy distributions (SED) and polarization angles, along and between lines of sight. As a consequence, the total Stokes parameters $Q$ and $U$ will have different distorted SEDs, so that the polarization angle becomes frequency dependent. In the present work, we show how this phenomenon similarly induces a different distorted SED for the three polarized angular power spectra $EE$, $BB$ and $EB$, implying a variation of the $EE/BB$ ratio with frequency. We demonstrate how the previously introduced spin-moment formalism provides a natural framework to grasp these effects, allowing us to derive analytical predictions for the spectral behaviors of the polarized spectra, focusing here on the example of thermal dust polarized emission. After a quantitative discussion based on a model combining emission from a filament with its background, we further reveal that the spectral complexity implemented in the dust models commonly used by the cosmic microwave background (CMB) community produce such effects. This new understanding is crucial for CMB component separation, in which an extreme accuracy is required in the modeling of the dust signal to allow for the search of the primordial imprints of inflation or cosmic birefringence. For the latter, as long as the dust $EB$ signal is not measured accurately, great caution is required about the assumptions made to model its spectral behavior, as it may not simply follow from the other dust angular power spectra. △ Less

Submitted 2 May, 2023; v1 submitted 26 October, 2022; originally announced October 2022.

Comments: Published in A&A

Journal ref: A&A 672, A146 (2023)

Dust polarization spectral dependence from Planck HFI data. Turning point on CMB polarization foregrounds modelling

Authors: Alessia Ritacco, François Boulanger, Vincent Guillet, Jean-Marc Delouis, Jean-Loup Puget, Jonathan Aumont, Léo Vacher

Abstract: The search for the primordial B-modes of the cosmic microwave background (CMB) relies on the separation from the brighter foreground dust signal. In this context, the characterisation of the spectral energy distribution (SED) of thermal dust in polarization has become a critical subject of study. We present a power-spectra analysis of Planck data, which improves on previous studies by using the ne… ▽ More The search for the primordial B-modes of the cosmic microwave background (CMB) relies on the separation from the brighter foreground dust signal. In this context, the characterisation of the spectral energy distribution (SED) of thermal dust in polarization has become a critical subject of study. We present a power-spectra analysis of Planck data, which improves on previous studies by using the newly released SRoll2 maps that correct residual data systematics, and by extending the analysis to regions near the Galactic plane. Our analysis focuses on the lowest multipoles between l=4 and 32, and three sky areas with sky fractions of fsky = 80%, 90%, and 97%. The mean dust SED for polarization and the 353 GHz Q and U maps are used to compute residual maps at 100, 143 and 217 GHz, highlighting spatial variations of the dust polarization SED. Residuals are detected at the three frequencies for the three sky areas. We show that models based on total intensity data are underestimating by a significant factor the complexity of dust polarized CMB foreground. Our analysis emphasizes the need to include variations of polarization angles of the dust polarized CMB foreground. The frequency dependence of the EE and BB power spectra of the residual maps yields further insight. We find that the moments expansion to the first order of the modified black-body (MBB) spectrum provides a good fit to the EE power-spectra. This result suggests that the residuals could follow mainly from variations of dust MBB spectral parameters. However, this conclusion is challenged by cross-spectra showing that the residuals maps at the three frequencies are not fully correlated, and the fact that the BB power-spectra do not match the first order moment expansion of a MBB SED. This work sets new requirements for simulations of the dust polarized foreground and component separation methods (abridged) △ Less

Submitted 20 December, 2022; v1 submitted 15 June, 2022; originally announced June 2022.

Comments: 13 pages, 13 figures, Accepted for publication in the Cosmology section of A&A journal

Journal ref: A&A 670, A163 (2023)

Performance of the polarization leakage correction in the PILOT data

Authors: J-Ph. Bernard, A. Bernard, H. Roussel, I. Choubani, D. Alina, J. Aumont, A. Hughes, I. Ristorcelli, S. Stever, T. Matsumura S. Sugiyama, K. Komatsu, G. de Gasperis, K. Ferriere, V. Guillet, N. Ysard, P. Ade, P. de Bernardis, N. Bray, B. Crane, J. P. Dubois, M. Griffin, P. Hargrave, Y. Longval, S. Louvel, B. Maffei , et al. (11 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 7 May, 2022; originally announced May 2022.

Fast methods to track grain coagulation and ionization. II. Extension to thermal ionization

Authors: Pierre Marchand, Vincent Guillet, Ugo Lebreuilly, Mordecai-Mark Mac Low

Abstract: Thermal ionization is a critical process at temperatures T > 10 3 K, particularly during star formation. An increase in ionization leads to a decrease in nonideal magnetohydrodynamics (MHD) resistivities, which has a significant impact on protoplanetary disks and protostar formation. We developed an extension of the fast computational ionization method presented in our recent paper to include ther… ▽ More Thermal ionization is a critical process at temperatures T > 10 3 K, particularly during star formation. An increase in ionization leads to a decrease in nonideal magnetohydrodynamics (MHD) resistivities, which has a significant impact on protoplanetary disks and protostar formation. We developed an extension of the fast computational ionization method presented in our recent paper to include thermal ionization. The model can be used to inexpensively calculate the density of ions and electrons and the electric charge of each size of grains for an arbitrary size distribution. This tool should be particularly useful for the self-consistent calculation of nonideal MHD resistivities in multidimensional simulations, especially of protostellar collapse and protoplanetary disks. △ Less

Submitted 30 June, 2022; v1 submitted 23 February, 2022; originally announced February 2022.

Comments: Accepted in A&A. 7 pages, 6 figures

Journal ref: A&A 666, A27 (2022)

Fast methods for tracking grain coagulation and ionization. I. Analytic derivation

Authors: Pierre Marchand, Vincent Guillet, Ugo Lebreuilly, Mordecai-Mark Mac Low

Abstract: Dust grains play a major role in many astrophysical contexts. They affect the chemical, magnetic, dynamical, and optical properties of their environment, from galaxies down to the interstellar medium, star-forming regions, and protoplanetary disks. Their coagulation leads to shifts in their size distribution and ultimately to the formation of planets. However, although the coagulation process is r… ▽ More Dust grains play a major role in many astrophysical contexts. They affect the chemical, magnetic, dynamical, and optical properties of their environment, from galaxies down to the interstellar medium, star-forming regions, and protoplanetary disks. Their coagulation leads to shifts in their size distribution and ultimately to the formation of planets. However, although the coagulation process is reasonably uncomplicated to numerically implement by itself, it is difficult to couple it with multidimensional hydrodynamics numerical simulations because of its high computational cost. We propose here a simple method for tracking the coagulation of grains at far lower cost. Given an initial grain size distribution, the state of the distribution at time t is solely determined by the value of a single variable integrated along the trajectory, independently of the specific path taken by the grains. Although this method cannot account for processes other than coagulation, it is mathematically exact, fast, inexpensive, and can be used to evaluate the effect of grain coagulation in most astrophysical contexts. It is applicable to all coagulation kernels in which local physical conditions and grain properties can be separated. We also describe another method for calculating the average electric charge of grains and the density of ions and electrons in environments that are shielded from radiation fields, given the density and temperature of the gas, the cosmic-ray ionization rate, and the average mass of the ions. The equations we provide are fast to integrate numerically and can be used in multidimensional numerical simulations to self-consistently calculate on the fly the local resistivities that are required to model nonideal magnetohydrodynamics. △ Less

Submitted 17 May, 2021; v1 submitted 25 February, 2021; originally announced March 2021.

Comments: 11 pages, 5 figures, accepted in A&A

Journal ref: A&A 649, A50 (2021)

A statistical analysis of dust polarization properties in ALMA observations of Class 0 protostellar cores

Authors: V. J. M. Le Gouellec, A. J. Maury, V. Guillet, C. L. H. Hull, J. M. Girart, A. Verliat, R. Mignon-Risse, V. Valdivia, P. Hennebelle, M. González, F. Louvet

Abstract: Recent observational progress has challenged the dust grain-alignment theories used to explain the polarized dust emission routinely observed in star-forming cores. In an effort to improve our understanding of the dust grain alignment mechanism(s), we have gathered a dozen ALMA maps of (sub)millimeter-wavelength polarized dust emission from Class 0 protostars, and carried out a comprehensive stati… ▽ More Recent observational progress has challenged the dust grain-alignment theories used to explain the polarized dust emission routinely observed in star-forming cores. In an effort to improve our understanding of the dust grain alignment mechanism(s), we have gathered a dozen ALMA maps of (sub)millimeter-wavelength polarized dust emission from Class 0 protostars, and carried out a comprehensive statistical analysis of dust polarization quantities. We analyze the statistical properties of the polarization fraction P_frac and dispersion of polarization position angles S. More specifically, we investigate the relationship between S and P_frac as well as the evolution of the product S*P_frac as a function of the column density of the gas in the protostellar envelopes. We find a significant correlation in the polarized dust emission from protostellar envelopes seen with ALMA; the power-law index differs significantly from the one observed by Planck in star-forming clouds. The product S*P_frac, which is sensitive to the dust grain alignment efficiency, is approximately constant across three orders of magnitude in envelope column density. This suggests that the grain alignment mechanism producing the bulk of the polarized dust emission in star-forming cores may not depend systematically on the local conditions such as local gas density. Ultimately, our results suggest dust alignment mechanism(s) are efficient at producing dust polarized emission in the various local conditions typical of Class 0 protostars. The grain alignment efficiency found in these objects seems to be higher than the efficiency produced by the standard RAT alignment of paramagnetic grains. Further study will be needed to understand how more efficient grain alignment via, e.g., different irradiation conditions, dust grain characteristics, or additional grain alignment mechanisms can reproduce the observations. △ Less

Submitted 15 September, 2020; originally announced September 2020.

Journal ref: A&A 644, A11 (2020)

Dust coagulation feedback on magnetohydrodynamic resistivities in protostellar collapse

Authors: V. Guillet, P. Hennebelle, G. Pineau des Forêts, A. Marcowith, B. Commerçon, P. Marchand

Abstract: The degree of coupling between the gas and the magnetic field during the collapse of a core and the subsequent formation of a disk depends on the assumed dust size distribution. We study the impact of grain-grain coagulation on the evolution of magnetohydrodynamic (MHD) resistivities during the collapse of a prestellar core. We use a 1-D model to follow the evolution of the dust size distribution,… ▽ More The degree of coupling between the gas and the magnetic field during the collapse of a core and the subsequent formation of a disk depends on the assumed dust size distribution. We study the impact of grain-grain coagulation on the evolution of magnetohydrodynamic (MHD) resistivities during the collapse of a prestellar core. We use a 1-D model to follow the evolution of the dust size distribution, out-of-equilibrium ionization state and gas chemistry during the collapse of a prestellar core. To compute the grain-grain collisional rate, we consider models for both random and systematic, size-dependent, velocities. We include grain growth through grain-grain coagulation and ice accretion, but ignore grain fragmentation. Starting with a MRN (Mathis et al. 1977) size distribution, we find that coagulation in grain-grain collisions generated by hydrodynamical turbulence is not efficient at removing the smallest grains, and as a consequence does not affect much the evolution of the Hall and ambipolar diffusion MHD resistivities which still severly drop during the collapse like in models without coagulation. The inclusion of systematic velocities, possibly induced by the presence of ambipolar diffusion, increases the coagulation rate between small and large grains, removing small grains earlier in the collapse and therefore limiting the drop in the Hall and ambipolar diffusion resistivities. At intermediate densities ($n_{\rm H} \sim 10^8\,$cm$^{-3}$), the Hall and ambipolar diffusion resistivities are found to be higher by 1 to 2 orders of magnitude in models with coagulation than in models where coagulation is ignored, and also higher than in a toy model without coagulation where all grains smaller than $0.1\,μ$m would have been removed in the parent cloud before the collapse. When grain drift velocities induced by ambipolar diffusion are included, dust coagulation ... (abridged) △ Less

Submitted 8 July, 2020; originally announced July 2020.

Journal ref: A&A 643, A17 (2020)

A systematic study of radiative torque grain alignment in the diffuse interstellar medium

Authors: Stefan Reissl, Vincent Guillet, Robert Brauer, François Levrier, François Boulanger, Ralf S. Klessen

Abstract: Context. Planck observations demonstrated that the grain alignment efficiency is almost constant in the diffuse ISM. Aims. We test if the Radiative Torque (RAT) theory is compatible with observational constraints on grain alignment. Methods. We combine a numerical simulation with the radiative transfer code POLARIS that incorporates a physical dust model and the detailed grain alignment physics of… ▽ More Context. Planck observations demonstrated that the grain alignment efficiency is almost constant in the diffuse ISM. Aims. We test if the Radiative Torque (RAT) theory is compatible with observational constraints on grain alignment. Methods. We combine a numerical simulation with the radiative transfer code POLARIS that incorporates a physical dust model and the detailed grain alignment physics of RATs. A dust model is designed to reproduce the spectral dependence of extinction of the ISM. From a RAMSES simulation of interstellar turbulence, we extract a cube representative of the diffuse ISM. We post-process the cube with POLARIS to get the grain temperature and RATs to simulate synthetic dust polarization maps. Results. In our simulation the grain alignment efficiency is correlated with gas pressure, but not with the RAT intensity. Because of the low dust extinction, the magnitude of RATs varies little, decreasing only for high column densities $N_H$. Comparing our maps with a uniform alignment efficiency, we find no systematic difference. The dependence of polarization fraction $p$ with $N_H$ or polarization dispersion $S$ is similar. The drop of RATs in dense regions barely affects the polarization pattern, the signal being dominated by the LOS and magnetic field geometry. If a star is inserted, the polarization increases, with no specific pattern around the star. The angle-dependence of RATs is not observed in the maps, and is weak using a uniform magnetic field. Conclusions. RATs are compatible with Planck data for the diffuse ISM such that both uniform alignment and RAT alignment lead to similar observations. To further test the predictions of RATs where an important drop of grain alignment is expected, polarization observations of dense regions must be confronted to numerical simulations sampling high column densities through dense clouds, with enough statistics. △ Less

Submitted 26 February, 2020; originally announced February 2020.

Comments: 22 pages, 22 figures

Journal ref: A&A 640, A118 (2020)

Polarized emission by aligned grains in the Mie regime : application to protoplanetary disks observed by ALMA

Authors: V. Guillet, J. M. Girart, A. J. Maury, F. O. Alves

Abstract: The azimuthal polarization patterns observed in some protoplanetary disks by ALMA at millimeter wavelength have raised doubts about their being produced by dust grains aligned with the magnetic field lines. These conclusions were based on the calculations of dust polarized emission in the Rayleigh regime, i.e. for grain sizes much smaller than the wavelength. However, the grain size in such disks… ▽ More The azimuthal polarization patterns observed in some protoplanetary disks by ALMA at millimeter wavelength have raised doubts about their being produced by dust grains aligned with the magnetic field lines. These conclusions were based on the calculations of dust polarized emission in the Rayleigh regime, i.e. for grain sizes much smaller than the wavelength. However, the grain size in such disks is estimated to be typically in the range 0.1 - 1 mm from independent observations. We study the dust polarization properties of aligned grains in emission in the Mie regime, i.e. when the mean grain size approches the wavelength. Using the T-MATRIX and DustEM codes, we compute the spectral dependence of the polarization fraction in emission for grains in perfect spinning alignment, for various grain size distributions of weakly-elongated oblate and prolate grains of astrosilicate composition, with a mean size ranging from 10 μm to 1 mm. In the submillimeter and millimeter wavelength range, the polarization by B-field aligned grains becomes negative for grains larger than ~ 250 μm, meaning that the polarization vector becomes parallel to the B-field. The transition from the positive to the negative polarization occurs at a wavelength λ ~ 1 mm. The regime of negative polarization does not exist for grains smaller than ~ 100 μm. When using realistic grain size distributions for disks with grains up to the submillimeter sizes, the polarization direction of thermal emission by aligned grains is shown to be parallel to the direction of the magnetic field over a significant fraction of the wavelengths typically used to observe young protoplanetary disks. This property may explain the peculiar azimuthal orientation of the polarization vectors in some of the disks observed with ALMA and attest of the conserved ability of dust polarized emission to trace the magnetic field in disks. △ Less

Submitted 23 January, 2020; originally announced January 2020.

Journal ref: A&A 634, L15 (2020)

Indirect evidence of significant grain growth in young protostellar envelopes from polarized dust emission

Authors: Valeska Valdivia, Anaelle Maury, Robert Brauer, Patrick Hennebelle, Maud Galametz, Vincent Guillet, Stefan Reissl

Abstract: How and when in the star formation sequence do dust grains start to grow into pebbles is a cornerstone question to both star and planet formation. We compute the polarized radiative transfer from a model solar-type protostellar core, using the POLARIS code, aligning the dust grains with the local magnetic field, following the radiative torques (RATs) theory. We test the dependency of the resulting… ▽ More How and when in the star formation sequence do dust grains start to grow into pebbles is a cornerstone question to both star and planet formation. We compute the polarized radiative transfer from a model solar-type protostellar core, using the POLARIS code, aligning the dust grains with the local magnetic field, following the radiative torques (RATs) theory. We test the dependency of the resulting dust polarized emission with the maximum grain size of the dust size distribution at the envelope scale, from amax = 1 micron to 50 micron. Our work shows that, in the framework of RAT alignment, large dust grains are required to produce polarized dust emission at levels similar to those currently observed in solar-type protostellar envelopes at millimeter wavelengths. Considering the current theoretical dificulties to align a large fraction of small ISM-like grains in the conditions typical of protostellar envelopes, our results suggest that grain growth (typically > 10 micron) might have already significantly progressed at scales 100-1000 au in the youngest objects, observed less than 10^5 years after the onset of collapse. Observations of dust polarized emission might open a new avenue to explore dust pristine properties and describe, for example, the initial conditions for the formation of planetesimals. △ Less

Submitted 25 July, 2019; originally announced July 2019.

Comments: Accepted for publication in MNRAS. Accepted 2019 July 23. Received 2019 July 23; in original form 2019 April 05

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

Authors: Ph. André, A. Hughes, V. Guillet, F. Boulanger, A. Bracco, E. Ntormousi, D. Arzoumanian, A. J. Maury, J. -Ph. Bernard, S. Bontemps, I. Ristorcelli, J. M. Girart, F. Motte, K. Tassis, E. Pantin, T. Montmerle, D. Johnstone, S. Gabici, A. Efstathiou, Shantanu Basu, M. Béthermin, H. Beuther, J. Braine, J. Di Francesco, E. Falgarone , et al. (31 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 9 May, 2019; originally announced May 2019.

Comments: Accepted for publication in PASA; 37 pages, 16 figures

Journal ref: PASA, 36, e029 (2019)

Dust spectrum and polarisation at 850 um in the massive IRDC G035.39-00.33

Authors: Mika Juvela, Vincent Guillet, Tie Liu, Isabelle Ristorcelli, Veli-Matti Pelkonen, Dana Alina, Leonardo Bronfman, David J. Eden, Kee Tae Kim, Patrick M. Koch, Woojin Kwon, Chang Won Lee, Johanna Malinen, Elisabetta Micelotta, Julien Montillaud, Mark G. Rawlings, Patricio Sanhueza, Archana Soam, Alessio Traficante, Nathalie Ysard, Chuan-Peng Zhang

Abstract: 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… ▽ More 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. △ Less

Submitted 6 September, 2018; v1 submitted 4 September, 2018; originally announced September 2018.

Comments: accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 620, A26 (2018)

Planck 2018 results. XII. Galactic astrophysics using polarized dust emission

Authors: Planck Collaboration, N. Aghanim, Y. Akrami, M. I. R. Alves, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, A. Bracco, M. Bucher, C. Burigana, E. Calabrese , et al. (138 additional authors not shown)

Abstract: We present 353 GHz full-sky maps of the polarization fraction $p$, angle $ψ$, and dispersion of angles $S$ of Galactic dust thermal emission produced from the 2018 release of Planck data. We confirm that the mean and maximum of $p$ decrease with increasing $N_H$. The uncertainty on the maximum polarization fraction, $p_\mathrm{max}=22.0$% at 80 arcmin resolution, is dominated by the uncertainty on… ▽ More We present 353 GHz full-sky maps of the polarization fraction $p$, angle $ψ$, and dispersion of angles $S$ of Galactic dust thermal emission produced from the 2018 release of Planck data. We confirm that the mean and maximum of $p$ decrease with increasing $N_H$. The uncertainty on the maximum polarization fraction, $p_\mathrm{max}=22.0$% at 80 arcmin resolution, is dominated by the uncertainty on the zero level in total intensity. The observed inverse behaviour between $p$ and $S$ is interpreted with models of the polarized sky that include effects from only the topology of the turbulent Galactic magnetic field. Thus, the statistical properties of $p$, $ψ$, and $S$ mostly reflect the structure of the magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map $S \times p$, looking for residual trends. While $p$ decreases by a factor of 3--4 between $N_H=10^{20}$ cm$^{-2}$ and $N_H=2\times 10^{22}$ cm$^{-2}$, $S \times p$ decreases by only about 25%, a systematic trend observed in both the diffuse ISM and molecular clouds. Second, we find no systematic trend of $S \times p$ with the dust temperature, even though in the diffuse ISM lines of sight with high $p$ and low $S$ tend to have colder dust. We also compare Planck data with starlight polarization in the visible at high latitudes. The agreement in polarization angles is remarkable. Two polarization emission-to-extinction ratios that characterize dust optical properties depend only weakly on $N_H$ and converge towards the values previously determined for translucent lines of sight. We determine an upper limit for the polarization fraction in extinction of 13%, compatible with the $p_\mathrm{max}$ observed in emission. These results provide strong constraints for models of Galactic dust in diffuse gas. △ Less

Submitted 7 March, 2019; v1 submitted 17 July, 2018; originally announced July 2018.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 641, A12 (2020)

Planck 2018 results. XI. Polarized dust foregrounds

Authors: Planck Collaboration, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, A. Bracco, M. Bucher, C. Burigana, E. Calabrese, J. -F. Cardoso, J. Carron, H. C. Chiang , et al. (109 additional authors not shown)

Abstract: The study of polarized dust emission has become entwined with the analysis of the cosmic microwave background (CMB) polarization. We use new Planck maps to characterize Galactic dust emission as a foreground to the CMB polarization. We present Planck EE, BB, and TE power spectra of dust polarization at 353 GHz for six nested sky regions covering from 24 to 71 % of the sky. We present power-law fit… ▽ More The study of polarized dust emission has become entwined with the analysis of the cosmic microwave background (CMB) polarization. We use new Planck maps to characterize Galactic dust emission as a foreground to the CMB polarization. We present Planck EE, BB, and TE power spectra of dust polarization at 353 GHz for six nested sky regions covering from 24 to 71 % of the sky. We present power-law fits to the angular power spectra, yielding evidence for statistically significant variations of the exponents over sky regions and a difference between the values for the EE and BB spectra. The TE correlation and E/B power asymmetry extend to low multipoles that were not included in earlier Planck polarization papers. We also report evidence for a positive TB dust signal. Combining data from Planck and WMAP, we determine the amplitudes and spectral energy distributions (SEDs) of polarized foregrounds, including the correlation between dust and synchrotron polarized emission, for the six sky regions as a function of multipole. This quantifies the challenge of the component separation procedure required for detecting the reionization and recombination peaks of primordial CMB B modes. The SED of polarized dust emission is fit well by a single-temperature modified blackbody emission law from 353 GHz to below 70 GHz. For a dust temperature of 19.6 K, the mean spectral index for dust polarization is $β_{\rm d}^{P} = 1.53\pm0.02 $. By fitting multi-frequency cross-spectra, we examine the correlation of the dust polarization maps across frequency. We find no evidence for decorrelation. If the Planck limit for the largest sky region applies to the smaller sky regions observed by sub-orbital experiments, then decorrelation might not be a problem for CMB experiments aiming at a primordial B-mode detection limit on the tensor-to-scalar ratio $r\simeq0.01$ at the recombination peak. △ Less

Submitted 12 November, 2018; v1 submitted 15 January, 2018; originally announced January 2018.

Comments: Final version to appear in A&A

Journal ref: A&A 641, A11 (2020)

Dust models compatible with Planck intensity and polarization data in translucent lines of sight

Authors: V. Guillet, L. Fanciullo, L. Verstraete, F. Boulanger, A. P. Jones, M. -A. Miville-Deschênes, N. Ysard, F. Levrier, M. Alves

Abstract: 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… ▽ More 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. △ Less

Submitted 12 October, 2017; originally announced October 2017.

Journal ref: A&A 610, A16 (2018)

Interplay of dust alignment, grain growth and magnetic fields in polarization: lessons from the emission-to-extinction ratio

Authors: Lapo Fanciullo, Vincent Guillet, François Boulanger, Anthony Jones

Abstract: Polarized extinction and emission from dust in the interstellar medium (ISM) are hard to interpret, as they have a complex dependence on dust optical properties, grain alignment and magnetic field orientation. This is particularly true in molecular clouds. The data available today are not yet used to their full potential. The combination of emission and extinction, in particular, provides inform… ▽ More Polarized extinction and emission from dust in the interstellar medium (ISM) are hard to interpret, as they have a complex dependence on dust optical properties, grain alignment and magnetic field orientation. This is particularly true in molecular clouds. The data available today are not yet used to their full potential. The combination of emission and extinction, in particular, provides information not available from either of them alone. We combine data from the scientific literature on polarized dust extinction with Planck data on polarized emission and we use them to constrain the possible variations in dust and environmental conditions inside molecular clouds, and especially translucent lines of sight, taking into account magnetic field orientation. We focus on the dependence between λ_max -- the wavelength of maximum polarization in extinction -- and other observables such as the extinction polarization, the emission polarization and the ratio of the two. We set out to reproduce these correlations using Monte-Carlo simulations where the relevant quantities in a dust model -- grain alignment, size distribution and magnetic field orientation -- vary to mimic the diverse conditions expected inside molecular clouds. None of the quantities chosen can explain the observational data on its own: the best results are obtained when all quantities vary significantly across and within clouds. However, some of the data -- most notably the stars with low emission-to-extinction polarization ratio -- are not reproduced by our simulation. Our results suggest not only that dust evolution is necessary to explain polarization in molecular clouds, but that a simple change in size distribution is not sufficient to explain the data, and point the way for future and more sophisticated models. △ Less

Submitted 8 March, 2017; v1 submitted 27 February, 2017; originally announced February 2017.

Journal ref: A&A 602, A7 (2017)

Magnetic field morphology in nearby molecular clouds as revealed by starlight and submillimetre polarization

Authors: J. D. Soler, F. Alves, F. Boulanger, A. Bracco, E. Falgarone, G. A. P. Franco, V. Guillet, P. Hennebelle, F. Levrier, P. G. Martin, M. -A. Miville-Deschênes

Abstract: Within four nearby (d < 160 pc) molecular clouds, we statistically evaluate the structure of the interstellar magnetic field, projected on the plane of the sky and integrated along the line of sight, as inferred from the polarized thermal emission of Galactic dust observed by Planck at 353 GHz and from the optical and NIR polarization of background starlight. We compare the dispersion of the field… ▽ More Within four nearby (d < 160 pc) molecular clouds, we statistically evaluate the structure of the interstellar magnetic field, projected on the plane of the sky and integrated along the line of sight, as inferred from the polarized thermal emission of Galactic dust observed by Planck at 353 GHz and from the optical and NIR polarization of background starlight. We compare the dispersion of the field orientation directly in vicinities with an area equivalent to that subtended by the Planck effective beam at 353 GHz (10') and using the second-order structure functions of the field orientation angles. We find that the average dispersion of the starlight-inferred field orientations within 10'-diameter vicinities is less than 20 deg, and that at these scales the mean field orientation is on average within 5 deg of that inferred from the submillimetre polarization observations in the considered regions. We also find that the dispersion of starlight polarization orientations and the polarization fractions within these vicinities are well reproduced by a Gaussian model of the turbulent structure of the magnetic field, in agreement with the findings reported by the Planck collaboration at scales greater than 10' and for comparable column densities. At scales greater than 10', we find differences of up to 14.7 deg between the second-order structure functions obtained from starlight and submillimetre polarization observations in the same positions in the plane of the sky, but comparison with a Gaussian model of the turbulent structure of the magnetic field indicates that these differences are small and are consistent with the difference in angular resolution between both techniques. △ Less

Submitted 30 May, 2016; originally announced May 2016.

Comments: 15 pages, 10 figures, submitted to A&A

Report number: AA/2016/28996

Journal ref: A&A 596, A93 (2016)

Planck intermediate results. XLIV. The structure of the Galactic magnetic field from dust polarization maps of the southern Galactic cap

Authors: Planck Collaboration, N. Aghanim, M. I. R. Alves, D. Arzoumanian, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, A. Bracco, M. Bucher, C. Burigana, E. Calabrese, J. -F. Cardoso , et al. (143 additional authors not shown)

Abstract: We study the statistical properties of interstellar dust polarization at high Galactic latitude, using the Stokes parameter Planck maps at 353 GHz. Our aim is to advance the understanding of the magnetized interstellar medium (ISM), and to provide a model of the polarized dust foreground for cosmic microwave background component-separation procedures. Focusing on the southern Galactic cap, we exam… ▽ More We study the statistical properties of interstellar dust polarization at high Galactic latitude, using the Stokes parameter Planck maps at 353 GHz. Our aim is to advance the understanding of the magnetized interstellar medium (ISM), and to provide a model of the polarized dust foreground for cosmic microwave background component-separation procedures. Focusing on the southern Galactic cap, we examine the statistical distributions of the polarization fraction ($p$) and angle ($ψ$) to characterize the ordered and turbulent components of the Galactic magnetic field (GMF) in the solar neighbourhood. We relate patterns at large angular scales in polarization to the orientation of the mean (ordered) GMF towards Galactic coordinates $(l_0,b_0)=(70^\circ \pm 5^\circ,24^\circ \pm 5^\circ)$. The histogram of $p$ shows a wide dispersion up to 25 %. The histogram of $ψ$ has a standard deviation of $12^\circ$ about the regular pattern expected from the ordered GMF. We use these histograms to build a phenomenological model of the turbulent component of the GMF, assuming a uniform effective polarization fraction ($p_0$) of dust emission. To model the Stokes parameters, we approximate the integration along the line of sight (LOS) as a sum over a set of $N$ independent polarization layers, in each of which the turbulent component of the GMF is obtained from Gaussian realizations of a power-law power spectrum. We are able to reproduce the observed $p$ and $ψ$ distributions using: a $p_0$ value of (26 $\pm$ 3)%; a ratio of 0.9 $\pm$ 0.1 between the strengths of the turbulent and mean components of the GMF; and a small value of $N$. We relate the polarization layers to the density structure and to the correlation length of the GMF along the LOS. △ Less

Submitted 4 April, 2016; originally announced April 2016.

Journal ref: A&A 596, A105 (2016)

Dust models post-Planck: constraining the far-infrared opacity of dust in the diffuse interstellar medium

Authors: Lapo Fanciullo, Vincent Guillet, Gonzalo Aniano, Anthony P. Jones, Nathalie Ysard, Marc-Antoine Miville-Deschênes, François Boulanger, M. Köhler

Abstract: 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… ▽ More 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. △ Less

Submitted 23 June, 2015; originally announced June 2015.

Comments: A&A, in press

Planck intermediate results. XXXVIII. E- and B-modes of dust polarization from the magnetized filamentary structure of the interstellar medium

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, N. Bartolo, E. Battaner, K. Benabed, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, A. Bonaldi, L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, A. Bracco, C. Burigana, E. Calabrese , et al. (170 additional authors not shown)

Abstract: The quest for a B-mode imprint from primordial gravity waves on the polarization of the cosmic microwave background (CMB) requires the characterization of foreground polarization from Galactic dust. We present a statistical study of the filamentary structure of the 353 GHz Planck Stokes maps at high Galactic latitude, relevant to the study of dust emission as a polarized foreground to the CMB. We… ▽ More The quest for a B-mode imprint from primordial gravity waves on the polarization of the cosmic microwave background (CMB) requires the characterization of foreground polarization from Galactic dust. We present a statistical study of the filamentary structure of the 353 GHz Planck Stokes maps at high Galactic latitude, relevant to the study of dust emission as a polarized foreground to the CMB. We filter the intensity and polarization maps to isolate filaments in the range of angular scales where the power asymmetry between E-modes and B-modes is observed. Using the Smoothed Hessian Major Axis Filament Finder, we identify 259 filaments at high Galactic latitude, with lengths larger or equal to 2° (corresponding to 3.5 pc in length for a typical distance of 100 pc). These filaments show a preferred orientation parallel to the magnetic field projected onto the plane of the sky, derived from their polarization angles. We present mean maps of the filaments in Stokes I, Q, U, E, and B, computed by stacking individual images rotated to align the orientations of the filaments. Combining the stacked images and the histogram of relative orientations, we estimate the mean polarization fraction of the filaments to be 11 %. Furthermore, we show that the correlation between the filaments and the magnetic field orientations may account for the E and B asymmetry and the $C_{\ell}^{TE}/C_{\ell}^{EE}$ ratio, reported in the power spectra analysis of the Planck 353 GHz polarization maps. Future models of the dust foreground for CMB polarization studies will need to take into account the observed correlation between the dust polarization and the structure of interstellar matter. △ Less

Submitted 5 October, 2015; v1 submitted 11 May, 2015; originally announced May 2015.

Comments: 18 pages, 15 figures. A&A accepted. Corresponding author: T. Ghosh

Journal ref: A&A 586, A141 (2016)

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

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, M. I. R. Alves, M. Arnaud, D. Arzoumanian, M. Ashdown, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, N. Bartolo, E. Battaner, K. Benabed, A. Benoît, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger , et al. (173 additional authors not shown)

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

Submitted 8 September, 2015; v1 submitted 13 February, 2015; originally announced February 2015.

Comments: 28 pages, 20 figures. A&A accepted. Corresponding author: J.D. Soler

Report number: AA/2015/25896

Journal ref: A&A, 586 (2016) A138

Planck intermediate results. XXXIV. The magnetic field structure in the Rosette Nebula

Authors: Planck Collaboration, N. Aghanim, M. I. R. Alves, M. Arnaud, D. Arzoumanian, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, N. Bartolo, E. Battaner, K. Benabed, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, A. Bonaldi, L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, A. Bracco, C. Burigana, E. Calabrese , et al. (175 additional authors not shown)

Abstract: Planck has mapped the polarized dust emission over the whole sky, making it possible to trace the Galactic magnetic field structure that pervades the interstellar medium (ISM). We combine polarization data from Planck with rotation measure (RM) observations towards a massive star-forming region, the Rosette Nebula in the Monoceros molecular cloud, to study its magnetic field structure and the impa… ▽ More Planck has mapped the polarized dust emission over the whole sky, making it possible to trace the Galactic magnetic field structure that pervades the interstellar medium (ISM). We combine polarization data from Planck with rotation measure (RM) observations towards a massive star-forming region, the Rosette Nebula in the Monoceros molecular cloud, to study its magnetic field structure and the impact of an expanding HII region on the morphology of the field. We derive an analytical solution for the magnetic field, assumed to evolve from an initially uniform configuration following the expansion of ionized gas and the formation of a shell of swept-up ISM. From the RM data we estimate a mean value of the line-of-sight component of the magnetic field of about 3microG (towards the observer) in the Rosette Nebula, for a uniform electron density of about 12cm-3. The dust shell that surrounds the Rosette HII region is clearly observed in the Planck intensity map at 353 GHz, with a polarization signal significantly different from that of the local background when considered as a whole. The Planck observations constrain the plane-of-the-sky orientation of the magnetic field in the Rosette's parent molecular cloud to be mostly aligned with the large-scale field along the Galactic plane. The Planck data are compared with the analytical model, which predicts the mean polarization properties of a spherical and uniform dust shell for a given orientation of the field. This comparison leads to an upper limit of about 45degr on the angle between the line of sight and the magnetic field in the Rosette complex, for an assumed intrinsic dust polarization fraction of 4%. This field direction can reproduce the RM values detected in the ionized region if the magnetic field strength in the Monoceros molecular cloud is in the range 6.5--9microG. △ Less

Submitted 15 September, 2015; v1 submitted 5 January, 2015; originally announced January 2015.

Comments: 16 pages, 11 figures. A&A accepted. Corresponding author: M.I.R. Alves

Journal ref: A&A 586, A137 (2016)

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

Authors: A. P. Jones, L. Fanciullo, M. Koehler, L. Verstraete, V. Guillet, M. Bocchio, N. Ysard

Abstract: 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 studie… ▽ More 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. △ Less

Submitted 23 November, 2014; originally announced November 2014.

Comments: 22 pages, 20 figures

Journal ref: published; A&A 558, A62 (2013)

The cycling of carbon into and out of dust

Authors: Anthony P. Jones, Nathalie Ysard, Melanie Koehler, Lapo Fanciullo, Marco Bocchio, Elisabetta Micelotta, Laurent Verstraete, Vincent Guillet

Abstract: 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 hydroc… ▽ More 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. △ Less

Submitted 21 November, 2014; originally announced November 2014.

Comments: 14 pages, 6 figures

Journal ref: Faraday Discussions (2014), Volume 168, Astrochemistry of Dust, Ice and Gas, p. 313

Planck intermediate results. XXXIII. Signature of the magnetic field geometry of interstellar filaments in dust polarization maps

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, M. I. R. Alves, M. Arnaud, D. Arzoumanian, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, N. Bartolo, E. Battaner, K. Benabed, A. Benoit-Lévy, J. -P. Bernard, O. Berné, M. Bersanelli, P. Bielewicz, A. Bonaldi, L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, A. Bracco , et al. (165 additional authors not shown)

Abstract: Planck observations at 353GHz provide the first fully-sampled maps of the polarized dust emission towards interstellar filaments and their backgrounds. The polarization data provide insight on the structure of their magnetic field (B). We present the polarization maps of three nearby star forming filament of moderate column density (NH~10^22cm^-2): Musca, B211, and L1506. We use the spatial inform… ▽ More Planck observations at 353GHz provide the first fully-sampled maps of the polarized dust emission towards interstellar filaments and their backgrounds. The polarization data provide insight on the structure of their magnetic field (B). We present the polarization maps of three nearby star forming filament of moderate column density (NH~10^22cm^-2): Musca, B211, and L1506. We use the spatial information to separate Stokes I, Q, and U of the filaments from those of their backgrounds, an essential step to measure the intrinsic polarization fraction (p) and angle (psi) of each emission component. We find that the polarization angles in the three filaments (psi_fil) are coherent along their lengths and not the same as in their backgrounds (psi_bg). The differences between psi_fil and psi_bg are 12deg, 6deg, and 54deg for Musca, B211, and L1506, respectively. These differences for Musca and L1506 are larger than the dispersions of psi, both along the filaments and in their backgrounds. The observed changes of psi are direct evidence for variations of the orientation of the plane of the sky (POS) projection of the B-field. As in previous studies, we find a decrease of several percent of p with NH. We show that the drop in p cannot be explained by random fluctuations of the orientation of B within the filaments because they are too small (sigma_psi<10deg). We recognize the degeneracy between dust alignment efficiency and the structure of B in causing variations in p, but we argue that the decrease of p from the backgrounds to the filaments results in part from depolarization associated with the 3D structure of B: both its orientation in the POS and with respect to the POS. We do not resolve the inner structure of the filaments, but at the smallest scales accessible with Planck (~0.2pc), the observed changes of psi and p hold information on the B-field structure within filaments. △ Less

Submitted 1 July, 2015; v1 submitted 9 November, 2014; originally announced November 2014.

Comments: 16 pages, 14 figures, Submitted to Astronomy & Astrophysics

Journal ref: A&A 586, A136 (2016)

Planck intermediate results. XXXII. The relative orientation between the magnetic field and structures traced by interstellar dust

Authors: Planck Collaboration, R. Adam, P. A. R. Ade, N. Aghanim, M. I. R. Alves, M. Arnaud, D. Arzoumanian, M. Ashdown, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, N. Bartolo, E. Battaner, K. Benabed, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, A. Bonaldi, L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger , et al. (177 additional authors not shown)

Abstract: The role of the magnetic field in the formation of the filamentary structures observed in the interstellar medium (ISM) is a debated topic. The Planck all-sky maps of linearly polarized emission from dust at 353GHz provide the required combination of imaging and statistics to study the correlation between the structures of the Galactic magnetic field and of interstellar matter, both in the diffuse… ▽ More The role of the magnetic field in the formation of the filamentary structures observed in the interstellar medium (ISM) is a debated topic. The Planck all-sky maps of linearly polarized emission from dust at 353GHz provide the required combination of imaging and statistics to study the correlation between the structures of the Galactic magnetic field and of interstellar matter, both in the diffuse ISM and in molecular clouds. The data reveal structures, or ridges, in the intensity map with counterparts in the Stokes Q and/or U maps. We focus on structures at intermediate and high Galactic latitudes with column density from $10^{20}$ to $10^{22}$ cm$^{-2}$. We measure the magnetic field orientation on the plane of the sky from the polarization data, and present an algorithm to estimate the orientation of the ridges from the dust intensity map. We use analytical models to account for projection effects. Comparing polarization angles on and off the structures, we estimate the mean ratio between the strengths of the turbulent and mean components of the magnetic field to be between 0.6 and 1.0, with a preferred value of 0.8. We find that the ridges are preferentially aligned with the magnetic field measured on the structures. This trend becomes more striking for increasing polarization fraction and decreasing column density. We interpret the increase of alignment with polarization fraction as a consequence of projections effects. The decrease of alignment for high column density is not due to a loss of correlation between the structures and the geometry of the magnetic field. In molecular complexes, we observe structures perpendicular to the magnetic field, which cannot be accounted for by projection effects. We discuss our results in the context of models and MHD simulations, which describe the formation of structures in the magnetized ISM. △ Less

Submitted 31 March, 2015; v1 submitted 23 September, 2014; originally announced September 2014.

Comments: Submitted to Astronomy & Astrophysics

Journal ref: A&A 586, A135 (2016)

Planck intermediate results. XXX. The angular power spectrum of polarized dust emission at intermediate and high Galactic latitudes

Authors: Planck Collaboration, R. Adam, P. A. R. Ade, N. Aghanim, M. Arnaud, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, J. G. Bartlett, N. Bartolo, E. Battaner, K. Benabed, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, A. Bonaldi, L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, A. Bracco, M. Bucher , et al. (208 additional authors not shown)

Abstract: The polarized thermal emission from Galactic dust is the main foreground present in measurements of the polarization of the cosmic microwave background (CMB) at frequencies above 100GHz. We exploit the Planck HFI polarization data from 100 to 353GHz to measure the dust angular power spectra $C_\ell^{EE,BB}$ over the range $40<\ell<600$ well away from the Galactic plane. These will bring new insigh… ▽ More The polarized thermal emission from Galactic dust is the main foreground present in measurements of the polarization of the cosmic microwave background (CMB) at frequencies above 100GHz. We exploit the Planck HFI polarization data from 100 to 353GHz to measure the dust angular power spectra $C_\ell^{EE,BB}$ over the range $40<\ell<600$ well away from the Galactic plane. These will bring new insights into interstellar dust physics and a precise determination of the level of contamination for CMB polarization experiments. We show that statistical properties of the emission can be characterized over large fractions of the sky using $C_\ell$. For the dust, they are well described by power laws in $\ell$ with exponents $α^{EE,BB}=-2.42\pm0.02$. The amplitudes of the polarization $C_\ell$ vary with the average brightness in a way similar to the intensity ones. The dust polarization frequency dependence is consistent with modified blackbody emission with $β_d=1.59$ and $T_d=19.6$K. We find a systematic ratio between the amplitudes of the Galactic $B$- and $E$-modes of 0.5. We show that even in the faintest dust-emitting regions there are no "clean" windows where primordial CMB $B$-mode polarization could be measured without subtraction of dust emission. Finally, we investigate the level of dust polarization in the BICEP2 experiment field. Extrapolation of the Planck 353GHz data to 150GHz gives a dust power $\ell(\ell+1)C_\ell^{BB}/(2π)$ of $1.32\times10^{-2}μ$K$_{CMB}^2$ over the $40<\ell<120$ range; the statistical uncertainty is $\pm0.29$ and there is an additional uncertainty (+0.28,-0.24) from the extrapolation, both in the same units. This is the same magnitude as reported by BICEP2 over this $\ell$ range, which highlights the need for assessment of the polarized dust signal even in the cleanest windows of the sky. △ Less

Submitted 8 December, 2014; v1 submitted 19 September, 2014; originally announced September 2014.

Comments: 25 pages, 25 figures, accepted in A&A

Journal ref: A&A 586, A133 (2016)

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

Authors: Plank Collaboration, P. A. R. Ade, N. Aghanim, M. I. R. Alves, G. Aniano, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, N. Bartolo, E. Battaner, K. Benabed, A. Benoit-Levy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, A. Bonaldi, L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, C. Burigana , et al. (173 additional authors not shown)

Abstract: 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 d… ▽ More 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. △ Less

Submitted 15 February, 2016; v1 submitted 8 September, 2014; originally announced September 2014.

Comments: Final version that has appeared in A&A

Journal ref: A&A 586, A132 (2016)

Shocks in dense clouds. IV. Effects of grain-grain processing on molecular line emission

Authors: S. Anderl, V. Guillet, G. Pineau des Forêts, D. R. Flower

Abstract: Grain-grain processing has been shown to be an indispensable ingredient of shock modelling in high density environments. For densities higher than \sim10^5 cm-3, shattering becomes a self-enhanced process that imposes severe chemical and dynamical consequences on the shock characteristics. Shattering is accompanied by the vaporization of grains, which can directly release SiO to the gas phase. Giv… ▽ More Grain-grain processing has been shown to be an indispensable ingredient of shock modelling in high density environments. For densities higher than \sim10^5 cm-3, shattering becomes a self-enhanced process that imposes severe chemical and dynamical consequences on the shock characteristics. Shattering is accompanied by the vaporization of grains, which can directly release SiO to the gas phase. Given that SiO rotational line radiation is used as a major tracer of shocks in dense clouds, it is crucial to understand the influence of vaporization on SiO line emission. We have developed a recipe for implementing the effects of shattering and vaporization into a 2-fluid shock model, resulting in a reduction of computation time by a factor \sim100 compared to a multi-fluid modelling approach. This implementation was combined with an LVG-based modelling of molecular line radiation transport. Using this model we calculated grids of shock models to explore the consequences of different dust-processing scenarios. Grain-grain processing is shown to have a strong influence on C-type shocks for a broad range of magnetic fields: they become hotter and thinner. The reduction in column density of shocked gas lowers the intensity of molecular lines, at the same time as higher peak temperatures increase the intensity of highly excited transitions compared to shocks without grain-grain processing. For OH the net effect is an increase in line intensities, while for CO and H2O it is the contrary. The intensity of H2 emission is decreased in low transitions and increased for highly excited lines. For all molecules, the highly excited lines become sensitive to the value of the magnetic field. Although vaporization increases the intensity of SiO rotational lines, this effect is weakened by the reduced shock width. The release of SiO early in the hot shock changes the excitation characteristics of SiO radiation. △ Less

Submitted 1 August, 2014; originally announced August 2014.

Comments: Published in Astronomy and Astrophysics (2013). 26 pages, 16 figures, 14 tables

Journal ref: A&A 556, A69 (2013)

Polarization measurements analysis II. Best estimators of polarization fraction and angle

Authors: L. Montier, S. Plaszczynski, F. Levrier, M. Tristram, D. Alina, I. Ristorcelli, J. -P. Bernard, V. Guillet

Abstract: With the forthcoming release of high precision polarization measurements, such as from the Planck satellite, it becomes critical to evaluate the performance of estimators for the polarization fraction and angle. These two physical quantities suffer from a well-known bias in the presence of measurement noise, as has been described in part I of this series. In this paper, part II of the series, we e… ▽ More With the forthcoming release of high precision polarization measurements, such as from the Planck satellite, it becomes critical to evaluate the performance of estimators for the polarization fraction and angle. These two physical quantities suffer from a well-known bias in the presence of measurement noise, as has been described in part I of this series. In this paper, part II of the series, we explore the extent to which various estimators may correct the bias. Traditional frequentist estimators of the polarization fraction are compared with two recent estimators: one inspired by a Bayesian analysis and a second following an asymptotic method. We investigate the sensitivity of these estimators to the asymmetry of the covariance matrix which may vary over large datasets. We present for the first time a comparison among polarization angle estimators, and evaluate the statistical bias on the angle that appears when the covariance matrix exhibits effective ellipticity. We also address the question of the accuracy of the polarization fraction and angle uncertainty estimators. The methods linked to the credible intervals and to the variance estimates are tested against the robust confidence interval method. From this pool of estimators, we build recipes adapted to different use-cases: build a mask, compute large maps, and deal with low S/N data. More generally, we show that the traditional estimators suffer from discontinuous distributions at low S/N, while the asymptotic and Bayesian methods do not. Attention is given to the shape of the output distribution of the estimators, and is compared with a Gaussian. In this regard, the new asymptotic method presents the best performance, while the Bayesian output distribution is shown to be strongly asymmetric with a sharp cut at low S/N.Finally, we present an optimization of the estimator derived from the Bayesian analysis using adapted priors. △ Less

Submitted 1 July, 2014; originally announced July 2014.

Journal ref: A&A 574, A136 (2015)

Planck intermediate results. XXII. Frequency dependence of thermal emission from Galactic dust in intensity and polarization

Authors: Planck Collaboration, P. A. R. Ade, M. I. R. Alves, G. Aniano, C. Armitage-Caplan, M. Arnaud, F. Atrio-Barandela, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, E. Battaner, K. Benabed, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, C. Burigana, J. -F. Cardoso, A. Catalano , et al. (155 additional authors not shown)

Abstract: Planck has mapped the intensity and polarization of the sky at microwave frequencies with unprecedented sensitivity. We make use of the Planck 353 GHz I, Q, and U Stokes maps as dust templates, and cross-correlate them with the Planck and WMAP data at 12 frequencies from 23 to 353 GHz, over circular patches with 10 degree radius. The cross-correlation analysis is performed for both intensity and p… ▽ More Planck has mapped the intensity and polarization of the sky at microwave frequencies with unprecedented sensitivity. We make use of the Planck 353 GHz I, Q, and U Stokes maps as dust templates, and cross-correlate them with the Planck and WMAP data at 12 frequencies from 23 to 353 GHz, over circular patches with 10 degree radius. The cross-correlation analysis is performed for both intensity and polarization data in a consistent manner. We use a mask that focuses our analysis on the diffuse interstellar medium at intermediate Galactic latitudes. We determine the spectral indices of dust emission in intensity and polarization between 100 and 353 GHz, for each sky-patch. The mean values, $1.59\pm0.02$ for polarization and $1.51\pm0.01$ for intensity, for a mean dust temperature of 19.6 K, are close, but significantly different ($3.6\,σ$). We determine the mean spectral energy distribution (SED) of the microwave emission, correlated with the 353 GHz dust templates, by averaging the results of the correlation over all sky-patches. We find that the mean SED increases for decreasing frequencies at $ν< 60$ GHz, for both intensity and polarization. The rise of the polarization SED towards low frequencies may be accounted for by a synchrotron component correlated with dust, with no need for any polarization of the anomalous microwave emission. We use a spectral model to separate the synchrotron and dust polarization and to characterize the spectral dependence of the dust polarization fraction. The polarization fraction ($p$) of the dust emission decreases by $(21\pm6)$ % from 353 to 70 GHz. The decrease of $p$ could indicate differences in polarization efficiency among components of interstellar dust (e.g., carbon versus silicate grains). Our observational results provide inputs to quantify and optimize the separation between Galactic and cosmological polarization. △ Less

Submitted 11 May, 2015; v1 submitted 5 May, 2014; originally announced May 2014.

Comments: 26 pages, 15 figures, 9 tables; updated version; small changes to match the published version

Journal ref: A&A 576, A107 (2015)

Planck intermediate results. XXI. Comparison of polarized thermal emission from Galactic dust at 353 GHz with interstellar polarization in the visible

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, D. Alina, G. Aniano, C. Armitage-Caplan, M. Arnaud, M. Ashdown, F. Atrio-Barandela, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, E. Battaner, C. Beichman, K. Benabed, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger , et al. (170 additional authors not shown)

Abstract: The Planck survey provides unprecedented full-sky coverage of the submillimetre polarized emission from Galactic dust, bringing new constraints on the properties of dust. The dust grains that emit the radiation seen by Planck in the submillimetre also extinguish and polarize starlight in the visible. Comparison of the polarization of the emission and of the interstellar polarization on selected li… ▽ More The Planck survey provides unprecedented full-sky coverage of the submillimetre polarized emission from Galactic dust, bringing new constraints on the properties of dust. The dust grains that emit the radiation seen by Planck in the submillimetre also extinguish and polarize starlight in the visible. Comparison of the polarization of the emission and of the interstellar polarization on selected lines of sight probed by stars provides unique new diagnostics of the emission and light scattering properties of dust. Using ancillary catalogues of interstellar polarization and extinction of starlight, we obtain the degree of polarization, $p_V$ , and the optical depth in the V band to the star, $τ_V$. Toward these stars we measure the submillimetre polarized intensity, $P_S$, and total intensity, $I_S$, in the Planck 353 GHz channel. For those lines of sight through the diffuse interstellar medium with comparable values of the estimated column density and polarization directions close to orthogonal, we correlate properties in the submillimetre and visible to find two ratios, $R_{S/V} = (P_S/I_S)/(p_V/τ_V)$ and $R_{P/p} = P_S/p_V$ , the latter focusing directly on the polarization properties of the aligned grain population alone. We find $R_{S/V}$ = 4.2, with statistical and systematic uncertainties 0.2 and 0.3, respectively, and $R_{P/p}$ = 5.4 MJy sr$^{-1}$, with uncertainties 0.2 and 0.3 MJy sr$^{-1}$, respectively. Our estimate of $R_{S/V}$ is compatible with predictions based on a range of polarizing dust models that have been developed for the diffuse interstellar medium. However, our estimate of $R_{P/p}$ is not compatible with predictions, which are too low by a factor of about 2.5. This more discriminating diagnostic, $R_{P/p}$, indicates that changes to the optical properties in the models of the aligned grain population are required. △ Less

Submitted 12 March, 2015; v1 submitted 5 May, 2014; originally announced May 2014.

Journal ref: A&A 576, A106 (2015)

Planck intermediate results. XX. Comparison of polarized thermal emission from Galactic dust with simulations of MHD turbulence

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, D. Alina, M. I. R. Alves, G. Aniano, C. Armitage-Caplan, M. Arnaud, D. Arzoumanian, M. Ashdown, F. Atrio-Barandela, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, E. Battaner, K. Benabed, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger , et al. (164 additional authors not shown)

Abstract: Polarized emission observed by Planck HFI at 353 GHz towards a sample of nearby fields is presented, focusing on the statistics of polarization fractions $p$ and angles $ψ$. The polarization fractions and column densities in these nearby fields are representative of the range of values obtained over the whole sky. We find that: (i) the largest polarization fractions are reached in the most diffuse… ▽ More Polarized emission observed by Planck HFI at 353 GHz towards a sample of nearby fields is presented, focusing on the statistics of polarization fractions $p$ and angles $ψ$. The polarization fractions and column densities in these nearby fields are representative of the range of values obtained over the whole sky. We find that: (i) the largest polarization fractions are reached in the most diffuse fields; (ii) the maximum polarization fraction $p_\mathrm{max}$ decreases with column density $N_\mathrm{H}$ in the more opaque fields with $N_\mathrm{H} > 10^{21}\,\mathrm{cm}^{-2}$; and (iii) the polarization fraction along a given line of sight is correlated with the local spatial coherence of the polarization angle. These observations are compared to polarized emission maps computed in simulations of anisotropic magnetohydrodynamical (MHD) turbulence in which we assume a uniform intrinsic polarization fraction of the dust grains. We find that an estimate of this parameter may be recovered from the maximum polarization fraction $p_\mathrm{max}$ in diffuse regions where the magnetic field is ordered on large scales and perpendicular to the line of sight. This emphasizes the impact of anisotropies of the magnetic field on the emerging polarization signal. The decrease of the polarization fraction with column density in nearby molecular clouds is well reproduced in the simulations, indicating that it is essentially due to the turbulent structure of the magnetic field: an accumulation of variously polarized structures along the line of sight leads to such an anti-correlation. In the simulations, polarization fractions are also found to anti-correlate with the angle dispersion function $\mathcal{S}$. [abridged] △ Less

Submitted 25 September, 2014; v1 submitted 5 May, 2014; originally announced May 2014.

Comments: Part of a series of four papers, accepted for publication in A&A

Journal ref: A&A 576, A105 (2015)

Planck intermediate results. XIX. An overview of the polarized thermal emission from Galactic dust

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, D. Alina, M. I. R. Alves, C. Armitage-Caplan, M. Arnaud, D. Arzoumanian, M. Ashdown, F. Atrio-Barandela, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, E. Battaner, K. Benabed, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger , et al. (177 additional authors not shown)

Abstract: This paper presents the large-scale polarized sky as seen by Planck HFI at 353 GHz, which is the most sensitive Planck channel for dust polarization. We construct and analyse large-scale maps of dust polarization fraction and polarization direction, while taking account of noise bias and possible systematic effects. We find that the maximum observed dust polarization fraction is high (pmax > 18%),… ▽ More This paper presents the large-scale polarized sky as seen by Planck HFI at 353 GHz, which is the most sensitive Planck channel for dust polarization. We construct and analyse large-scale maps of dust polarization fraction and polarization direction, while taking account of noise bias and possible systematic effects. We find that the maximum observed dust polarization fraction is high (pmax > 18%), in particular in some of the intermediate dust column density (AV < 1mag) regions. There is a systematic decrease in the dust polarization fraction with increasing dust column density, and we interpret the features of this correlation in light of both radiative grain alignment predictions and fluctuations in the magnetic field orientation. We also characterize the spatial structure of the polarization angle using the angle dispersion function and find that, in nearby fields at intermediate latitudes, the polarization angle is ordered over extended areas that are separated by filamentary structures, which appear as interfaces where the magnetic field sky projection rotates abruptly without apparent variations in the dust column density. The polarization fraction is found to be anti-correlated with the dispersion of the polarization angle, implying that the variations are likely due to fluctuations in the 3D magnetic field orientation along the line of sight sampling the diffuse interstellar medium.We also compare the dust emission with the polarized synchrotron emission measured with the Planck LFI, with low-frequency radio data, and with Faraday rotation measurements of extragalactic sources. The two polarized components are globally similar in structure along the plane and notably in the Fan and North Polar Spur regions. A detailed comparison of these three tracers shows, however, that dust and cosmic rays generally sample different parts of the line of sight and confirms that much of the variation observed in the Planck data is due to the 3D structure of the magnetic field. △ Less

Submitted 5 May, 2014; originally announced May 2014.

Journal ref: A&A 576, A104 (2015)

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

Authors: Planck Collaboration, A. Abergel, P. A. R. Ade, N. Aghanim, M. I. R. Alves, G. Aniano, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, J. G. Bartlett, E. Battaner, K. Benabed, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. Bobin, A. Bonaldi, J. R. Bond, F. R. Bouchet, F. Boulanger, C. Burigana , et al. (174 additional authors not shown)

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

Submitted 7 May, 2014; v1 submitted 19 December, 2013; originally announced December 2013.

Comments: Final Version to appear in A&A

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

Authors: Planck Collaboration, A. Abergel, P. A. R. Ade, N. Aghanim, M. I. R. Alves, G. Aniano, C. Armitage-Caplan, M. Arnaud, M. Ashdown, F. Atrio-Barandela, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, J. G. Bartlett, E. Battaner, K. Benabed, A. Benoît, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. Bobin, J. J. Bock, A. Bonaldi , et al. (222 additional authors not shown)

Abstract: 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 wide… ▽ More 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. △ Less

Submitted 23 September, 2014; v1 submitted 4 December, 2013; originally announced December 2013.

Comments: 37 pages. accepted version. This paper is one of a set associated with the 2013 data release from the Planck mission

Journal ref: A&A 571, A11 (2014)

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

Authors: N. Ysard, M. Juvela, K. Demyk, V. Guillet, A. Abergel, J. -P. Bernard, J. Malinen, C. Mény, L. Montier, D. Paradis, I. Ristorcelli, L. Verstraete

Abstract: 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 dus… ▽ More 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. △ Less

Submitted 25 April, 2012; v1 submitted 27 February, 2012; originally announced February 2012.

Comments: accepted for publication by A&A, abstract modified to match arXiv format, version 2 => correction of numerical values throughout the text/figures

Planck Early Results: Thermal dust in Nearby Molecular Clouds

Authors: Planck Collaboration, A. Abergel, P. A. R. Ade, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. Balbi, A. J. Banday, R. B. Barreiro, J. G. Bartlett, E. Battaner, K. Benabed, A. Benoît, J. -P. Bernard, M. Bersanelli, R. Bhatia, J. J. Bock, A. Bonaldi, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher , et al. (175 additional authors not shown)

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

Submitted 25 August, 2011; v1 submitted 11 January, 2011; originally announced January 2011.

Comments: Planck early paper accepted for publication in Astronomy & Astrophysics, 18 pages