It follows from the fact that the eignevalues of the hermitian operator [itex]\hat{n} = \hat{a}^{\dagger} \, \hat{a}[/itex] are non-negative! Namely, let:
[tex]
\hat{n} \, \vert \nu \rangle = \nu \, \vert \nu \rangle[/tex]
Assuming the eigenket [itex]\vert \nu \rangle[/itex] is normalized to unity ([itex]\langle \nu \vert \nu \rangle = 1[/itex]), we have:
[tex]
\nu = \langle \nu \vert \hat{n} \vert \nu \rangle = \langle \nu \vert \hat{a}^{\dagger} \, \hat{a} \vert \nu \rangle[/tex]
Now, insert a compete orthonormal basis [itex]\lbrace k' \rbrace[/itex] between [itex]\hat{a}[/itex] and [itex]\hat{a}^\dagger[/itex]:
[tex]
\nu = \sum_{k'}{\langle \nu \vert \hat{a}^{\dagger} \vert k' \rangle \, \langle k' \vert \hat{a} \vert \nu \rangle} = \sum_{k'}{\left\vert \langle k' \vert \hat{a} \vert \nu \rangle \right\vert^2}[/tex]
Each of the summands is a non-zero real number, being a modulus of a complex number. Therefore, [itex]\nu \ge 0[/itex].
Now, from the commutation relation [itex]\left[ \hat{a}, \hat{a}^\dagger \right] = 1[/itex], it follows that the ket:
[tex]
\hat{a} \, \vert \nu \rangle[/tex]
is also an eigenket of [itex]\hat{n}[/itex], corresponding to an eigenvalue [itex]\nu - 1[/itex]. By inductive reasoning, we conclude that [itex]\nu - n, \forall n \in \mathbb{N}[/itex] are also eigenvalues of the number operator.
So, let us suppose that [itex]N + 1 > \nu \ge N \ge 0, \ N \in \mathbb{N}[/itex]. This means that [itex]\nu - N - 1 < 0[/itex], which is not allowed (because all the eigenvalues are non-negative!). The only way out is if [itex]\left( \hat{a} \right)^{N + 1} \, \vert \nu \rangle = 0[/itex] (not the ground state [itex]\vert 0 \rangle[/itex], but the zero ket in the Hilbert space!). But, then apply [itex]\left( \hat{a}^\dagger \right)^{N + 1}[/itex] from the left, and use:
===================================================================================
[tex]
\hat{X}_m \equiv \left[ \hat{a}^\dagger, \hat{a}^m \right][/tex]
[tex]
\hat{X}_{m + 1} = \left[ \hat{a}^\dagger, \hat{a} \, \hat{a}^m \right] = \hat{a} \, \left[ \hat{a}^\dagger, \hat{a}^m \right] + \left[ \hat{a}^\dagger, \hat{a} \right] \, \hat{a}^m = \hat{a} \, \hat{X}_m - \hat{a}^m[/tex]
[tex]
\hat{X}_0 = 0[/tex]
[tex]
\hat{X}_1 = \hat{a} \, 0 - \hat{1} = -\hat{1}[/tex]
[tex]
\hat{X}_2 = \hat{a} \, (-\hat{1}) - \hat{a} = -2 \, \hat{a}[/tex]
...
Lemma:
[tex]
\boxed{\hat{X}_m = -m \, \hat{a}^{m - 1}, \ m \ge 1}[/tex]
Proof:
[tex]
\hat{X}_{m + 1} = \hat{a} \, \left(-m \, \hat{a}^{m - 1} \right) - \hat{a}^{m} = -(m + 1) \, \hat{a}^{m}[/tex]
By Principle of Mathematical Induction, Q.E.D.
===================================================================================
Then, we have:
[tex]
\begin{array}{lcl}<br />
\left(\hat{a}^\dagger\right)^{N + 1} \, \hat{a}^{N + 1} \, \vert \nu \rangle & = & 0 \\<br />
<br />
\left( \hat{a}^\dagger \right)^N \, \hat{a}^\dagger \, \hat{a}^{N + 1} \, \vert \nu \rangle & = & 0 \\<br />
<br />
\left( \hat{a}^\dagger \right)^N \, \left( \hat{a}^{N + 1} \, \hat{a}^\dagger + X_{N + 1} \right) \, \vert \nu \rangle & = & 0 \\<br />
<br />
\left[ \left( \hat{a}^\dagger \right)^N \, \hat{a}^{N} \, \hat{a} \, \hat{a}^\dagger - (N + 1) \, \left( \hat{a}^\dagger\right )^{N} \, \hat{a}^{N} \right] \, \vert \nu \rangle & = & 0 \\<br />
<br />
\left[ \left( \hat{a}^\dagger \right)^N \, \hat{a}^{N} \left( \hat{n} + 1 \right) - (N + 1) \, \left( \hat{a}^\dagger \right)^N \, \hat{a}^{N} \right] \, \vert \nu \rangle & = & 0 \\<br />
<br />
\left( \nu - N \right) \, \left( \hat{a}^\dagger \right)^N \, \hat{a}^{N} \, \vert \nu \rangle & = & 0 \\<br />
<br />
\ldots & & \\<br />
<br />
\left( \nu - N \right) \, \left( \nu - N + 1\right) \, \ldots \nu \, \vert \nu \rangle & = & 0<br />
\end{array}[/tex]
Now, from the initial assumption [itex]\vert \nu \rangle \neq 0[/itex] (not a trivial eigenket), and [itex]N \le \nu < N + 1[/itex], we conclude that [itex]\nu = N[/itex].
This means that the only eigenvalues of the number operator [itex]\hat{n} = \hat{a}^\dagger \, \hat{a}[/itex] are the non-negative integers. Q.E.D.