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)

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 Marchand et al. (2018). 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 are 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 Tsukamoto et al. (2015) with consistent resistivities, the three non-ideal MHD 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.

We present an implementation of the Hall term in the non-ideal magnetohydrodynamics equations into the adaptive-mesh-refinement code RAMSES to study its impact on star formation. Recent works show that the Hall effect heavily influences the regulation of the angular momentum in collapsing dense cores, strengthening or weakening the magnetic braking. Our method consists of a modification of the two-dimensional constrained transport scheme. Our scheme shows convergence of second-order in space and the frequency of the propagation of whistler waves is accurate. We confirm previous results, namely that during the collapse, the Hall effect generates a rotation of the fluid with a direction in the mid-plane that depends on the sign of the Hall resistivity, while counter-rotating envelopes develop on each side of the mid-plane. However, we find that the predictability of our numerical results is severely limited. The angular momentum is not conserved in any of our dense core-collapse simulations with the Hall effect: a large amount of angular momentum is generated within the first Larson core, a few hundred years after its formation, without compensation by the surrounding gas. This issue is not mentioned in previous studies and may be correlated to the formation of the accretion shock on the Larson core. We expect that this numerical effect could be a serious issue in star formation simulations.

The formation of protoplanetary discs during the collapse of molecular dense cores is significantly influenced by angular momentum transport, notably by the magnetic torque. In turn, the evolution of the magnetic field is determined by dynamical processes and non-ideal MHD effects such as ambipolar diffusion. Considering simple relations between various timescales characteristic of the magnetized collapse, we derive an expression for the early disc radius, $ r \simeq 18 \, {\rm AU} \, \left({\eta_{\rm AD} / 0.1 \, {\rm s}} \right)^{2/9} \left({B_z / 0.1\, {\rm G}} \right) ^{-4/9} \left({M / 0.1 \msol} \right) ^{1/3},$ where $M$ is the total disc plus protostar mass, $\eta_\mathrm{AD}$ is the ambipolar diffusion coefficient and $B_z$ is the magnetic field in the inner part of the core. This is about significantly smaller than the discs that would form if angular momentum was conserved. The analytical predictions are confronted against a large sample of 3D, non-ideal MHD collapse calculations covering variations of a factor 100 in core mass, a factor 10 in the level of turbulence, a factor 5 in rotation, and magnetic mass-to-flux over critical mass-to-flux ratios 2 and 5. The disc radius estimates are found to agree with the numerical simulations within less than a factor 2. A striking prediction of our analysis is the weak dependence of circumstellar disc radii upon the various relevant quantities, suggesting weak variations among class-0 disc sizes. In some cases, we note the onset of large spiral arms beyond this radius.

We develop a detailed chemical network relevant to the conditions characteristic of prestellar core collapse. We solve the system of time-dependent differential equations to calculate the equilibrium abundances of molecules and dust grains, with a size distribution given by size-bins for these latter. These abundances are used to compute the different non-ideal magneto-hydrodynamics resistivities (ambipolar, Ohmic and Hall), needed to carry out simulations of protostellar collapse. For the first time in this context, we take into account the evaporation of the grains, the thermal ionisation of Potassium, Sodium and Hydrogen at high temperature, and the thermionic emission of grains in the chemical network, and we explore the impact of various cosmic ray ionisation rates. All these processes significantly affect the non-ideal magneto-hydrodynamics resistivities, which will modify the dynamics of the collapse. Ambipolar diffusion and Hall effect dominate at low densities, up to n_H = 10^12 cm^-3, after which Ohmic diffusion takes over. We find that the time-scale needed to reach chemical equilibrium is always shorter than the typical dynamical (free fall) one. This allows us to build a large, multi-dimensional multi-species equilibrium abundance table over a large temperature, density and ionisation rate ranges. This table, which we make accessible to the community, is used during first and second prestellar core collapse calculations to compute the non-ideal magneto-hydrodynamics resistivities, yielding a consistent dynamical-chemical description of this process.