where S is the spin operator with eigenvalue s and L is the angular momentum operator with eigenvalue l. g J is the Landé g-factor, calculated as g J = 1 + j ( j + 1 ) + s ( s + 1 ) − − l ( l + 1 ) 2 j ( j + 1 )