Consider a spin-1/2 particle in a magnetic field (say in z direction) and in a harmonic potential. For the 3D harmonic oscillator component, The Hamiltonian $H_1= \frac{p^2}{2m}+\frac{1}{2}m\omega ^2r^2.$ For the spin component, the Hamiltonian $H_2=-\gamma B_z S_z$, where $\gamma$ is the gyroscopic ratio.
Questions:
Is it possible to represent the an eigenstate of the system as tensor product of eigenstate of each of the two Hamiltonian? I.e. Is the eigenstate $\left|n_x,n_y,n_z\right>\otimes \left|1/2,m_s\right>$, where $\left|n_x,n_y,n_z\right>$ is an eigenstate of $H_1$, and $\left|1/2,m_s\right>$ is an eigenstate of $H_2$? Is the explicit form of a state, for example, $\left|0,0,0\right>\otimes \left|1/2,+1/2\right>=(\frac{m\omega}{\pi\hbar})^{1/4}exp(-\frac{m\omega}{2\hbar} r^2)\otimes \begin{pmatrix}1\\0\end{pmatrix}$?
If answer to 2 is yes, are the states $\left|0,0,0\right>\otimes \left|1/2,+1/2\right>$ and $\left|0,0,0\right>\otimes \left|1/2,-1/2\right>$ orthogonal? My guess would be that the inner product of the two tensor states is the inner product of each of the component, and since for the second component $\left<1/2,+1/2\middle| 1/2,-1/2\right>=0$, the two tensor states should be orthogonal?
I have been trying to answer the last question, but different reasoning (some might be too naive) seems to give contradicting answers, so I want to see if the above reasoning is valid.