We have the usual Hamiltonian for the 1D Harmonic Oscillator: $\hat{H_{0}}=\frac{\hat{P^2}}{2m} + \frac{1}{2}m \omega \hat{X^2}$
Now a new term has been added to the Hamiltonian, $\hat{H} = \hat{H_0} + \mu B\hat{S_z}$
The system has a spin degree of freedom with $s=\frac{1}{2}$
What are the new eigenstates of the Hamiltonian and what are the energy eigenvalues? We have been denoting the stationary states as $|n,s_z\rangle$ and the non-spin eigenvalues as usual $E = h \omega (n+1/2)$
Can anyone give me some help with this question? I'm guessing that the eigenstates don't change and then it's straightforward to find the eigenvalues but I'm not sure if this is correct (or why it would be correct).