I am trying to calculate the trace of $4$ generators of $U(N)$, i.e. $\text{Tr}\ T^a T^b T^c T^d$? I found a plausible result, but I would also like to show that the result is cyclic with respect to the indices, as the trace is. Here is first my calculation of $\text{Tr}\ T^a T^b T^c T^d$:
We start with the following:
$$T^a T^b = \frac{1}{2} \left( [T^a,T^b]+\lbrace T^a,T^b\rbrace \right) \tag{1}$$
with
$$[T^a,T^b] := if^{abc} T^c \tag{2}$$
and
$$\lbrace T^a, T^b \rbrace = d^{abc} T^c \tag{3}$$
which can be shown using the completeness relation for $U(N)$:
$$T^a_{ij}T^a_{lk} = \frac{1}{2}\delta_{ik}\delta_{jl} \tag{4}$$
and the definition of $d^{abc}$:
$$d^{abc}:= 2\ \text{Tr} (\lbrace T^a,T^b\rbrace T^c) \tag{5}$$
Putting everything together, I find:
$$T^a T^b T^c T^d = \frac{1}{4} (if^{abe}T^e + d^{abe}T^e)(if^{cdf}T^f + d^{cdf}T^f) \tag{6}$$
and thus
$$\text{Tr}\ T^a T^b T^c T^d = \frac{1}{8} (-f^{abe}f^{cde} + if^{abe}d^{cde} + id^{abe}f^{cde} + d^{abe}d^{cde}) \tag{7}$$
Now the issue is the following: if I change the indices cyclically, I get:
$$\text{Tr}\ T^b T^c T^d T^a = \frac{1}{8} (-f^{bce}f^{dae} + if^{bce}d^{dae} + id^{bce}f^{dae} + d^{bce}d^{dae}) \tag{8}$$
which should be equal to $(7)$ because of the cyclicity of the trace. How can I show that?