Conformal Killing Tensors and covariant Hamiltonian Dynamics
  A covariant algorithm for deriving the conserved quantities for natural Hamiltonian systems is combined with the non-relativistic framework of Eisenhart, and of Duval, in which the classical trajectories arise as geodesics in a higher dimensional space-time, realized by Brinkmann manifolds. Conserved quantities which are polynomial in the momenta can be built using time-dependent conformal Killing tensors with flux. The latter are associated with terms proportional to the Hamiltonian in the lower dimensional theory and with spectrum generating algebras for higher dimensional quantities of order $1$ and $2$ in the momenta. Illustrations of the general theory include the Runge-Lenz vector for planetary motion with a time-dependent gravitational constant $G(t)$, motion in a time-dependent electromagnetic field of a certain form, quantum dots, the H\'enon-Heiles and Holt systems, respectively, providing us with Killing tensors of rank that ranges from one to six. 
M. Cariglia, G.W. Gibbons, J.-W. van Holten, P.A. Horvathy, P.-M.   Zhang
math-ph hep-th math.MP