If you set a displaced Gaussian initial state
##\Psi (x,t_0 ) = Ae^{-b(x+\Delta x)^2}##
in motion in a harmonic oscillator potential ##V(x) = \frac{1}{2}kx^2##, the expectation values ##\left<x\right>## and ##\left<p\right>## will evolve with cosine and sine time dependence just like the classical harmonic oscillator set in motion from rest at a point with distance ##\Delta x## from equilibrium. However, if the function ##V(x)## has any powers of ##x## higher than ##x^2##, this will not work exactly like that.
Some ions can have so called "Kepler wavepackets" orbiting around them, if I can remember correctly. This means that there's a relatively localized electron wavepacket at a quite large distance from the nucleus, and it moves somewhat similarly to a planet orbiting a star, except that it probably comes down gradually because of photon emission. But I don't think that the ##m\frac{d\left<x\right>}{dt} = \left<p\right>## correspondence is true in that case.