2025MNRAS.543..565F

Modelling the gravitational interaction between an eccentric perturber and a differentially shearing gas disc is a longstanding problem with various astrophysical applications, ranging from the evolution of planetary systems to the migration of black holes embedded in AGN discs. Recent work has advanced a global, linear, modal approach for calculating the excited wake and the resulting feedback on the perturber's orbital evolution. In this work, we perform a complementary suite of targeted hydrodynamic simulations to test this linear framework across a range of disc temperature and density profiles. In particular, we push from circular orbits to highly eccentric trajectories for which the perturber moves supersonically with respect to the background gas. We find remarkable agreement between our simulations and the linear methodology across a range of diagnostics - lending support to the predicted wake morphologies, complex radial torque density profiles, and torque reversal phenomena, which occur when the eccentricity exceeds the local aspect ratio. In contrast, comparison with previous fitting functions exposes noticeable discrepancies, cautioning against their indiscriminate use in studies which explore a wide range of perturber eccentricities, in varied disc structures. Our simulations also probe the fundamentally nonlinear effects of shock induced angular momentum deposition and coorbital horseshoe drags, which exhibit clear variations with eccentricity. Finally, this careful comparison between linear theory and numerics provides a detailed benchmark for planet-disc interaction problems and therefore we have provided a repository of our linear calculations for use as a rigorous test of future numerical investigations.