As I'd mentioned in the email, you shouldn't use the `.doit()` here, since symbolic `j` values will raise errors when evaluating the Wigner-D matrix element, e.g. `qapply(Rotation(alpha, beta, gamma) * JzKet(j, m))`.