Testing the locality of transport in self-gravitating accretion discs

Giuseppe Lodato,W.K.M. Rice

DOI: https://doi.org/10.1063/1.1718465

2003-09-22

Abstract:In this paper we examine the issue of characterizing the transport associated with gravitational instabilities in relatively cold discs, discussing in particular under which condition 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 parameterization for the cooling function. Our simulations show that the disc settles in a ``self-regulated'' state, where the axisymmetric stability parameter $Q\approx 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, for disc masses smaller than $0.25M_{\star}$ and aspect ratio $H/r\lesssim 0.1$, transport is determined locally, thus allowing for a viscous treatment of the disc evolution.

Astrophysics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: **How to characterize the angular momentum transfer and energy dissipation caused by gravitational instability, and determine whether this transfer can be described in a local, viscous framework**.

Specifically, the author focuses on whether gravitational instability will provide an effective angular momentum transfer mechanism in self - gravitating accretion discs when the discs are cold. The key to this problem lies in:

1. **Angular momentum transfer mechanism**: In accretion disc theory, a fundamental unknown is the specific physical mechanism responsible for angular momentum transfer and energy dissipation. It is usually assumed that this transfer is dominated by some viscous process, and the α - model (as shown in equation (1)) is used to describe the radial - tangential component of the stress tensor:
\[
T_{r\phi} = \left| \frac{d \ln \Omega}{d \ln r} \right| \alpha \Sigma c_s^2
\]
where \(\Omega\) is the angular velocity of the disc, \(c_s\) is the speed of sound, \(\Sigma\) is the surface density, and \(\alpha\) is a free parameter.

2. **Gravitational instability condition**: For a sufficiently cold disc, the condition for axisymmetric stability (as shown in equation (2)) is:
\[
Q = \frac{c_s \kappa}{\pi G \Sigma} > \bar{Q} \approx 1
\]
where \(\kappa\) is the epicyclic frequency. When \(Q < 1\), the disc may become unstable, causing gravitational instability to generate spiral structures, which may then become the source of angular momentum transfer.

3. **Applicability of local viscous description**: The author hopes to test whether the transfer caused by gravitational instability can be described by a local viscous framework through global three - dimensional SPH simulations (smoothed particle hydrodynamics). In particular, they want to verify the hypothesis proposed by Lin and Pringle, that is, the transfer caused by gravitational instability can be described in a local viscous framework.

Through these simulations, the author hopes to answer the following questions:
- Can gravitational instability effectively heat the disc and bring it to a self - regulated state?
- Under what conditions can the transfer caused by gravitational instability be described by a local viscous model?

Finally, the conclusion of the paper shows that when the disc mass is less than 0.25 times the central star mass and the disc thickness - to - radius ratio \(H/r < 0.1\), the transfer caused by gravitational instability can be well described by a local viscous framework.