A numerical method for solution of the generalized Liouville equation (Q2564467)

A numerical method for the time evolution of systems described by Liouville-type equations is derived. The algorithm uses a lattice of numerical markers, which follow exactly Hamiltonian trajectories, to represent the operator \(d/dt\) in moving (i.e. Lagrangian) coordinates. However, nonconservative effects such as particle drag, creation, and annihilation are allowed in the evolution of the physical distribution function. The method is suited to the study of a general class of systems involving the resonant interaction of energetic particles with plasma waves.