Planet Migration in Protoplanetary Disks with Rims

Complex structures, including sharp edges, rings and gaps, have been commonly observed in protoplanetary disks with or without planetary candidates. Here we consider the possibility that they are the intrinsic consequences of angular momentum transfer mechanisms, and investigate how they may influence the dynamical evolution of embedded planets. With the aid of numerical hydrodynamic simulations, we show that gas giants have a tendency to migrate away from sharp edges, whereas super-Earths embedded in the annuli tend to be retained. This implies that, observationally, Jupiters are preferentially detected in dark rings (gaps), whereas super-Earths tend to be found in bright rings (density bumps). Moreover, planets’ tidal torque provide, not necessarily predominant, feedback on the surface density profile. This tendency implies that Jupiter’s gap-opening process deepens and widens the density gap associated with the dark ring, while super-Earths can be halted by steep surface density gradient near the disk or ring boundaries. 13Hence, we expect there would be a desert for super-Earths in the surface density gap.

💡 Research Summary

This paper investigates how sharp radial variations in viscosity—commonly referred to as “disk rims” or dead‑zone boundaries—affect the migration of embedded planets in protoplanetary disks. Using two‑dimensional hydrodynamic simulations performed with the Athena++ code, the authors construct a disk model in which the α‑viscosity parameter transitions from a low value (α_dead) inside a dead zone to a higher value (α_active) in the surrounding active zone. The transition is described by a hyperbolic‑tangent function that defines inner and outer rim locations (r_dz_in and r_dz_out) and a steepness parameter δ. For the Jupiter‑mass case α_dead=10⁻⁴ (α_active=10⁻²); for the super‑Earth case α_dead=10⁻³ (α_active=10⁻²). This viscosity profile naturally produces a surface‑density bump (bright ring) inside the dead zone and a corresponding gap (dark ring) where viscosity is higher.

Two planetary masses are examined: a Jupiter‑mass planet (M_p=10⁻³ M_⊙) and a super‑Earth (M_p=3×10⁻⁵ M_⊙). Both are placed on a fixed circular orbit at r₀, with softening lengths ε=0.8 H (Jupiter) and 0.6 H (super‑Earth) to mimic three‑dimensional effects. In some runs an accretion sink (radius 0.1 r_H, removal rate η=10 Ω₀) is added to test the influence of planetary growth. The computational domain spans 0.3 r₀–5 r₀ with 240 radial and 512 azimuthal cells, refined to Δr≈0.001 r₀ near the planet. Each simulation runs for 1,000–2,000 planetary orbits until the torque on the planet reaches a quasi‑steady state.

The total torque is obtained by direct summation of the softened gravitational force from each grid cell, averaged over the final 160 orbits. The authors map torque values onto a plane defined by Δ_dz_in (distance from the planet to the inner rim) and Δ_dz_out (distance to the outer rim). Blue points denote negative torque (inward migration), red points positive torque (outward migration).

Key results:

1. Jupiter‑mass planets experience a strong positive corotation torque when located just inside the density bump (∂Σ/∂r>0). This torque overwhelms the negative Lindblad torque, driving outward migration. When the planet sits outside the bump, the Lindblad torque dominates and the planet migrates inward. A narrow region near Δ_dz_in≈0.8 r₀ and Δ_dz_out≈1.3 r₀ yields near‑zero net torque, acting as a potential trapping site.

2. Super‑Earths have much weaker corotation torques; their migration is governed primarily by Lindblad torques. Inside the bump they tend to migrate inward, while outside the bump the net torque is close to zero, allowing the planet to stall. The torque reversal is much less pronounced than for the Jupiter case, leading to a “parking” of low‑mass planets at the rim edges.

3. Feedback on the disk: The massive planet’s tidal torques deepen the pre‑existing gap, widening the dark ring and slightly shifting the rim location. The super‑Earth’s influence is modest but can locally steepen the surface‑density gradient, reinforcing the trapping mechanism.

4. Dependence on rim steepness: Smaller δ (sharper viscosity transition) produces a more abrupt density jump, which amplifies both corotation and Lindblad torques, enlarging the region where migration can be halted or reversed.

Observational implications follow directly. High‑resolution ALMA images often reveal concentric bright rings and dark gaps. According to the simulations, giant planets should preferentially be found within the dark gaps, whereas super‑Earths should accumulate in the bright rings. This offers a physical explanation for the apparent “Jupiter desert” inside bright rings and the “super‑Earth pile‑up” in high‑density regions.

In summary, the study adds a crucial “rim effect” to the classic Type I/II migration framework. Sharp viscosity transitions, by reshaping the surface‑density profile, can halt, reverse, or accelerate planetary migration depending on planet mass. This mechanism helps reconcile theoretical migration rates with the observed distribution of exoplanets and disk substructures, highlighting the importance of detailed disk viscosity structure in shaping planetary system architectures.