We present the latest development of the disk gravitational instability and fragmentation model, originally introduced by us to explain episodic accretion bursts in the early stages of star formation. Using our numerical hydrodynamics model with impr
oved disk thermal balance and star-disk interaction, we computed the evolution of protostellar disks formed from the gravitational collapse of prestellar cores. In agreement with our previous studies, we find that cores of higher initial mass and angular momentum produce disks that are more favorable to gravitational instability and fragmentation, while a higher background irradiation and magnetic fields moderate the disk tendency to fragment. The protostellar accretion in our models is time-variable, thanks to the nonlinear interaction between different spiral modes in the gravitationally unstable disk, and can undergo episodic bursts when fragments migrate onto the star owing to the gravitational interaction with other fragments or spiral arms. Most bursts occur in the partly embedded Class I phase, with a smaller fraction taking place in the deeply embedded Class 0 phase and a few possible bursts in the optically visible Class II phase. The average burst duration and mean luminosity are found to be in good agreement with those inferred from observations of FU-Orionis-type eruptions. The model predicts the existence of two types of bursts: the isolated ones, showing well-defined luminosity peaks separated with prolonged periods (~ 10^4 yr) of quiescent accretion, and clustered ones, demonstrating several bursts occurring one after another during just a few hundred years. Finally, we estimate that 40%--70% of the star-forming cores can display bursts after forming a star-disk system.
We demonstrate that the formation of collapsing cores in subcritical clouds is accelerated by nonlinear flows, by performing three-dimensional non-ideal MHD simulations. An initial random supersonic (and trans-Alfvenic) turbulent-like flow is input i
nto a self-gravitating gas layer that is threaded by a uniform magnetic field (perpendicular to the layer) such that the initial mass-to-flux ratio is subcritical. Magnetic ambipolar diffusion occurs very rapidly initially due to the sharp gradients introduced by the turbulent flow. It subsequently occurs more slowly in the traditional near-quasistatic manner, but in regions of greater mean density than present in the initial state. The overall timescale for runaway growth of the first core(s) is several times, 10^6 yr, even though previous studies have found a timescale of several times, 10^7 yr when starting with linear perturbations and similar physical parameters. Large-scale supersonic flows exist in the cloud and provide an observationally testable distinguishing characteristic from core formation due to linear initial perturbations. However, the nonlinear flows have decayed sufficiently that the relative infall motions onto the first core are subsonic, as in the case of starting from linear initial perturbations. The ion infall motions are very similar to those of neutrals; however, they lag the neutral infall in directions perpendicular to the mean magnetic field direction and lead the neutral infall in the direction parallel to the mean magnetic field.
We have studied numerically the evolution of protostellar disks around intermediate and upper mass T Tauri stars (0.25 M_sun < M_st < 3.0 M_sun) that have formed self-consistently from the collapse of molecular cloud cores. In the T Tauri phase, disk
s settle into a self-regulated state, with low-amplitude nonaxisymmetric density perturbations persisting for at least several million years. Our main finding is that the global effect of gravitational torques due to these perturbations is to produce disk accretion rates that are of the correct magnitude to explain observed accretion onto T Tauri stars. Our models yield a correlation between accretion rate M_dot and stellar mass M_st that has a best fit M_dot propto M_st^{1.7}, in good agreement with recent observations. We also predict a near-linear correlation between the disk accretion rate and the disk mass.
We employ the first fully three-dimensional simulation to study the role of magnetic fields and ion-neutral friction in regulating gravitationally-driven fragmentation of molecular clouds. The cores in an initially subcritical cloud develop gradually
over an ambipolar diffusion time while the cores in an initially supercritical cloud develop in a dynamical time. The infall speeds on to cores are subsonic in the case of an initially subcritical cloud, while an extended (ga 0.1 pc) region of supersonic infall exists in the case of an initially supercritical cloud. These results are consistent with previous two-dimensional simulations. We also found that a snapshot of the relation between density (rho) and the strength of the magnetic field (B) at different spatial points of the cloud coincides with the evolutionary track of an individual core. When the density becomes large, both relations tend to B propto rho^{0.5}.
