We tackled these issues as follows. The symmetry of one-and two-electron parts of the spin-orbit operator is the same as the symmetry of ˆ L • ˆ S. Because dot product is rotationally invariant, we compute matrix elements in different molecular orientations and trace how the states change. Using properties of Wigner's D-matrices and the resolution-of-the-identity, one can establish the following relation: