No Arabic abstract
In this paper we examine the issue of characterising the transport associated with gravitational instabilities in relatively cold discs, discussing in particular the conditions under which it can be described within a local, viscous framework. We present the results of global, three-dimensional, SPH simulations of self-gravitating accretion discs, in which the disc is cooled using a simple parametrisation for the cooling function. Our simulations show that the disc settles in a ``self-regulated state, where the axisymmetric stability parameter $Qapprox 1$ and where transport and energy dissipation are dominated by self-gravity. We have computed the gravitational stress tensor and compared our results with expectations based on a local theory of transport. We find that, as long as the disc mass is smaller than $0.25M_{star}$ and the aspect ratio $H/Rlesssim 0.1$, transport is determined locally, thus allowing for a viscous treatment of the disc evolution.
In this paper, we extend our previous analysis (Lodato & Rice 2004) of the transport properties induced by gravitational instabilities in cooling, gaseous accretion discs to the case where the disc mass is comparable to the central object. In order to do so, we have performed global, three-dimensional smoothed particle hydrodynamics simulations of massive discs. These new simulations show a much more complex temporal evolution with respect to the less massive case. Whereas in the low disc mass case a self-regulated, marginally stable state (characterized by an approximately constant radial profile of the stability parameter $Q$) is easily established, in the high disc mass case we observe the development of an initial transient and subsequent settling down in a self-regulated state in some simulations, or a series or recurrent spiral episodes, with low azimuthal wave number $m$, in others. Accretion in this last case can therefore be a highly variable process. On the other hand, we find that the secular evolution of the disc is relatively slow. In fact, the time-average of the stress induced by self-gravity results in accretion time-scales much longer than the dynamical timescale, in contrast with previous isothermal simulations of massive accretion discs. We have also compared the resulting stress tensor with the expectations based on a local theory of transport, finding no significant evidence for global wave energy transport.
I review recent progresses in the dynamics and the evolution of self-gravitating accretion discs. Accretion discs are a fundamental component of several astrophysical systems on very diverse scales, and can be found around supermassive black holes in Active Galactic Nuclei (AGN), and also in our Galaxy around stellar mass compact objects and around young stars. Notwithstanding the specific differences arising from such diversity in physical extent, all these systems share a common feature where a central object is fed from the accretion disc, due to the effect of turbulence and disc instabilities, which are able to remove the angular momentum from the gas and allow its accretion. In recent years, it has become increasingly apparent that the gravitational field produced by the disc itself (the discs self-gravity) is an important ingredient in the models, especially in the context of protostellar discs and of AGN discs. Indeed, it appears that in many cases (and especially in the colder outer parts of the disc) the development of gravitational instabilities can be one of the main agents in the redistribution of angular momentum. In some cases, the instability can be strong enough to lead to the formation of gravitationally bound clumps within the disc, and thus to determine the disc fragmentation. As a result, progress in our understanding of the dynamics of self-gravitating discs is essential to understand the processes that lead to the feeding of both young stars and of supermassive black holes in AGN. At the same time, understanding the fragmentation conditions is important to determine under which conditions AGN discs would fragment and form stars and whether protostellar discs might form giant gaseous planets through disc fragmentation.
It is quite likely that self-gravity will play an important role in the evolution of accretion discs, in particular those around young stars, and those around supermassive black holes. We summarise, here, our current understanding of the evolution of such discs, focussing more on discs in young stellar system, than on discs in active galactic nuclei. We consider the conditions under which such discs may fragment to form bound objects, and when they might, instead, be expected to settle into a quasi-steady, self-regulated state. We also discuss how this understanding may depend on the mass of the disc relative to the mass of the central object, and how it might depend on the presence of external irradiation. Additionally, we consider whether or not fragmentation might be stochastic, where we might expect it to occur in an actual protostellar disc, and if there is any evidence for fragmentation actually playing a role in the formation of planetary-mass bodies. Although there are still a number of outstanding issue, such as the convergence of simulations of self-gravitating discs, whether or not there is more than one mode of fragmentation, and quite what role self-gravitating discs may play in the planet formation process, our general understanding of these systems seems quite robust.
Large-scale, dust-trapping vortices may account for observations of asymmetric protoplanetary discs. Disc vortices are also potential sites for accelerated planetesimal formation by concentrating dust grains. However, in 3D discs vortices are subject to destructive `elliptic instabilities, which reduces their viability as dust traps. The survival of vortices in 3D accretion discs is thus an important issue to address. In this work, we perform shearing box simulations to show that disc self-gravity enhances the survival of 3D vortices, even when self-gravity is weak in the classic sense (e.g. with a Toomre $Qsimeq5$). We find a 3D, self-gravitating vortex can grow on secular timescales in spite of the elliptic instability. The vortex aspect-ratio decreases as it strengthens, which feeds the elliptic instability. The result is a 3D vortex with a turbulent core that persists for $sim 10^{3}$ orbits. We find when gravitational and hydrodynamic stresses become comparable, the vortex may undergo episodic bursts, which we interpret as interaction between elliptic and gravitational instabilities. We estimate the distribution of dust particles in self-gravitating, turbulent vortices. Our results suggest large-scale vortices in protoplanetary discs are more easily observed at large radii.
The steady-state structure of a disc with a corona is analyzed when the vertical component of the gravitational force due to the self-gravity of the disc is considered. For the energy exchange between the disc and the corona, we assume a fraction f of the dissipated energy inside the accretion disc is transported to the corona via the magnetic tubes. Analytical solutions corresponding to a prescription for f (in which this parameter directly depends on the ratio of the gas pressure to the total pressure) or free f are presented and their physical properties are studied in detail. We show that the existence of the corona not only decreases the temperature of the disc, but also increases the surface density.The vertical component of the gravitational force due to the self-gravity of the disc decreases the self-gravitating radius and the mass of the fragments at this radius. However, as more energy is transported from the disc to the corona, the effect of the vertical component of the gravitational force due to the self-gravity of the disc on the self-gravitating radius becomes weaker, though the mass of the fragments is reduced irrespective of the amount of the energy exchange from the disc to the corona.