We study the formation and long-term evolution of primordial protostellar disks harbored by first stars using numerical hydrodynamics simulations in the thin-disk limit. The initial conditions are specified by pre-stellar cores with distinct mass, angular momentum, and temperature. This allows us to probe several tens of thousand years of the disk's initial evolution, during which we observe multiple episodes of fragmentation leading to the formation of gravitationally bound gaseous clumps within spiral arms. These fragments are torqued inward due to gravitational interaction with the spiral arms on timescales of 10 3 -10 4 yr and accreted onto the growing protostar, giving rise to accretion and luminosity bursts. The burst phenomenon is fueled by continuing accretion of material falling onto the disk from the collapsing parent core, which replenishes the mass lost by the disk due to accretion, and triggers repetitive episodes of disk fragmentation. We show that the burst phenomenon is expected to occur for a wide spectrum of initial conditions in primordial pre-stellar cores and speculate on how the intense luminosities (∼10 7 L ⊙ ) produced by this mechanism may have important consequences for the disk evolution and subsequent growth of the protostar.

We analyze the time evolution of the luminosity of a cluster of Population III protostars formed in the early universe. We argue from the Jeans criterion that primordial gas can collapse to form a cluster of first stars that evolve relatively independently of one another (i.e., with negligible gravitational interaction). We model the collapse of individual protostellar clumps using 2+1D nonaxisymmetric numerical hydrodynamics simulations. Each collapse produces a protostar surrounded by a massive disk (i.e., M disk /M * 0.1), whose evolution we follow for a further 30-40 kyr. Gravitational instabilities result in the fragmentation and the formation of gravitationally bound clumps within the disk. The accretion of these fragments by the host protostar produces accretion and luminosity bursts on the order of 10 6 L ⊙ . Within the cluster, we show that a simultaneity of such events across several protostellar cluster members can elevate the cluster luminosity to 5-10× greater than expected, and that the cluster spends ∼ 15% of its star-forming history at these levels. This enhanced luminosity effect is particularly enabled in clusters of modest size with ≃ 10-20 members. In one such instance, we identify a confluence of burst events that raise the luminosity to nearly 1000× greater than the cluster mean luminosity, resulting in L > 10 8 L ⊙ . This phenomenon arises solely through the gravitational-instability-driven episodic fragmentation and accretion that characterizes this early stage of protostellar evolution.

We model the mass accretion rateṀ to stellar mass M * correlation that has been inferred from observations of intermediate to upper mass T Tauri stars-that iṡ M ∝ M 1.3±0.3 * . We explain this correlation within the framework of quiescent disk evolution, in which accretion is driven largely by gravitational torques acting in the bulk of the mass and volume of the disk. Stresses within the disk arise from the action of gravitationally driven torques parameterized in our 1D model in terms of Toomre's Q criterion. We do not model the hot inner sub-AU scale region of the disk that is likely stable according to this criterion, and appeal to other mechanisms to remove or redistribute angular momentum and allow accretion onto the star. Our model has the advantage of agreeing with large-scale angle-averaged values from more complex nonaxisymmetric calculations. The model disk transitions from an early phase (dominated by initial conditions inherited from the burst mode of accretion) into a later self-similar mode characterized by a steeper temporal decline inṀ . The models effectively reproduce the spread in mass accretion rates that have been observed for protostellar objects of 0.2 M ⊙ ≤ M * ≤ 3.0 M ⊙ , such as those found in the ρ Ophiuchus and Taurus star forming regions. We then compare realistically sampled populations of young stellar objects produced by our model to their observational counterparts. We find these populations to be statistically coincident, which we argue is evidence for the role of gravitational torques in the late time evolution of quiescent protostellar disks.

We present simulation results for the formation and long-term evolution of a primordial protostellar disk harbored by a first star. Using a 2+1D nonaxisymmetric thin disk numerical simulation, together with a barotropic relation for the gas, we are able to probe ∼20 kyr of the disk's evolution. During this time period we observe fragmentation leading to loosely bound gaseous clumps within the disk. These are then torqued inward and accreted onto the growing protostar, giving rise to a burst phenomenon. The luminous feedback produced by this mechanism may have important consequences for the subsequent growth of the protostar.

Abstract.Modeling of the self-consistent formation and evolution of disks as a result of prestellar core collapse reveals an intense early phase of recurrent gravitational instability and clump formation. These clumps generally migrate inward due to gravitational interaction with trailing spiral arms, and can be absorbed into the central object. However, in situations of multiple clump formation, gravitational scattering of clumps can result in the ejection of a low mass clump. These clumps can then give rise to free-floating low mass stars, brown dwarfs, or even giant planets. Detailed modeling of this process in the context of present-day star formation reveals that these clumps start out essentially as Larson first cores and grow subsequently by accretion. In the context of Pop III star formation, preliminary indications are that the disk clumps may also be of low mass. This mechanism of clump formation and possible ejection provides a channel for the formation of low mass objects in the first generation of stars.