Global change in action due to trapping: How to derive it whatever the rate of variation of the dynamics

Phys Rev E Stat Nonlin Soft Matter Phys. 2015 Apr;91(4):042915. doi: 10.1103/PhysRevE.91.042915. Epub 2015 Apr 29.

Abstract

In this paper, we investigate the motion of a set of charged particles acted upon by a growing electrostatic wave in the limit when the initial wave amplitude is vanishingly small and when all the particles have the same initial action, I(0). We show, both theoretically and numerically, that when all the particles have been trapped in the wave potential, the distribution in action exhibits a very sharp peak about the smallest action. Moreover, as the wave keeps growing, the most probable action tends toward a constant, I(f), which we estimate theoretically. In particular, we show that I(f) may be calculated very accurately when the particles' motion before trapping is far from adiabatic by making use of a perturbation analysis in the wave amplitude. This fills a gap regarding the computation of the action change, which, in the past, has only been addressed for slowly varying dynamics. Moreover, when the variations of the dynamics are fast enough, we show that the Fourier components of the particles' distribution function can be calculated by connecting estimates from our perturbation analysis with those obtained by assuming that all the particles have the same constant action, I=I(f). This result is used to compute theoretically the imaginary part of the electron susceptibility of an electrostatic wave in a plasma. Moreover, using our formula for the electron susceptibility, we can extend the range in ε(a) (the parameter that quantifies the slowness of the dynamics) for our perturbative estimate of I(f)-I(0). This range can actually be pushed down to values of ε(a) allowing the use of neoadiabatic techniques to compute the jump in action. Hence, this paper shows that the action change due to trapping can be calculated theoretically, regardless of the rate of variation of the dynamics, by connecting perturbative results with neoadiabatic ones.