I'm struggling with the following problem:
Consider a system of an arbitrary number of indistinguishable bosonic particles. The system has two sites and $a_i^{\dagger}$ and $a_i$ are the corresponding creating and annihiliation operators. If $$\displaystyle \hat{H}=-a_1^{\dagger}a_2-a_2^{\dagger}a_1+\sum_{i=1,2}a_i^{\dagger}a_i(a_i^{\dagger}a_i-1),$$ show that applying this Hamiltonian to state does not affect the overall number of particles in that state.
I use occupation number representation and that
$$a_i|n_1,\dots,n_i,\dots ⟩=\sqrt{n_i}|n_1,\dots,n_i-1,\dots⟩ $$ and $$ a_i^{\dagger}|n_1,\dots,n_i,\dots ⟩=\sqrt{n_i+1}|n_1,\dots,n_i+1,\dots⟩$$
So, \begin{align} \hat{H}|n_1,n_2⟩ &= \left[ -a_1^{\dagger}a_2-a_2^{\dagger}a_1+\sum_{i=1,2}a_i^{\dagger}a_i(a_i^{\dagger}a_i-1) \right] |n_1,n_2⟩ \\ &= -a_1^{\dagger}a_2 |n_1,n_2⟩ -a_2^{\dagger}a_1 |n_1,n_2⟩ +\sum_{i=1,2}a_i^{\dagger}a_i(a_i^{\dagger}a_i-1) |n_1,n_2⟩ \\&= -\sqrt{n_1+1}\sqrt{n_2}|n_1+1,n_2-1⟩ -\sqrt{n_1}\sqrt{n_2+1} |n_1-1,n_2+1⟩ \\&\quad + \frac{1}{2}n_1^2|n_1,n_2⟩ +\frac{1}{2}n_2^2|n_1,n_2⟩-\frac{1}{2}n_1|n_1,n_2⟩-\frac{1}{2}n_2|n_1,n_2⟩ \end{align}
Am I able to just say $|n_1+1,n_2-1⟩ + |n_1-1,n_2+1⟩ =|n_1,n_2⟩$? Im not sure how to interpret the answer I found.
Is anyone able to do this by computing the action of $\hat{H}$ on an arbitrary occupation number representation state?