U=g_s\mu_B\vec{s}\cdot\vec{B}