Dust coagulation feedback on magnetohydrodynamic resistivities in protostellar collapse

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 (nH ~ 108 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.

Protostellar collapse: the conditions to form dust-rich protoplanetary disks

Context. Dust plays a key role during star, disk, and planet formation. Yet, its dynamics during the protostellar collapse remain a poorly investigated field. Recent studies seem to indicate that dust may decouple efficiently from the gas during these early stages. Aims. We aim to understand how much and in which regions dust grains concentrate during the early phases of the protostellar collapse, and to see how this depends on the properties of the initial cloud and of the solid particles. Methods. We used the multiple species dust dynamics MULTIGRAIN solver of the grid-based code RAMSES to perform various simulations of dusty collapses. We performed hydrodynamical and magnetohydrodynamical simulations where we varied the maximum size of the dust distribution, the thermal-to-gravitational energy ratio, and the magnetic properties of the cloud. We simulated the simultaneous evolution of ten neutral dust grain species with grain sizes varying from a few nanometers to a few hundreds of microns. Results. We obtain a significant decoupling between the gas and the dust for grains of typical sizes of a few tens of microns. This decoupling strongly depends on the thermal-to-gravitational energy ratio, the grain sizes, and the inclusion of a magnetic field. With a semi-analytic model calibrated on our results, we show that the dust ratio mostly varies exponentially with the initial Stokes number at a rate that depends on the local cloud properties. Conclusions. We find that larger grains tend to settle and drift efficiently in the first-core and in the newly formed disk. This can produce dust-to-gas ratios of several times the initial value. Dust concentrates in high-density regions (cores, disk, and pseudo-disk) and is depleted in low-density regions (envelope and outflows). The size at which grains decouple from the gas depends on the initial properties of the clouds. Since dust cannot necessarily be used as a proxy for gas during the collapse, we emphasize the necessity of including the treatment of its dynamics in protostellar collapse simulations.

Ionization: a possible explanation for the difference of mean disk sizes in star-forming regions

Context. Surveys of protoplanetary disks in star-forming regions of similar age revealed significant variations in average disk mass in some regions. For instance, disks in the Orion Nebular Cluster (ONC) and Corona Australis (CrA) are on average smaller than disks observed in Lupus, Taurus, Chamaeleon I, or Ophiuchus. Aims. In contrast to previous models that studied the truncation of disks at a late stage of their evolution, we investigate whether disks may already be born with systematically smaller disk sizes in more massive star-forming regions as a consequence of higher ionization rates. Methods. Assuming various cosmic-ray ionization rates, we computed the resistivities for ambipolar diffusion and Ohmic dissipation with a chemical network, and performed 2D nonideal magnetohydrodynamical protostellar collapse simulations. Results. A higher ionization rate leads to stronger magnetic braking, and hence to the formation of smaller disks. Accounting for recent findings that protostars act as forges of cosmic rays and considering only mild attenuation during the collapse phase, we show that a high average cosmic-ray ionization rate in star-forming regions such as the ONC or CrA can explain the detection of smaller disks in these regions. Conclusions. Our results show that on average, a higher ionization rate leads to the formation of smaller disks. Smaller disks in regions of similar age can therefore be the consequence of different levels of ionization, and may not exclusively be caused by disk truncation through external photoevaporation. We strongly encourage observations that allow measuring the cosmic-ray ionization degrees in different star-forming regions to test our hypothesis.

Episodic accretion constrained by a rich cluster of outflows

Context. The accretion history of protostars remains widely mysterious, even though it represents one of the best ways to understand the protostellar collapse that leads to the formation of stars. Aims. Molecular outflows, which are easier to detect than the direct accretion onto the prostellar embryo, are here used to characterize the protostellar accretion phase in W43-MM1. Methods. The W43-MM1 protocluster hosts a sufficient number of protostars to statistically investigate molecular outflows in a single, homogeneous region. We used the CO(2–1) and SiO(5–4) line datacubes, taken as part of an ALMA mosaic with a 2000 AU resolution, to search for protostellar outflows, evaluate the influence that the environment has on these outflows’ characteristics and put constraints on outflow variability in W43-MM1. Results. We discovered a rich cluster of 46 outflow lobes, driven by 27 protostars with masses of 1−100 M⊙. The complex environment inside which these outflow lobes develop has a definite influence on their length, limiting the validity of using outflows’ dynamical timescale as a proxy of the ejection timescale in clouds with high dynamics and varying conditions. We performed a detailed study of Position–Velocity diagrams of outflows that revealed clear events of episodic ejection. The time variability of W43-MM1 outflows is a general trend and is more generally observed than in nearby, low- to intermediate-mass star-forming regions. The typical timescale found between two ejecta, ~500 yr, is consistent with that found in nearby protostars. Conclusions. If ejection episodicity reflects variability in the accretion process, either protostellar accretion is more variable, or episodicity is easier to detect in high-mass star-forming regions than in nearby clouds. The timescale found between accretion events could result from instabilities associated with bursts of inflowing gas arising from the close dynamical environment of high-mass star-forming cores.

Impact of the Hall effect in star formation, improving the angular momentum conservation

We present here a minor modification of our numerical implementation of the Hall effect for the 2D Riemann solver used in Constrained Transport schemes, as described in a former paper. In the previous work, the tests showed that the angular momentum was not conserved during protostellar collapse simulations, with significant impact. By removing the whistler waves speed from the characteristic speeds of non-magnetic variables in the 1D Riemann solver, we were able to improve the angular momentum conservation in our test-case by one order of magnitude, while keeping the second-order numerical convergence of the scheme. We also reproduce the simulations of a previous study with consistent resistivities, the three non-ideal magnetohydrodynamic effects and initial rotation, and agree with their results. In this case, the violation of angular momentum conservation is negligible in regard to the total angular momentum and the angular momentum of the disk.

Small dust grain dynamics on adaptive mesh refinement grids

Context. Small dust grains are essential ingredients of star, disk and planet formation. Aims. We present an Eulerian numerical approach to study small dust grain dynamics in the context of star and protoplanetary disk formation. It is designed for finite volume codes. We use it to investigate dust dynamics during the protostellar collapse. Methods. We present a method to solve the monofluid equations of gas and dust mixtures with several dust species in the diffusion approximation implemented in the adaptive-mesh-refinement code RAMSES. It uses a finite volume second-order Godunov method with a predictor-corrector MUSCL scheme to estimate the fluxes between the grid cells. Results. We benchmark our method against six distinct tests, DUSTYADVECT, DUSTYDIFFUSE, DUSTYSHOCK, DUSTYWAVE, SETTLING, and DUSTYCOLLAPSE. We show that the scheme is second-order accurate in space on uniform grids and intermediate between second- and first-order on non-uniform grids. We apply our method on various DUSTYCOLLAPSE simulations of 1 M⊙ cores composed of gas and dust. Conclusions. We developed an efficient approach to treat gas and dust dynamics in the diffusion regime on grid-based codes. The canonical tests were successfully passed. In the context of protostellar collapse, we show that dust is less coupled to the gas in the outer regions of the collapse where grains larger than ≃100 μm fall significantly faster than the gas.

Characterizing young protostellar disks with the CALYPSO IRAM-PdBI survey: large Class 0 disks are rare

Context. Understanding the formation mechanisms of protoplanetary disks and multiple systems and also their pristine properties are key questions for modern astrophysics. The properties of the youngest disks, embedded in rotating infalling protostellar envelopes, have largely remained unconstrained up to now. Aims. We aim to observe the youngest protostars with a spatial resolution that is high enough to resolve and characterize the progenitors of protoplanetary disks. This can only be achieved using submillimeter and millimeter interferometric facilities. In the framework of the IRAM Plateau de Bure Interferometer survey CALYPSO, we have obtained subarcsecond observations of the dust continuum emission at 231 and 94 GHz for a sample of 16 solar-type Class 0 protostars. Methods. In an attempt to identify disk-like structures embedded at small scales in the protostellar envelopes, we modeled the dust continuum emission visibility profiles using Plummer-like envelope models and envelope models that include additional Gaussian disk-like components. Results. Our analysis shows that in the CALYPSO sample, 11 of the 16 Class 0 protostars are better reproduced by models including a disk-like dust continuum component contributing to the flux at small scales, but less than 25% of these candidate protostellar disks are resolved at radii >60 au. Including all available literature constraints on Class 0 disks at subarcsecond scales, we show that our results are representative: most (>72% in a sample of 26 protostars) Class 0 protostellar disks are small and emerge only at radii <60 au. We find a multiplicity fraction of the CALYPSO protostars ≲57% ± 10% at the scales 100–5000 au, which generally agrees with the multiplicity properties of Class I protostars at similar scales. Conclusions. We compare our observational constraints on the disk size distribution in Class 0 protostars to the typical disk properties from protostellar formation models. If Class 0 protostars contain similar rotational energy as is currently estimated for prestellar cores, then hydrodynamical models of protostellar collapse systematically predict a high occurrence of large disks. Our observations suggest that these are rarely observed, however. Because they reduce the centrifugal radius and produce a disk size distribution that peaks at radii <100 au during the main accretion phase, magnetized models of rotating protostellar collapse are favored by our observations.