We normally observe classical behaviour due to the time dependent schrodinger equation in simple quantum systems when we introduce 'Gaussian wavepackets' which have bell shaped uncertainty in energy, momentum and position. A wavefunction which is often constructed by tapering the energy eigenstates. This can be seen in the case of a particle in a step function, harmonic oscillator, or under gravity (exhibits bouncing). This is especially true when $dL>>\hbar$ and the mean of quantum number $n$ is large.
Can this be modelled for a hydrogen atom? For example let us construct an atom with a small Rydberg constant, and have our initial state as $G(n,l,m)\psi_{nlm}(r)$ summed over all n,l,m where G is a tapered bell curve in all 3 variables with a large mean, i.e n = 50, l = 25 and m = 0. Then we plot this |psi|^2 as probability cloud and allow it to evolve over time. What does the shape of this probability cloud look like? Does it approximate classical behaviour such as orbital motion?
EDIT: My question asks in particular for the evolution of simple Gaussian coefficients in n,l,m. I'm not sure what the other question appears to model, or whether they have used the same approach.
Just to be clear I'm not asking for a necessarily Gaussian wavepacket. The distribution in n,l,m should be bell shaped. This is much simpler