Not Just Gas: How Solid-Driven Torques Shaped the Migration of the Galilean Moons

Surviving rapid inward orbital migration is a crucial aspect of formation models for the Jupiter’s Galilean moons. The primary aim of this study is to investigate the orbital migration of the Galilean moons by incorporating self-consistent solid dynamics in circumjovian disk models. We perform two-fluid simulations using the FARGO3D code on a 2D polar grid. The simulations model a satellite with the mass of a proto-moon, Europa, or Ganymede interacting with a circumjovian disk. The dust component, coupled to the gas via a drag force, is characterized by the dust-to-gas mass ratio ($ε$) and the Stokes number ($T_s$). The effect of solids fundamentally alter the satellites’ evolution. We identify a vast parameter space where migration is slowed, halted, robustly reversed -leading to outward migration-, or significantly accelerated inward. The migration rate is dependent on satellite mass, providing a natural source of differential migration. Solid dynamics provides a robust and self-consistent mechanism that fundamentally alters the migration of the Galilean moons, potentially addressing the long-standing migration catastrophe. This mechanism critically affects the survival of satellites and could offer a viable physical process to explain the establishment of resonances through differential migration. These findings establish that solid torques are a critical, non-negligible factor in shaping the final architecture of satellite systems.

💡 Research Summary

This paper presents the first systematic investigation of how solid particles embedded in a circumjovian disk modify the orbital migration of the Galilean satellites. Using the two‑fluid version of the FARGO3D code, the authors simulate a 2‑D polar grid (768 × 4096 cells) that resolves both the gas dynamics and the dust component coupled through Epstein/Stokes drag. The gas is modeled as locally isothermal (aspect ratio h = 0.1) with an α‑viscosity of 10⁻³ and a surface density profile Σ ∝ r⁻¹/₂. The dust is characterized by a constant dust‑to‑gas mass ratio ε (ranging from 0 to 0.5) and a dimensionless stopping time (Stokes number) Tₛ (0.01–5). Three satellite masses are considered: a proto‑moon (Mₛ≈1.25×10⁻⁵ M_J), Europa (Mₛ≈2.53×10⁻⁵ M_J), and Ganymede (Mₛ≈7.82×10⁻⁵ M_J). The satellite’s gravitational potential is softened over its Hill radius, and back‑reaction of dust on gas is included.

The simulations are run until the net torque (gas + dust) on the fixed‑orbit satellite reaches a quasi‑steady value (typically 28–56 orbital periods). The results are summarized in a torque map (Figure 1) that shows the normalized net torque Γ/Γ₀,gas as a function of (ε, Tₛ) for each satellite mass. The key findings are:

1. Positive (outward) torques are common when solids are abundant or large. For the proto‑moon and Europa, ε < 0.01 requires Tₛ ≳ 2 to produce outward migration, but increasing ε relaxes the Tₛ threshold. Ganymede, being more massive, needs ε ≳ 0.025 and Tₛ ≳ 1. In these regimes the solid torque dominates the gas torque, reversing the direction of migration.

2. Intermediate Stokes numbers (0.1 ≲ Tₛ ≲ 1) can strongly amplify inward migration. In this regime the net torque magnitude can reach 3–7 times the gas‑only value, especially for the low‑mass proto‑moon. This “danger zone” could cause rapid loss of a nascent satellite unless it grows quickly enough to exit the regime.

3. Very high dust‑to‑gas ratios (ε ≈ 0.5) trigger dramatic torque enhancements. Dust feedback generates vortices in the horseshoe region, leading to stochastic torque fluctuations and an average torque >100 × the gas‑only case. This suggests that in gas‑starved disk models, where ε may approach unity, outward migration of the regular satellites is almost inevitable once they exceed Mₛ ≈ 10⁻⁵ M_J.

4. The torque response is largely independent of the gas surface density Σ. Because the torque is measured on a fixed orbit, the sign and magnitude depend primarily on ε, Tₛ, and satellite mass, making the results applicable across a wide range of circumjovian disk models (gas‑starved, minimum‑mass subnebula, etc.).

5. Differential migration naturally emerges. The mass‑dependent torque curves mean that, for a given (ε, Tₛ), a more massive satellite (e.g., Ganymede) may experience reduced inward drift or even outward drift, while a lighter one (e.g., Europa) may stall or migrate outward at different thresholds. This provides a plausible mechanism for establishing the observed 1:2:4 mean‑motion resonance among Io, Europa, and Ganymede.

6. Limitations and future directions. The study adopts a locally isothermal equation of state, neglecting non‑isothermal effects such as cold thermal torques and heating torques, which can modify gas torques by orders of magnitude. Moreover, the satellites are held on fixed circular orbits and do not accrete solids during the runs. The authors discuss that solid accretion would reshape the dust distribution (enhancing asymmetries) and introduce an “accretion torque” that can further boost outward migration, especially for low‑mass bodies. Incorporating realistic thermodynamics, particle growth, and self‑consistent accretion in future 3‑D simulations will be essential to confirm and extend these findings.

In summary, the paper demonstrates that solid‑driven torques are not a secondary correction but a dominant factor that can halt, reverse, or accelerate the migration of Galilean moons. By mapping the parameter space of dust‑to‑gas ratio and Stokes number, the authors provide a robust physical solution to the long‑standing “migration catastrophe” problem and offer a natural pathway for the formation of the Galilean satellite system’s resonant architecture.