No estoy muy puesto con esto de los Casimires de orden superior, pero yo lo que haría es pasar a la base de \( \mathfrak{so}(4) \cong \mathfrak{so}(3)\oplus \mathfrak{so}(3) \) dada por los elementos \( B_{+,i} = \frac{J_i+K_i}{2} \) y \( B_{-,i} = \frac{J_i-K_i}{2} \) de las que hablábamos en otro hilo. La ventaja de esta base es que cada conjunto por separado forma base de \( \mathfrak{so}(3) \) y además los conmutadores de elementos de conjuntos distintos son cero.
Entonces, tienes un casimir (en el sentido de un elemento del centro del álgebra universal envolvente) en cada factor, dados por \( B_+^2 = B_{+,x}^2+B_{+,y}^2+B_{+,z}^2 \) y \( B_-^2 = B_{-,x}^2+B_{-,y}^2+B_{-,z}^2 \). Si no me equivoco, estos dos elementos te generan todo el centro del álgebra universal envolvente de \( \mathfrak{so}(4) \). Si los sumas, obtienes (esencialmente) el Casimir cuadrático que mencionas, \( B_+^2+B_-^2 = \frac{J^2+K^2}{2} \) (aparece un factor \( 2 \) por ahí), y el Casimir de orden superior en esta base yo diría que es \( B^2_+B^2_- = \sum_{i,j} B^2_{+,i}B^2_{-,j} \).