Accretion across scales: streamers, surface-layer transport, and rapid replenishment in young protoplanetary discs

Abstract

Protoplanetary discs evolve around newly-formed stars through an interplay of infall from surrounding turbulent cloud material, accretion towards the young star, and outflow driven mass-loss. It has been challenging to determine if discs are fed predominantly through infall along the disc midplane, or along the poles, and if accretion occurs in a steady or burst-like fashion. Here, we present a suite of 3D ideal magnetohydrodynamical simulations of protoplanetary disc formation and evolution in a dynamic, large-scale molecular cloud environment using the adaptive mesh refinement framework DISPATCH. We focus on nine stellar systems, where we resolve discs down to a scale of 0.8 au. Across the sample, stellar accretion proceeds at rates of $\sim$10${}^{-5}$ M${}_{\odot }$ yr${}^{-1}$ over 10${}^5$ yr, with significant variability. Discs grow to 100 au scales and remain gravitationally stable in time, with disc-to-star mass ratios below 10 %. Transient high-density streamers, with 10 kyr infall times, can drive anisotropic mass delivery at rates comparable to the background accretion flow. Their interaction with discs typically results in a temporary reduction of the disc size by half, and disc mass by 40 %. During later quiescent disc evolution stages ($t\gtrsim$50 kyr), accretion predominantly occurs through the midplane and disc surface layers. This is associated with the development of a toroidal magnetic field morphology, which includes field reversals across both disc surfaces. In this way, the full disc mass reservoir is replenished on 10 kyr-timescales. These findings support that the outer parts of very young discs, when well-ionised and close to the ideal MHD regime, are not yet conducive to planet formation, due to high replenishment rates, strong turbulence, and disruptive streamer infall events.

21 pages, 19 figures