# N=1 SUSY Casimir operator

I followed a book on SUSY where the SUSY Casimir operator, modified from the Pauli-Lubanski vector, is stated and then proven to commute with all SUSY generators. However, I'm wondering whether there's a systematic approach of deriving the Casimir operator, other than trial-and-error? Is there a general method of deriving Casimir operators for an arbitrary graded Lie algebra?

Generally, the idea is to start with [T_a, T_b] = f_ab^c T_c, form adjoint rep ad_T_a = J_a, then form the Killing form g_ab = Tr(J_aJ_b), find its inverse g^ab such that g^ab g_bc = del^a_c, form the dual basis elements J^a = g^ab J_b, and then finally, the Casimir invariant is given by C = J_a J^a = g^ab J_a J_b. Looking at the definition of g_ab, we see that this is a symmetric metric, as one would hope for. Other representations of C can be formed as well; it'd simply be C = g^ab T_a T_b.

More details here: http://en.wikipedia.org/wiki/Casimir_invariant

Edit: What I said applies to algebras. I assume you would need to modify this some for superalgebras; Tr -> Str and so on. I don't think anything out of the ordinary happens beyond that, but I'm not positive. I'm less familiar with superalgebras.