Fundamental limits to orbit reconstruction due to non-conservation of stellar actions in a Milky Way-like simulation

ABSTRACT The conservation of stellar actions is a fundamental assumption in orbit reconstruction studies in the Milky Way. However, the disc is highly dynamic, with time-dependent, non-axisymmetric features like transient spiral arms and giant molecular clouds (GMCs) driving local fluctuations in the gravitational potential on top of the near-axisymmetric background. Using high-resolution magnetohydrodynamic simulations that incorporate gas dynamics and star formation, we quantify the rate at which these effects drive non-conservation of the actions of young stars from Myr to Gyr time-scales. We find that action evolution is well described as a logarithmic random walk, with vertical action evolving more rapidly than radial action; the diffusion rate associated with this random walk is weakly dependent on the stellar birth environment and scales approximately linearly with the galactic orbital frequency at a star’s position. The diffusion rates we measure imply a fundamental limit of ∼100 Myr as the time-scale over which stellar orbits can be reliably reconstructed using methods that assume action conservation. By comparing diffusion rates for younger stars to those measured for an older and more vertically extended control population, we conclude that radial action evolution is driven primarily by transient spiral arms, while vertical action evolution is driven by gravitational scattering off gaseous structures. Our results have significant implications for galactic archaeology and disc dynamics studies, necessitating a closer look at the time-scales over which actions are assumed to be conserved in the disc.