Beyond solar metallicity: How enhanced solid content in disks re-shape low-mass planet torques

The migration of low-mass planets is tightly controlled by the torques exerted by both gas and solids in their natal disks. While canonical models assume a solar solid-to-gas mass ratio (epsilon=0.01) and neglect the back-reaction of solid component of the disk, recent work suggests that enhanced metallicity can radically alter these torques. We quantify how elevated metallicities (epsilon=0.03 and epsilon=0.1) modify the gas and solid torques, test widely used linear scaling prescriptions, and identify the regimes where solid back-reaction becomes decisive. We performed global, 2D hydrodynamic simulations that treat solid material as a pressureless fluid fully coupled to the gas through drag and include the reciprocal back-reaction force. The planet was maintained on a fixed circular orbit, thus we computed static torques. The Stokes number was varied from 0.01 to 10, three surface-density slopes (p=0.5, 1.0, and 1.5) and three accretion efficiencies (eta=0, 10, and 100%) were explored. Torques, obtained by rescaling canonical epsilon=0.01 results, were compared with direct simulations. Solid torques scale linearly with epsilon, but gas torques deviate by 50-100% and can even reverse sign for St<=1 in epsilon=0.1 disks. These are due to strong, feedback-driven, asymmetric gas perturbations in the co-orbital region, amplified by rapid planetary accretion. Solid back-reaction in high-metallicity environments can dominate the migration torque budget of low-mass planets. Simple metallicity rescalings are therefore unreliable for St<=2, implying that precise migration tracks - particularly in metal-rich disks – require simulations that fully couple solid and gas dynamics. These results highlight metallicity as a key parameter in shaping the early orbital architecture of planetary systems.

💡 Research Summary

This paper investigates how elevated solid-to-gas mass ratios (metallicities) in protoplanetary disks affect the torques exerted on low‑mass planets (≤10 M⊕). While traditional planet‑disk interaction models assume a solar metallicity (ε≈0.01) and neglect the back‑reaction of solids on the gas, recent observations suggest that many disks may be significantly richer in solids (ε=0.03–0.1). The authors perform a suite of global, two‑dimensional hydrodynamic simulations using the GPU‑accelerated GFARGO2 code. The solid component is treated as a pressureless fluid coupled to the gas through aerodynamic drag, and the reciprocal back‑reaction force is fully included. An Earth‑mass planet (1 M⊕) is held on a fixed circular orbit at R=1 for 200 orbital periods, allowing torques to reach a saturated state.

Key parameters explored are: metallicity ε (0.01, 0.03, 0.1), Stokes number St (0.01, 0.1, 1, 2, 3, 4, 5, 10), surface‑density slope p (0.5, 1.0, 1.5), and planetary accretion efficiency η (0 %, 10 %, 100 %). Torques are measured separately for the solid component (Γ_d) and the gas component (Γ_g), and then combined to obtain the total torque.

The main findings are:

1. Solid torque scaling – The torque contributed by solids scales almost linearly with metallicity across all Stokes numbers. This confirms the widely used assumption that solid torque ∝ ε when back‑reaction is weak.

2. Gas torque non‑linearity – Gas torques deviate strongly from a simple linear scaling when ε is increased. For ε=0.03 the gas torque magnitude can change by up to ~50 % relative to the canonical case, but its sign generally remains negative. At ε=0.1, however, the gas torque can be reduced or even reversed for St ≤ 1, with deviations up to 100 % in magnitude. The reversal is caused by strong, asymmetric gas density perturbations in the co‑orbital region, driven by the back‑reaction of a dense solid layer.

3. Role of back‑reaction – The back‑reaction of solids on the gas creates a pressure dip that is displaced asymmetrically ahead of or behind the planet. This asymmetry is amplified when the planet accretes solids efficiently (η=100 %). Rapid accretion empties the planet’s Hill sphere, steepening solid density gradients and thereby strengthening the feedback on the gas.

4. Stokes number dependence – For large particles (St ≥ 3) the coupling between solids and gas is weak, and both solid and gas torques follow the simple linear metallicity scaling. In contrast, for tightly coupled particles (St ≤ 2) the back‑reaction becomes dominant, leading to non‑linear torque behavior and, in some cases, sign changes.

5. Comparison with analytic prescriptions – The authors test a common analytic prescription that assumes solid torque ∝ ε and gas torque corrected linearly with ε. Direct simulation results show that this prescription overestimates total torques by up to a factor of two for 1 ≤ St ≤ 5 at ε=0.03, and can even predict the wrong sign for St ≤ 2 at ε=0.1. The discrepancy originates mainly from the gas torque, while solid torque predictions remain accurate.

6. Astrophysical implications – In metal‑rich disks, solid torques can become comparable to or exceed gas torques, potentially halting inward migration or driving outward migration, especially near regions of enhanced solid concentration such as the water‑ice line. This suggests that metallicity is a crucial parameter shaping the early orbital architecture of planetary systems.

7. Conclusions and future work – The study concludes that simple metallicity rescaling is unreliable for St ≤ 2, and that accurate migration tracks in high‑metallicity environments require fully coupled gas‑solid simulations. The authors recommend extending the work to three dimensions, incorporating realistic thermodynamics, and modeling dust growth and fragmentation to capture the full complexity of planet‑disk interactions in metal‑rich disks.