Why are the energy levels of a simple harmonic oscillator equally spaced? The energy level of a simple harmonic oscillator is $E_n=(n+\frac{1}{2})\hbar\omega$. Is there any physical explanation why these levels are equally spaced ($= \hbar\omega$)?
Maybe this link can be helpful.
 A: I might as well condense the neat phase-space argument of Royer's 1996 pedagogical paper, since it is so short and standard for the phase-space behavior of the QHO, and the other answers did not utilize it.  
The essential point is that, both classically and quantum mechanically, the time evolution of the oscillator is rigid rotation in phase space, as discovered by Groenewold in 1946.
That is to say that, just as in classical evolution, any quantum configuration is uniformly rotated in phase space by the hamiltonian. Absorbing all constants into the variables to rescale elliptic trajectories to circles, the hamiltonian (classical and quantum)  presents as 
$$
H=\tfrac{1}{2} (p^2+x^2).
$$
Consequently, remarkably, the action of the hamiltonian on any phase space configuration, classical Liouville density or Wigner function, amounts to 
$$
\{\{ H,f\}\}= \{ H,f\}= (x\partial_p- p\partial_x)~ f(x,p),
$$
readily recognized as rigid rotation in phase space, since, exceptionally  here, the quantum Moyal Brackets 
$$\{\{H,f\}\}=\left (H\left (x+{i\hbar\over 2} \overrightarrow{\partial }_{p} ~,~ 
   p-{i\hbar\over 2} \overrightarrow{\partial }_{x}\right ) ~ f(x,p) -H\left (x-{i\hbar\over 2} \overrightarrow{\partial }_{p} ~,~ 
   p+{i\hbar\over 2} \overrightarrow{\partial }_{x}\right ) ~ f(x,p)\right )/i\hbar
   $$
collapse to classical Poisson Brackets (!). 
Rigid rotation around a circle, of course, has integral eigenvalues, as its angular eigenfunctions are single-valued exponentials. This is a feature of both the Liouville operator and the QM Hamiltonian. 
The Fourier modes of the classical propagator also involve integrally spaced energy levels, but in classical mechanics this type of separation of variables is not necessary, or useful. Standard phase-space QM reviews further solve the radial $\star$-genvalue stationary equation (a Laguerre equation, not a Hermite one!) to confirm the eigenvalues are necessarily integrally spaced, as  demonstrated here.
Note the integral eigenvalue spectrum is not a mere consequence of phase space U(1) symmetry: it is a feature of the Hamiltonian being the very generator of rotations itself!  A Hamiltonian such as $(p^2+x^2)^a_\star$ corresponding to the a-th power of the Hilbert space oscillator Hamiltonian would still be symmetric under x-p rotations, but would have a spectrum of $(n+1/2)^a$, of course.
The argument, sadly, does not prove uniqueness of this hamiltonian for that spectrum, but a reasonable presumption given such a spectrum is the existence of a map proving equivalence of the fortunate hamiltonian to that of the oscillator, above. 
A: Yes. There is a simple explanation.
The classical harmonic oscillator has a well defined frequency $\omega$, independent of initial conditions. This can only happen if the quantum system has precisely equally spaced energies with gap $\hbar \omega$.
The reasons for this is that motion in a quantum system can only happen if more than one energy level is occupied. The frequency of the motion is then set by the energy difference of the different occupied energy levels.
Thus if every classical trajectory has the same frequency $\omega$, then any combination of quantum states must also have frequency differences that are multiples of $\hbar\omega$.
The argument can be made even more refined by thinking about the existence of coherent states and semi-classical superpositions in the quantum system, but the basic argument is as above.
The reverse argument can also be used to conclude that the hydrogen atom will not have equally spaced energy levels, since the frequency for motion in a Coloumb potential depends on the initial conditions.
A: This is probably not a physical explanation but it's worth to track the speciality of the Harmonic oscillator that leads to equispaced energy levels in quantum mechanics.
$\bullet$ A Hamiltonian of the form $H=\hbar(a^\dagger a+1/2)$ leads to equally space energy levels.
$\bullet$ The Hamiltonian of a simple harmonic oscillator is given by $$H=\frac{1}{2}(p^2+x^2)$$ in units where $m=\omega=1$. Both the kinetic and the potential energies are quadratic. Whenever the potential is quadratic, the Hamiltonian can be reduced to $$H=\hbar(a^\dagger a+1/2)$$ which is not possible if the potential $V(x)$ is not quadratic.
