2023MNRAS.521.5917F

The distribution of stars in the Milky Way disc shows a spiral structure - the snail - in the space of velocity and position normal to the Galactic mid-plane. The snail appears as straight lines in the vertical frequency-vertical phase plane when effects from sample selection are removed. Their slope has the dimension of inverse time, with the simplest interpretation being the inverse age of the snail. Here, we devise and fit a simple model in which the spiral starts as a lopsided perturbation from steady state, that winds up into the present-day morphology. The winding occurs because the vertical frequency decreases with vertical action. We use data from stars in Gaia EDR3 that have measured radial velocities, pruned by simple distance and photometric selection functions. We divide the data into boxels of dynamical invariants (radial action, angular momentum); our model fits the data well in many of the boxels. The model parameters have physical interpretations: one, A, is a perturbation amplitude, and one, t, is interpretable in the simplest models as the time since the event that caused the snail. We find trends relating the strength and age to angular momentum: (i) the amplitude A is small at low angular momentum (${\lt 1\, 600{\,}kpc\ km\ s}^{-1}$ or guiding-centre radius $\lt 7.3\,$kpc) and over a factor of three larger, with strong variations, in the outer disc; (ii) there is no single well-defined perturbation time, with t varying between 0.2 and 0.6 Gyr. Residuals between the data and the model display systematic trends, implying that the data call for more complex models.