On the horseshoe drag of a low-mass planet. I - Migration in isothermal disks

We investigate the unsaturated horseshoe drag exerted on a low-mass planet by an isothermal gaseous disk. In the globally isothermal case, we use a formal- ism, based on the use of a Bernoulli invariant, that takes into account pressure effects, and that extends the torque estimate to a region wider than the horse- shoe region. We find a result that is strictly identical to the standard horseshoe drag. This shows that the horseshoe drag accounts for the torque of the whole corotation region, and not only of the horseshoe region, thereby deserving to be called corotation torque. We find that evanescent waves launched downstream of the horseshoe U-turns by the perturbations of vortensity exert a feed-back on the upstream region, that render the horseshoe region asymmetric. This asymmetry scales with the vortensity gradient and with the disk’s aspect ratio. It does not depend on the planetary mass, and it does not have any impact on the horseshoe drag. Since the horseshoe drag has a steep dependence on the width of the horseshoe region, we provide an adequate definition of the width that needs to be used in horseshoe drag estimates. We then consider the case of locally isothermal disks, in which the tempera- ture is constant in time but depends on the distance to the star. The horseshoe drag appears to be different from the case of a globally isothermal disk. The difference, which is due to the driving of vortensity in the vicinity of the planet, is intimately linked to the topology of the flow. We provide a descriptive inter- pretation of these effects, as well as a crude estimate of the dependency of the excess on the temperature gradient.

💡 Research Summary

The paper presents a thorough analytical study of the corotation (horseshoe) torque exerted on a low‑mass planet embedded in an isothermal protoplanetary disk. The authors first treat the case of a globally isothermal disk, where the temperature is uniform throughout the disk. By introducing a Bernoulli invariant that simultaneously conserves the total specific energy and the specific angular momentum along streamlines, they are able to incorporate pressure gradients into the torque calculation. This formalism extends the torque estimate beyond the traditional horseshoe region, yet the resulting expression for the torque is exactly identical to the classic horseshoe‑drag formula. Consequently, the horseshoe drag is shown to represent the torque of the entire corotation region, justifying the use of the term “corotation torque.”

A key finding is that evanescent (damped) waves generated downstream of the horseshoe U‑turns feed back on the upstream flow, producing a modest asymmetry of the horseshoe region. The magnitude of this asymmetry scales linearly with the vortensity gradient (d ln ζ/d ln r) and with the disk aspect ratio h = H/r, but it is independent of the planet mass. Because the horseshoe drag depends on the fourth power of the horseshoe half‑width, the authors argue that a precise definition of this width is essential. They propose a definition based on the Bernoulli invariant and the actual stream‑tube boundaries, which remains valid when pressure effects deform the flow.

The second part of the paper addresses locally isothermal disks, where the temperature varies radially as T(r) ∝ r^β while remaining constant in time. In this situation the radial temperature gradient creates a source term for vortensity in the vicinity of the planet. This “vortensity generation” modifies the topology of the horseshoe flow, leading to a torque that differs from the globally isothermal case. The authors provide a qualitative description of the altered streamlines and derive a crude scaling for the excess torque: ΔΓ ≈ C β h² Γ₀, where Γ₀ is the standard horseshoe torque for β = 0 and C is an order‑unity constant.

The analytical results are validated with two‑dimensional hydrodynamic simulations. In the globally isothermal runs the measured torque matches the analytic horseshoe drag to within a few percent, and the predicted asymmetry of the horseshoe region is observed in the density and velocity fields. In locally isothermal runs, the torque grows with the magnitude of the temperature gradient, and the flow exhibits the predicted distortion of the horseshoe U‑turns.

Overall, the paper makes three major contributions: (1) it rigorously proves that the horseshoe drag accounts for the full corotation torque in an isothermal disk; (2) it clarifies how pressure, vortensity gradients, and temperature gradients each affect the torque, including the subtle role of evanescent wave feedback; and (3) it supplies a robust definition of the horseshoe half‑width and a simple scaling law for the additional torque arising in locally isothermal disks. These results provide a solid theoretical foundation for future migration studies, especially those that aim to incorporate realistic temperature structures and vortensity evolution in planet‑disk interaction models.