PiPER Arms in Protoplanetary Disks

• PiPER Arms are spiral density wave structures in protoplanetary disks generated by constructive interference of multiple azimuthal modes from a single embedded planet.
• Their number, profile, and pitch angles depend sensitively on planet mass and disk thermodynamics, offering a practical tool to infer disk conditions.
• Observational fits to theoretical PiPER models have accurately inferred planet properties in systems like MWC 758 and TW Hya, resolving previous misconceptions about spiral arm origin.

Planet-driven spiral arms (abbreviated "PiPER Arms") in protoplanetary disks refer to the spiral density wave structures excited by a single embedded planet. Contrary to the earlier understanding that a planet produces only a single primary spiral wake, recent simulations and analysis demonstrate that a single planet can excite multiple spiral arms due to constructive interference of density waves launched at different azimuthal modes. The number, profile, and location of these arms depend sensitively on planet mass and the physical state of the disk, notably the temperature and aspect ratio. PiPER Arms provide a robust physical mechanism explaining high-resolution observations of multiple spiral arms in transition disks and form a diagnostic tool for inferring planet properties from disk morphologies (Bae et al., 2017).

1. Governing Equations and Linear Excitation

The dynamics of PiPER Arms are governed by the equations of mass and momentum conservation in a two-dimensional, isothermal disk:

• Mass conservation: $$\partial \Sigma/\partial t + \nabla\cdot (\Sigma \mathbf{v}) = 0$$
• Momentum conservation: $$\Sigma[\partial\mathbf{v}/\partial t + (\mathbf{v}\cdot\nabla)\mathbf{v}] = -\nabla P - \Sigma\nabla(\Phi_* + \Phi_p)$$

Where $P = \Sigma c_s^2$, $\Phi_* = -GM_*/r$, and $\Phi_p$ is the planetary potential. The epicyclic frequency $$\kappa^2(r) \equiv \frac{1}{r^3}\frac{d}{dr}[(r^2\Omega)^2]$$ and the sound speed $c_s(r) \sim r^{-q/2}$ for temperature $T\propto r^{-q}$ set the disk response.

The planetary potential is decomposed as $$\Phi_p(r,\phi,t) = \sum_m \Phi_m(r)\exp[im(\phi - \Omega_p t)]$$. The resulting density perturbations are described by an $m$-armed spatial mode, obeying the WKBJ dispersion relation:

$$\Sigma[\partial\mathbf{v}/\partial t + (\mathbf{v}\cdot\nabla)\mathbf{v}] = -\nabla P - \Sigma\nabla(\Phi_* + \Phi_p)$$0

Spiral waves are launched at Lindblad resonances, with resonance radii $$\Sigma[\partial\mathbf{v}/\partial t + (\mathbf{v}\cdot\nabla)\mathbf{v}] = -\nabla P - \Sigma\nabla(\Phi_* + \Phi_p)$$1 for Keplerian disks, where $$\Sigma[\partial\mathbf{v}/\partial t + (\mathbf{v}\cdot\nabla)\mathbf{v}] = -\nabla P - \Sigma\nabla(\Phi_* + \Phi_p)$$2 is the planetary orbital radius.

2. Constructive Interference and Emergence of Multiple Arms

Each azimuthal mode $$\Sigma[\partial\mathbf{v}/\partial t + (\mathbf{v}\cdot\nabla)\mathbf{v}] = -\nabla P - \Sigma\nabla(\Phi_* + \Phi_p)$$3 excites $$\Sigma[\partial\mathbf{v}/\partial t + (\mathbf{v}\cdot\nabla)\mathbf{v}] = -\nabla P - \Sigma\nabla(\Phi_* + \Phi_p)$$4 evenly spaced wave crests at $$\Sigma[\partial\mathbf{v}/\partial t + (\mathbf{v}\cdot\nabla)\mathbf{v}] = -\nabla P - \Sigma\nabla(\Phi_* + \Phi_p)$$5. The phase of the $$\Sigma[\partial\mathbf{v}/\partial t + (\mathbf{v}\cdot\nabla)\mathbf{v}] = -\nabla P - \Sigma\nabla(\Phi_* + \Phi_p)$$6 crest for azimuthal mode $$\Sigma[\partial\mathbf{v}/\partial t + (\mathbf{v}\cdot\nabla)\mathbf{v}] = -\nabla P - \Sigma\nabla(\Phi_* + \Phi_p)$$7 at radius $$\Sigma[\partial\mathbf{v}/\partial t + (\mathbf{v}\cdot\nabla)\mathbf{v}] = -\nabla P - \Sigma\nabla(\Phi_* + \Phi_p)$$8 is given by:

$$\Sigma[\partial\mathbf{v}/\partial t + (\mathbf{v}\cdot\nabla)\mathbf{v}] = -\nabla P - \Sigma\nabla(\Phi_* + \Phi_p)$$9

Constructive interference occurs when crests from a range of $P = \Sigma c_s^2$0 align in azimuth; the degree of phasing is quantified by the coherence criterion:

$P = \Sigma c_s^2$1

• Primary Arm: $P = \Sigma c_s^2$2 crests from several $P = \Sigma c_s^2$3 pile up coherently near the resonance, forming the primary arm.
• Secondary and Tertiary Arms: For $P = \Sigma c_s^2$4 (secondary) and $P = \Sigma c_s^2$5 (tertiary), crests phase up further inward, at smaller radii. The secondary arm in a typical disk with $P = \Sigma c_s^2$6 forms near $P = \Sigma c_s^2$7, the tertiary near $P = \Sigma c_s^2$8.

The emergence of multiple arms is thus a linear phenomenon for sufficiently low planet mass, governed by the disk's ability to maintain phase coherence among the excited wave modes (Bae et al., 2017).

3. Dependence on Planet Mass and Disk Thermodynamics

The critical parameter delineating the linear regime is the planet's thermal mass:

$P = \Sigma c_s^2$9

• For $\Phi_* = -GM_*/r$0, arm formation is linear and multiple arms can form interior to the planet's orbit. The number of interior arms $\Phi_* = -GM_*/r$1 is dictated by how many distinct $\Phi_* = -GM_*/r$2 crests align before reaching the disk's inner edge: up to four for $\Phi_* = -GM_*/r$3, and up to three for $\Phi_* = -GM_*/r$4 and $\Phi_* = -GM_*/r$5.
• Colder disks (smaller $\Phi_* = -GM_*/r$6) allow for more arms because the coherence criterion is easier to satisfy.

The pitch angle along an arm is:

$\Phi_* = -GM_*/r$7

and, to first order in the linear regime,

$\Phi_* = -GM_*/r$8

The arm-to-arm separation and pitch angle deviations scale with planet mass:

• Separation $\Phi_* = -GM_*/r$9
• Pitch angle deviations from linear theory grow as $\Phi_p$0 increases.

4. Nonlinear Evolution and Arm Mergers for Massive Planets

As $\Phi_p$1 increases $\Phi_p$2 for $\Phi_p$3, non-linear steepening of density waves becomes significant. Waves transition into shocks beyond the steepening distance

$\Phi_p$4

where $\Phi_p$5 a few. Shocked arms propagate faster than predicted linearly, resulting in pitch angles $\Phi_p$6; the excess scales as $\Phi_p$7. In simulations with $\Phi_p$8 at $\Phi_p$9, the pitch angle at $$\kappa^2(r) \equiv \frac{1}{r^3}\frac{d}{dr}[(r^2\Omega)^2]$$0 exceeds linear predictions by $$\kappa^2(r) \equiv \frac{1}{r^3}\frac{d}{dr}[(r^2\Omega)^2]$$1 respectively.

For $$\kappa^2(r) \equiv \frac{1}{r^3}\frac{d}{dr}[(r^2\Omega)^2]$$2, the tertiary crest merges with the broadened primary shock by $$\kappa^2(r) \equiv \frac{1}{r^3}\frac{d}{dr}[(r^2\Omega)^2]$$3, so only two distinct interior arms remain. This sets an upper limit on the number of observable arms for massive planets (Bae et al., 2017).

5. Inference of Planet Properties from Observational Data

PiPER Arms provide powerful diagnostics for planet detection and characterization in protoplanetary disks. Observable quantities include:

• Number of interior arms $$\kappa^2(r) \equiv \frac{1}{r^3}\frac{d}{dr}[(r^2\Omega)^2]$$4 (typically 1–3 for single planets),
• Radial launch points $$\kappa^2(r) \equiv \frac{1}{r^3}\frac{d}{dr}[(r^2\Omega)^2]$$5,
• Pitch angle $$\kappa^2(r) \equiv \frac{1}{r^3}\frac{d}{dr}[(r^2\Omega)^2]$$6,
• Azimuthal arm separations $$\kappa^2(r) \equiv \frac{1}{r^3}\frac{d}{dr}[(r^2\Omega)^2]$$7.

The observed arm morphology is fit to the theoretical phase model:

$$\kappa^2(r) \equiv \frac{1}{r^3}\frac{d}{dr}[(r^2\Omega)^2]$$8

Fitting procedures follow a data-driven recipe:

1. Measure $$\kappa^2(r) \equiv \frac{1}{r^3}\frac{d}{dr}[(r^2\Omega)^2]$$9 of the primary arm; use $c_s(r) \sim r^{-q/2}$0 to infer $c_s(r) \sim r^{-q/2}$1 (temperature).
2. Measure secondary separation $c_s(r) \sim r^{-q/2}$2 at $c_s(r) \sim r^{-q/2}$3; invert the scaling law to estimate $c_s(r) \sim r^{-q/2}$4.
3. Confirm $c_s(r) \sim r^{-q/2}$5 is consistent with the inferred parameters.

Empirical scalings aid in quantitative inference:

$c_s(r) \sim r^{-q/2}$6

with $c_s(r) \sim r^{-q/2}$7, $c_s(r) \sim r^{-q/2}$8 typical. The PiPER Arms framework has successfully yielded planet masses and locations in multi-arm systems such as MWC 758, HD 100453, and TW Hya, in agreement with other planet-disk diagnostics (Bae et al., 2017).

6. Observational Significance and Applicability

PiPER Arms provide a physical mechanism explaining the multiplicity of spiral arms in protoplanetary disks with a single embedded planet. The analytical and numerical results indicate all observed arms interior to the planet, as well as the sole outer arm, result from constructive interference of planet-driven density waves. Inference of disk temperature structures and embedded planet properties is made possible by fitting observed spiral morphologies to the theoretical PiPER model.

A plausible implication is that variations in $c_s(r) \sim r^{-q/2}$9, pitch angle, or arm separation can be traced directly to changing disk thermodynamics or to the planet's mass entering the non-linear regime. This framework resolves prior misconceptions that additional arms require multiple planets or fundamentally nonlinear mechanisms, establishing linear wave interference as the primary generative process in the low-mass regime.

The linear theory’s predictive accuracy for low $T\propto r^{-q}$0 and cold disks supports its application to current high-resolution imaging datasets, facilitating planet detection and characterization where other techniques are infeasible (Bae et al., 2017).