Ambipolar diffusion in low-mass star formation. I. General comparison with the ideal MHD case

In this paper, we provide a more accurate description of the evolution of the magnetic flux redistribution during prestellar core collapse by including resistive terms in the magnetohydrodynamics (MHD) equations. We focus more particularly on the impact of ambipolar diffusion. We use the adaptive mesh refinement code RAMSES to carry out such calculations. The resistivities required to calculate the ambipolar diffusion terms were computed using a reduced chemical network of charged, neutral and grain species. The inclusion of ambipolar diffusion leads to the formation of a magnetic diffusion barrier in the vicinity of the core, preventing accumulation of magnetic flux in and around the core and amplification of the field above 0.1G. The mass and radius of the first Larson core remain similar between ideal and non-ideal MHD models. This diffusion plateau has crucial consequences on magnetic braking processes, allowing the formation of disk structures. Magnetically supported outflows launched in ideal MHD models are weakened when using non-ideal MHD. Contrary to ideal MHD misalignment between the initial rotation axis and the magnetic field direction does not significantly affect the results for a given mu, showing that the physical dissipation truly dominate over numerical diffusion. We demonstrate severe limits of the ideal MHD formalism, which yield unphysical behaviours in the long-term evolution of the system. This includes counter rotation inside the outflow, interchange instabilities, and flux redistribution triggered by numerical diffusion, none observed in non-ideal MHD. Disks with Keplerian velocity profiles form in all our non-ideal MHD simulations, with final mass and size which depend on the initial magnetisation. This ranges from a few 0.01 solar masses and 20-30 au for the most magnetised case (mu=2) to 0.2 solar masses and 40-80 au for a lower magnetisation (mu=5).Comment: Accepted in A&A section 7 (on Wednesday, september the 16th, year 2015

Simulations of early structure formation: Properties of halos that host primordial star formation

International audiencePopulation III (pop III) stars were born in halos characterised by a pristine gas composition. In such a halo, once the gas density reaches n$_{\mathrm{H}} \sim$ 1 cm$^{-3}$, molecular cooling leads to the collapse of the gas and the birth of pop III stars. Halo properties, such as the chemical abundances, mass, and angular momentum can affect the collapse of the gas, thereby leading to the pop III initial mass function (IMF) of star formation. We want to study the properties of primordial halos and how halos that host early star formation differ from other types of halos. The aim of this study is to obtain a representative population of halos at a given redshift hosting a cold and massive gas cloud that enables the birth of the first stars. We investigated the growth of primordial halos in a $\Lambda$CDM Universe in a large cosmological simulation. We used the hydrodynamic code RAMSES and the chemical solver KROME to study halo formation with non-equilibrium thermochemistry. We then identified structures in the dark and baryonic matter fields, thereby linking the presence or absence of dense gas clouds to the mass and the physical properties of the hosting halos. In our simulations, the mass threshold for a halo for hosting a cold dense gas cloud is $\simeq 7 \times 10^5 M_{\odot}$ and the threshold in the H$_2$ mass fraction is found to be $\simeq 2 \times 10^{-4}$. This is in agreement with previous works. We find that the halo history and accretion rate play a minor role. Here, we present halos with higher HD abundances, which are shown to be colder, as the temperature in the range between $10^2 - 10^4 \, \mathrm{cm^{-3}}$ depends on the HD abundance to a large extent. The higher fraction of HD is linked to the higher spin parameter that is seen for the dense gas

What determines the formation and characteristics of protoplanetary discs?

International audienceContext. Planets form in protoplanetary discs. Their masses, distribution, and orbits sensitively depend on the structure of the protoplanetary discs. However, what sets the initial structure of the discs in terms of mass, radius and accretion rate is still unknown.Aims. It is therefore of great importance to understand exactly how protoplanetary discs form and what determines their physical properties. We aim to quantify the role of the initial dense core magnetisation, rotation, turbulence, and misalignment between rotation and magnetic field axis as well as the role of the accretion scheme onto the central object.Methods. We performed non-ideal magnetohydrodynamics numerical simulations using the adaptive mesh refinement code Ramses of a collapsing, one solar mass molecular core to study the disc formation and early, up to 100 kyr, evolution. We paid particular attention to the impact of numerical resolution and accretion scheme.Results. We found that the mass of the central object is almost independent of the numerical parameters such as the resolution and the accretion scheme onto the sink particle. The disc mass and to a lower extent its size, however heavily depend on the accretion scheme, which we found is itself resolution dependent. This implies that the accretion onto the star and through the disc are largely decoupled. For a relatively large domain of initial conditions (except at low magnetisation), we found that the properties of the disc do not change too significantly. In particular both the level of initial rotation and turbulence do not influence the disc properties provide the core is sufficiently magnetised. After a short relaxation phase, the disc settles in a stationary state. It then slowly grows in size but not in mass. The disc itself is weakly magnetised but its immediate surrounding on the contrary is highly magnetised.Conclusions. Our results show that the disc properties directly depend on the inner boundary condition, i.e. the accretion scheme onto the central object. This suggests that the disc mass is eventually controlled by a small-scale accretion process, possibly the star-disc interaction. Because of ambipolar diffusion and its significant resistivity, the disc diversity remains limited and except for low magnetisation, their properties are weakly sensitive to initial conditions such as rotation and turbulence

Fragmentation of massive dense cores down to <~ 1000 AU: relation between fragmentation and density structure

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 μm (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. Even though the estimated fragmentation level is strictly speaking a lower limit, its relative value is significant and several trends could be explored with our data. 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

Fragmentation of massive dense cores down to ~1000 AU: Relation between fragmentation and density structure

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