It is well-known in quantum mechanics that the harmonic oscillator Hamiltonian given by $\mathcal{H} = -\frac{1}{2}\frac{d^2}{dx^2} + \frac{1}{2}x^2 - \frac{1}{2}$ admits a basis of eigenfunctions on $L^2(\Bbb R,dx)$. There are many proofs of this ranging from hard analysis (brute force proving density) to Bargmann transform techniques to showing that the resolvent is compact. An important part of the analysis is that the Gaussian $e^{-\frac{x^2}{2}}$ is an eigenfunction with eigenvalue $0$.
Associated to this system is a very nice Lie algebra. Define the operator $a$ by
$$a = \frac{1}{\sqrt{2}}\left(\frac{d}{dx}+x\right)$$
and its (formal - I'm not going through the awful nitty gritty details of domains of definition, etc) adjoint
$$a^* = \frac{1}{\sqrt{2}}\left(-\frac{d}{dx}+x\right).$$
$\mathcal{H}$ has the following representation:
$$\mathcal{H} = a^*a.$$
Moreover, this triple of operators satisfies (on sufficiently nice functions):
$$[\mathcal{H},a^*] = a^*, \qquad [\mathcal{H},a] = -a, \qquad [a,a^*] = -1.$$
The quadruplet $\mathcal{H},a^*,a,1$ generates a Lie algebra. Because of these relations, applying $a^*$ repeatedly to the Gaussian leads to successive eigenfunctions (that they are in $L^2$ needs to be checked, but it is straightforward by induction) and indeed these eigenfunctions form a basis.
My question is this: does a representation of the Lie algebra
$$[a^*a, a^*] = a^*, \qquad [a^*a, a] = -a, \qquad [a,a^*] = -1$$
on an abstract separable Hilbert space $\mathfrak{H}$ (not necessarily $L^2(\Bbb R,dx)$) necessarily admit a basis of eigenfunctions for $a^*a$ assuming that $\ker a\neq \{0\}$? That is, if $\{\psi_i:i\in I\}$ denotes a basis for the kernel of $a$, does $\{(a^*)^m\psi_i:m\in\Bbb N_0,i\in I\}$ form a basis for $\mathfrak{H}$? Or is this somehow a happy accident of working on $\Bbb R$ (or, as we know, more generally $\Bbb R^n$)?