Dynamics of tidal spiral arms: Machine learning-assisted identification of equations and application to the Milky Way

TL;DR

This paper uses machine learning to identify equations governing tidally induced spiral arms in the Milky Way, combining simulations and data to develop models that match observations and can be extended to complex galactic phenomena.

Contribution

It introduces a novel machine learning approach to derive equations for tidal spiral arms, including non-linear dynamics, validated with simulations and applied to Gaia data.

Findings

Linear wave patterns match small perturbation simulations

Non-linear equations describe large impact spiral dynamics

Fitted models to Gaia data suggest complex origins for observed features

Abstract

Understanding the spiral arms of the Milky Way (MW) remains a key open question in galactic dynamics. Tidal perturbations, such as the recent passage of the Sagittarius dwarf galaxy (Sgr), could play a significant role in exciting them. We aim to analytically characterize the dynamics of tidally induced spiral arms, including their phase-space signatures. We ran idealized test-particle simulations resembling impulsive satellite impacts, and used the Sparse Identification of Non-linear Dynamics (SINDy) method to infer their governing Partial Differential Equations (PDEs). We validated the method with analytical derivations and a realistic $N$-body simulation of a MW-Sgr encounter analogue. For small perturbations, a linear system of equations was recovered with SINDy, consistent with predictions from linearised collisionless dynamics. In this case, two distinct waves wrapping at pattern speeds $\Omega \pm \kappa/m$ emerge. For large impacts, we empirically discovered a non-linear system of equations, representing a novel formulation for the dynamics of tidally induced spiral arms. For both cases, these equations describe wave properties like amplitude and pattern speed, and their shape and temporal evolution in different phase-space projections. We fit the Gaia $L_Z-V_R$ waves with the linear model, providing a reasonable fit and plausible parameters for the Sgr passage. However, the predicted amplitude ratio of the two waves is inconsistent with observations, supporting a more complex origin for this feature (e.g. multiple passages, bar, spiral arms). We merge data-driven discovery with theory to create simple, accurate models of tidal spiral arms that match simulations and provide a simple tool to fit Gaia and external galaxy data. This methodology could be extended to model complex phenomena like self-gravity and dynamical friction. (ABR)