If $R$ and $S$ are symmetric covariant tensor fields of rank $r$ and $s$, respectively, then T = KillingBracket(g, R, S) is a symmetric covariant tensor field of rank $r\+s-1$. If $R$ and $S$ correspond to Killing tensors for the metric $g$, then $T$ is also a Killing tensor.

KillingBracket(g, R, S) can be defined in terms of the Schouten bracket for symmetric contravariant tensors by using the metric $g$ to raise all the indices on the tensors $R$ and $S$ to obtain contravariant tensors $U$ and $V$, then computing the Schouten bracket W = TensorBrackets(g, U, V, "Schouten") and then lowering all the indices of $W$. An explicit formula for the KillingBracket can be found in the article of N. M. J. Woodhouse.

This command is part of the DifferentialGeometry:-Tensor package, and so can be used in the form KillingBracket(...) only after executing the commands with(DifferentialGeometry), with(Tensor) in that order.  It can always be used in the long form DifferentialGeometry:-Tensor:-KillingBracket.

$\mathrm{with}\left(\mathrm{DifferentialGeometry}\right)\:$$\mathrm{with}\left(\mathrm{Tensor}\right)\:$

$\mathrm{DGsetup}\left(\left[x\,y\,z\right]\,M\right)\:$

$g≔\mathrm{evalDG}\left(x\mathrm{dx}\&s\mathrm{dy}+x\mathrm{dz}\&t\mathrm{dz}\right)$

$g:=\frac{x}{2}\mathrm{dx}\mathrm{dy}+\frac{x}{2}\mathrm{dy}\mathrm{dx}+x\mathrm{dz}\mathrm{dz}$

$R≔\mathrm{evalDG}\left(z{x}^2\mathrm{dx}\&t\mathrm{dx}-\frac{1}{2}{x}^3\mathrm{dx}\&t\mathrm{dz}-\frac{1}{2}{x}^3\mathrm{dz}\&t\mathrm{dx}\right)$

$R:=z{x}^2\mathrm{dx}\mathrm{dx}-\frac{x^3}{2}\mathrm{dx}\mathrm{dz}-\frac{x^3}{2}\mathrm{dz}\mathrm{dx}$

$S≔\mathrm{evalDG}\left(\frac{1}{2}{x}^2\mathrm{dx}\&t\mathrm{dz}+\frac{1}{2}{x}^2\mathrm{dz}\&t\mathrm{dx}\right)$

$S:=\frac{x^2}{2}\mathrm{dx}\mathrm{dz}+\frac{x^2}{2}\mathrm{dz}\mathrm{dx}$

$\mathrm{CheckKillingTensor}\left(g\,\left[R\,S\right]\right)$

$\left[0\mathrm{dx}\mathrm{dx}\mathrm{dx}\,0\mathrm{dx}\mathrm{dx}\mathrm{dx}\right]$

$T≔\mathrm{KillingBracket}\left(g\,R\,S\right)$

$T:=-x^3\mathrm{dx}\mathrm{dx}\mathrm{dx}$

$\mathrm{CheckKillingTensor}\left(g\,T\right)$

$0\mathrm{dx}\mathrm{dx}\mathrm{dx}\mathrm{dx}$