Search

The polarized dust emission observed in Class 0 protostellar cores at high angular resolution with ALMA has raised several concerns about the grain alignment conditions in these regions. We aim to study the role of the radiation field on the grain alignment mechanisms occurring in the interior (<1000 au) of Class 0 protostars. We produce synthetic observations of the polarized dust emission from a MHD model of protostellar formation, using the POLARIS dust radiative transfer tool, which includes dust alignment with Radiative Torques Alignment (RATs). We test how the polarized dust emission from the model core depends on the irradiation conditions in the protostellar envelope, by varying the radiation due to accretion luminosity propagating from the central protostellar embryo throughout the envelope. The level of grain alignment efficiency obtained in the radiative transfer models is then compared to (sub-) millimeter ALMA dust polarization observations of Class 0 protostars. Our radiative transfer calculations have a central irradiation that reproduces the protostellar luminosities typically observed towards low- to intermediate-mass protostars, as well as super-paramagnetic grains, and grains >10 micron, which are required to bring the dust grain alignment efficiencies of the synthetic observations up to observed levels. Our radiative transfer calculations show that irradiation plays an important role in the mechanisms that dictate the size range of aligned grains in Class 0 protostars. Regions of the envelope that are preferentially irradiated harbor strong polarized dust emission but can be affected by the rotational disruption of dust grains. Episodes of high luminosity could affect grain alignment and trigger grain disruption mechanisms. [abridged]

Aims: The large abundances of CH$^+$ in the diffuse interstellar medium (ISM) are a long standing issue of our understanding of the thermodynamical and chemical states of the gas. We investigate, here, the formation of CH+ in turbulent and multiphase environments, where the heating of the gas is almost solely driven by the photoelectric effect. Methods: The diffuse ISM is simulated using the magnetohydrodynamic (MHD) code RAMSES which self-consistently computes the dynamical and thermal evolution of the gas along with the time-dependent evolutions of the abundances of H$^+$, H, and H$_2$. The rest of the chemistry, including the abundance of CH$^+$, is computed in post-processing, at equilibrium, under the constraint of out-ofequilibrium of H$^+$, H, and H$_2$. The comparison with the observations is performed taking into account an often neglected, yet paramount, piece of information, namely the length of the intercepted diffuse matter along the observed lines of sight. Results: The quasi totality of the mass of CH$^+$ originates from the unstable gas, in environments where the kinetic temperature is larger than 600 K, the density ranges between 0.6 and 10 cm$^{-3}$, the electronic fraction ranges between 3 x 10$^{-4}$ and 6 x 10$^{-3}$, and the molecular fraction is smaller than 0.4. Its formation is driven by warm and out-of-equilibrium H$_2$ initially formed in the cold neutral medium (CNM) and injected in more diffuse environments and even the warm neutral medium (WNM) through a combination of advection and thermal instability. The simulation which displays the tightest agreement with the HI-to-H$_2$ transition and the thermal pressure distribution observed in the Solar Neighborhood is found to naturally reproduce the observed abundances of CH$^+$, the dispersion of observations, the probability of occurrence of most of the lines of sight, the fraction of non-detections of CH$^+$, and the distribution of its line profiles. The amount of CH$^+$ and the statistical properties of the simulated lines of sight are set by the fraction of unstable gas rich in H$_2$ which is controlled, on Galactic scales, by the mean density of the diffuse ISM (or, equivalently, its total mass), the amplitude of the mean UV radiation field, and the strength of the turbulent forcing. Conclusions: This work offers a new and natural solution to an 80 years old chemical riddle. The almost ubiquitous presence of CH$^+$ in the diffuse ISM likely results from the exchanges of matter between the CNM and the WNM induced by the combination of turbulent advection and thermal instability, without the need to invoke ambipolar diffusion or regions of intermittent turbulent dissipation. Through two phase turbulent mixing, CH$^+$ might thus be a tracer of the H$_2$ mass loss rate of CNM clouds., Comment: Astronomy and Astrophysics, EDP Sciences, In press

High resolution (sub-)millimeter polarization observations have opened a new era in the understanding of how B-fields are organized in star forming regions, unveiling an intricate interplay between the B-fields and the gas in protostellar cores. However, to assess the role of the B-field in the process of solar-type star formation, it is key to be able to understand to what extent these polarized dust emissions are good tracers of the B-field in the youngest protostellar objects. We present a thorough investigation of the fidelity and limitations of using dust polarized emission to map the B-field topologies in low-mass protostars. To assess the importance of these effects, we performed the analysis of B-field properties in 27 realizations of MHD models of star-forming cores. Assuming a uniform population of dust grains whose sizes follow the standard MRN, we analyze the synthetic polarized dust emission maps produced if these grains align with the local B-field thanks to B-RATs. We find that (sub-)millimeter polarized dust emission is a robust tracer of the B-field topologies in inner protostellar envelopes and is successful at capturing the details of the B-field spatial distribution down to radii ~100 au. Measurements of the los averaged B-field orientation using the polarized dust emission are precise to < 15{\deg} in about 75 - 95% of the independent lines of sight peering through protostellar envelopes. Large discrepancies between the integrated B-field mean orientation and the orientation reconstructed from the polarized dust emission are mostly observed in (i) lines of sight where the B-field is highly disorganized and (ii) lines of sight probing large column densities. Our analysis shows that high opacity of the thermal dust emission and low polarization fractions could be used to avoid utilizing the small fraction of measurements affected by large errors., Comment: Accepted for publication in A&A

The Davis-Chandrasekhar-Fermi (DCF) method is widely used to indirectly estimate the magnetic field strength from the plane-of-sky field orientation. In this work, we present a set of 3D MHD simulations and synthetic polarization images using radiative transfer of clustered massive star-forming regions. We apply the DCF method on the synthetic polarization maps to investigate its reliability in high-density molecular clumps and dense cores where self-gravity is significant. We investigate the validity of the assumptions of the DCF method step by step and compare the model and estimated field strength to derive the correction factors for the estimated uniform and total (rms) magnetic field strength at clump and core scales. The correction factors in different situations are catalogued. We find the DCF method works well in strong field cases. However, the magnetic field strength in weak field cases could be significantly overestimated by the DCF method when the turbulent magnetic energy is smaller than the turbulent kinetic energy. We investigate the accuracy of the angular dispersion function (ADF, a modified DCF method) method on the effects that may affect the measured angular dispersion and find that the ADF method correctly accounts for the ordered field structure, the beam-smoothing, and the interferometric filtering, but may not be applicable to account for the signal integration along the line of sight in most cases. Our results suggest that the DCF methods should be avoided to be applied below $\sim$0.1 pc scales if the effect of line-of-sight signal integration is not properly addressed., Comment: 40 pages, 24 figures, accepted by ApJ

To assess the potential role of magnetic fields in regulating the envelope rotation and the fragmentation of Class 0 protostars, we carried out observations of the dust polarized emission at 0.87 mm with the SMA, in the envelopes of a large sample of 20 Class 0 protostars. We estimate the mean magnetic field orientation over the central 1000 au envelope scales and compared it to that of the protostellar outflow in order to study the relation between their misalignment and the kinematics of the circumstellar gas. We discover a strong relationship between the misalignment of the magnetic field orientation with the outflow and the amount of angular momentum observed at similar scales in the protostellar envelope, revealing a potential link between the kinetic and the magnetic energy at envelope scales. The relation could be driven by favored B misalignments in more dynamical envelopes or a dependence of the envelope dynamics with the large-scale B initial configuration. Comparing the trend with the presence of fragmentation, we observe that single sources are mostly associated with conditions of low angular momentum in the inner envelope and good alignment of the magnetic field with protostellar outflows, at intermediate scales. Our results suggest that the properties of the magnetic field in protostellar envelopes bear a tight relationship with the rotating-infalling gas directly involved in the star and disk formation: we find that it may not only influence the fragmentation of protostellar cores into multiple stellar systems, but also set the conditions establishing the pristine properties of planet-forming disks., Comment: 21 pages, 7 figures, accepted for publication in A&A

We studied the statistical properties of the HI-to-H$_2$ transition observed in absorption in the local diffuse and multiphase ISM to identify the physical processes controlling the probability of occurrence of any line of sight. The turbulent diffuse ISM is modeled using the RAMSES code, which includes detailed treatments of the magnetohydrodynamics, the thermal evolution of the gas, and the chemistry of H$_2$. The impacts of the UV radiation field, the mean density, the turbulent forcing, the integral scale, the magnetic field, and the gravity on the molecular content of the gas are explored through a parametric study covering a wide range of physical conditions. The statistics of the HI-to-H$_2$ transition are interpreted through analytical prescriptions and compared with the observations using a modified and robust version of the Kolmogorov-Smirnov test. The results of one simulation, convolved with the distribution of distances of the observational sample, are able to explain most of the statistical properties of the HI-to-H$_2$ transition observed in the local ISM. The tightest agreement is obtained for a neutral diffuse gas modeled over ~200 pc, with a mean density of $1-2$ cm$^{-3}$, illuminated by the standard interstellar UV radiation field, and stirred up by a large-scale compressive turbulent forcing. Within this configuration, the 2D probability histogram of the column densities of H and H$_2$ is remarkably stable and is almost unaltered by gravity, the strength of the turbulent forcing, the resolution of the simulation, or the strength of the magnetic field $B_x$. The weak effect of the resolution and our analytical prescription suggest that the column densities of HI are likely built up in large-scale WNM and CNM structures correlated in density over ~20 pc and ~10 pc, respectively, while those of H$_2$ are built up in CNM structures between ~3 pc and ~10 pc.

Analyzing the properties of dust and its evolution in the early phases of star formation is crucial to put constraints on the collapse and accretion processes as well as on the pristine properties of planet-forming seeds. We analyze PdBI observations at 1.3 and 3.2 mm for 12 Class 0 protostars obtained as part of the CALYPSO survey and derive dust emissivity index (beta) profiles as a function of the envelope radius at 200-2000 au scales. Most of the protostellar envelopes show low dust emissivity indices decreasing toward the central regions. The decreasing trend remains after correction of the (potentially optically-thick) central region emission, with surprisingly low beta (lower than 1) values across most of the envelope radii of NGC1333-IRAS4A, NGC1333-IRAS4B, SVS13B, and Serpens-SMM4. We discuss the various processes that could explain such low and varying dust emissivity indices at envelope radii 200-2000 au. Our observations of extremely low dust emissivity indices could trace the presence of large (mm-size) grains in Class 0 envelopes and a radial increase of the dust grain size toward the inner envelope regions. While it is expected that large grains in young protostellar envelopes could be built via grain growth and coagulation, we stress that the typical timescales required to build mm grains in current coagulation models are at odds with the youth of our Class 0 protostars. Additional variations in the dust composition could also partly contribute to the low beta we observe. The steepness of the beta radial gradient depends strongly on the envelope mass, which might favor a scenario where large grains are built in high density protostellar disks and transported to the intermediate envelope radii, for example with the help of outflows and winds., Comment: 17 pages, 11 figures, accepted for publication in A&A

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., Comment: Accepted for publication in MNRAS. Accepted 2019 July 23. Received 2019 July 23; in original form 2019 April 05

Although, from a theoretical point of view, magnetic fields are believed to have a significant role during the early stages of star formation, especially during the main accretion phase, the magnetic field strength and topology is poorly constrained in the youngest accreting Class 0 protostars that lead to the formation of solar-type stars. We carried out observations of the polarized dust continuum emission with the SMA interferometer at 0.87mm, in order to probe the structure of the magnetic field in a sample of 12 low-mass Class 0 envelopes, including both single protostars and multiple systems, in nearby clouds. Our SMA observations probe the envelope emission at scales of 600-5000 au with a spatial resolution ranging from 600 to 1500 au depending on the source distance. We report the detection of linearly polarized dust continuum emission in all of our targets, with average polarization fractions ranging from 2% to 10% in these protostellar envelopes. The polarization fraction decreases with the continuum flux density, which translates into a decrease with the H2 column density within an individual envelope. Our analysis show that the envelope-scale magnetic field is preferentially observed either aligned or perpendicular to the outflow direction. Interestingly, our results suggest for the first time a relation between the orientation of the magnetic field and the rotational energy of envelopes, with a larger occurrence of misalignment in sources where strong rotational motions are detected at hundreds to thousands of au scales. We also show that the best agreement between the magnetic field and outflow orientation is found in sources showing no small-scale multiplicity and no large disks at about 100 au scales., Comment: 22 pages, 9 figures, accepted for publication in A&A

Molecular hydrogen is the most abundant molecule in the universe. It is the first one to form and survive photo-dissociation in tenuous environments. Its formation involves catalytic reactions on the surface of interstellar grains. The micro-physics of the formation process has been investigated intensively in the last 20 years, in parallel of new astrophysical observational and modeling progresses. In the perspectives of the probable revolution brought by the future satellite JWST, this article has been written to present what we think we know about the H$_2$ formation in a variety of interstellar environments., Comment: Accepted for publication in Molecular Astrophysics

Molecular clouds are essentially made up of atomic and molecular hydrogen, which in spite of being the simplest molecule in the ISM plays a key role in the chemical evolution of molecular clouds. Since its formation time is very long, the H2 molecules can be transported by the turbulent motions within the cloud toward low density and warm regions, where its enhanced abundance can boost the abundances of molecules with high endothermicities. We present high resolution simulations where we include the evolution of the molecular gas under the effect of the dynamics, and we analyze its impact on the abundance of CH+., Comment: To appear in: Proceedings IAU Symposium No. 332, 2017

This paper assesses the roles of the presence of warm H2, and the increased formation rate due to the ion-neutral drift. We performed ideal MHD simulations that include the heating and cooling of the multiphase ISM, and where we treat dynamically the formation of H2. In a post-processing step we compute the abundances of species at chemical equilibrium. We show that CH+ is efficiently formed at the edge of clumps, in regions where the H2 fraction is low, but nevertheless higher than its equilibrium value, and where the gas temperature is high. We show that warm and out of equilibrium H2 increases the integrated column densities of CH+ by one order of magnitude, up to values still 3-10 times lower than those observed in the diffuse ISM. We balance the Lorentz force with the ion-neutral drag to estimate the ion-drift velocities (vd). We find that the vd distribution peaks around 0.04 km s-1, and that high vd are too rare to have a significant statistical impact on the abundances of CH+. Compared to previous works, our multiphase simulations reduce the spread in vd, and our self-consistent treatment of the ionisation leads to much reduced vd. Nevertheless, our resolution study shows that this velocity distribution is not converged: the ion-neutral drift has a higher impact on CH+ at higher resolution. On the other hand, our ideal MHD simulations do not include ambipolar diffusion, which would yield lower drift velocities. Within these limitations, we conclude that warm H2 is a key ingredient in the efficient formation of CH+ and that the ambipolar diffusion has very little influence on the abundance of CH+, mainly due to the small drift velocities obtained. However, we point out that small-scale processes and other non-thermal processes not included in our MHD simulation may be of crucial importance, and higher resolution studies with better controlled dissipation processes are needed., Comment: Accepted for publication in section 6. Interstellar and circumstellar matter of Astronomy and Astrophysics. Date of acceptance: 08/12/2016. (13 pages, 12 figures.)

H2 is the simplest and the most abundant molecule in the ISM, and its formation precedes the formation of other molecules. Understanding the dynamical influence of the environment and the interplay between the thermal processes related to the formation and destruction of H2 and the structure of the cloud is mandatory to understand correctly the observations of H2. We perform high resolution MHD colliding flow simulations with the AMR code RAMSES in which the physics of H2 has been included. We compare the simulation results with various observations including the column densities of excited rotational levels. Due to a combination of thermal pressure, ram pressure and gravity, the clouds produced at the converging point of HI streams are highly inhomogeneous. H2 molecules quickly form in relatively dense clumps and spread into the diffuse interclump gas. This in particular leads to the existence of significant abundances of H2 in the diffuse and warm gas that lies in between clumps. Simulations and observations show similar trends, specially for the HI-to-H2 transition. The abundances of excited rotational levels, calculated at equilibrium in the simulations are very similar to the observed abundances inferred from FUSE results. This is a direct consequence of the presence of the H2 enriched diffuse and warm gas. Our simulations show that H2 rapidly forms in the dense clumps and, due to the complex structure of molecular clouds, quickly spreads at lower densities. Consequently a significant fraction of warm H2 exists in the low density gas. This warm H2 leads to column densities of excited rotational levels close to the observed ones likely revealing the complex intermix between the warm and the cold gas in molecular clouds. This suggests that the 2-phase structure of molecular clouds is an essential ingredient to fully understand molecular hydrogen in these objects., Comment: 16 pages, 19 figures. Accepted for publication in A&A

Context. Ultraviolet radiation plays a crucial role in molecular clouds. Radiation and matter are tightly coupled and their interplay influences the physical and chemical properties of gas. In particular, modeling the radiation propagation requires calculating column densities, which can be numerically expensive in high-resolution multidimensional simulations. Aims. Developing fast methods for estimating column densities is mandatory if we are interested in the dynamical influence of the radiative transfer. In particular, we focus on the effect of the UV screening on the dynamics and on the statistical properties of molecular clouds. Methods. We have developed a tree-based method for a fast estimate of column densities, implemented in the adaptive mesh refinement code RAMSES. We performed numerical simulations using this method in order to analyze the influence of the screening on the clump formation. Results. We find that the accuracy for the extinction of the tree-based method is better than 10%, while the relative error for the column density can be much more. We describe the implementation of a method based on precalculating the geometrical terms that noticeably reduces the calculation time. To study the influence of the screening on the statistical properties of molecular clouds we present the probability distribution function (PDF) of gas and the associated temperature per density bin and the mass spectra for different density thresholds. Conclusions. The tree-based method is fast and accurate enough to be used during numerical simulations since no communication is needed between CPUs when using a fully threaded tree. It is then suitable to parallel computing. We show that the screening for far UV radiation mainly affects the dense gas, thereby favoring low temperatures and affecting the fragmentation., Comment: 14 pages, 20 figures