You can use the second syntax of [`Transpose`](https://reference.wolfram.com/language/ref/Transpose.html) and it should work like the code below, but I am not sure for the each 4-length `List` as the second argument of `Transpose` ``` n = {a, b}; id = IdentityMatrix[2]; id4 = TensorProduct[id, n, n]; mat = TensorProduct[n, n, n, n] - (id4 + Transpose[id4, {1, 3, 2, 4}] + Transpose[id4, {1, 4, 3, 2}] + Transpose[id4, {3, 1, 4, 2}] + Transpose[id4, {3, 1, 2, 4}] + Transpose[id4, {3, 4, 1, 2}])/6 ```