Search

The first hydrostatic core (FHSC) phase is a brief stage in the protostellar evolution that is difficult to detect. Our goal is to characterize the chemical evolution of gas and dust during the formation of the FHSC. Moreover, we are interested in analyzing, for the first time with 3D magnetohydrodynamic (MHD) simulations, the role of grain growth in its chemistry. We postprocessed $2\times10^{5}$ tracer particles from a $\texttt{RAMSES}$ non-ideal MHD simulation using the codes $\texttt{NAUTILUS}$ and $\texttt{SHARK}$ to follow the chemistry and grain growth throughout the simulation. A great chemical inheritance is seen, as gas-phase abundances of most of the C, O, N, and S reservoirs in the hot corino at the end of the simulation match the ice-phase abundances from the prestellar phase. Additionally, interstellar complex organic molecules (iCOMs) such as methyl formate, acetaldehyde, and formamide are formed during the warm-up process. The typical grain size in the hot corino $(n_{\rm H}>10^{11}\ {\rm cm^{-3}})$ increases forty-fold during the last 30 kyr, with negligible effects on its chemical composition. At moderate densities $(10^{10}

Protostellar disks are the product of angular momentum conservation during the protostellar collapse. Understanding their formation is crucial because they are the birthplace of planets and because their formation is tightly related to star formation. Unfortunately, the initial properties of Class 0 disks and their evolution are still poorly constrained observationally and theoretically. We aim to better understand the mechanisms that set the statistics of disk properties as well as to study their formation in massive protostellar clumps. We also want to provide the community with synthetic disk populations to better interpret young disk observations. We use the ramses code to model star and disk formation in massive protostellar clumps with MHD including the effect of ambipolar diffusion and RT including the stellar radiative feedback. Those simulations, resolved up to the astronomical unit scale, allow to investigate the formation of disk populations. Magnetic fields play a crucial role in disk formation. A weaker initial field leads to larger and massive disks and weakens the stellar radiative feedback by increasing fragmentation. We find that ambipolar diffusion impacts disk and star formation and leads to very different disk magnetic properties. The stellar radiative feedback also have a strong influence, increasing the temperature and reducing fragmentation. Comparing our disk populations with observations reveals that our models with a mass-to-flux ratio of 10 seems to better reproduce observed disk sizes. This also sheds light on a tension between models and observations for the disk masses. The clump properties and physical modeling impact disk populations significantly. The tension between observations and models for disk mass estimates is critical to solve with synthetic observations in future years, in particular for our comprehension of planet formation., Comment: 20 pages, 15 figures, accepted for publication in Astronomy & Astrophysics

Context: Chondrules originate from the reprocessing of dust grains. They are key building blocks of telluric planets, yet their formation, which must happen in strongly localized regions of high temperature, remains poorly understood. Aims: We examine the dust spatial distribution near regions of strong local heating produced by current sheets, as a step toward exploring a potential path for chondrule formation. We further aim to investigate current sheet formation under various conditions in protoplanetary disks in the presence of ambipolar diffusion and Ohmic resistivity and the effect of current sheet morphology on dust dynamics in their vicinity. Methods: We use the RAMSES code including modules for non-ideal magnetohydrodynamics and solution of the dynamics of multiple sizes of dust grains to compute unstratified shearing box simulations of current sheet formation. We investigate, through seven models the effect of the ambipolar diffusion and Ohmic resistivity strength, the initial density, and magnetic field, as well as the resolution and box size. Results: We find that current sheets form in all our models, with typical widths of 0.001-0.01 AU and that strong dust fraction variations occur for millimeter-sized grains. These variations are typically of an order of magnitude and up to two orders of magnitude for the most favorable cases. We also show that the box size and resolution has a strong impact on the current sheet distribution and intensity. Conclusions: The formation of current sheets that can intensely heat their surroundings near strong dynamical dust fraction variations could have important implications for chondrule formation, as it appears likely to happen in regions of large dust fraction., Comment: 11 pages, 8 figures, accepted for publication in A&A

(Abridged) Context. Most massive stars are located in multiple stellar systems. Magnetic fields are believed to be essential in the accretion and ejection processes around single massive protostars. Aims. Our aim is to unveil the influence of magnetic fields in the formation of multiple massive stars, in particular on the fragmentation modes and properties of the multiple protostellar system. Methods. Using RAMSES, we follow the collapse of a massive pre-stellar core with (non-ideal) radiation-(magneto-)hydrodynamics. We choose a setup which promotes multiple stellar system formation. Results. In the purely hydrodynamical models, we always obtain (at least) binary systems. When more than two stars are present, their gravitational interaction triggers mergers until there are two stars left. The following gas accretion increases their orbital separation and hierarchical fragmentation occurs so that both stars host a comparable disk and stellar system which then form similar disks as well. We identify several modes of fragmentation: Toomre-unstable disk fragmentation, arm-arm collision and arm-filament collision. Disks grow in size until they fragment and become truncated as the newly-formed companion gains mass. When including magnetic fields, the picture evolves: the primary disk produces less fragments, arm-filament collision is absent. Magnetic fields reduce the initial orbital separation but do not affect its further evolution, which is mainly driven by gas accretion. With magnetic fields, the growth of individual disks is regulated even in the absence of fragmentation or truncation. Conclusions. Hierarchical fragmentation is seen in unmagnetized and magnetized models. Magnetic fields, including non-ideal effects, are important because they remove certain fragmentation modes and limit the growth of disks, which is otherwise only limited through fragmentation., Comment: 16 pages, 11 figures. Accepted for publication in Astronomy & Astrophysics

(Abridged) Most massive protostars exhibit bipolar outflows. Nonetheless, there is no consensus regarding the mechanism at the origin of these outflows, nor on the cause of the less-frequently observed monopolar outflows. We aim to identify the origin of early massive protostellar outflows, focusing on the combined effects of radiative transfer and magnetic fields in a turbulent medium. We use four state-of-the-art radiation-magnetohydrodynamical simulations following the collapse of massive 100 Msun pre-stellar cores with the Ramses code. Turbulence is taken into account via initial velocity dispersion. We use a hybrid radiative transfer method and include ambipolar diffusion. We find that turbulence delays the launching of outflows, which appear to be mainly driven by magnetohydrodynamical processes. Magnetic tower flow and the magneto-centrifugal acceleration contribute to the acceleration and the former operates on larger volumes than the latter. Our finest resolution, 5 AU, does not allow us to get converged results on magneto-centrifugally accelerated outflows. Radiative acceleration takes place as well, dominates in the star vicinity, enlarges the outflow extent, and has no negative impact on the launching of magnetic outflows (up to M~17 Msun, L~1e5 Lsun). The associated opening angles (20-30 deg when magnetic fields dominate) suggest additional (de-)collimating effects to meet observational constraints. Outflows are launched nearly perpendicular to the disk and are misaligned with the initial core-scale magnetic fields, in agreement with several observational studies. In the most turbulent run, the outflow is monopolar. We conclude that magnetic processes dominate the acceleration of massive protostellar outflows up to ~17 Msun, against radiative processes. Turbulence perturbs the outflow launching and is a possible explanation for monopolar outflows., Comment: Accepted for publication in A&A, 20 pages, 15 figures

(Abridged) Context. Massive stars form in magnetized and turbulent environments, and are often located in stellar clusters. Their accretion mechanism, as well as the origin of their system's stellar multiplicity are poorly understood. Aims. We study the influence of both magnetic fields and turbulence on the accretion mechanism of massive protostars and their multiplicity. Methods. We present a series of four Radiation-MHD simulations of the collapse of a massive magnetized, turbulent core of 100 $M_\odot$ with the AMR code Ramses, including a hybrid radiative transfer method for stellar irradiation and ambipolar diffusion. We vary the Mach and Alfvenic Mach numbers to probe sub- and superalfvenic turbulence as well as sub- and supersonic turbulence regimes. Results. Subalfvenic turbulence leads to single stellar systems while superalfvenic turbulence leads to binary formation from disk fragmentation following spiral arm collision, with mass ratios of 1.1-1.6. In those runs, infalling gas reaches the individual disks via a transient circumbinary structure. Magnetically-regulated, thermally-dominated (plasma beta $\beta>1$), Keplerian disks form in all runs, with sizes 100-200 AU and masses 1-8 $M_\odot$. The disks around primary and secondary sink particles share similar properties. We observe higher accretion rates onto the secondary stars than onto their primary star companion. The primary disk orientation is found to be set by the initial angular momentum carried by turbulence. Conclusions. Small (300 AU) massive protostellar disks as those frequently observed nowadays can only be reproduced so far in the presence of (moderate) magnetic fields with ambipolar diffusion, even in a turbulent medium. The interplay between magnetic fields and turbulence sets the multiplicity of stellar clusters. A plasma beta $\beta>1$ is a good indicator of streamers and disks., Comment: Accepted in A&A. Main body: 20 pages, 17 figures

Gas phase Elemental abundances in Molecular CloudS (GEMS) is an IRAM 30m Large Program designed to estimate the S, C, N, and O depletions and gas ionization degree, X(e-), in a set of star-forming filaments of Taurus, Perseus and Orion. Our immediate goal is to build up a complete database of molecular abundances that can serve as an observational basis for estimating X(e-) and the C, O, N, and S depletions through chemical modeling. We observed and derived the abundances of 14 species (13CO, C18O, HCO+, H13CO+, HC18O+, HCN, H13CN, HNC, HCS+, CS, SO, 34SO, H2S, and OCS) in 244 positions, covering the AV 3 to 100 mag, n(H2) a few 10$^{3}$ to 10$^6$ cm$^{-3}$, and Tk 10 to 30 K ranges in these clouds, avoiding protostars, HII regions, and outflows. A statistical analysis is carried out to identify general trends between different species and with physical parameters. Relations between molecules reveal strong linear correlations which define three different families: (1) 13CO and C18O; (2) H13CO+, HC18O+, H13CN, and HNC; and (3) the S-bearing molecules. The abundances of the CO isotopologs increase with the gas kinetic temperature until TK 15 K. For higher temperatures, the abundance remains constant with a scatter of a factor of 3. The abundances of H13CO+, HC18O+, H13CN, and HNC are well correlated with each other, and all of them decrease with molecular hydrogen density, following the law n(H2)$^{-0.8\pm0.2}$. The abundances of S-bearing species also decrease with n(H2) at a rate of (S-bearing/H)gas n(H2)$^{-0.6\pm0.1}$. The abundances of molecules belonging to groups 2 and 3 do not present any clear trend with gas temperature. At scales of molecular clouds, the C18O abundance is the quantity that better correlates with the cloud mass. We discuss the utility of the 13CO/C18O, HCO+/H13CO+, and H13CO+/H13CN abundance ratios as chemical diagnostics of star formation in external galaxies., Comment: 28 pages, 22 figures. Accepted for publication in A&A

CS is among the most abundant gas-phase S-bearing molecules in cold dark molecular clouds. It is easily observable with several transitions in the millimeter wavelength range, and has been widely used as a tracer of the gas density in the interstellar medium in our Galaxy and external galaxies. Chemical models fail to account for the observed CS abundances when assuming the cosmic value for the elemental abundance of sulfur. The CS+O -> CO + S reaction has been proposed as a relevant CS destruction mechanism at low temperatures, and could explain the discrepancy between models and observations. Its reaction rate has been experimentally measured at temperatures of 150-400 K, but the extrapolation to lower temperatures is doubtful. Here we calculate the CS+O reaction rate at temperatures <150 K which are prevailing in the interstellar medium. We performed ab initio calculations to obtain the three lowest PES of the CS+O system. These PESs are used to study the reaction dynamics, using several methods to eventually calculate the CS+O thermal reaction rates. We compare the results of our theoretical calculations for 150-400 K with those obtained in the laboratory. Our detailed theoretical study on the CS+O reaction, which is in agreement with the experimental data obtained at 150-400 K, demonstrates the reliability of our approach. After a careful analysis at lower temperatures, we find that the rate constant at 10 K is negligible, which is consistent with the extrapolation of experimental data using the Arrhenius expression. We use the updated chemical network to model the sulfur chemistry in TMC1 based on molecular abundances determined from GEMS project observations. In our model, we take into account the expected decrease of the cosmic ray ionization rate along the cloud. The abundance of CS is still overestimated when assuming the cosmic value for the sulfur abundance., Comment: 11 pages, 10 figures, accepted in Astronomy & Astrophysics

(Abridged) The aim of this study is to investigate the chemical evolution from the prestellar phase to the formation of the disk, and to determine the impact that the chemical composition of the cold and dense core has on the final composition of the disk. We performed 3D nonideal magneto-hydrodynamic (MHD) simulations of a dense core collapse using the adaptive-mesh-refinement RAMSES code. For each particle ending in the young rotationally supported disk, we ran chemical simulations with the three-phase gas-grain chemistry code Nautilus. Two different sets of initial abundances, which are characteristic of cold cores, were considered. The final distributions of the abundances of common species were compared to each other, as well as with the initial abundances of the cold core. We find that the spatial distributions of molecules reflect their sensitivity to the temperature distribution. The main carriers of the chemical elements in the disk are usually the same as the ones in the cold core, except for the S-bearing species, where HS is replaced by H$_2$S$_3$, and the P-bearing species, where atomic P leads to the formation of PO, PN, HCP, and CP. However, the abundances of less abundant species change over time. This is especially the case for "large" complex organic molecules (COMs) such as CH$_3$CHO, CH$_3$NH$_2$, CH$_3$OCH$_3$, and HCOOCH$_3$ which see their abundances significantly increase during the collapse. These COMs often present similar abundances in the disk despite significantly different abundances in the cold core. In contrast, the abundances of many radicals decrease with time. A significant number of species still show the same abundances in the cold core and the disk, which indicates efficient formation of these molecules in the cold core. This includes H$_2$O, H$_2$CO, HNCO, and "small" COMs such as CH$_3$OH, CH$_3$CN, and NH$_2$CHO., Comment: Accepted in Astronomy & Astrophysics, 44 pages, 36 figures

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\,{\mu}$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)

Sulphur is one of the most abundant elements in the Universe. Surprisingly, sulphuretted molecules are not as abundant as expected in the interstellar medium, and the identity of the main sulphur reservoir is still an open question. Our goal is to investigate the H$_{2}$S chemistry in dark clouds, as this stable molecule is a potential sulphur reservoir. Using millimeter observations of CS, SO, H$_{2}$S, and their isotopologues, we determine the physical conditions and H$_{2}$S abundances along the cores TMC 1-C, TMC 1-CP, and Barnard 1b. The gas-grain model Nautilus is then used to model the sulphur chemistry and explore the impact of photo-desorption and chemical desorption on the H$_2$S abundance. Our model shows that chemical desorption is the main source of gas-phase H$_2$S in dark cores. The measured H$_{2}$S abundance can only be fitted if we assume that the chemical desorption rate decreases by more than a factor of 10 when $n_{\rm H}>2\times10^{4}$. This change in the desorption rate is consistent with the formation of thick H$_2$O and CO ice mantles on grain surfaces. The observed SO and H$_2$S abundances are in good agreement with our predictions adopting an undepleted value of the sulphur abundance. However, the CS abundance is overestimated by a factor of $5-10$. Along the three cores, atomic S is predicted to be the main sulphur reservoir. We conclude that the gaseous H$_2$S abundance is well reproduced, assuming undepleted sulphur abundance and chemical desorption as the main source of H$_2$S. The behavior of the observed H$_{2}$S abundance suggests a changing desorption efficiency, which would probe the snowline in these cores. Our model, however, overestimates the observed gas-phase CS abundance. Given the uncertainty in the sulphur chemistry, our data are consistent with a cosmic elemental S abundance with an uncertainty of a factor of 10., Comment: 23 pages, 16 figures

Frequency-dependent/hybrid approaches for stellar irradiation are of primary importance in numerical simulations of massive star formation. We seek to compare outflow and accretion mechanisms in star formation simulations. We investigate the accuracy of a hybrid radiative transfer method using the gray M1 closure relation for proto-stellar irradiation and gray flux-limited diffusion (FLD) for photons emitted everywhere else. We have coupled the FLD module of the adaptive-mesh refinement code Ramses with Ramses-RT, which is based on the M1 closure relation. Our hybrid (M1+FLD) method takes an average opacity at the stellar temperature for the M1 module, instead of the local environmental radiation field. We have tested this approach in radiative transfer tests of disks irradiated by a star for three levels of optical thickness and compared the temperature structure with RADMC-3D and MCFOST. We applied it to a radiation-hydrodynamical simulation of massive star formation. Our tests validate our hybrid approach for determining the temperature structure of an irradiated disk in the optically-thin and moderately optically-thick regimes and the most optically-thick test shows the limitation of our approach. The optically-thick setups highlight the ability of the hybrid method to partially capture the self-shielding in the disk while the FLD cannot. The radiative acceleration is 100 times greater with the hybrid method. It consistently leads to about +50% more extended and wider-angle radiative outflows in the massive star formation simulation. We obtain a $17.6 M_\odot$ at $t{\simeq}0.7 \tau_\mathrm{ff}$, while the accretion phase is ongoing. Finally, despite the use of refinement to resolve the radiative cavities, no Rayleigh-Taylor instability appears in our simulations, and we justify their absence by physical arguments based on the entropy gradient. (abridged), Comment: 17 pages, 17 figures, accepted for publication in A&A

On the basis of its low luminosity, its chemical composition, and the absence of a large-scale outflow, the dense core Cha-MMS1 located in the Chamaeleon I molecular cloud was proposed as a first hydrostatic core (FHSC) candidate a decade ago. Our goal is to test this hypothesis by searching for a slow, compact outflow driven by Cha-MMS1 that would match the predictions of MHD simulations for this short phase of star formation. We use the Atacama Large Millimetre/submillimetre Array (ALMA) to map Cha-MMS1 at high angular resolution in CO 3-2 and 13CO 3-2 as well as in continuum emission. We report the detection of a bipolar outflow emanating from the central core, along a (projected) direction roughly parallel to the filament in which Cha-MMS1 is embedded and perpendicular to the large-scale magnetic field. The morphology of the outflow indicates that its axis lies close to the plane of the sky. We measure velocities corrected for inclination of more than 90km/s which is clearly incompatible with the expected properties of a FHSC outflow. Several properties of the outflow are determined and compared to previous studies of Class 0 and Class I protostars. The outflow of Cha-MMS1 has a much smaller momentum force than the outflows of other Class 0 protostars. In addition, we find a dynamical age of 200-3000yr indicating that Cha-MMS1 might be one of the youngest ever observed Class 0 protostars. While the existence of the outflow suggests the presence of a disk, no disk is detected in continuum emission and we derive an upper limit of 55au to its radius. We conclude that Cha-MMS1 has already gone through the FHSC phase and is a young Class 0 protostar, but it has not brought its outflow to full power yet., Comment: Accepted for publication in A&A

The Barnard 1b core shows signatures of being at the earliest stages of low-mass star formation, with two extremely young and deeply embedded protostellar objects. Hence, this core is an ideal target to study the structure and chemistry of the first objects formed in the collapse of prestellar cores. We present ALMA Band 6 spectral line observations at ~0.6'' of angular resolution towards Barnard 1b. We have extracted the spectra towards both protostars, and used a Local Thermodynamic Equilibrium (LTE) model to reproduce the observed line profiles. B1b-S shows rich and complex spectra, with emission from high energy transitions of complex molecules, such as CH3OCOH and CH3CHO, including vibrational level transitions. We have tentatively detected for the first time in this source emission from NH2CN, NH2CHO, CH3CH2OH, CH2OHCHO, CH3CH2OCOH and both aGg' and gGg' conformers of (CH2OH)2. This is the first detection of ethyl formate (CH3CH2OCOH) towards a low-mass star forming region. On the other hand, the spectra of the FHSC candidate B1b-N are free of COMs emission. In order to fit the observed line profiles in B1b-S, we used a source model with two components: an inner hot and compact component (200 K, 0.35'') and an outer and colder one (60 K, 0.6''). The resulting COM abundances in B1b-S range from 1e-13 for NH2CN and NH2CHO, up to 1e-9 for CH3OCOH. Our ALMA Band 6 observations reveal the presence of a compact and hot component in B1b-S, with moderate abundances of complex organics. These results indicate that a hot corino is being formed in this very young Class 0 source., Comment: Accepted for publication in A&A

GEMS is an IRAM 30m Large Program whose aim is determining the elemental depletions and the ionization fraction in a set of prototypical star-forming regions. This paper presents the first results from the prototypical dark cloud TMC 1. Extensive millimeter observations have been carried out with the IRAM 30m telescope (3mm and 2mm) and the 40m Yebes telescope (1.3cm and 7mm) to determine the fractional abundances of CO, HCO+, HCN, CS, SO, HCS+, and N2H+ in three cuts which intersect the dense filament at the well-known positions TMC 1-CP, TMC 1-NH3, and TMC 1-C, covering a visual extinction range from Av~3 to ~20mag. Two phases with differentiated chemistry can be distinguished: i) the translucent envelope with molecular hydrogen densities of (1-5)x10$^3$ cm$^{-3}$; and ii) the dense phase, located at Av>10mag, with molecular hydrogen densities >10$^4$ cm$^{-3}$. Observations and modeling show that the gas phase abundances of C and O progressively decrease along the C+/C/CO transition zone where C/H~8x10$^{-5}$ and C/O~0.8-1, until the beginning of the dense phase at Av~10 mag. This is consistent with the grain temperatures being below the CO evaporation temperature in this region. In the case of sulfur, a strong depletion should occur before the translucent phase where we estimate a S/H~(0.4 - 2.2) x10$^{-6}$, an abundance ~7-40 times lower than the solar value. A second strong depletion must be present during the formation of the thick icy mantles to achieve the values of S/H measured in the dense cold cores (S/H~8x10$^{-8}$). Based on our chemical modeling, we constrain the value of $\zeta_{\rm H_2}$ to ~(0.5 - 1.8) x10$^{-16}$ s$^{-1}$ in the translucent cloud., Comment: 26 pages, 20 figures, A&A in press

Fragmentation of massive dense molecular clouds is the starting point in the formation of rich clusters and massive stars. Theory and numerical simulations indicate that the population of the fragments (number, mass, diameter, separation) resulting from the gravitational collapse of such clumps is probably regulated by the balance between the magnetic field and the other competitors of self-gravity, in particular turbulence and protostellar feedback. We have observed 11 massive, dense and young star-forming clumps with the Atacama Large Millimeter Array (ALMA) in the thermal dust continuum emission at $\sim 1$~mm with an angular resolution of 0.25 arcseconds with the aim of determining their population of fragments. We find fragments on sub-arcsecond scales in 8 out of the 11 sources. The ALMA images indicate two different fragmentation modes: a dominant fragment surrounded by companions with much smaller mass and size, and many ($\geq 8$) fragments with a gradual change in masses and sizes. On average, the largest number of fragments is found towards the warmer and more massive clumps. Also, the warmer clumps tend to form fragments with larger mass and size. To understand the role of the different physical parameters to regulate the final population of the fragments, we have simulated the collapse of a massive clump of $100$ and $300$ M$_{\odot}$ having different magnetic support. The simulations indicate that: (1) fragmentation is inhibited when the initial turbulence is low, independent of the other physical parameters. (2) a filamentary distribution of the fragments is favoured in a highly magnetised clump. We conclude that the clumps that show many fragments distributed in a filamentary-like structure are likely characterised by a strong magnetic field, while the others are possible also in a weaker magnetic field., Comment: 20 pages, 11 figures, 2 Appendices, accepted for publication in A/&A

The formation epoch of protostellar disks is debated because of the competing roles of rotation, turbulence, and magnetic fields in the early stages of low-mass star formation. Magnetohydrodynamics simulations of collapsing cores predict that rotationally supported disks may form in strongly magnetized cores through ambipolar diffusion or misalignment between the rotation axis and the magnetic field orientation. Detailed studies of individual sources are needed to cross check the theoretical predictions. We present 0.06-0.1" resolution images at 350 GHz toward B1b-N and B1b-S, which are young class 0 protostars, possibly first hydrostatic cores. The images have been obtained with ALMA, and we compare these data with magnetohydrodynamics simulations of a collapsing turbulent and magnetized core. The submillimeter continuum emission is spatially resolved by ALMA. Compact structures with optically thick 350 GHz emission are detected toward both B1b-N and B1b-S, with 0.2 and 0.35" radii (46 and 80 au at the Perseus distance of 230 pc), within a more extended envelope. The flux ratio between the compact structure and the envelope is lower in B1b-N than in B1b-S, in agreement with its earlier evolutionary status. The size and orientation of the compact structure are consistent with 0.2" resolution 32 GHz observations obtained with the Very Large Array as a part of the VANDAM survey, suggesting that grains have grown through coagulation. The morphology, temperature, and densities of the compact structures are consistent with those of disks formed in numerical simulations of collapsing cores. Moreover, the properties of B1b-N are consistent with those of a very young protostar, possibly a first hydrostatic core. These observations provide support for the early formation of disks around low-mass protostars., Comment: 8 pages, 10 PDF figures. Accepted for publication in A&A. Uses aa and natbib latex macros

The early evolution of massive cluster progenitors is poorly understood. We investigate the fragmentation properties from 0.3 pc to 0.06 pc scales of a homogenous sample of infrared-quiet massive clumps within 4.5 kpc selected from the ATLASGAL survey. Using the ALMA 7m array we detect compact dust continuum emission towards all targets, and find that fragmentation, at these scales, is limited. The mass distribution of the fragments uncovers a large fraction of cores above 40 $M_\odot$, corresponding to massive dense cores (MDCs) with masses up to ~400 $M_\odot$. 77 % of the clumps contain at most 3 MDCs per clump, and we also reveal single clumps/MDCs. The most massive cores are formed within the more massive clumps, and a high concentration of mass on small scales reveals a high core formation efficiency. The mass of MDCs highly exceeds the local thermal Jeans-mass, and observational evidence is lacking for a sufficiently high level of turbulence or strong enough magnetic fields to keep the most massive MDCs in equilibrium. If already collapsing, the observed fragmentation properties with a high core formation efficiency are consistent with the collapse setting in at parsec scales., Comment: accepted by A&A Letters

Massive stars, multiple stellar systems and clusters are born from the gravitational collapse of massive dense gaseous clumps, and the way these systems form strongly depends on how the parent clump fragments into cores during collapse. Numerical simulations show that magnetic fields may be the key ingredient in regulating fragmentation. Here we present ALMA observations at ~0.25'' resolution of the thermal dust continuum emission at ~278 GHz towards a turbulent, dense, and massive clump, IRAS 16061-5048c1, in a very early evolutionary stage. The ALMA image shows that the clump has fragmented into many cores along a filamentary structure. We find that the number, the total mass and the spatial distribution of the fragments are consistent with fragmentation dominated by a strong magnetic field. Our observations support the theoretical prediction that the magnetic field plays a dominant role in the fragmentation process of massive turbulent clump., Comment: 7 pages, 5 figures, 2 appendices (4 pages and 2 figures in the main text). Accepted for publication in A&A Letters

The first hydrostatic core, also called the first Larson core, is one of the first steps in low-mass star formation, as predicted by theory. With recent and future high performance telescopes, details of these first phases become accessible, and observations may confirm theory and even bring new challenges for theoreticians. In this context, we study from a theoretical point of view the chemical and physical evolution of the collapse of prestellar cores until the formation of the first Larson core, in order to better characterize this early phase in the star formation process. We couple a state-of-the-art hydrodynamical model with full gas-grain chemistry, using different assumptions on the magnetic field strength and orientation. We extract the different components of each collapsing core (i.e., the central core, the outflow, the disk, the pseudodisk, and the envelope) to highlight their specific physical and chemical characteristics. Each component often presents a specific physical history, as well as a specific chemical evolution. From some species, the components can clearly be differentiated. The different core models can also be chemically differentiated. Our simulation suggests some chemical species as tracers of the different components of a collapsing prestellar dense core, and as tracers of the magnetic field characteristics of the core. From this result, we pinpoint promising key chemical species to be observed., Comment: Accepted for publication in ApJ

The origin of massive stars is a fundamental open issue in modern astrophysics. Pre-ALMA interferometric studies reveal precursors to early B to late O type stars with collapsing envelopes of 15-20 M$_\odot$ on 1000-3000 AU size-scales. To search for more massive envelopes we selected the most massive nearby young clumps from the ATLASGAL survey to study their protostellar content with ALMA. Our first results using the intermediate scales revealed by the ALMA ACA array providing 3-5" angular resolution, corresponding to $\sim$0.05-0.1 pc size-scales, reveals a sample of compact objects. These massive dense cores are on average two-times more massive than previous studies of similar types of objects. We expect that once the full survey is completed, it will provide a comprehensive view on the origin of the most massive stars., Comment: 4 pages, 1 figure, to appear in proceedings of The 6th Zermatt ISM Symposium: Conditions and Impact of Star Formation From Lab to Space, eds. R. Simon, M. Rollig

In the theory of star formation, the first hydrostatic core (FHSC) phase is a critical step in which a condensed object emerges from a prestellar core. This step lasts about one thousand years, a very short time compared with the lifetime of prestellar cores, and therefore is hard to detect unambiguously. We present IRAM Plateau de Bure observations of the Barnard 1b dense molecular core, combining detections of H2CO and CH3OH spectral lines and dust continuum at 2.3" resolution (~ 500 AU). The two compact cores B1b-N and B1b-S are detected in the dust continuum at 2mm, with fluxes that agree with their spectral energy distribution. Molecular outflows associated with both cores are detected. They are inclined relative to the direction of the magnetic field, in agreement with predictions of collapse in turbulent and magnetized gas with a ratio of mass to magnetic flux somewhat higher than the critical value, \mu ~ 2 - 7. The outflow associated with B1b-S presents sharp spatial structures, with ejection velocities of up to ~ 7 kms from the mean velocity. Its dynamical age is estimated to be ~2000 yrs. The B1b-N outflow is smaller and slower, with a short dynamical age of ~1000 yrs. The B1b-N outflow mass, mass-loss rate, and mechanical luminosity agree well with theoretical predictions of FHSC. These observations confirm the early evolutionary stage of B1b-N and the slightly more evolved stage of B1b-S., Comment: 6 pages, 3 figures

In order to shed light on the main physical processes controlling fragmentation of massive dense cores, we present a uniform study of the density structure of 19 massive dense cores, selected to be at similar evolutionary stages, for which their relative fragmentation level was assessed in a previous work. We inferred the density structure of the 19 cores through a simultaneous fit of the radial intensity profiles at 450 and 850 micron (or 1.2 mm in two cases) and the Spectral Energy Distribution, assuming spherical symmetry and that the density and temperature of the cores decrease with radius following power-laws. We find a weak (inverse) trend of fragmentation level and density power-law index, with steeper density profiles tending to show lower fragmentation, and vice versa. In addition, we find a trend of fragmentation increasing with density within a given radius, which arises from a combination of flat density profile and high central density and is consistent with Jeans fragmentation. We considered the effects of rotational-to-gravitational energy ratio, non-thermal velocity dispersion, and turbulence mode on the density structure of the cores, and found that compressive turbulence seems to yield higher central densities. Finally, a possible explanation for the origin of cores with concentrated density profiles, which are the cores showing no fragmentation, could be related with a strong magnetic field, consistent with the outcome of radiation magnetohydrodynamic simulations., Comment: accepted for publication in ApJ; a version with all (high resolution) figures is available at: http://www.ice.csic.es/files/palau/Palau14.pdf

Our aim is to study the thermal and dynamical evolution of protoplanetary disks in global simulations, including the physics of radiation transfer and magneto-hydrodynamic (MHD) turbulence caused by the magneto-rotational instability. We develop a radiative transfer method based on the flux-limited diffusion approximation that includes frequency dependent irradiation by the central star. This hybrid scheme is implemented in the PLUTO code. The focus of our implementation is on the performance of the radiative transfer method. Using an optimized Jacobi preconditioned BiCGSTAB solver, the radiative module is three times faster than the MHD step for the disk setup we consider. We obtain weak scaling efficiencies of 70% up to 1024 cores. We present the first global 3D radiation MHD simulations of a stratified protoplanetary disk. The disk model parameters are chosen to approximate those of the system AS 209 in the star-forming region Ophiuchus. Starting the simulation from a disk in radiative and hydrostatic equilibrium, the magnetorotational instability quickly causes MHD turbulence and heating in the disk. For the disk parameters we use, turbulent dissipation heats the disk midplane and raises the temperature by about 15% compared to passive disk models. A roughly flat vertical temperature profile establishes in the disk optically thick region close to the midplane. We reproduce the vertical temperature profile with a viscous disk models for which the stress tensor vertical profile is flat in the bulk of the disk and vanishes in the disk corona. The present paper demonstrates for the first time that global radiation MHD simulations of turbulent protoplanetary disks are feasible with current computational facilities. This opens up the windows to a wide range of studies of the dynamics of protoplanetary disks inner parts, for which there are significant observational constraints., Comment: Accepted to A&A

Star formation begins with the gravitational collapse of a dense core inside a molecular cloud. As the collapse progresses, the centre of the core begins to heat up as it becomes optically thick. The temperature and density in the centre eventually reach high enough values where fusion reactions can ignite; the protostar is born. This sequence of events entail many physical processes, of which radiative transfer is of paramount importance. Many simulations of protostellar collapse make use of a grey treatment of radiative transfer coupled to the hydrodynamics. However, interstellar gas and dust opacities present large variations as a function of frequency. In this paper, we follow-up on a previous paper on the collapse and formation of Larson's first core using multigroup radiation hydrodynamics (Paper I) by extending the calculations to the second phase of the collapse and the formation of Larson's second core. We have made the use of a non-ideal gas equation of state as well as an extensive set of spectral opacities in a spherically symmetric fully implicit Godunov code to model all the phases of the collapse of a 0.1, 1 and 10 solar mass cloud cores. We find that, for a same central density, there are only small differences between the grey and multigroup simulations. The first core accretion shock remains supercritical while the shock at the second core border is found to be strongly subcritical with all the accreted energy being transfered to the core. The size of the first core was found to vary somewhat in the different simulations (more unstable clouds form smaller first cores) while the size, mass and temperature of the second cores are independent of initial cloud mass, size and temperature. Our simulations support the idea of a standard (universal) initial second core size of 0.003 AU and mass 0.0014 solar masses., Comment: 15 pages, 11 figures, accepted for publication in A&A

Observations of First Hydrostatic Core candidates, a theoretically predicted evolutionary link between the prestellar and protostellar phases, are vital for probing the earliest phases of star formation. We aim to determine the dynamical state of the First Hydrostatic Core candidate Cha-MMS1. We observed Cha-MMS1 in various transitions with the APEX and Mopra telescopes. The molecular emission was modeled with a radiative transfer code to derive constraints on the envelope kinematics. We derive an internal luminosity of 0.08 - 0.18 Lsol. An average velocity gradient of 3.1(0.1) km/s/pc over 0.08 pc is found perpendicular to the filament in which Cha-MMS1 is embedded. The gradient is flatter in the outer parts and at the innermost 2000 - 4000 AU. These features suggest solid-body rotation beyond 4000 AU and slower, differential rotation beyond 8000 AU. The origin of the flatter gradient in the innermost parts is unclear. The classical infall signature is detected in HCO+ 3-2 and CS 2-1. The radiative transfer modeling indicates a uniform infall velocity in the outer parts of the envelope. An infall velocity field scaling with r^(-0.5) is consistent with the data for r < 9000 AU. The infall velocities are 0.1 - 0.2 km/s at r > 3300 AU and 0.04 - 0.6 km/s at r < 3300 AU. Both the internal luminosity of Cha-MMS1 and the infall velocity field in its envelope are consistent with predictions of MHD simulations for the first core phase. There is no evidence for a fast, large-scale outflow stemming from Cha-MMS1 but excess emission from the high-density tracers CS 5-4, CO 6-5, and CO 7-6 suggests the presence of higher-velocity material at the inner core. Its internal luminosity excludes Cha-MMS1 being a prestellar core. The kinematical properties of its envelope are consistent with Cha-MMS1 being a first core candidate or a very young Class 0 protostar.(abridged)., Comment: Accepted for publication in A&A. 27 pages, 22 figures, 13 tables. A version with high-resolution figures is available on request to the first author

We present a possible identification strategy for first hydrostatic core (FHSC) candidates and make predictions of ALMA dust continuum emission maps from these objects. We analyze the results given by the different bands and array configurations and identify which combinations of the two represent our best chance of solving the fragmentation issue in these objects. If the magnetic field is playing a role, the emission pattern will show evidence of a pseudo-disk and even of a magnetically driven outflow, which pure hydrodynamical calculations cannot reproduce., Comment: 2 pages, to appear in the proceedings of the conference "New Trends in Radio Astronomy in the ALMA Era", Hakone, Japan, December 3-8, 2012

Stars and more particularly massive stars, have a drastic impact on galaxy evolution. Yet the conditions in which they form and collapse are still not fully understood. In particular, the influence of the magnetic field on the collapse of massive clumps is relatively unexplored, it is thus of great relevance in the context of the formation of massive stars to investigate its impact. We perform high resolution, MHD simulations of the collapse of hundred solar masses, turbulent and magnetized clouds, using the adaptive mesh refinement code RAMSES. We compute various quantities such as mass distribution, magnetic field and angular momentum within the collapsing core and study the episodic outflows and the fragmentation that occurs during the collapse. The magnetic field has a drastic impact on the cloud evolution. We find that magnetic braking is able to substantially reduce the angular momentum in the inner part of the collapsing cloud. Fast and episodic outflows are being launched with typical velocities of the order of 3-5 km s$^{-1}$ although the highest velocities can be as high as 30-40 km s$^{-1}$. The fragmentation in several objects, is reduced in substantially magnetized clouds with respect to hydrodynamical ones by a factor of the order of 1.5-2. We conclude that magnetic fields have a significant impact on the evolution of massive clumps. In combination with radiation, magnetic fields largely determine the outcome of massive core collapse. We stress that numerical convergence of MHD collapse is a challenging issue. In particular, numerical diffusion appears to be important at high density therefore possibly leading to an over-estimation of the number of fragments., Comment: accepted for publication in A&A

It is established that both radiative transfer and magnetic field have a strong impact on the collapse and the fragmentation of prestellar dense cores, but no consistent calculation exists yet at such scales. We present original AMR calculations including magnetic field (in the ideal MHD limit) and radiative transfer, within the Flux Limited Diffusion approximation, of the collapse of a 1 solar mass dense core. We compare the results with calculations performed with a barotropic EOS. We show that radiative transfer has an important impact on the collapse and the fragmentation, through the cooling or heating of the gas, and is complementary of the magnetic field. A larger field yields a stronger magnetic braking, increasing the accretion rate and thus the effect of the radiative feedback. Even for a strongly magnetized core, where the dynamics of the collapse is dominated by the magnetic field, radiative transfer is crucial to determine the temperature and optical depth distributions, two potentially accessible observational diagnostics. A barotropic EOS cannot account for realistic fragmentation. The diffusivity of the numerical scheme, however, is found to strongly affect the output of the collapse, leading eventually to spurious fragmentation. Both radiative transfer and magnetic field must be included in numerical calculations of star formation to obtain realistic collapse configurations and observable signatures. Nevertheless, the numerical resolution and the robustness of the solver are of prime importance to obtain reliable results. When using an accurate solver, the fragmentation is found to always remain inhibited by the magnetic field, at least in the ideal MHD limit, even when radiative transfer is included., Comment: 6 pages, 7 figures, accepted for publication in Astronomy and Astrophysics

The development of parallel supercomputers allows today the detailed study of the collapse and the fragmentation of prestellar cores with increasingly accurate numerical simulations. Thanks to the advances in sub-millimeter observations, a wide range of observed initial conditions enable us to study the different modes of low-mass star formation. The challenge for the simulations is to reproduce the observational results. Two main numerical methods, namely AMR and SPH, are widely used to simulate the collapse and the fragmentation of prestellar cores. We compare here thoroughly these two methods with numerical resolution requirements deduced from previous studies. Our physical model is as simple as possible, and consists of an isothermal sphere rotating around the z-axis. We first study the conservation of angular momentum as a function of the resolution. Then, we explore a wide range of simulation parameters to study the fragmentation of prestellar cores. There seems to be a convergence between the two methods, provided resolution in each case is sufficient. Resolution criteria adapted to our physical cases, in terms of resolution per Jeans mass, for an accurate description of the formation of protostellar cores are deduced from the present study. This convergence is encouraging for future work in simulations of low-mass star formation, providing the aforementioned criteria are fulfilled., Comment: 6 pages, 2 figures, To appear in the proceedings of "SF2A-2007: Semaine de l'Astrophysique Francaise", (J. Bouvier, A. Chalabaev, and C. Charbonnel eds)