ABSTRACT
In slowly evolving spherical potentials, Φ(r, t), radial actions are typically assumed to remain constant. Here, we construct dynamical invariants that allow us to derive the evolution of radial actions in spherical central potentials with an arbitrary time dependence. We show that to linear order, radial actions oscillate around a constant value with an amplitude $\propto \dot{\Phi }/\Phi \, P(E,L)$. Using this result, we develop a diffusion theory that describes the evolution of the radial action distributions of ensembles of tracer particles orbiting in generic time-dependent spherical potentials. Tests against restricted N-body simulations in a varying Kepler potential indicate that our linear theory is accurate in regions of phase-space in which the diffusion coefficient $\tilde {D}(J_r) \lt 0.01\, J_r^2$. For illustration, we apply our theory to two astrophysical processes. We show that the median mass accretion rate of a Milky Way (MW) dark matter (DM) halo leads to slow global time-variation of the gravitational potential, in which the evolution of radial actions is linear (i.e. either adiabatic or diffusive) for ∼84 per cent of the DM halo at redshift z = 0. This fraction grows considerably with look-back time, suggesting that diffusion may be relevant to the modelling of several Gyr old tidal streams in action-angle space. As a second application, we show that dynamical tracers in a dwarf-size self-interacting DM halo (with $\sigma /m_\chi = 1\, {\rm cm^2g^{-1}}$) have invariant radial actions during the formation of a cored density profile.
Dynamical invariants and diffusion of merger substructures in time-dependent gravitational potentials
