---
title: Massive Self-Gravitating Accretion Discs
url: https://www.emergentmind.com/topics/massive-self-gravitating-accretion-discs-smds
type: topic
---

# Massive Self-Gravitating Accretion Discs

Searching arXiv for the cited SMD papers to ground the article in the literature.
Massive self-gravitating accretion discs are accretion discs in which the disc’s own gravity is dynamically important, so that self-gravity competes with shear, pressure, cooling, and, in some settings, irradiation, magnetic transport, or binary torques. A common benchmark for “massive” is a disc-to-central-object mass ratio of order \(q \equiv M_d/M_* \approx 0.1\), although several regimes discussed in the literature extend to \(q \sim 0.5\) or higher [1101.2448; 1602.08390]. In such discs, gravitational instability is commonly diagnosed with the Toomre parameter \(Q = c_s \kappa / (\pi G \Sigma)\), with \(\kappa \simeq \Omega\) in Keplerian flows; when \(Q\) approaches unity, spiral structure, gravito-turbulence, non-local torques, and, under sufficiently rapid cooling, fragmentation into bound objects may occur [1108.1194; 1602.08390]. Across protoplanetary discs, discs around massive young stars, active galactic nuclei, and circumbinary black-hole discs, the central problem is how self-gravity modifies transport, thermal balance, morphology, and collapse thresholds [1108.1194; 1202.6063; 1503.05099; 1608.05539].

## 1. Dynamical definition and stability criteria

The standard local diagnostic of gravitational instability in an accretion disc is the Toomre parameter,
\[
Q = \frac{c_s \kappa}{\pi G \Sigma},
\]
where \(c_s\) is the sound speed, \(\kappa\) the epicyclic frequency, \(G\) the gravitational constant, and \(\Sigma\) the surface density. For razor-thin discs, axisymmetric linear instability occurs if \(Q < 1\); in Keplerian discs, \(\kappa = \Omega\) [1108.1194]. The same criterion is used broadly across the SMD literature, with non-axisymmetric structure typically appearing when \(Q \lesssim 1.5\)–\(1.7\) or, more generally, when the disc approaches a marginally self-regulated state with \(Q \sim 1\) [1101.2448; 1602.08390].

The distinction between local linear instability and nonlinear saturated states is important. In two-dimensional local shearing-sheet calculations with irradiation, quasi-steady self-gravitating states saturate at \(Q \simeq 1.8\)–\(1.9\), whereas comparable three-dimensional discs typically saturate near \(Q \simeq 1.1\) because finite thickness dilutes self-gravity by \(\approx (1+kH)\) for the most unstable modes with \(k \simeq 1/H\) [1108.1194]. This implies that quoted threshold values are model-dependent, especially with respect to dimensionality and vertical structure.

Cooling introduces a second control parameter. In the standard \(\beta\)-cooling notation,
\[
\tau_c = \frac{\beta}{\Omega},
\]
and fragmentation occurs if cooling is sufficiently rapid, conventionally written as
\[
\tau_c < \frac{\beta_{\rm crit}}{\Omega}
\quad \text{or} \quad
\beta < \beta_{\rm crit}.
\]
The review literature emphasizes that the fragmentation boundary is more physically interpreted as a maximum sustainable self-gravitating stress, often near \(\alpha \simeq 0.06\), rather than as a single universal \(\beta_{\rm crit}\) [1602.08390]. In two-dimensional irradiated local simulations with \(\gamma = 1.6\), however, explicit thresholds are reported: \(\beta_{\rm crit} \simeq 8\) without irradiation, declining to \(\beta_{\rm crit} \simeq 4\) at the strongest irradiation studied [1108.1194].

A further distinction concerns whether self-gravity remains local. Global three-dimensional radiative SPH calculations show that low-mass, geometrically thin discs can be represented by a local \(\alpha\)-parametrization when \(H/R < 0.1\), generally the case for \(q < 0.5\), but this approximation breaks down as \(q\) approaches unity and low-\(m\) global spirals dominate [1008.1547]. This suggests that “massive” is not a single regime but a continuum from local gravito-turbulence to strongly non-local, wave-mediated transport.

## 2. Thermal balance, effective stress, and transport laws

In non-fragmenting SMDs, transport is often described by an effective viscous stress whose magnitude is set by local thermal balance. Without irradiation, local equilibrium yields
\[
\alpha = \frac{4}{9 \gamma (\gamma - 1)\tau_c \Omega}.
\]
With irradiation, the self-gravitating component of the stress is reduced by the fraction of pressure support provided by the irradiation floor,
\[
\alpha \approx \frac{4}{9 \gamma (\gamma - 1)\Omega \tau_c}
\left(1 - \frac{\langle \Sigma \rangle c_{so}^2}{\langle \Sigma c_s^2\rangle}\right)
\approx
\frac{4}{9 \gamma (\gamma - 1)\Omega \tau_c}
\left(1 - \frac{Q_{\rm irr}^2}{Q_{\rm sat}^2}\right),
\]
with \(Q_{\rm sat} \simeq 1.9\) used in the irradiated shearing-sheet study [1108.1194].

That same study measures the transport coefficient directly from Reynolds and gravitational stresses as
\[
\alpha = \frac{2}{3 \langle \Sigma c_s^2 \rangle}
\left(\langle G_{xy}\rangle + \langle H_{xy}\rangle\right),
\]
where
\[
\langle H_{xy}\rangle = \langle \Sigma u_x u_y \rangle,
\qquad
\langle G_{xy}\rangle = \sum_k \frac{\pi G k_x k_y |\Sigma_k|^2}{|\mathbf{k}|^3}.
\]
The main result is that, in non-fragmenting runs, \(\alpha\) remains set by local thermal equilibrium even for very long cooling times, up to \(\tau_c = 240\,\Omega^{-1}\) or \(\beta = 240\) [1108.1194]. This is one of the clearest statements of the gravito-turbulent closure in irradiated local models.

Other works frame the same physics in alternative transport prescriptions. One-dimensional self-similar models for geometrically thin, viscous self-gravitating discs reduce the full vertically integrated system to a nonlinear advection–diffusion equation for \(\Omega(r,t)\),
\[
-r^4 \Omega^2 \frac{\partial \Omega}{\partial t}
=
\frac{\partial}{\partial r}
\left[\nu r^3 \Omega^3 x (2x+3)\right],
\qquad
x \equiv \frac{\partial \ln \Omega}{\partial \ln r},
\]
with the surface density coupled through
\[
2\pi G \Sigma = r \Omega^2 (2x+3).
\]
Three viscosity prescriptions are considered there: a self-regulated “LP” form, a DSB \(\beta\)-viscosity \(\nu = \beta r^2 \Omega\), and an RZ form \(\nu = \beta |x| r^2 \Omega\) [1503.05099]. In these similarity solutions, the outer rotation-law exponent \(n\) is the key control parameter; flatter rotation laws at large radii yield higher accretion rates, and fully self-gravitating discs evolve faster than nearly Keplerian discs [1503.05099].

The transport locality question remains central. Global radiative SPH simulations find that for \(q \lesssim 0.5\), \(\alpha_{\rm total}(r)\) is consistent with the cooling-based \(\alpha_{\rm cool}(r)\), whereas for \(q \sim 1\)–\(1.5\), global \(m=2\) spirals produce non-local energy flux and systematic departures from local thermodynamic equilibrium [1008.1547]. A plausible implication is that the \(\alpha\)-closure is robust in the thin-disc, modest-\(q\) limit, but should be interpreted as a time-averaged phenomenology in the most massive discs.

## 3. Fragmentation, stochasticity, and numerical convergence

Fragmentation in SMDs is usually defined by the sustained formation of overdense, long-lived clumps. In the irradiated local simulations, fragmenting runs are identified by clumps with densities \(>100\times\) the mean that survive for many cooling times, while the disc-averaged \(Q\) rises rather than settling [1108.1194]. In protostellar disc reviews, fragmentation is located predominantly in the outer disc, typically beyond \(\sim 30\)–\(50\) au under standard opacities, and rarely inside \(\sim 10\)–\(20\) au even when opacity is reduced [1602.08390].

Irradiation weakens instability but does not generally remove fragmentation. As irradiation is increased toward the quasi-linear stability threshold, the critical cooling time decreases only by about a factor of two, from \(\beta_{\rm crit} \simeq 8\) to \(\simeq 4\) for \(\gamma = 1.6\), and the fragmentation boundary in the \((Q_{\rm irr},\beta)\) plane does not follow contours of constant \(\alpha\) [1108.1194]. The interpretation advanced there is that stronger irradiation raises \(c_s\) and \(Q\), but if other transport channels are weak, mass builds up until self-gravity turns on and fragmentation ensues anyway [1108.1194].

A long-standing controversy concerns convergence of the fragmentation threshold. SPH resolution studies argue that the measured threshold depends on \(h/H\), with finer resolution needed to fragment at larger \(\beta\). Interpreting the Meru–Bate results as resolution-driven trends, Lodato and Clarke infer tentative convergence at \(N \approx\) a few \(\times 10^6\)–\(10^7\) particles for \(q \approx 0.1\), with a converged threshold around \(\beta_{\rm crit} \approx 10\)–\(15\), and a preferred estimate \(\beta_{\rm crit} \approx 15\Omega^{-1}\) [1101.2448]. Two numerical origins are proposed: artificial-viscosity heating, quantified by
\[
\alpha_{\rm art} = \frac{1}{10}\alpha_{\rm SPH}\left(\frac{h}{H}\right),
\]
and smoothing of density peaks over finite \(h\), leading to a saturating relation
\[
\beta_{\rm res} = \frac{\beta_0}{(1+a h/H)^2},
\]
with best-fit \(\beta_0 \approx 14.7\) and \(a \approx 1.77\) [1101.2448].

The review literature treats this issue cautiously. It notes reports of non-convergence and proposed fragmentation at \(\beta \sim 30\), but also cites high-resolution calculations showing no fragmentation for \(\beta = 10\) and analyses indicating that simulations with \(\alpha < 0.06\) typically do not fragment [1602.08390]. The robust consensus is narrower than the numerical debate: transport relations and the perturbation–stress connection are comparatively stable, whereas exact thermodynamic thresholds are sensitive to resolution and implementation [1101.2448; 1602.08390].

Stochastic fragmentation has also been examined explicitly. Two-dimensional SPH studies of the waiting-time distribution between strong shocks find an exponential law,
\[
P(\Delta t) \approx \lambda e^{-\lambda \Delta t},
\qquad
S(\Delta t)=e^{-\lambda \Delta t},
\]
with \(\lambda \approx 0.2\Omega\), most probable waiting times of \(\approx 4\)–\(5\,t_{\rm dyn}\), and negligible probability of shock-free intervals longer than \(\gtrsim 100\,t_{\rm dyn}\) [1510.05290]. Combining this with \(t_{\rm cool}=\beta \Omega^{-1}\), the survival probability of a contracting clump scales as \(e^{-A\beta}\) with \(A \approx 0.2\), leading to the conclusion that stochastic fragmentation cannot move the fragmentation radius inward by more than \(\sim 20\%\) [1510.05290]. This suggests that stochasticity does not qualitatively alter the standard view that direct gravitational collapse is largely confined to the outer disc.

## 4. Morphology, spectra, and the universality of gravito-turbulent structure

In non-fragmenting SMDs, the turbulent structure is not arbitrary. The irradiated local simulations show that the power spectrum of surface-density perturbations is uniquely set by \(\alpha\), not by \(\beta\) or irradiation level. All spectra share the same shape and peak at
\[
kL/(2\pi) \simeq 7,
\]
implying that most of the stress arises from wavelengths much smaller than the local box size, and simulations with the same \(\alpha\) but different \((\beta,Q_{\rm irr})\) have nearly identical spectra [1108.1194]. The perturbation amplitude scales as
\[
\alpha \propto \left\langle (\delta \Sigma/\Sigma)^2 \right\rangle,
\]
generalizing earlier results that \(\langle \delta\Sigma/\Sigma\rangle \propto \beta^{-1/2}\) [1108.1194].

Global radiative SPH calculations further distinguish between local, high-\(m\) structure and global low-\(m\) spirals. Thin low-\(q\) discs distribute power over higher azimuthal mode number, with low variability in temperature and \(Q\), while massive thick discs are dominated by \(m=2\) spirals, transient bursts, and outward wave-mediated energy transport [1008.1547]. In this sense, morphology is itself a diagnostic of whether the disc is in the local gravito-turbulent regime or the non-local global regime.

In massive protostellar systems, disc asymmetry feeds back dynamically on the central star. Radiation-hydrodynamic simulations including the indirect potential from stellar wobbling find that this backreaction makes discs smaller and rounder, delays and reduces fragmentation, changes angular momentum redistribution, and suppresses gaseous clump ejection [2405.19905]. Without wobbling, fragmentation begins between \(R \simeq 300\)–\(500\) au by \(t \simeq 19.85\) kyr and becomes violent by \(t \simeq 23.25\) kyr, with many clumps out to \(\gtrsim 800\)–\(1000\) au and later to \(\simeq 1500\) au. With wobbling, fragmentation is delayed until \(\simeq 23.25\) kyr and remains milder, with only one migrating clump at \(27.55\) kyr and no clump ejections reported [2405.19905]. This suggests that non-axisymmetric backreaction can act as an additional self-regulation channel in massive protostellar SMDs.

The same work links morphology to observability through synthetic millimetre imaging. Post-processing with RADMC-3D and CASA shows that wobbling models yield images in better agreement with ALMA observations of AFGL 4176 mm1, G17.64+0.16, and G353.273 than fixed-star models [2405.19905]. The significance is not that self-gravity simply produces spirals, but that the detailed nonlinear morphology of SMDs depends on feedbacks internal to the star–disc system.

## 5. Astrophysical settings

SMDs occur in several astrophysical environments, with the same basic instability criteria but different thermal and dynamical consequences.

In protoplanetary discs, self-gravity can dominate angular momentum transport at large radii. Under self-luminous conditions, fragmentation is expected beyond \(r \gtrsim 70\)–\(100\) au; in more general protostellar-disc discussions, a practical fragmentation zone of \(\gtrsim 30\)–\(50\) au is emphasized under realistic opacities [1108.1194; 1602.08390]. Early, massive protostellar discs with \(M_{\rm disc}\sim 0.1\,M_\odot\) around a solar-mass star and \(Q \gtrsim 1.5\) are globally stable but still self-gravity-modified. Radiation-hydrodynamic SPH calculations of \(1\,M_J\) seeds inserted into such discs show an initial rapid inward migration phase with \(\tau_{\rm mig}\sim (1\)–\(2)\times 10^4\) yr, followed by gap opening and either slower inward Type II–like migration or outward migration if the gap edges become gravitationally unstable [1804.00583]. Without radiative feedback from the protoplanet, gap edges reach \(Q \sim 1.5\)–2 and drive outward migration; with feedback, edges remain at \(Q \gtrsim 3\) and migration remains inward [1804.00583]. This is not fragmentation of the disc itself, but it shows how proximity to self-gravity changes the migration and mass-growth pathways of embedded objects.

The solids budget is also altered in self-gravitating protostellar discs. Quasi-steady \(Q \simeq 1\) models with radiative cooling and no irradiation indicate that mm-sized pebbles are fragmentation-limited to \(St \lesssim 0.06\), that midplane pebble-to-gas density ratios are typically \(\sim 10^{-3}\), and that the streaming instability is therefore generally suppressed [1902.05385]. By contrast, GI fragments with initial masses \(\gtrsim 3\,M_{\rm Jup}\) can accrete pebbles efficiently while continuing to migrate, since they open gaps in the pebble component but generally fail to open gas gaps [1902.05385]. A plausible implication is that early self-gravitating phases redistribute solids in a way that imprints later core-accretion conditions.

Around massive young stars, SMDs are both a transport engine and an observational challenge. Semi-analytic \(Q=2\) models of candidate discs show that continuum-based masses can underpredict true masses by factors of \(2\)–\(5\) because high optical depth suppresses mm emission [1608.05539]. In G11.92−0.61 MM1 and NGC 6334 I(N) SMA1b, self-gravitating models match observed continuum masses within a factor \(<1.5\) and predict outer-disc fragmentation into low-mass stellar companions, with fragment masses typically \(\ge 0.08\,M_\odot\) [1608.05539]. By contrast, AFGL 4176 mm1 and IRAS 16547−4247 are interpreted there as gravitationally stable because irradiation keeps \(Q\) above the instability regime [1608.05539]. The same theme appears in radiation-hydrodynamic collapse calculations, where self-gravity-generated torques in massive circumstellar discs sustain mean stellar accretion rates of \(\sim 10^{-3}\,M_\odot\,{\rm yr}^{-1}\), comparable to explicit \(\alpha\)-disc models, even though the three-dimensional accretion is episodic [1102.4090].

In active galactic nuclei and related massive black-hole environments, SMDs can be both transport layers and star-forming reservoirs. One-dimensional similarity solutions applied to AGN conclude that flatter outer rotation laws yield higher accretion rates and that fully self-gravitating discs evolve faster than nearly Keplerian discs, a point linked there to supermassive black-hole growth and quasar evolution [1503.05099]. The review literature likewise notes that AGN discs may self-regulate or fragment to form stars depending on cooling conditions [1602.08390].

Circumbinary massive black-hole discs provide a distinct dynamical realization. Three-dimensional SPH simulations of binaries in self-gravitating circumbinary discs with \(M_d/M = 0.2\), \(\beta = 10\), and \(Q \approx 1\)–2 show a leaky cavity, edge overdensities, and spiral arms at \(3a\)–\(5a\), with the net gravitational torque changing sign across corotation at \(r \approx a\) [1202.6063]. The dominant net torque is described there as kinematic and non-resonant, arising from cavity streams on eccentric orbits rather than from a simple resonant tidal picture [1202.6063]. In retrograde circumbinary discs, where outer Lindblad resonances are absent, the cavity is smaller, \(R_{\rm cav,retro}\approx a(1+e)\), yet circular binaries shrink at essentially the same rate as in prograde discs; for eccentric retrograde binaries, both \(a\) decay and \(e\) growth are approximately exponential while the binary remains coplanar [1307.6283]. These results situate SMDs within binary hardening theory and show that disc self-gravity can influence not only transport but also the orbital architecture of compact-object binaries.

## 6. Caveats, controversies, and theoretical limits

Several caveats recur across the SMD literature. Dimensionality matters: two-dimensional razor-thin calculations can overestimate the effective strength of self-gravity, and finite thickness lowers the saturated \(Q\) in three dimensions [1108.1194]. Local shearing-sheet models omit global modes, non-local torques, and infall, whereas global models with realistic radiation often have limited resolution in the inner disc or in collapsing clumps [1108.1194; 1102.4090; 2405.19905].

Cooling prescriptions are another major uncertainty. Idealized \(\beta\)-cooling is analytically and numerically convenient, but realistic discs experience opacity transitions, irradiation, and radiative diffusion. This matters acutely for fragmentation thresholds. Resolution studies explicitly argue that thermodynamic thresholds are less robust than transport diagnostics [1101.2448], and the review literature continues to regard convergence, multiple fragmentation modes, and stochastic fragmentation as open issues even while judging the overall physical picture robust [1602.08390].

The validity of local transport closures has a similarly conditional status. For \(M_d/M_* \lesssim 0.25\), local \(\alpha\)-models can be a good approximation; by \(M_d/M_* \sim 0.5\), transient low-\(m\) spirals become important; and by \(M_d/M_* \gtrsim 1\), only time-averaged local descriptions remain plausible because global wave transport is significant [1602.08390]. Global radiative simulations sharpen this by associating the transition with \(H/R \gtrsim 0.1\) and \(q \gtrsim 0.5\) [1008.1547].

Observational interpretation is also nontrivial. Massive discs may appear light in the continuum because they are optically thick [1608.05539]; spirals may indicate self-gravity only if the outer disc midplane is not overly stabilized by irradiation [2001.06225]; and compact structure in massive protostellar discs may depend sensitively on whether stellar wobbling is included [2405.19905]. The literature therefore repeatedly warns against direct inference of disc mass or GI state from continuum data alone.

Taken together, these caveats do not negate the core framework. Rather, they delimit it. Self-gravitating discs are consistently described as systems that hover near \(Q \sim 1\), transport angular momentum through gravitationally driven stresses or global torques, fragment when cooling overwhelms self-regulation, and shift from local to non-local behaviour as mass ratio and thickness increase [1108.1194; 1602.08390; 1008.1547]. What remains under active refinement is the precise placement of boundaries between these regimes, and the degree to which irradiation, resolution, inflow, magnetic transport, or stellar backreaction modify them.

Source: https://www.emergentmind.com/topics/massive-self-gravitating-accretion-discs-smds