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.

(abridged) The far-UV wavelength range (912-2000A) provides access to atomic and molecular transitions of many species the interstellar medium (ISM), circumgalactic medium (CGM), and intergalactic medium, within phases spanning a wide range of ionization, density, temperature, and molecular gas fraction. Far-UV space telescopes have enabled detailed studies of the ISM in the Milky Way thanks to absorption features appearing in the UV spectra of hot stars and yielding fundamental insights into the composition and physical characteristics of all phases of the ISM along with the processes that influence them. However, we have yet to design a spectrometer able to observe the full UV domain at resolving power R>10^5 with a signal-to-noise ratio SNR>500. Such a resolution is necessary to resolve lines from both the cold molecular hydrogen and the warm metal ions with a turbulent velocity of about 1 km s-1, and to differentiate distinct velocity components. Future UV spectroscopic studies of the Milky Way ISM must revolutionize our understanding of the ISM as a dynamical, unstable, and magnetized medium, and rise to the challenge brought forward by current theories. Another interesting prospect is to transpose the same level of details that has been reached for the Milky Way to the ISM in external galaxies, in particular in metal-poor galaxies, where the ISM chemical composition, physical conditions, and topology change dramatically, with significant consequences on the star-formation properties. Finally, we need to be able to perform statistical analyses of background quasar lines of sight intersecting the CGM of galaxies at various redshifts and to comprehend the role of gas exchanges and flows for galaxy evolution.

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.

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.

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.

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

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.

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.

We performed high-resolution, 3D MHD simulations and we compared to observations of translucent molecular clouds. We show that the observed populations of rotational levels of H2 can arise as a consequence of the multi-phase structure of the ISM.

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.