 Solution

On a suitable dense domain of smooth rapidly decreasing wave functions, the canonical commutator is $[x,p]=i\hslash$. Expanding the two ladder operators givesTheir sum and difference giveDirect multiplication, retaining the operator order, givesUsing ${\beta }^2=m\omega /\hslash$ identifies the Hamiltonian operator:For any admissible wave function, adjointness givesFor the conventional normalized state $\Vert \Psi \Vert =1$, this is the stated lower bound $E\ge \hslash \omega /2$. For an unnormalized function the norm factor must be retained; rescaling a wave function cannot leave its unnormalized energy integral bounded below by a fixed positive number.

Equality holds exactly when $a{\Psi }_0=0$. In position representation this equation is $(\beta x+{\beta }^{-1}{\partial }_x){\Psi }_0=0$, whose normalized square-integrable solution isThis construction establishes existence of a state attaining the bound, and hence proves that it is the ground state.

The commutator $[a,a^{†}]=1$ implies $[H,a^{†}]=\hslash \omega a^{†}$. Repeatedly commuting through the product givesThese vectors are nonzero: the identity $a(a^{†}{)}^n{\Psi }_0=n(a^{†}{)}^{n-1}{\Psi }_0$ gives their norm squared $n!$ by induction. ThereforeFor time-dependent stationary states multiply each spatial eigenfunction by $e^{-i{E}_{n}t/\hslash }$.