Vortex stretching in self-gravitating protoplanetary discs

Horseshoe-shaped brightness asymmetries of several transitional discs are thought to be caused by large-scale vortices. Anticyclonic vortices are efficiently collect dust particles, therefore they can play a major role in planet formation. Former studies suggest that the disc self-gravity weakens vortices formed at the edge of the gap opened by a massive planet in discs whose masses are in the range of 0.01<=M_disc/M_*<=0.1. Here we present an investigation on the long-term evolution of the large-scale vortices formed at the viscosity transition of the discs' dead zone outer edge by means of two-dimensional hydrodynamic simulations taking disc self-gravity into account. We perform a numerical study of low mass, 0.001<=M_disc/M_*<=0.01, discs, for which cases disc self-gravity was previously neglected. The large-scale vortices are found to be stretched due to disc self-gravity even for low-mass discs with M_disc/M_*>=0.005 where initially the Toomre Q-parameter was <=50 at the vortex distance. As a result of stretching, the vortex aspect ratio increases and a weaker azimuthal density contrast develops. The strength of the vortex stretching is proportional to the disc mass. The vortex stretching can be explained by a combined action of a non-vanishing gravitational torque caused by the vortex, and the Keplerian shear of the disc. Self-gravitating vortices are subject to significantly faster decay than non-self-gravitating ones. We found that vortices developed at sharp viscosity transitions of self-gravitating discs can be described by a GNG model as long as the disc viscosity is low, i.e. alpha_dz<=10^-5.Comment: 13 pages, 8 figures, appear in MNRA

Mass accretion rates in self-regulated disks of T Tauri stars

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, disks 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.Comment: Accepted for publication in ApJ Letter

Self-regulated gravitational accretion in protostellar discs

We present a numerical model for the evolution of a protostellar disc that has formed self-consistently from the collapse of a molecular cloud core. The global evolution of the disc is followed for several million years after its formation. The capture of a wide range of spatial and temporal scales is made possible by use of the thin-disc approximation. We focus on the role of gravitational torques in transporting mass inward and angular momentum outward during different evolutionary phases of a protostellar disc with disc-to-star mass ratio of order 0.1. In the early phase, when the infall of matter from the surrounding envelope is substantial, mass is transported inward by the gravitational torques from spiral arms that are a manifestation of the envelope-induced gravitational instability in the disc. In the late phase, when the gas reservoir of the envelope is depleted, the distinct spiral structure is replaced by ongoing irregular nonaxisymmetric density perturbations. The amplitude of these density perturbations decreases with time, though this process is moderated by swing amplification aided by the existence of the disc's sharp outer edge. Our global modelling of the protostellar disc reveals that there is typically a residual nonzero gravitational torque from these density perturbations, i.e. their effects do not exactly cancel out in each region. In particular, the net gravitational torque in the inner disc tends to be negative during first several million years of the evolution, while the outer disc has a net positive gravitational torque. Our global model of a self-consistently formed disc shows that it is also self-regulated in the late phase, so that it is near the Toomre stability limit, with a near-uniform Toomre parameter Q\approx 1.5-2.0. (Abstract abridged).Comment: 9 pages, 9 figures, accepted for publication in MNRA

Self-regulated gravitational accretion in protostellar discs

We present a numerical model for the evolution of a protostellar disc that has formed self-consistently from the collapse of a molecular cloud core. The global evolution of the disc is followed for several million years after its formation. The capture of a wide range of spatial and temporal scales is made possible by use of the thin-disc approximation. We focus on the role of gravitational torques in transporting mass inward and angular momentum outward during different evolutionary phases of a protostellar disc with disc-to-star mass ratio of order 0.1. In the early phase, when the infall of matter from the surrounding envelope is substantial, mass is transported inward by the gravitational torques from spiral arms that are a manifestation of the envelope-induced gravitational instability in the disc. In the late phase, when the gas reservoir of the envelope is depleted, the distinct spiral structure is replaced by ongoing irregular nonaxisymmetric density perturbations. The amplitude of these density perturbations decreases with time, though this process is moderated by swing amplification aided by the existence of the disc's sharp outer edge. Our global modelling of the protostellar disc reveals that there is typically a residual nonzero gravitational torque from these density perturbations, i.e. their effects do not exactly cancel out in each region. In particular, the net gravitational torque in the inner disc tends to be negative during first several million years of the evolution, while the outer disc has a net positive gravitational torque. Our global model of a self-consistently formed disc shows that it is also self-regulated in the late phase, so that it is near the Toomre stability limit, with a near-uniform Toomre parameter Q\approx 1.5-2.0. (Abstract abridged).Comment: 9 pages, 9 figures, accepted for publication in MNRA

Observed luminosity spread in young clusters and Fu Ori stars: a unified picture

The idea that non steady accretion during the embedded phase of protostar evolution can produce the observed luminosity spread in the Herzsprung-Russell diagram (HRD) of young clusters has recently been called into question. Observations of Fu Ori, for instance, suggest an expansion of the star during strong accretion events whereas the luminosity spread implies a contraction of the accreting objects, decreasing their radiating surface. In this paper, we present a global scenario based on calculations coupling episodic accretion histories derived from numerical simulations of collapsing cloud prestellar cores of various masses and subsequent protostar evolution. Our calculations show that, assuming an initial protostar mass \mi \sim 1\,\mjup, typical of the second Larson's core, both the luminosity spread in the HRD and the inferred properties of Fu Ori events (mass, radius, accretion rate) can be explained by this scenario, providing two conditions. First, there must be some variation within the fraction of accretion energy absorbed by the protostar during the accretion process. Second the range of this variation should increase with increasing accretion burst intensity, and thus with the initial core mass and final star mass. The numerical hydrodynamics simulations of collapsing cloud prestellar cores indeed show that the intensity of the accretion bursts correlates with the mass and initial angular momentum of the prestellar core. Massive prestellar cores with high initial angular momentum are found to produce intense bursts characteristic of Fu Ori like events. Our results thus suggest a link between the burst intensities and the fraction of accretion energy absorbed by the protostar, with some threshold in the accretion rate, of the order of 10^{-5}\msolyr, delimitating the transition from "cold" to "hot" accretion. [Abridged]Comment: 23 pages, 5 figures, ApJ accepte

The Burst Mode of Accretion in Primordial Star Formation

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.Comment: 3 pages, 2 figures, to appear in proceedings of First Stars IV meeting (Kyoto, Japan; 2012