manyspin

◼ manyspin[{a, b, c}] returns X, Y, and Z components of the total spin, as well as the total spin operator squared, for levels corresponding to a list of operators given as the argument.

◼ Operators must be expressed in a form which excludes the first index (creation/annihilation) and the last index (spin).

◼ See also: spinx, spiny, spinz, spinss, manyisospin.

Further examples

Load sneg and declared c and d to be fermionic operators

<< sneg`sneg`;
snegfermionoperators[c, d];

{(Subscript[S, x]), (Subscript[S, y]), (Subscript[S, z]), SS} = 
      manyspin[{c[], d[]}]
komutator[(Subscript[S, x]), (Subscript[S, y])] == I (Subscript[S, z]) // 
      ExpandAll
komutator[(Subscript[S, y]), (Subscript[S, z])] == I (Subscript[S, x]) // 
      ExpandAll
komutator[(Subscript[S, z]), (Subscript[S, x])] == I (Subscript[S, y]) // 
      ExpandAll
komutator[SS, (Subscript[S, x])] == 0 // ExpandAll
komutator[SS, (Subscript[S, y])] == 0 // ExpandAll
komutator[SS, (Subscript[S, z])] == 0 // ExpandAll