---
title: Gravitoturbulence in Astrophysical Disks
url: https://www.emergentmind.com/topics/gravitoturbulence
type: topic
---

# Gravitoturbulence in Astrophysical Disks

Gravitoturbulence refers to the self-regulated turbulent state arising in sufficiently massive, differentially rotating astrophysical disks when gravitational instability (GI) is present but rapid fragmentation is averted by sufficiently slow radiative cooling. This phenomenon is fundamental to the angular momentum transport, disk evolution, planetesimal formation, and dynamo activity in a range of systems, notably protoplanetary and circumplanetary disks, as well as in outer regions of AGN disks. At its core, gravitoturbulence results from the interplay between disk self-gravity, differential rotation, and the cooling timescale, producing persistent spiral density waves, large-scale and small-scale turbulence, and a rich set of secondary phenomena.

## 1. Instability Criteria, Onset, and Turbulent State

The onset of gravitoturbulence is governed by the Toomre parameter:
$$
Q \equiv \frac{c_s \kappa}{\pi G \Sigma}
$$
where $c_s$ is the local sound speed, $\kappa$ the epicyclic (orbital) frequency (typically $\Omega$ for Keplerian disks), $G$ the gravitational constant, and $\Sigma$ the surface density. Disks are unstable to axisymmetric self-gravity when $Q \lesssim 1$ [1812.05644][1909.08883][1405.3291][1605.01873][1402.1180]. Non-axisymmetric gravitational instability persists for $Q\lesssim1.4$–1.7 depending on disk thickness and other effects.

Instability alone is not sufficient to sustain gravitoturbulence. For this, radiative cooling ($t_{\mathrm{cool}}$) must be slow enough to prevent immediate fragmentation into bound clumps. The key dimensionless parameter is the cooling time in orbital units, $\beta \equiv t_{\mathrm{cool}}\Omega$. If $\beta$ is too short, fragmentation dominates; for moderate or long $\beta$, the disk self-regulates into a quasi-steady, turbulent state with $Q\sim1$ [1812.05644][1909.08883].

Gravitoturbulence self-consistently balances shock heating from spiral arms and compressional work against imposed cooling, determining both the amplitude of turbulence and effective angular momentum transport.

## 2. Statistical, Morphological, and Spectral Properties

Gravitoturbulent disks are characterized by strong, large-scale, transient spiral density waves and a multiscale turbulent cascade.

**Key features:**
- **Spiral Wakes:** Dominant, non-axisymmetric spiral density features with measured pitch angles $i \sim 12^\circ$–$16^\circ$ and azimuthal wavelengths $\lambda_y\sim 60 H$ ($H$ is disk scale height) [1812.05644][2102.00775].
- **Spatial Power:** On large scales ($\lambda \gtrsim H$), turbulence exhibits quasi-2D Kolmogorov-like spectra $E(k) \propto k^{-5/3}$ in the velocity field, with in-plane motions dominating. On smaller scales ($\lambda \lesssim 0.3 H$; $kH \gtrsim 20$), the flow becomes nearly isotropic [1812.05644][1706.06537].
- **Solenoidal Dominance:** Decomposition of the turbulent velocity shows $>$70% of the power is in solenoidal (incompressible) modes at all resolved scales, consistent with a turbulent cascade driven by inertial waves [1812.05644][1706.06537].
- **Locality:** In thin disks ($H/R\lesssim 0.1$), angular-momentum and energy transport is largely local in both time and space when time-averaged over several orbits or scale heights. The Shakura–Sunyaev $\alpha$ parameter, measured as $\alpha=\langle$stress$\rangle/\langle P\rangle$, closely tracks the inverse of the effective cooling time: $\alpha \approx 0.8\,\beta_{\mathrm{heat}}^{-1}$ [2504.18751][2102.00775].
- **Clumpy Morphology:** In saturated gravito-turbulence, spiral structure is dominated by recurrent, high-contrast clumps at corotation, stemming from nonlinear mode coupling across a range of azimuthal wavenumbers. These structures have typical scales $\sim H$ in both $R$ and $\phi$, with lifetimes of $1$–$2$ orbits [2504.18751].

## 3. Fragmentation, Convergence, and Domain Size

Fragmentation occurs when cooling is sufficiently rapid ($\beta \lesssim 3$ in 3D; higher in some 2D studies), so that heating via shocks cannot balance energy loss, and overdensities collapse gravitationally on a dynamical timescale [1812.05644][1909.08883]. However, for marginally longer cooling ($\beta\sim 4$–$5$), stochastically triggered fragmentation may occur, typically at highest numerical resolutions. For $\beta \gtrsim 10$, disks are robustly stable against fragmentation for simulation times of hundreds of orbits [1812.05644][1909.08883].

Numerically, convergence of the fragmentation boundary and turbulent properties requires large simulation domains and high resolution. In 3D shearing-box setups, the azimuthal box size must satisfy $L_y \gtrsim 60\,H$ in order to resolve the dominant spiral wavelength and avoid artificial burstiness or suppressed turbulence. Similarly, in 2D, box sizes $L\gtrsim40\,H$ are required [1812.05644][1706.06537]. Adequate resolution, generally $>40$ zones per $H$, is necessary to properly capture the critical scale-height and Toomre length turbulent components that mediate fragmentation [1909.08883].

## 4. Vertical Structure, Parametric Instabilities, and Energy Transport

Gravitoturbulent velocity fluctuations are nearly vertically uniform, rising by only a factor of $\sim 2$ from the midplane to several scale heights above, in stark contrast with MRI-driven turbulence where vertical gradients in velocity amplitude can exceed factors of $15$ [1405.3291][1802.06620]. Such vertical uniformity arises because even non-self-gravitating material at $|z|\gg0$ is strongly forced by midplane overdensities through the long-range nature of gravity.

Three-dimensional simulations reveal the existence of poloidal, baroclinically-generated rolls associated with spiral wakes, with vertical velocities $|v_z|\sim 0.2$–$0.5\,c_s$ and typical widths $\sim 1$–$2 H$ [1802.06620]. Parametric instabilities in the upper layers, triggered by large-scale epicyclic oscillations, drive nearly isotropic inertial-wave turbulence at $z \sim 0.5$–$1\,H$, seen as fine-scale, high-wavenumber kinetic motions [1706.06537]. These features enhance vertical mixing, support the disk against collapse, and can hinder inward solid sedimentation.

## 5. Gravitoturbulence and Dust Evolution

Gravitoturbulent structure critically alters the collision velocities, spatial distribution, and growth conditions of solid particles:
- **Diffusion and Forcing:** Small dust grains ($\mathrm{St}\ll1$; Stokes number) couple to the turbulent gas and diffuse with $D_p \sim D_g\sim0.01$–$0.02\,c_s^2/\Omega$, with relative velocities scaling as $\Delta v \sim 0.4\,\sqrt{\mathrm{St}}\,c_s$ [1812.05644][1603.06575].
- **Large Solids:** Particles with $\mathrm{St}\gtrsim 1$ undergo strong stochastic gravitational stirring, leading to elevated radial diffusion and high eccentricities, raising destructive collision rates [1603.06575].
- **Filament Formation:** Intermediate-size particles ($\mathrm{St}\sim 0.1$–$1$) are uniquely susceptible to being swept into narrow, dense filaments by the interaction of spiral gravitational forcing and gas drag. Within these filaments, local dust-to-gas ratios can reach $10$–$10^3$ [1603.06575][2204.13310].
- **Collapse to Planetesimals:** When the local dust density in filaments exceeds the Hill/Roche limit, direct gravitational collapse into planetesimals or planetary embryos is triggered, with bound objects up to several $M_\oplus$ for $\mathrm{St}=1$ at distances $\sim50$~AU [2204.13310]. Dust-gas backreaction modulates clump mass and multiplicity but does not prevent embryo formation.
- **Vertical Mixing:** Poloidal rolls and inertial-wave turbulence loft small grains into the disk atmosphere ($H_d\sim 0.5 H$ for $D_z\sim0.2 c_s H$), impeding efficient midplane sedimentation and thereby potentially delaying or modifying planetesimal formation [1802.06620].

## 6. Interaction with Other Angular Momentum Transport Mechanisms and Dynamos

Gravitoturbulence coexists and interacts with magnetorotational instability (MRI) and large-scale magnetic fields.
- **Angular Momentum Transport:** In pure GI, spiral wakes and correlated overdensities drive gravitational $+$ Reynolds stresses, measured via the effective $\alpha$-parameter, typically $\alpha\sim 0.01$–$0.1$ for $\beta=10$ [1812.05644][2102.00775][1605.01873].
- **MRI Suppression:** Strong gravito-turbulence suppresses canonical zero-net-flux MRI for $\tau_c\lesssim100\,\Omega^{-1}$, yet can sustain its own spiral-wave dynamo, amplifying magnetic fields to near-equipartition levels via midplane vortex stretching and $\Omega$-effects [1709.06845][2206.03917].
- **Dynamo Activity:** The GI-driven dynamo, distinct from classic mean-field $\alpha$-$\Omega$ mechanisms, relies on turbulent electromotive force arising from velocity-magnetic field correlations induced by spiral wakes and poloidal rolls. The process is efficient at moderate cooling ($\tau_{\mathrm{cool}}\sim 10$–$30$) and moderate magnetic Reynolds numbers ($\mathrm{Rm} \sim 3$–$10$), saturating when Maxwell stresses and energy reach those of gravitational stresses, at which point the background $Q$ and disk scale height increase beyond the instability threshold [2206.03917][1709.06845].

## 7. Observational and Astrophysical Implications

Gravitoturbulence imprints distinctive signatures on disk morphology, kinematics, and solid body evolution:
- **Observational Diagnostics:** Turbulent velocity fields with nearly flat vertical profiles may be revealed via multi-line molecular spectroscopy, distinguishing gravito-turbulence from MRI-driven turbulence [1405.3291]. Filamentary, clumpy spiral structure can masquerade as planets or companions in continuum imaging when telescope resolution is comparable to $H$ [2504.18751].
- **Variability and Accretion:** Stochastic dissipation in transient spiral wakes can cause significant accretion-rate variability on orbital timescales [2102.00775][2504.18751].
- **Planet and Core Formation:** Direct embryo formation in filaments or clumps—implicated for the earliest stages of core and giant planet formation at large orbital separation—proceeds rapidly, providing seeds for later growth via pebble or gas accretion [2204.13310][1603.06575].
- **Circumplanetary and AGN Disks:** In circumplanetary disks, gravitoturbulence dominates transport in the outer disk ($r\gtrsim200R_J$), sets the effective viscosity, prevents fragmentation for high cooling times, and determines the disk mass and temperature structure [1402.1180]. Analogous mechanisms operate in the outer regions of AGN disks, with potentially similar implications for black hole feeding and star formation.

---

**References:**

- [1812.05644] Booth & Clarke, "Characterizing gravito-turbulence in 3D..."  
- [1909.08883] Klee et al., "Closing the gap to convergence of gravitoturbulence in local simulations"  
- [1405.3291] Shi & Chiang, "Gravito-Turbulent Disks in 3D: Turbulent Velocities vs. Depth"  
- [1706.06537] Riols et al., "Gravito-turbulence and the excitation of small-scale parametric instability..."  
- [1802.06620] Riols & Latter, "Spiral density waves and vertical circulation in protoplanetary discs"  
- [2102.00775] Riols et al., "Spiral structures in gravito-turbulent gaseous disks"  
- [2504.18751] Zhu et al., "Global Simulations of Gravitational Instability in Protostellar Disks..."  
- [1603.06575] Shi et al., "Dust dynamics in 2D gravito-turbulent disks"  
- [2204.13310] Baehr et al., "Direct Formation of Planetary Embryos in Self-Gravitating Disks"  
- [1703.10292] Hirose & Shi, "Gravito-turbulence in irradiated protoplanetary discs"  
- [1402.1180] Keith & Wardle, "Accretion in giant planet circumplanetary disks"  
- [1605.01873] Riols & Latter, "Gravitoturbulence in magnetised protostellar discs"  
- [2206.03917] Riols et al., "Gravitoturbulent dynamo in global simulations of gaseous disks"  
- [1709.06845] Riols & Latter, "Magnetorotational instability and dynamo action in gravitoturbulent astrophysical discs"

Source: https://www.emergentmind.com/topics/gravitoturbulence