Disc formation in magnetized dense cores with turbulence and ambipolar diffusion

ABSTRACT Discs are essential to the formation of both stars and planets, but how they form in magnetized molecular cloud cores remains debated. This work focuses on how the disc formation is affected by turbulence and ambipolar diffusion (AD), both separately and in combination, with an emphasis on the protostellar mass accretion phase of star formation. We find that a relatively strong, sonic turbulence on the core scale strongly warps but does not completely disrupt the well-known magnetically induced flattened pseudo-disc that dominates the inner protostellar accretion flow in the laminar case, in agreement with previous work. The turbulence enables the formation of a relatively large disc at early times with or without AD, but such a disc remains strongly magnetized and does not persist to the end of our simulation unless a relatively strong AD is also present. The AD-enabled discs in laminar simulations tend to fragment gravitationally. The disc fragmentation is suppressed by initial turbulence. The AD facilitates the disc formation and survival by reducing the field strength in the circumstellar region through magnetic flux redistribution and by making the field lines there less pinched azimuthally, especially at late times. We conclude that turbulence and AD complement each other in promoting disc formation. The discs formed in our simulations inherit a rather strong magnetic field from its parental core, with a typical plasma-β of order a few tens or smaller, which is 2–3 orders of magnitude lower than the values commonly adopted in magnetohydrodynamic simulations of protoplanetary discs. To resolve this potential tension, longer term simulations of disc formation and evolution with increasingly more realistic physics are needed.

Stratified shear instability in a field of pre-existing turbulence

Turbulent mixing of heat and momentum in the stably-stratified ocean interior occurs in discrete events driven by vertical variations of the horizontal velocity. Typically, these events have been modelled assuming an initially laminar stratified shear flow which develops wavelike instabilities, becomes fully turbulent, and then relaminarizes into a stable state. However, in the real ocean there is always some level of turbulence left over from previous events. Using direct numerical simulations, we show that the evolution of a stably-stratified shear layer may be significantly modified by pre-existing turbulence. The classical billow structure associated with Kelvin–Helmholtz instability is suppressed and eventually eliminated as the strength of the initial turbulence is increased. A corresponding energetics analysis shows that potential energy changes and dissipation of kinetic energy depend non-monotonically on initial turbulence strength, with the largest effects when initial turbulence is present but insufficient to prevent billow formation. The mixing efficiency decreases with increasing initial turbulence amplitude as the development of the Kelvin–Helmholtz billow, with its large pre-turbulent mixing efficiency, is arrested.

Fragmentation properties of massive protocluster gas clumps: an ALMA study

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, and 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 ~1 mm with an angular resolution of 0.′′25 with the aim of determining their population of fragments. The targets have been selected from a sample of massive molecular clumps with limited or absent star formation activity and hence limited feedback. 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 lower mass and smaller size, and many (≥8) fragments with a gradual change in masses and sizes. The morphologies are very different, with three sources that show filament-like distributions of the fragments, while the others have irregular geometry. On average, the largest number of fragments is found towards the warmer and more massive clumps. The warmer clumps also tend to form fragments with higher mass and larger size. To understand the role of the different physical parameters in regulating the final population of the fragments, we simulated the collapse of a massive clump of 100 and 300M⊙ with different magnetic support. The 300 M⊙ case was also run for different initial temperatures and Mach numbers M to evaluate the separate role of each of these parameters. The simulations indicate that (1) fragmentation is inhibited when the initial turbulence is low (M ~ 3), independent of the other physical parameters. This would indicate that the number of fragments in our clumps can be explained assuming a high (M ~ 6) initial turbulence, although an initial density profile different to that assumed can play a relevant role. (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 filament-like structure are likely characterised by a strong magnetic field, while the other morphologies are also possible in a weaker magnetic field.