Planet-Driven Spirals in Protoplanetary Disks

• Planet-driven spirals are non-axisymmetric density waves in protoplanetary disks, generated by the gravitational torques of an embedded planet.
• Analytic theory and high-resolution simulations link spiral arm features—such as number, pitch angle, and separation—to variations in planet mass and disk thermodynamics.
• Multi-wavelength observations, including ALMA and NIR imaging, confirm that spiral kinematics and morphology serve as precise diagnostics of unseen protoplanets and disk conditions.

Planet-driven spirals are non-axisymmetric density structures in protoplanetary disks, generated by the gravitational interaction between an embedded planet and the circumstellar gas. These spirals encode dynamical information about the planet’s mass, location, and the physical conditions of the disk. Both analytic theory and high-resolution hydrodynamic simulations show that planet-driven spiral morphology—number of arms, pitch angle, contrast, and kinematics—serves as a powerful probe of unseen protoplanets and disk properties. Recent multi-epoch imaging and ALMA kinematic studies have provided direct empirical evidence that such spiral arms generally co-rotate with the planet, and their properties are mutually determined by the planet–disk mass ratio and local thermodynamic structure.

1. Physical Mechanisms and Theoretical Framework

An embedded planet excites spiral density waves via gravitational torques at Lindblad resonances. The classical linear theory decomposes the planet's potential into Fourier modes (azimuthal number $m$). Each $m$ excites density waves near its corresponding Lindblad resonance at radius $r_m^\pm = \left(1\pm 1/m\right)^{2/3} r_p$. The superposition and phase coherence of wave modes—especially the $n=0,1,2$ "wavelets" for each $m$—result in primary and secondary spiral arms. Linear density wave theory predicts the arm phase as

$$\phi_{m,n}(r) = -\mathrm{sgn}(r-r_p)\frac{\pi}{4m} + 2\pi\frac{n}{m} - \int_{r_m^{\pm}}^{r} \frac{\Omega(r')}{c_s(r')}\sqrt{[1-(r'/r_p)^{3/2}]^2 - 1/m^2}\, dr'$$

where $c_s(r)$ is the sound speed and $\Omega(r)$ the Keplerian frequency (Bae et al., 2017, Bae et al., 2017).

The number of arms ($N$) and their morphology are set by the ratio of planet mass to disk thermal mass $M_{\rm th} \equiv (h/r)_p^3 M_*$, where $m$0 is the disk aspect ratio at $m$1. Low-mass planets ($m$2) in cold disks ($m$3) launch multiple (up to 5) arms; higher-mass planets ($m$4) excite only two dominant arms as non-linear wave steepening causes the secondary/tertiary arms to merge with the primary (Bae et al., 2017, Bae et al., 2017).

Nonlinear effects, including spiral shock formation and shock heating, set in for $m$5. The local pitch angle $m$6 and the arm-to-arm separation both increase with planet mass, whereas the number of arms decreases (Bae et al., 2017, Bae et al., 2017, Richert et al., 2015).

2. Spiral Morphology, Scaling Relations, and Diagnostics

Quantitative relations derived from simulations and analytic theory form the basis for “morphological planetology”:

• Number of Arms ($m$7): $m$8 increases as $m$9 and $r_m^\pm = \left(1\pm 1/m\right)^{2/3} r_p$0 decrease; up to five arms possible for small $r_m^\pm = \left(1\pm 1/m\right)^{2/3} r_p$1 and $r_m^\pm = \left(1\pm 1/m\right)^{2/3} r_p$2 (Bae et al., 2017).
• Arm Separations: Primary-to-secondary separation $r_m^\pm = \left(1\pm 1/m\right)^{2/3} r_p$3 (with $r_m^\pm = \left(1\pm 1/m\right)^{2/3} r_p$4), steeper dependencies for colder disks or smaller radii (Bae et al., 2017).
• Pitch Angle ($r_m^\pm = \left(1\pm 1/m\right)^{2/3} r_p$5): Increases with $r_m^\pm = \left(1\pm 1/m\right)^{2/3} r_p$6 and $r_m^\pm = \left(1\pm 1/m\right)^{2/3} r_p$7; always decreases monotonically outwards, i.e., arms "unwind" with increasing $r_m^\pm = \left(1\pm 1/m\right)^{2/3} r_p$8 (Bae et al., 2017, Bae et al., 2017).
• Contrast: In the linear regime, density and brightness contrasts scale as $r_m^\pm = \left(1\pm 1/m\right)^{2/3} r_p$9; in the shock regime, the scaling becomes steeper (Speedie et al., 2022).
• Arm Convergence: In thermally stratified disks, arm loci measured at different heights (e.g., NIR scattered light vs. mm continuum) "converge" on the planet, providing an unambiguous outside/inside diagnostic (Juhasz et al., 2017).

Multiple arms interior to the planet are generic for sub-thermal-mass planets, while exterior multi-armed structures demand colder disks. For each observed disk, joint fitting of arm count, separation, and pitch angle breaks degeneracies in $n=0,1,2$0, $n=0,1,2$1, and planet location (Bae et al., 2017, Bae et al., 2017).

3. Kinematic and Multi-Wavelength Observational Signatures

Observational diagnostics span scattered light, dust continuum, and molecular line emission:

High-Contrast Imaging and Spiral Motion

Long baseline (multi-epoch) polarimetric imaging measures the rotation of spiral arms. In planet-driven models, the arms share a single rigid pattern speed, directly revealing the driver’s orbital radius via Kepler’s law. Examples:

• MWC 758: spiral arm pattern speed $n=0,1,2$2 yr$n=0,1,2$3, corresponding to a planet at $n=0,1,2$4 au; rules out gravitational instability (GI) (Ren et al., 2020, Ren et al., 2018).
• SAO 206462: measured $n=0,1,2$5 deg yr$n=0,1,2$6, implying a planet at $n=0,1,2$7 au precisely coincident with an ALMA dust continuum gap and co-orbital filament—direct linkage between spiral and gap (Xie et al., 2024).

ALMA Continuum and Kinematics

ALMA Band 7 continuum surveys show that spirals driven by planets with $n=0,1,2$8 (Neptune mass at 50 au in typical parameters) are detectable in hours of integration; spiral contrast is maximized in adiabatic disks with slow cooling ($n=0,1,2$9) (Speedie et al., 2022). Dust spirals are non-dust-trapping; their amplitude fades with Stokes number $m$0 as grains decouple from gas (Sturm et al., 2020).

Gas kinematic signatures (e.g., CO channel-map “kinks” or velocity residuals) trace the gas spiral even where dust is not visible. Spiral wakes cause characteristic $m$1–$m$2 m/s line-of-sight velocity shifts, best detected in intermediate-altitude CO isotopologues at inclinations $m$3 (Muley et al., 27 Nov 2025, Wölfer et al., 15 Dec 2025). Buoyancy-driven “vertical” spirals (distinct from Lindblad arms) exhibit tightly wound, predominantly vertical velocity features at the CO emission surface, detectable as $m$4 m/s velocity residuals in ALMA data (Bae et al., 2021).

Thermal Stratification Effects

Observed pitch angles are systematically larger in NIR (upper disk) than in sub-mm (midplane), in agreement with 3D models incorporating a positive vertical temperature gradient. The difference in spiral loci between surface and midplane converges toward the planet’s location (Juhasz et al., 2017).

4. Limitations and Model Uncertainties

The semi-analytic framework, combining linear density wave theory and Burgers’ equation for nonlinear wave propagation, is predictive for $m$5 but underestimates azimuthal width, overpredicts amplitude decay, and lags in shock-front advance at higher masses. For $m$6 and above, fitting procedures based solely on analytic wake shapes and residual minimization do not reliably estimate planet mass, due to pitch-angle mismatch and failure to reproduce secondary/tertiary arms (Fasano et al., 2024). Sophisticated parametric models or machine-learned emulators from high-mass hydrodynamic grids are motivated for future analyses.

Observationally, the prevalence of rings/gaps as compared to spirals depends sensitively on disk viscosity, dust properties, and thermodynamics. In mm continuum, co-located gaps and rings can mask spirals even when present (Speedie et al., 2022). In scattered light, gap-edge illumination effects, i.e., radiative shadowing from non-axisymmetric gap rims, can drive prominent $m$7 arms that mimic classical Lindblad spirals unless distinguished by kinematics or polarization (Muley et al., 2024, Muley et al., 27 Nov 2025).

5. Applications to Specific Systems and Constraints on Planet Populations

Multi-wavelength and kinematic analyses of spiral arms now provide robust constraints on individual planets:

System	Arm Pattern Speed	Implied Planet Radius	Mass Estimate	Key Evidence
MWC 758	$m$8 deg/yr	$m$9 au	$$\phi_{m,n}(r) = -\mathrm{sgn}(r-r_p)\frac{\pi}{4m} + 2\pi\frac{n}{m} - \int_{r_m^{\pm}}^{r} \frac{\Omega(r')}{c_s(r')}\sqrt{[1-(r'/r_p)^{3/2}]^2 - 1/m^2}\, dr'$$0 $$\phi_{m,n}(r) = -\mathrm{sgn}(r-r_p)\frac{\pi}{4m} + 2\pi\frac{n}{m} - \int_{r_m^{\pm}}^{r} \frac{\Omega(r')}{c_s(r')}\sqrt{[1-(r'/r_p)^{3/2}]^2 - 1/m^2}\, dr'$$1	Rigid-body rotation, non-GI (Ren et al., 2020, Ren et al., 2018)
SAO 206462	$$\phi_{m,n}(r) = -\mathrm{sgn}(r-r_p)\frac{\pi}{4m} + 2\pi\frac{n}{m} - \int_{r_m^{\pm}}^{r} \frac{\Omega(r')}{c_s(r')}\sqrt{[1-(r'/r_p)^{3/2}]^2 - 1/m^2}\, dr'$$2 deg/yr	$$\phi_{m,n}(r) = -\mathrm{sgn}(r-r_p)\frac{\pi}{4m} + 2\pi\frac{n}{m} - \int_{r_m^{\pm}}^{r} \frac{\Omega(r')}{c_s(r')}\sqrt{[1-(r'/r_p)^{3/2}]^2 - 1/m^2}\, dr'$$3 au	$$\phi_{m,n}(r) = -\mathrm{sgn}(r-r_p)\frac{\pi}{4m} + 2\pi\frac{n}{m} - \int_{r_m^{\pm}}^{r} \frac{\Omega(r')}{c_s(r')}\sqrt{[1-(r'/r_p)^{3/2}]^2 - 1/m^2}\, dr'$$4–$$\phi_{m,n}(r) = -\mathrm{sgn}(r-r_p)\frac{\pi}{4m} + 2\pi\frac{n}{m} - \int_{r_m^{\pm}}^{r} \frac{\Omega(r')}{c_s(r')}\sqrt{[1-(r'/r_p)^{3/2}]^2 - 1/m^2}\, dr'$$5 $$\phi_{m,n}(r) = -\mathrm{sgn}(r-r_p)\frac{\pi}{4m} + 2\pi\frac{n}{m} - \int_{r_m^{\pm}}^{r} \frac{\Omega(r')}{c_s(r')}\sqrt{[1-(r'/r_p)^{3/2}]^2 - 1/m^2}\, dr'$$6	Spiral/gap/filament match (Xie et al., 2024)
HD 100546	Spiral pitch $$\phi_{m,n}(r) = -\mathrm{sgn}(r-r_p)\frac{\pi}{4m} + 2\pi\frac{n}{m} - \int_{r_m^{\pm}}^{r} \frac{\Omega(r')}{c_s(r')}\sqrt{[1-(r'/r_p)^{3/2}]^2 - 1/m^2}\, dr'$$7	$$\phi_{m,n}(r) = -\mathrm{sgn}(r-r_p)\frac{\pi}{4m} + 2\pi\frac{n}{m} - \int_{r_m^{\pm}}^{r} \frac{\Omega(r')}{c_s(r')}\sqrt{[1-(r'/r_p)^{3/2}]^2 - 1/m^2}\, dr'$$850 au, $$\phi_{m,n}(r) = -\mathrm{sgn}(r-r_p)\frac{\pi}{4m} + 2\pi\frac{n}{m} - \int_{r_m^{\pm}}^{r} \frac{\Omega(r')}{c_s(r')}\sqrt{[1-(r'/r_p)^{3/2}]^2 - 1/m^2}\, dr'$$9100 au	—	Three arms: kinematics, vertical flows, gas gaps (Wölfer et al., 15 Dec 2025)
SR 21	—	$c_s(r)$0 au	$c_s(r)$1 $c_s(r)$2	Scattered light spirals and mm ring match hydro models (Muro-Arena et al., 2020)

The absence of detectable spiral structure sets exclusion limits on planet masses, e.g., in Taurus disks, $c_s(r)$3 over 20–60 au can be ruled out if no spiral is seen (Stevenson et al., 2024).

6. Impact on Disk Evolution and Planet Formation

Spiral shocks not only sculpt gas and dust distributions but also promote dust fragmentation inside the arms. Numerical simulations with particle tracking find collision velocities in spirals exceeding fragmentation thresholds ($c_s(r)$4–$c_s(r)$5 m/s), implying rapid downsizing of solids and potential leakage of small grains through planetary gaps, challenging the idealized model of pebble isolation (Eriksson et al., 2024). This feedback between dynamically excited spirals and dust evolution influences gap clearing, ring morphology, planetesimal formation, and accretion efficiency.

Over time, spiral-induced gaps and rings develop alongside or outlive the arms; detection of spirals in continuum and gas is temporally limited compared to more persistent rings/gaps (Bae et al., 2017).

7. Distinguishing Competing Mechanisms and Future Directions

While planet-driven spirals (Lindblad wakes) are now robustly observed in several disks—anchored by dynamical constraints—alternate phenomena can also generate spiral structure:

• Gravitational Instabilities (GI): Produce m=2 spirals with nearly constant pitch, typically at higher disk masses, and show radially invariant temporal shearing (Bae et al., 2017, Ren et al., 2020).
• Gap-Edge Illumination Spirals: m=2 arms in the disk surface launched by time-varying irradiation due to gap-edge asymmetries or RWI, distinct by high surface polarization, vanishing in midplane tracers, and discrepant kinematic signatures (Muley et al., 2024, Muley et al., 27 Nov 2025).
• Buoyancy-Resonant Spirals: Vertical, tightly-wound arms traceable in CO line surface emission, distinguishable by their vertical, not planar, perturbed velocities (Bae et al., 2021).

Discriminating among these requires coordinated use of pattern speed measurements, multi-tracer spatial mapping, polarization, and velocity diagnostics. The combination of future high-sensitivity ALMA, JWST coronagraphy, and high-cadence imaging will enable increasingly precise measurements, refine planet occurrence statistics, and further elucidate the role of spiral arms in shaping planetary systems (Xie et al., 2024).