We prove an abstract result giving a t ǫ upper bound on the growth of the Sobolev norms of a time dependent Schrödinger equation of the form i ψ = H 0 ψ + V (t)ψ. H 0 is assumed to be the Hamiltonian of a steep quantum integrable system and to be a pseudodifferential operator of order d > 1; V (t) is a time dependent family of pseudodifferential operators, unbounded, but of order b < d. The abstract theorem is then applied to perturbations of the quantum anharmonic oscillators in dimension 2 and to perturbations of the Laplacian on a manifold with integrable geodesic flow, and in particular Zoll manifolds, rotation invariant surfaces and Lie groups. We also cover the case of several particles on a Zoll manifold or on a Lie group, possibly obeying some restrictions due to the Fermionic or Bosonic nature of the particles. The proof is based a on quantum version of the proof of the classical Nekhoroshev theorem.