Disk-Jet Lense-Thirring Precession Model

Updated 3 July 2026

The Disk–Jet Lense–Thirring Precession Model is an astrophysical framework describing the global precession of a misaligned accretion disk and its coupled jet due to frame dragging by a spinning Kerr black hole.
It integrates fluid dynamics, general relativity, and jet geometry to explain observable jet wobble and disk warping in systems like TDEs, microquasars, and AGN.
By comparing theoretical precession rates with quasi-periodic oscillations, the model constrains key parameters such as black-hole spin, disk structure, and alignment dynamics.

The Disk–Jet Lense–Thirring Precession Model describes the global, coupled precession of a misaligned accretion disk and its associated relativistic jet under the influence of the Lense–Thirring torque from a spinning compact object, typically a Kerr black hole. This framework provides a rigorous relativistic treatment of jet “wobble” and disk warping, with broad applications from tidal disruption events (TDEs) to microquasars and active galactic nuclei (AGN). The model integrates axisymmetric fluid dynamics, general relativistic torque theory, disk–jet coupling geometry, and observable consequences for multi-wavelength time-domain astrophysics.

1. Lense–Thirring Torque and Precession Frequency

The fundamental physical effect is frame dragging in the Kerr metric, which imparts a nodal precession to the orbital angular momentum of matter on orbits inclined with respect to the black-hole spin. For a test ring at radius $R$ , the local nodal precession (Lense–Thirring) frequency is

$\Omega_{\rm LT}(R) = \frac{2GJ}{c^2 R^3} = 2a_{*}\,\frac{G^2M^2}{c^3 R^3},$

where $J = a_{*}GM^2/c$ is the angular momentum (dimensionless spin $a_{*} \in [-1, 1]$ ). For a differentially rotating disk with surface density $\Sigma(R)$ and Keplerian angular velocity $\Omega(R)$ ,

$\mathbf{L}(R) = \Sigma R^2 \Omega(R) \mathbf{l}(R)$

is the local angular momentum per unit area and $\mathbf{l}(R)$ is the tilt unit vector.

The total Lense–Thirring torque acting on the disk is computed as a moment-weighted average: $T_{\rm tot} = 2\pi \int_{R_{\rm in}}^{R_{\rm out}} \Omega_{\rm LT}(R) \Sigma(R) R^3 \Omega(R) dR,$ with total disk angular momentum

$L_{\rm tot} = 2\pi \int_{R_{\rm in}}^{R_{\rm out}} \Sigma(R) R^3 \Omega(R) dR.$

Under the rigid-body precession hypothesis, the global precession frequency is

$\Omega_{\rm LT}(R) = \frac{2GJ}{c^2 R^3} = 2a_{*}\,\frac{G^2M^2}{c^3 R^3},$ 0

The precession period is then $\Omega_{\rm LT}(R) = \frac{2GJ}{c^2 R^3} = 2a_{*}\,\frac{G^2M^2}{c^3 R^3},$ 1 (Franchini et al., 2015).

2. Disk Structure, Warp Propagation, and Rigidity Conditions

A solid-body precession solution is only valid if warp communication across the disk is rapid compared to differential precession. For geometrically thick disks, as arise during TDE super-Eddington flows, warps propagate as bending waves at $\Omega_{\rm LT}(R) = \frac{2GJ}{c^2 R^3} = 2a_{*}\,\frac{G^2M^2}{c^3 R^3},$ 2. The wave crossing time

$\Omega_{\rm LT}(R) = \frac{2GJ}{c^2 R^3} = 2a_{*}\,\frac{G^2M^2}{c^3 R^3},$ 3

must be less than the differential precession timescale. The solid-body regime is realized when

$\Omega_{\rm LT}(R) = \frac{2GJ}{c^2 R^3} = 2a_{*}\,\frac{G^2M^2}{c^3 R^3},$ 4

with $\Omega_{\rm LT}(R) = \frac{2GJ}{c^2 R^3} = 2a_{*}\,\frac{G^2M^2}{c^3 R^3},$ 5 the Shakura–Sunyaev viscosity parameter. If $\Omega_{\rm LT}(R) = \frac{2GJ}{c^2 R^3} = 2a_{*}\,\frac{G^2M^2}{c^3 R^3},$ 6, warp communication is diffusive; the classic Bardeen–Petterson alignment occurs and global precession is quickly suppressed (Franchini et al., 2015, Stone et al., 2011, Wang et al., 16 Nov 2025). In geometrically thin disks, disk “breaking” and multi-zone precession can occur for strong misalignment and sufficiently small $\Omega_{\rm LT}(R) = \frac{2GJ}{c^2 R^3} = 2a_{*}\,\frac{G^2M^2}{c^3 R^3},$ 7 (Dyda et al., 2020, Shen et al., 2024).

3. Alignment, Damping, and Bardeen–Petterson Transition

Misaligned disks subject to Lense–Thirring torque gradually align with the black-hole spin. Two alignment pathways are dominant:

Wave-driven (bending waves): For thick disks, viscous dissipation of bending waves damps the tilt on a timescale

$\Omega_{\rm LT}(R) = \frac{2GJ}{c^2 R^3} = 2a_{*}\,\frac{G^2M^2}{c^3 R^3},$ 8

with more precise expressions incorporating disk structure and surface density (Franchini et al., 2015).

Diffusive (Bardeen–Petterson, cooling-induced): When accretion cools and the disk becomes thin, the inner portion aligns rapidly over a timescale

$\Omega_{\rm LT}(R) = \frac{2GJ}{c^2 R^3} = 2a_{*}\,\frac{G^2M^2}{c^3 R^3},$ 9

scaling as $J = a_{*}GM^2/c$ 0 (Franchini et al., 2015). Once $J = a_{*}GM^2/c$ 1, warp propagation becomes diffusive, leading to classical Bardeen–Petterson alignment (Franchini et al., 2015, Lei et al., 2012).

4. Disk–Jet Coupling and Observable Precession

If a relativistic jet is tied to the instantaneous orientation of the inner disk angular momentum, the jet will precess with the same period and cone opening angle as the disk. The jet precession forms a conical motion with half-opening angle equal to the initial disk–spin misalignment $J = a_{*}GM^2/c$ 2. This produces quasi-periodic jet position angle swings or flux modulations: $J = a_{*}GM^2/c$ 3 The amplitude of the observed wobble decays exponentially as viscous or cooling-driven alignment proceeds (Franchini et al., 2015, Wang et al., 16 Nov 2025). The phenomenology of short-timescale ( $J = a_{*}GM^2/c$ 4days) jet swings in microquasars and daily to weekly periods in TDEs (e.g., Swift J1644+57) is quantitatively consistent with this framework (Franchini et al., 2015, Stone et al., 2011, Wang et al., 16 Nov 2025, Lei et al., 2012).

5. Generalized Models: No-Hair Extensions and Quadrupole Effects

Advanced models incorporate the general relativistic “no-hair” theorem, adding a black-hole quadrupole precession contribution: $J = a_{*}GM^2/c$ 5 where $J = a_{*}GM^2/c$ 6 is the Kerr quadrupole (Iorio, 24 Feb 2026). The total precession rate then reads

$J = a_{*}GM^2/c$ 7

Inclusion of the quadrupole term breaks the degeneracy between prograde and retrograde configurations, allowing observational constraints on the sign and value of the black-hole spin (Iorio, 24 Feb 2026, Iorio, 2024). For M87*, combined fits of Lense–Thirring and quadrupole dynamics reproduce both jet precession rates and observed sky-projected jet angles (Iorio, 2024, Iorio, 2024).

6. Magnetically Driven Retrograde Precession

In highly magnetized flows (MAD regime), large-scale poloidal magnetic fields can generate retrograde (opposing) precession torques that surpass the classical Lense–Thirring effect. The net precession rate is given by

$J = a_{*}GM^2/c$ 8

where

$J = a_{*}GM^2/c$ 9

and $a_{*} \in [-1, 1]$ 0 is the poloidal field, $a_{*} \in [-1, 1]$ 1 the surface density. Jiang et al. demonstrated that in the MAD state, the net precession can be both retrograde and several times larger in amplitude than the frame-dragging torque alone, modifying both jet direction and precession period. Observational signatures include reversed drift of jet position angle, as potentially seen in M87* (Jiang et al., 18 Jul 2025).

7. Application to Quasi-Periodic Variability and Parameter Constraints

Disk–jet Lense–Thirring precession provides a direct physical explanation for X-ray/radio quasi-periodic oscillations (QPOs) in TDEs, black-hole X-ray binaries, microquasars, and AGN. Comparison between theoretical precession periods,

$a_{*} \in [-1, 1]$ 2

and observed QPO periodicities allows constraints on black-hole spin, disk truncation radius, disk aspect ratio, and alignment properties (Franchini et al., 2015, Massi et al., 2010, Kubota et al., 2023, Ma et al., 23 Jun 2025, Wang et al., 16 Nov 2025, Iorio, 24 Feb 2026, Iorio, 2024). Integrating multi-band timing, spectral fitting of the inner disk radius, and relativistic precession theory yields self-consistent system parameter estimations and supports the global solid-body precession scenario.

References: