We investigate the formation of protoplanetary disks around nine solar mass stars formed in the context of a (40 pc) 3 Giant Molecular Cloud model, using ramses adaptive-mesh refinement simulations extending over a scale range of about 4 million, from an outer scale of 40 pc down to cell sizes of 2 AU. Our most important result is that the accretion process is heterogeneous in multiple ways; in time, in space, and among protostars of otherwise similar mass. Accretion is heterogeneous in time, in the sense that accretion rates vary during the evolution, with generally decreasing profiles, whose slopes vary over a wide range, and where accretion can increase again if a protostar enters a region with increased density and low speed. Accretion is heterogeneous in space, because of the mass distribution, with mass approaching the accreting star-disk system in filaments and sheets. Finally, accretion is heterogeneous among stars, since the detailed conditions and dynamics in the neighborhood of each star can vary widely. We also investigate the sensitivity of disk formation to physical conditions, and test their robustness by varying numerical parameters. We find that disk formation is robust even when choosing the least favorable sink particle parameters, and that turbulence cascading from larger scales is a decisive factor in disk formation. We also investigate the transport of angular momentum, finding that the net inward mechanical transport is compensated for mainly by an outward directed magnetic transport, with a contribution from gravitational torques usually subordinate to the magnetic transport.

After its formation, a young star spends some time traversing the molecular cloud complex in which it was born. It is therefore not unlikely that, well after the initial cloud collapse event which produced the star, it will encounter one or more low mass cloud fragments, which we call "cloudlets" to distinguish them from full-fledged molecular clouds. Some of this cloudlet material may accrete onto the star+disk system, while other material may fly by in a hyperbolic orbit. In contrast to the original cloud collapse event, this process will be a "cloudlet flyby" and/or "cloudlet capture" event: A Bondi-Hoyle-Lyttleton type accretion event, driven by the relative velocity between the star and the cloudlet. As we will show in this paper, if the cloudlet is small enough and has an impact parameter similar or less than GM * /v 2 ∞ (with v∞ being the approach velocity), such a flyby and/or capture event would lead to arc-shaped or tail-shaped reflection nebulosity near the star. Those shapes of reflection nebulosity can be seen around several transitional disks and FU Orionis stars. Although the masses in the those arcs appears to be much less than the disk masses in these sources, we speculate that higher-mass cloudlet capture events may also happen occasionally. If so, they may lead to the tilting of the outer disk, because the newly infalling matter will have an angular momentum orientation entirely unrelated to that of the disk. This may be one possible explanation for the highly warped/tilted inner/outer disk geometries found in several transitional disks. We also speculate that such events, if massive enough, may lead to FU Orionis outbursts.

Using zoom-simulations carried out with the adaptive mesh-refinement code ramses with a dynamic range of up to 2 27 ≈ 1.34 × 10 8 we investigate the accretion profiles around six stars embedded in different environments inside a (40 pc) 3 giant molecular cloud, the role of mass infall and disc instabilities on the accretion profile, and thus on the luminosity of the forming protostar. Our results show that the environment in which the protostar is embedded determines the overall accretion profile of the protostar. Infall on to the circumstellar disc may trigger gravitational disc instabilities in the disc at distances of around ∼10 to ∼50 au leading to rapid transport of angular momentum and strong accretion bursts. These bursts typically last for about ∼10 to a ∼100 yr, consistent with typical orbital times at the location of the instability, and enhance the luminosity of the protostar. Calculations with the stellar evolution code mesa show that the accretion bursts induce significant changes in the protostellar properties, such as the stellar temperature and radius. We apply the obtained protostellar properties to produce synthetic observables with radmc3d and predict that accretion bursts lead to observable enhancements around 20 to 200 µm in the spectral energy distribution of Class 0 type young stellar objects.