epsilon_to_delta

Replace a product of two epsilon tensors with a generalised delta according to the expression \begin{equation} \epsilon^{r_1\cdots r_{d}} \epsilon_{s_1\cdots s_{d}} = \frac{1}{\sqrt{|g|}}\varepsilon^{r_1 \cdots r_{d}} \sqrt{|g|}\varepsilon_{s_1\cdots s_{d}} = {\rm sign}(g)\, d!\, \delta^{r_1 \cdots r_{d}}_{s_1\cdots s_{d}}\, , \end{equation} where ${\rm sign}(g)$ denotes the signature of the metric $g$ used to raise/lower the indices (see EpsilonTensor for conventions on the epsilon tensor). When the indices are not ocurring up/down as in this expression, and the index position is not free, metric objects will be generated instead. Here is an example:

Remember that if the result is a generalised delta, you can expand it in terms of normal deltas using expand_delta,

In order for this algorithm to work, you need to make sure that the epsilon symbols in your expression are declared as EpsilonTensor and that these declarations involve a specification of the delta and/or metric symbols.

As you can see from this example, contracted indices inside the generalised delta are automatically eliminated, as the algorithm reduce_gendelta is called automatically; if you do not want this use the optional argument reduce=False.

Note that the results typically depend on the signature of the space-time, which can be introduced through the optional metric argument of the EpsilonTensor property. Compare the two examples below:

Note that you need to specify the full symbol for the metric, including the indices, whereas the Kronecker delta argument only requires the name without the indices (because a contraction can generate generalised Kronecker delta symbols with any number of indices).