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 dust polarized emission with the maximum grain size of the dust size distribution at the envelope scale, from $a_\mathrm{max}=1\,$ to $50 \,\hbox{$\mu $m}$. 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 millimetre wavelengths. Considering the current theoretical difficulties to align a large fraction of small ISM-like grains in the conditions typical of protostellar envelopes, our results suggest that grain growth (typically $\gt 10 \,\hbox{$\mu $m}$) might have already significantly progressed at scales 100–1000 au in the youngest objects, observed less than 105 yr 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.
Context. Analyses of Planck data have demonstrated that the grain alignment efficiency is almost constant in the diffuse and translucent interstellar medium (ISM). Aims. We aim to test whether the ...radiative torque (RAT) theory is compatible with these new observational constraints on grain alignment. Methods. We combine a numerical magnetohydrodynamical simulation with a state-of-the-art radiative transfer post-processing code POLARIS which incorporates a physical dust model and the detailed physics of grain alignment by RATs. A dust model based on two distinct power-law-sized distributions of spherical graphite grains and oblate silicate grains was designed to reproduce the mean spectral dependence of extinction and polarization observed in the diffuse ISM. From a simulation of interstellar turbulence obtained with the adaptive-mesh-refinement code RAMSES , we extracted a data cube with physical conditions representative of the diffuse ISM. We post-process the RAMSES cube with POLARIS to compute the grain temperature and alignment efficiency in each cell of the cube. Finally, we simulate synthetic dust emission and polarization observations. Results. In our simulation, the grain alignment efficiency is well-correlated with the gas pressure, but not with the radiative torque intensity. Because of the low dust extinction in our simulation, the magnitude of the radiative torque varies little, decreasing only for column densities larger than 10 22 cm −2 . In comparing our synthetic maps with those obtained assuming a uniform alignment efficiency, we find no systematic difference and very small random differences. The dependencies of the polarization fraction p with the column density N H or with the dispersion in polarization angle S are also similar in both cases. The drop of grain alignment produced by the RAT model in the denser cells of the data cube does not significantly affect the patterns of the synthetic polarization maps, the polarization signal being dominated by the line-of-sight and beam integration of the geometry of the magnetic field. If a star is artificially inserted at the center of the simulation, the polarization fraction is increased everywhere, with no specific pattern around the star. The angle-dependence of the RAT efficiency is not observed in simulated maps and where the magnetic field is artificially set to a uniform configuration in the plane of the sky, it is only seen to be very weak in the optimal configuration. Conclusions. The RAT alignment theory is found to be compatible with the Planck polarization data for the diffuse and translucent ISM in the sense that both uniform alignment and RAT alignment lead to very similar simulated maps. To further test the predictions of the RAT theory in an environment where an important drop of grain alignment is expected, high-resolution polarization observations of dense regions must be confronted with numerical simulations sampling high-column densities ( N H > 10 22 cm −2 ) through dense clouds, given a sufficient statistical basis.
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 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.
Context.
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.
Aims.
We ...study the impact of grain–grain coagulation on the evolution of magnetohydrodynamic (MHD) resistivities during the collapse of a prestellar core.
Methods.
We use a 1D model to follow the evolution of the dust size distribution, out-of-equilibrium ionisation 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.
Results.
Starting with a Mathis-Rumpl-Nordsieck (MRN) size distribution (Mathis et al. 1977, ApJ, 217, 425), 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 have a large effect on the evolution of the Hall and ambipolar diffusion MHD resistivities, which still drop significantly 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
H
~ 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.
Conclusions.
When grain drift velocities induced by ambipolar diffusion are included, dust coagulation happening during the collapse of a prestellar core starting from an initial MRN dust size distribution appears to be efficient enough to increase the MHD resistivities to the values necessary to strongly modify the magnetically regulated formation of a planet-forming disk. A consistent treatment of the competition between fragmentation and coagulation is, however, necessary before reaching firm conclusions.
Context.
The azimuthal polarization patterns observed in some protoplanetary disks by the Atacama Large Millimetre Array (ALMA) at millimeter wavelengths have raised doubts about whether they are ...truly produced by dust grains that are aligned with the magnetic field lines. These conclusions were based on the calculations of dust polarized emission in the Rayleigh regime, that is, for grain sizes that are much smaller than the wavelength. However, the grain size in such disks is typically estimated to be in the range of 0.1−1 mm from independent observations.
Aims.
We study the dust polarization properties of aligned grains in emission in the Mie regime, that is, when the mean grain size approaches the wavelength.
Methods.
By using the
T-MATRIX
and
DustEM
codes, we computed the spectral dependence of the polarization fraction in emission for grains in perfect spinning alignment for various grain size distributions. We restricted our study to weakly-elongated oblate and prolate grains of astrosilicate composition that have a mean size ranging from 10
μ
m to 1 mm.
Results.
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 of
λ
∼ 1 mm. The regime of negative polarization does not exist for grains that are smaller than ∼100
μ
m.
Conclusions.
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 to the conserved ability of dust polarized emission to trace the magnetic field in disks.
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.
The change of physical conditions across the turbulent and magnetized interstellar medium induces a 3D spatial variation of the properties of Galactic polarized emission. The observed signal results ...from the averaging of different spectral energy distributions (SEDs) and polarization angles along and between lines of sight. As a consequence, the total Stokes parameters
Q
and
U
will have different frequency dependencies, both departing from the canonical emission law, 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 and enables 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 includes different distortions for the three polarized power-spectra. This new understanding is crucial for CMB component separation, in which extreme accuracy is required for 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
E
B
signal is not measured accurately, great caution is required regarding the assumptions made to model its spectral behavior, as it may not be inferred from the other dust angular power spectra.
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)
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 characterization 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 upon previous studies by using the newly released
SRoll2
maps that include corrections on residual data systematics and by extending the analysis to regions near the Galactic plane. Our analysis focuses on the lowest multipoles between
ℓ
= 4 and 32, as well as three sky areas with sky fractions of
f
sky
= 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 variations of the dust polarization SED on the sky and along the line of sight. Residuals are detected at the three frequencies for the three sky areas. We show that models based on total-intensity data end up underestimating (by a significant factor) the complexity of dust polarized CMB foreground. Our analysis emphasizes the need to include variations of the 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 the 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, as well as 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, showing that a significant refinement to the dust modeling is necessary to ensure an unbiased detection of the CMB primordial
B
-modes at the precision required by future CMB experiments. Further works would also be required to theoretically model the impact of polarization-angle variations on the
EE
and
BB
power spectra of residual maps.
Context. The sub-millimetre polarisation of dust emission from star-forming clouds carries information on grain properties and on the effects that magnetic fields have on cloud evolution. Aims. Using ...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. Methods. JCMT SCUBA-2/POL-2 observations at 850 μm were combined with Planck 850 μm(353 GHz) data to map polarisation fraction at small and large scales. With previous total intensity SCUBA-2 observations (450 and 850 μm) and Herschel data, the column densities were determined via modified black-body fits and via radiative transfer modelling. Models were constructed to examine how the observed polarisation angles and fractions depend on potential magnetic field geometries and grain alignment processes. Results. POL-2 data show clear changes in the magnetic field orientation. These are not in contradiction with the uniform orientation and almost constant polarisation fraction seen by Planck, because of the difference in the beam sizes and the POL-2 data being affected by spatial filtering. The filament has a peak column density of N(H2) ~ 7 × 1022 cm−2, a minimum dust temperature of T ~ 12 K, and a mass of ~4300 M⊙ for the area N(H2) > 5 × 1021 cm−2. The estimated average value of the dust opacity spectral index is β ~ 1.9. The ratio of sub-millimetre and J-band optical depths is τ (250 μm)∕τ (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 × 1022 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 ~104 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. Conclusions. The G035.39-00.33 filament shows strong signs of dust evolution and the low polarisation fraction is suggestive of a loss of polarised emission from its densest parts.