Levi-Civita symbol

Short description: Antisymmetric permutation object acting on tensors

See Ricci calculus, Einstein notation, and Raising and lowering indices for the index notation used in the article.

In mathematics, particularly in linear algebra, tensor analysis, and differential geometry, the Levi-Civita symbol or Levi-Civita epsilon represents a collection of numbers; defined from the sign of a permutation of the natural numbers 1, 2, ..., n, for some positive integer n. It is named after the Italian mathematician and physicist Tullio Levi-Civita. Other names include the permutation symbol, antisymmetric symbol, or alternating symbol, which refer to its antisymmetric property and definition in terms of permutations.

The standard letters to denote the Levi-Civita symbol are the Greek lower case epsilon ε or ϵ, or the Latin lower case e for the e-matrix in the Cartesian reference system. Index notation allows one to display permutations in a way compatible with tensor analysis: [math]\displaystyle{ \varepsilon_{i_1 i_2 \dots i_n} }[/math] where each index i₁, i₂, ..., i_n takes values 1, 2, ..., n. There are nⁿ indexed values of ε_i1i2...in, which can be arranged into an n-dimensional array. The key defining property of the symbol is total antisymmetry in the indices. When any two indices are interchanged, equal or not, the symbol is negated: [math]\displaystyle{ \varepsilon_{\dots i_p \dots i_q \dots } = -\varepsilon_{\dots i_q \dots i_p \dots } . }[/math]

If any two indices are equal, the symbol is zero. When all indices are unequal, we have: [math]\displaystyle{ \varepsilon_{i_1 i_2 \dots i_n} = (-1)^p \varepsilon_{1 \, 2 \, \dots n} , }[/math] where p (called the parity of the permutation) is the number of pairwise interchanges of indices necessary to unscramble i₁, i₂, ..., i_n into the order 1, 2, ..., n, and the factor (−1)^p is called the sign, or signature of the permutation. The value ε_{1 2 ... n} must be defined, else the particular values of the symbol for all permutations are indeterminate. Most authors choose ε_{1 2 ... n} = +1, which means the Levi-Civita symbol equals the sign of a permutation when the indices are all unequal. This choice is used throughout this article.

The term "n-dimensional Levi-Civita symbol" refers to the fact that the number of indices on the symbol n matches the dimensionality of the vector space in question, which may be Euclidean or non-Euclidean, for example, [math]\displaystyle{ \mathbb{R}^3 }[/math] or Minkowski space. The values of the Levi-Civita symbol are independent of any metric tensor and coordinate system. Also, the specific term "symbol" emphasizes that it is not a tensor because of how it transforms between coordinate systems; however it can be interpreted as a tensor density.

The Levi-Civita symbol allows the determinant of a square matrix, and the cross product of two vectors in three-dimensional Euclidean space, to be expressed in Einstein index notation.

Definition

The Levi-Civita symbol is most often used in three and four dimensions, and to some extent in two dimensions, so these are given here before defining the general case.

Two dimensions

In two dimensions, the Levi-Civita symbol is defined by: [math]\displaystyle{ \varepsilon_{ij} = \begin{cases} +1 & \text{if } (i, j) = (1, 2) \\ -1 & \text{if } (i, j) = (2, 1) \\ \;\;\,0 & \text{if } i = j \end{cases} }[/math] The values can be arranged into a 2 × 2 antisymmetric matrix: [math]\displaystyle{ \begin{pmatrix} \varepsilon_{11} & \varepsilon_{12} \\ \varepsilon_{21} & \varepsilon_{22} \end{pmatrix} = \begin{pmatrix} 0 & 1 \\ -1 & 0 \end{pmatrix} }[/math]

Use of the two-dimensional symbol is common in condensed matter, and in certain specialized high-energy topics like supersymmetry and twistor theory, where it appears in the context of 2-spinors.

Three dimensions

For the indices (i, j, k) in ε_ijk, the values 1, 2, 3 occurring in the      cyclic order (1, 2, 3) correspond to ε = +1, while occurring in the      reverse cyclic order correspond to ε = −1, otherwise ε = 0.

In three dimensions, the Levi-Civita symbol is defined by: [math]\displaystyle{ \varepsilon_{ijk} = \begin{cases} +1 & \text{if } (i,j,k) \text{ is } (1,2,3), (2,3,1), \text{ or } (3,1,2), \\ -1 & \text{if } (i,j,k) \text{ is } (3,2,1), (1,3,2), \text{ or } (2,1,3), \\ \;\;\,0 & \text{if } i = j, \text{ or } j = k, \text{ or } k = i \end{cases} }[/math]

That is, ε_ijk is 1 if (i, j, k) is an even permutation of (1, 2, 3), −1 if it is an odd permutation, and 0 if any index is repeated. In three dimensions only, the cyclic permutations of (1, 2, 3) are all even permutations, similarly the anticyclic permutations are all odd permutations. This means in 3d it is sufficient to take cyclic or anticyclic permutations of (1, 2, 3) and easily obtain all the even or odd permutations.

Analogous to 2-dimensional matrices, the values of the 3-dimensional Levi-Civita symbol can be arranged into a 3 × 3 × 3 array:

where i is the depth (blue: i = 1; red: i = 2; green: i = 3), j is the row and k is the column.

Some examples: [math]\displaystyle{ \begin{align} \varepsilon_{\color{BrickRed}{1} \color{Violet}{3} \color{Orange}{2}} = -\varepsilon_{\color{BrickRed}{1} \color{Orange}{2} \color{Violet}{3}} &= - 1 \\ \varepsilon_{ \color{Violet}{3} \color{BrickRed}{1} \color{Orange}{2}} = -\varepsilon_{ \color{Orange}{2} \color{BrickRed}{1} \color{Violet}{3}} &= -(-\varepsilon_{\color{BrickRed}{1} \color{Orange}{2} \color{Violet}{3}}) = 1 \\ \varepsilon_{ \color{Orange}{2} \color{Violet}{3} \color{BrickRed}{1}} = -\varepsilon_{\color{BrickRed}{1} \color{Violet}{3} \color{Orange}{2}} &= -(-\varepsilon_{\color{BrickRed}{1} \color{Orange}{2} \color{Violet}{3}}) = 1 \\ \varepsilon_{ \color{Orange}{2} \color{Violet}{3} \color{Orange}{2}} = -\varepsilon_{ \color{Orange}{2} \color{Violet}{3} \color{Orange}{2}} &= 0 \end{align} }[/math]

Four dimensions

In four dimensions, the Levi-Civita symbol is defined by: [math]\displaystyle{ \varepsilon_{ijkl} = \begin{cases} +1 & \text{if }(i,j,k,l) \text{ is an even permutation of } (1,2,3,4) \\ -1 & \text{if }(i,j,k,l) \text{ is an odd permutation of } (1,2,3,4) \\ \;\;\,0 & \text{otherwise} \end{cases} }[/math]

These values can be arranged into a 4 × 4 × 4 × 4 array, although in 4 dimensions and higher this is difficult to draw.

Some examples: [math]\displaystyle{ \begin{align} \varepsilon_{\color{BrickRed}{1} \color{RedViolet}{4}\color{Violet}{3} \color{Orange}{\color{Orange}{2}}} = -\varepsilon_{\color{BrickRed}{1} \color{Orange}{\color{Orange}{2}} \color{Violet}{3} \color{RedViolet}{4}} &= - 1\\ \varepsilon_{\color{Orange}{\color{Orange}{2}} \color{BrickRed}{1} \color{Violet}{3} \color{RedViolet}{4}} = -\varepsilon_{\color{BrickRed}{1} \color{Orange}{\color{Orange}{2}} \color{Violet}{3} \color{RedViolet}{4}} &= -1\\ \varepsilon_{\color{RedViolet}{4} \color{Violet}{3} \color{Orange}{\color{Orange}{2}} \color{BrickRed}{1}} = -\varepsilon_{\color{BrickRed}{1} \color{Violet}{3} \color{Orange}{\color{Orange}{2}} \color{RedViolet}{4}} &= -(-\varepsilon_{\color{BrickRed}{1} \color{Orange}{\color{Orange}{2}} \color{Violet}{3} \color{RedViolet}{4}}) = 1\\ \varepsilon_{\color{Violet}{3} \color{Orange}{\color{Orange}{2}} \color{RedViolet}{4} \color{Violet}{3}} = -\varepsilon_{\color{Violet}{3} \color{Orange}{\color{Orange}{2}} \color{RedViolet}{4}\color{Violet}{3}} &= 0 \end{align} }[/math]

Generalization to n dimensions

More generally, in n dimensions, the Levi-Civita symbol is defined by: [math]\displaystyle{ \varepsilon_{a_1 a_2 a_3 \ldots a_n} = \begin{cases} +1 & \text{if }(a_1, a_2, a_3, \ldots, a_n) \text{ is an even permutation of } (1, 2, 3, \dots, n) \\ -1 & \text{if }(a_1, a_2, a_3, \ldots, a_n) \text{ is an odd permutation of } (1, 2, 3, \dots, n) \\ \;\;\,0 & \text{otherwise} \end{cases} }[/math]

Thus, it is the sign of the permutation in the case of a permutation, and zero otherwise.

Using the capital pi notation Π for ordinary multiplication of numbers, an explicit expression for the symbol is: [math]\displaystyle{ \begin{align} \varepsilon_{a_1 a_2 a_3 \ldots a_n} & = \prod_{1 \leq i \lt j \leq n} \sgn (a_j - a_i) \\ & = \sgn(a_2 - a_1) \sgn(a_3 - a_1) \dotsm \sgn(a_n - a_1) \sgn(a_3 - a_2) \sgn(a_4 - a_2) \dotsm \sgn(a_n - a_2) \dotsm \sgn(a_n - a_{n-1}) \end{align} }[/math] where the signum function (denoted sgn) returns the sign of its argument while discarding the absolute value if nonzero. The formula is valid for all index values, and for any n (when n = 0 or n = 1, this is the empty product). However, computing the formula above naively has a time complexity of O(n²), whereas the sign can be computed from the parity of the permutation from its disjoint cycles in only O(n log(n)) cost.

Properties

A tensor whose components in an orthonormal basis are given by the Levi-Civita symbol (a tensor of covariant rank n) is sometimes called a permutation tensor.

Under the ordinary transformation rules for tensors the Levi-Civita symbol is unchanged under pure rotations, consistent with that it is (by definition) the same in all coordinate systems related by orthogonal transformations. However, the Levi-Civita symbol is a pseudotensor because under an orthogonal transformation of Jacobian determinant −1, for example, a reflection in an odd number of dimensions, it should acquire a minus sign if it were a tensor. As it does not change at all, the Levi-Civita symbol is, by definition, a pseudotensor.

As the Levi-Civita symbol is a pseudotensor, the result of taking a cross product is a pseudovector, not a vector.

Under a general coordinate change, the components of the permutation tensor are multiplied by the Jacobian of the transformation matrix. This implies that in coordinate frames different from the one in which the tensor was defined, its components can differ from those of the Levi-Civita symbol by an overall factor. If the frame is orthonormal, the factor will be ±1 depending on whether the orientation of the frame is the same or not.

In index-free tensor notation, the Levi-Civita symbol is replaced by the concept of the Hodge dual.

Summation symbols can be eliminated by using Einstein notation, where an index repeated between two or more terms indicates summation over that index. For example,

[math]\displaystyle{ \varepsilon_{ijk} \varepsilon^{imn} \equiv \sum_{i=1,2,3} \varepsilon_{ijk} \varepsilon^{imn} }[/math].

In the following examples, Einstein notation is used.

Two dimensions

In two dimensions, when all i, j, m, n each take the values 1 and 2:

[math]\displaystyle{ \varepsilon_{ij} \varepsilon^{mn} = {\delta_i}^m {\delta_j}^n - {\delta_i}^n {\delta_j}^m }[/math]

( 1 )

[math]\displaystyle{ \varepsilon_{ij} \varepsilon^{in} = {\delta_j}^n }[/math]

( 2 )

[math]\displaystyle{ \varepsilon_{ij} \varepsilon^{ij} = 2. }[/math]

( 3 )

Three dimensions

Index and symbol values

In three dimensions, when all i, j, k, m, n each take values 1, 2, and 3:

[math]\displaystyle{ \varepsilon_{ijk} \varepsilon^{pqk}=\delta_i{}^{p}\delta_j{}^q - \delta_i{}^q\delta_j{}^p }[/math]

( 4 )

[math]\displaystyle{ \varepsilon_{jmn} \varepsilon^{imn}=2{\delta_j}^i }[/math]

( 5 )

[math]\displaystyle{ \varepsilon_{ijk} \varepsilon^{ijk}=6. }[/math]

( 6 )

Tensor Calculus

The Levi-Civita symbol [math]\displaystyle{ \varepsilon_{ijk} }[/math] resp. [math]\displaystyle{ \varepsilon^{ijk} }[/math] was already defined above in the identical Cartesian reference systems [math]\displaystyle{ \vec{e}{_i} = [\vec{e}{_0}\;\vec{e}{_1}\;\vec{e}{_2}] }[/math] resp. [math]\displaystyle{ \vec{e}{^k} = [\vec{e}{^0}\;\vec{e}{^1}\;\vec{e}{^2}] }[/math] as the permutation symbol through the skew-symmetric e-matrix [math]\displaystyle{ \varepsilon_{ijk} = e_{ijk} }[/math] resp. [math]\displaystyle{ \varepsilon^{ijk} = e^{ijk} }[/math] with numerical values [math]\displaystyle{ \pm1 }[/math]. In the Cartesian system, it should continue to be denoted with Latin indices.

However, in the following, it is now transferred for tensor calculus into two dual, skewed, non-normalized reference systems with covariant basis vectors[math]\displaystyle{ \vec{g}_\alpha = [\vec{g}{_0}\;\vec{g}_1\;\vec{g}_2] }[/math] on one hand, and contravariant basis vectors [math]\displaystyle{ \vec{g}_\alpha = [\vec{g}{_0}\;\vec{g}_1\;\vec{g}_2] }[/math] on the other hand. Both reference systems are orthogonally aligned to each other [math]\displaystyle{ \vec{g}_\alpha \cdot \vec{g}^\beta = \delta{_\alpha}^\beta }[/math] through the identity matrix [math]\displaystyle{ \delta{_\alpha}^\beta }[/math], and are always marked with Greek indices.

The transformation of tensor components between the Cartesian reference system and skewed reference systems is done through the components [math]\displaystyle{ g{_i}^\alpha \; bzw. \; g{^k}_\beta }[/math] of the contravariant basis vectors [math]\displaystyle{ \vec{g}{^\alpha}=g{_i}^{\alpha}\;\vec{e}{^i} }[/math] resp. the covariant basis vectors [math]\displaystyle{ \vec{g}{_\beta}=g{^k}{_\beta}\;{\vec{e}}{_k}, }[/math], where the orthogonality [math]\displaystyle{ g{_i}^{\alpha} * g{^i}{_\beta} = \delta{^\alpha}_\beta }[/math] arises from the above relationship.

The determinant of the covariant metric tensor [math]\displaystyle{ g_{\alpha\beta} = g{^i}{_\alpha} * g{^i}{_\beta} }[/math] is obtained from the volume form [math]\displaystyle{ g = det(g_{\alpha\beta})=det(g{^i}{_\alpha}) * det(g{^i}{_\beta}) }[/math] with [math]\displaystyle{ \sqrt{g} = det(g{^i}{_\alpha}) }[/math], and similarly, the determinant of the contravariant metric tensor [math]\displaystyle{ g^{\alpha\beta} = g{_i}{^\alpha} * g{_i}{^\beta} }[/math] is obtained from the volume form [math]\displaystyle{ 1/g = det(g^{\alpha\beta}) }[/math] with [math]\displaystyle{ 1/\sqrt{g} = det(g{_i}{^\alpha}) }[/math].

As shown above, the Levi-Civita symbol in the Cartesian reference system can be obtained through the following determinant formation from

the Kronecker delta: [math]\displaystyle{ \begin{align} \varepsilon_{ijk} = \begin{vmatrix} \delta{^0}{_i} & \delta{^0}{_j} & \delta{^0}{_k} \\ \delta{^1}{_i} & \delta{^1}{_j} & \delta{^1}{_k} \\ \delta{^2}{_i} & \delta{^2}{_j} & \delta{^2}{_k} \\ \end{vmatrix} = e_{ijk} \;\; resp. \;\; \varepsilon^{ijk} = \begin{vmatrix} \delta{_0}{^i} & \delta{_0}{^j} & \delta{_0}{^k} \\ \delta{_1}{^i} & \delta{_1}{^j} & \delta{_1}{^k} \\ \delta{_2}{^i} & \delta{_2}{^j} & \delta{_2}{^k} \\ \end{vmatrix} = e^{ijk}. \end{align} }[/math]

The [math]\displaystyle{ \varepsilon }[/math]-tensors in the skewed reference system are obtained through corresponding tensor transformations of the Levi-Civita symbol or alternatively by multiplying the e-matrix by the volume element [math]\displaystyle{ \sqrt{g} }[/math] for both

the covariant [math]\displaystyle{ \varepsilon }[/math]-tensor: [math]\displaystyle{ \;\;\;\;\; \varepsilon{_\alpha}{_\beta}{_\gamma} = \varepsilon{_i}{_j}{_k} * g{^i}{_\alpha}* g{^j}{_\beta} * g{^k}{_\gamma} = \begin{vmatrix} g{^0}_{\alpha} & g{^0}_{\beta} & g{^0}_{\gamma} \\ g{^1}_{\alpha} & g{^1}_{\beta} & g{^1}_{\gamma} \\ g{^2}_{\alpha} & g{^2}_{\beta} & g{^2}_{\gamma} \\ \end{vmatrix} = \begin{vmatrix} g{^s}_{\alpha} & g{^s}_{\beta} & g{^s}_{\gamma} \\ \end{vmatrix} = \;\;\;\sqrt{g} * e{_\alpha}{_\beta}{_\gamma} }[/math], and

the contravariant [math]\displaystyle{ \varepsilon }[/math]-tensor: [math]\displaystyle{ \varepsilon{^\kappa}{^\lambda}{^\mu} = \varepsilon{^i}{^j}{^k} * g{_i}^{\kappa}* g{_j}^{\lambda} * g{_k}^{\mu} = \begin{vmatrix} g{_0}^{\kappa} & g{_0}^{\lambda} & g{_0}^{\mu} \\ g{_1}^{\kappa} & g{_1}^{\lambda} & g{_1}^{\mu} \\ g{_2}^{\kappa} & g{_2}^{\lambda} & g{_2}^{\mu} \\ \end{vmatrix} = \begin{vmatrix} g{_s}^{\kappa} & g{_s}^{\lambda} & g{_s}^{\mu} \\ \end{vmatrix} = 1/\sqrt{g} * e{^\kappa}{^\lambda}{^\mu} }[/math].

Delta-Identities and Tensor Transformations

The Levi-Civita symbol resp. [math]\displaystyle{ \varepsilon }[/math]-tensor is related to the Kronecker delta by the Laplace expansion theorem. In three dimensions, this relationship is given by

[math]\displaystyle{ \begin{align} \delta{_{\alpha\beta\gamma}}^{\kappa\lambda\mu} = \varepsilon_{\alpha\beta\gamma} * \varepsilon^{\kappa\lambda\mu} = \begin{vmatrix} g{^s}{_\alpha} \\ g{^s}{_\beta} \\ g{^s}{_\gamma} \\ \end{vmatrix} * \begin{vmatrix} g{_s}{^\kappa} & g{_s}{^\lambda} & g{_s}{^\mu} \\ \end{vmatrix} = \begin{vmatrix} \delta{_\alpha}{^\kappa} & \delta{_\alpha}{^\lambda} & \delta{_\alpha}{^\mu} \\ \delta{_\beta} {^\kappa} & \delta{_\beta} {^\lambda} & \delta{_\beta} {^\mu} \\ \delta{_\gamma}{^\kappa} & \delta{_\gamma}{^\lambda} & \delta{_\gamma}{^\mu} \\ \end{vmatrix} = \delta{_\alpha}{^\kappa} * \begin{vmatrix} \delta{_\beta} {^\lambda} & \delta{_\beta}{^\mu} \\ \delta{_\gamma}{^\lambda} & \delta{_\gamma}{^\mu} \end{vmatrix} - \delta{_\alpha}{^\lambda} * \begin{vmatrix} \delta{_\beta} {^\kappa} & \delta{_\beta}{^\mu} \\ \delta{_\gamma}{^\kappa} & \delta{_\gamma}{^\mu} \end{vmatrix} + \delta{_\alpha}{^\mu} * \begin{vmatrix} \delta{_\beta} {^\kappa} & \delta{_\beta}{^\mu} \\ \delta{_\gamma}{^\kappa} & \delta{_\gamma}{^\mu} \end{vmatrix} \\ \end{align} }[/math]

[math]\displaystyle{ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; = \delta{_\alpha}{^\kappa} * ( \delta{_\beta}{^\lambda} * \delta{_\gamma}{^\mu} - \delta{_\beta}{^\mu} * \delta{_\gamma}{^\lambda}) - \delta{_\alpha}{^\lambda} * ( \delta{_\beta}{^\kappa} * \delta{_\gamma}{^\mu} - \delta{_\beta}{^\mu} * \delta{_\gamma}{^\kappa}) + \delta{_\alpha}{^\mu} * ( \delta{_\beta}{^\kappa} * \delta{_\gamma}{^\mu} - \delta{_\beta}{^\mu} * \delta{_\gamma}{^\kappa}) \textrm{} }[/math]

with the [math]\displaystyle{ \delta }[/math]-symbol of 6th order, representing the Lagrange identity.

By contraction of the last index, a special case follows from this identity as the Grassmann [math]\displaystyle{ \varepsilon }[/math]-identity with the Kronecker delta of 4th order

[math]\displaystyle{ \begin{align} \delta{_{\alpha\beta}}{^{\kappa\lambda}} = \varepsilon_{\alpha\beta\mu} * \varepsilon{^\kappa}^{\lambda\mu} = \begin{vmatrix} \delta{_\alpha}{^\kappa} & \delta{_\alpha}{^\lambda} \\ \delta{_\beta} {^\kappa} & \delta{_\beta} {^\lambda} \end{vmatrix} = \delta{_\alpha}{^\kappa} * \delta{_\beta}{^\lambda} - \delta{_\alpha}{^\lambda} * \delta{_\beta}{^\kappa}\\ \end{align} }[/math],

sometimes called as the "contracted epsilon identity".

This relationship is e.g. useful in deriving further identities for the cross product [math]\displaystyle{ \; \vec{c} = \vec{a}\times\vec{b} \; }[/math] with

[math]\displaystyle{ \; c^\mu = \varepsilon^{\kappa\lambda\mu} * a_\kappa * b_\lambda \;: }[/math]

like the Lagrange identity [math]\displaystyle{ \; (\vec{a} \times \vec{b}) \cdot (\vec{c} \times \vec{d}) = (\vec{a} \cdot \vec{c}) (\vec{b} \cdot \vec{d}) - (\vec{b} \cdot \vec{c}) (\vec{a} \cdot \vec{d}) \; }[/math] with

[math]\displaystyle{ \; \varepsilon_{\alpha\beta\mu} * a^\alpha * b^\beta * \varepsilon^{\kappa\lambda\mu} * c_\kappa * d_\lambda = a^\alpha * b^\beta * \delta{{_\alpha}{_\beta}{^\kappa}{^\lambda}} * c_\kappa * d_\lambda = a^\alpha * c_\alpha * b^\beta * d_\beta - b^\beta * c_\beta * a^\alpha * d_\alpha \; }[/math],

or the Grassmann identity [math]\displaystyle{ \; \vec{c}\times(\vec{a}\times\vec{b}) = (\vec{b} \cdot \vec{c} ) \vec{a} - (\vec{a} \cdot \vec{c}) \vec{b} \; }[/math] with

[math]\displaystyle{ \; c^\beta * \varepsilon_{\alpha\beta\mu} * \varepsilon^{\kappa\lambda\mu} * a_\kappa * b_\lambda = c^\beta * \delta{{_\alpha}{_\beta}{^\kappa}{^\lambda}} * a_\kappa * b_\lambda = b_\beta * c^\beta * a_\alpha - a_\beta * c^\beta * b_\alpha \; }[/math],

and the Jacobi identity [math]\displaystyle{ \; \vec{c}\times(\vec{a}\times\vec{b}) +\vec{a}\times (\vec{b}\times\vec{c}) +\vec{b}\times (\vec{c}\times\vec{a}) \equiv \vec{0} \; }[/math] with

[math]\displaystyle{ \; \varepsilon_{\alpha\beta\mu} * (c^\beta\ * a_\kappa * b_\lambda + a^\beta\ * b_\kappa * c_\lambda + b^\beta\ * c_\kappa * a_\lambda) * \varepsilon^{\kappa\lambda\mu} = (c^\beta\ * a_\kappa * b_\lambda + a^\beta\ * b_\kappa * c_\lambda + b^\beta\ * c_\kappa * a_\lambda) }[/math]

[math]\displaystyle{ \;\;\;\;\;\;\; =\; a_\alpha*(c^\beta\ * b_\beta - b^\beta\ * c_\beta) + b_\alpha * (a^\beta\ * c_\beta - c^\beta\ * a_\beta) + c_\alpha*(b^\beta\ * a_\beta - a^\beta\ * b_\beta) \equiv 0 \; }[/math].

Further contractions of the Grassmann [math]\displaystyle{ \varepsilon }[/math]-identität [math]\displaystyle{ \delta{_{\alpha\beta}}{^{\kappa\lambda}} }[/math] result in the [math]\displaystyle{ \delta }[/math]-identities

[math]\displaystyle{ in\; \mathbb R^3:\;\; \begin{align} \varepsilon_{\alpha\lambda\mu} * \varepsilon^{\beta\lambda\mu}\\ \end{align} \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; in\; \mathbb R^2:\;\; \delta{_\alpha}{^\beta} = \varepsilon_{\alpha\lambda} * \varepsilon^{\beta\lambda} }[/math]

[math]\displaystyle{ \begin{align} \;\;\;\;\;\;\;\;\; 2\;\delta{_\alpha}{^\alpha}= \varepsilon_{\alpha\beta\gamma} * \varepsilon^{\alpha\beta\gamma} & = 6 \end{align} \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; \delta{_\alpha}{^\alpha}= \varepsilon{_{\alpha\beta}} * \varepsilon^{\alpha\beta} = 2. }[/math]

By use of according tensor transformations, the upper [math]\displaystyle{ \delta }[/math]-identity may be used to determine the associated contravariant vector basis [math]\displaystyle{ \vec{g}{^\beta} }[/math] directly from a given covariant vector basis [math]\displaystyle{ \vec{g}{_\lambda} }[/math] - just by rotation and normalization of the basis vectors - without any matrix inversion!

For this purpose the [math]\displaystyle{ \delta }[/math]-identity of 2nd order is converted into the Meissner [math]\displaystyle{ \delta\varepsilon }[/math]-identities by several modifications as follow:

in [math]\displaystyle{ \mathbb R^3:\;\; }[/math] :[math]\displaystyle{ \delta{{_\alpha}{^\beta}} = 1/2 * \varepsilon{{_\alpha}{_j}{_k}} * g{^j}{_\lambda} * g{^k}{_\mu} * \varepsilon{^\beta}{^\lambda}{^\mu} \;\;\;\; }[/math] resp. [math]\displaystyle{ \;\;\;\; \delta{^{\alpha}{_\beta}} = 1/2 * \varepsilon{{^\alpha}{^\gamma}{^\kappa}} * g{_\gamma}{^\lambda} * g{_\kappa}{^\mu} * \varepsilon{{_\beta}{_\lambda}{_\mu}}, }[/math]

in [math]\displaystyle{ \mathbb R^2:\;\; }[/math] :[math]\displaystyle{ \delta{{_\alpha}{^\beta}} = \varepsilon{{_\alpha}{_j}} * g{{^j}{_\lambda}} * \varepsilon{{^\beta}{^\lambda}} \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; }[/math] resp. [math]\displaystyle{ \;\;\;\; \delta{{^\alpha}{_\beta}} = \varepsilon{{^\alpha}{^\gamma}} * g{_\gamma}{^\lambda} * \varepsilon{{_\beta}{_\lambda}} }[/math].

Thus, approriate tensor transformations into the skewed reference system lead then, by some realignments and renamings of indices, to the contravariant Meissner tensorbasis transformations for

the three-dimensional contravariant basis vectors [math]\displaystyle{ g{_i}{^\beta} = 1/2 * e{_i}{_j}{_k} * g{^j}{_\lambda} * g{^k}{_\mu} * \varepsilon{^\beta}{^\lambda}{^\mu} }[/math]

and the contravariant metric tensor [math]\displaystyle{ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; g{^\alpha}{^\beta} = 1/2 * \varepsilon{^\alpha}{^\gamma}{^\kappa} * g{_\gamma}{_\lambda} * g{_\kappa}{_\mu} * \varepsilon{^\beta}{^\lambda}{^\mu} }[/math],

respectively to

the two-dimensional contravariant basis vectors [math]\displaystyle{ \;\; g{{_i}{^\beta}} = e{_i}{_j} * g{^j}{_\lambda} * \varepsilon{{^\beta}{^\lambda}} }[/math]

and the contravariant metric tensor [math]\displaystyle{ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; g{{^\alpha}{^\beta}} = \varepsilon{{^\alpha}{^\gamma}} * g{{_\gamma}{_\lambda}} * \varepsilon{{^\beta}{^\lambda}} }[/math].

For the numerical evaluation of tensor equations, specific tensor/matrix algorithms were developed based on object-oriented principles with documentation in details. The associated C++ software with structured matrix and tensor classes implements the overloading of arithmetic and functional standard operators for the specific use of multi-dimensional tensor/matrix objects, incl. the generalized matrix multiplication [math]\displaystyle{ C_{IKLM} = A_{ISK} * B_{LSM} }[/math]. The program sources are available to the academic public and for private use over the WWW.

n dimensions

Index and symbol values

In n dimensions, when all i₁, ...,i_n, j₁, ..., j_n take values 1, 2, ..., n:

[math]\displaystyle{ \varepsilon_{i_1 \dots i_n} \varepsilon^{j_1 \dots j_n} = \delta^{j_1 \dots j_n}_{i_1 \dots i_n} }[/math]

( 7 )

[math]\displaystyle{ \varepsilon_{i_1 \dots i_k~i_{k+1} \dots i_n} \varepsilon^{i_1 \dots i_k~j_{k+1} \dots j_n} = \delta_{ i_1 \ldots i_k~i_{k+1} \ldots i_n}^{i_1 \dots i_k~j_{k+1}\ldots j_n} = k!~\delta^{j_{k+1} \dots j_n}_{i_{k+1} \dots i_n} }[/math]

( 8 )

[math]\displaystyle{ \varepsilon_{i_1 \dots i_n}\varepsilon^{i_1 \dots i_n} = n! }[/math]

( 9 )

where the exclamation mark (!) denotes the factorial, and δα...β... is the generalized Kronecker delta. For any n, the property

[math]\displaystyle{ \sum_{i, j, k, \dots = 1}^n \varepsilon_{ijk\dots}\varepsilon_{ijk\dots} = n! }[/math]

follows from the facts that

• every permutation is either even or odd,
• (+1)² = (−1)² = 1, and
• the number of permutations of any n-element set number is exactly n!.

The particular case of (8) with [math]\displaystyle{ k=n-2 }[/math] is [math]\displaystyle{ \varepsilon_{i_1\dots i_{n-2}jk}\varepsilon^{i_1\dots i_{n-2}lm} = (n-2)!(\delta_j^l\delta_k^m - \delta_j^m\delta_l^k)\,. }[/math]

Product

In general, for n dimensions, one can write the product of two Levi-Civita symbols as: [math]\displaystyle{ \varepsilon_{i_1 i_2 \dots i_n} \varepsilon_{j_1 j_2 \dots j_n} = \begin{vmatrix} \delta_{i_1 j_1} & \delta_{i_1 j_2} & \dots & \delta_{i_1 j_n} \\ \delta_{i_2 j_1} & \delta_{i_2 j_2} & \dots & \delta_{i_2 j_n} \\ \vdots & \vdots & \ddots & \vdots \\ \delta_{i_n j_1} & \delta_{i_n j_2} & \dots & \delta_{i_n j_n} \\ \end{vmatrix}. }[/math]Proof: Both sides change signs upon switching two indices, so without loss of generality assume [math]\displaystyle{ i_1 \leq \cdots\leq i_n, j_1 \leq \cdots \leq j_n }[/math]. If some [math]\displaystyle{ i_c = i_{c+1} }[/math] then left side is zero, and right side is also zero since two of its rows are equal. Similarly for [math]\displaystyle{ j_c = j_{c+1} }[/math]. Finally, if [math]\displaystyle{ i_1 \lt \cdots \lt i_n, j_1 \lt \cdots \lt j_n }[/math], then both sides are 1.

Proofs

For (1), both sides are antisymmetric with respect of ij and mn. We therefore only need to consider the case i ≠ j and m ≠ n. By substitution, we see that the equation holds for ε₁₂ε¹², that is, for i = m = 1 and j = n = 2. (Both sides are then one). Since the equation is antisymmetric in ij and mn, any set of values for these can be reduced to the above case (which holds). The equation thus holds for all values of ij and mn.

Using (1), we have for (2)

[math]\displaystyle{ \varepsilon_{ij} \varepsilon^{in} = \delta_i{}^i \delta_j{}^n - \delta_i{}^n \delta_j{}^i = 2 \delta_j{}^n - \delta_j{}^n = \delta_j{}^n \,. }[/math]

Here we used the Einstein summation convention with i going from 1 to 2. Next, (3) follows similarly from (2).

To establish (5), notice that both sides vanish when i ≠ j. Indeed, if i ≠ j, then one can not choose m and n such that both permutation symbols on the left are nonzero. Then, with i = j fixed, there are only two ways to choose m and n from the remaining two indices. For any such indices, we have

[math]\displaystyle{ \varepsilon_{jmn} \varepsilon^{imn} = \left(\varepsilon^{imn}\right)^2 = 1 }[/math]

(no summation), and the result follows.

Then (6) follows since 3! = 6 and for any distinct indices i, j, k taking values 1, 2, 3, we have

[math]\displaystyle{ \varepsilon_{ijk} \varepsilon^{ijk} = 1 }[/math] (no summation, distinct i, j, k)

Applications and examples

Determinants

In linear algebra, the determinant of a 3 × 3 square matrix A = [a_ij] can be written

[math]\displaystyle{ \det(\mathbf{A}) = \sum_{i=1}^3 \sum_{j=1}^3 \sum_{k=1}^3 \varepsilon_{ijk} a_{1i} a_{2j} a_{3k} }[/math]

Similarly the determinant of an n × n matrix A = [a_ij] can be written as

[math]\displaystyle{ \det(\mathbf{A}) = \varepsilon_{i_1\dots i_n} a_{1i_1} \dots a_{ni_n}, }[/math]

where each i_r should be summed over 1, ..., n, or equivalently:

[math]\displaystyle{ \det(\mathbf{A}) = \frac{1}{n!} \varepsilon_{i_1\dots i_n} \varepsilon_{j_1\dots j_n} a_{i_1 j_1} \dots a_{i_n j_n}, }[/math]

where now each i_r and each j_r should be summed over 1, ..., n. More generally, we have the identity

[math]\displaystyle{ \sum_{i_1,i_2,\dots}\varepsilon_{i_1\dots i_n} a_{i_1 \, j_1} \dots a_{i_n \, j_n} = \det(\mathbf{A}) \varepsilon_{j_1\dots j_n} }[/math]

Vector cross product

Cross product (two vectors)

Let [math]\displaystyle{ (\mathbf{e_1}, \mathbf{e_2}, \mathbf{e_3}) }[/math] a positively oriented orthonormal basis of a vector space. If (a¹, a², a³) and (b¹, b², b³) are the coordinates of the vectors a and b in this basis, then their cross product can be written as a determinant:

[math]\displaystyle{ \mathbf{a \times b} = \begin{vmatrix} \mathbf{e_1} & \mathbf{e_2} & \mathbf{e_3} \\ a^1 & a^2 & a^3 \\ b^1 & b^2 & b^3 \\ \end{vmatrix} = \sum_{i=1}^3 \sum_{j=1}^3 \sum_{k=1}^3 \varepsilon_{ijk} \mathbf{e}_i a^j b^k }[/math]

hence also using the Levi-Civita symbol, and more simply:

[math]\displaystyle{ (\mathbf{a \times b})^i = \sum_{j=1}^3 \sum_{k=1}^3 \varepsilon_{ijk} a^j b^k. }[/math]

In Einstein notation, the summation symbols may be omitted, and the ith component of their cross product equals

[math]\displaystyle{ (\mathbf{a\times b})^i = \varepsilon_{ijk} a^j b^k. }[/math]

The first component is

[math]\displaystyle{ (\mathbf{a\times b})^1 = a^2 b^3-a^3 b^2\,, }[/math]

then by cyclic permutations of 1, 2, 3 the others can be derived immediately, without explicitly calculating them from the above formulae:

[math]\displaystyle{ \begin{align} (\mathbf{a\times b})^2 &= a^3 b^1-a^1 b^3\,, \\ (\mathbf{a\times b})^3 &= a^1 b^2-a^2 b^1\,. \end{align} }[/math]

Triple scalar product (three vectors)

From the above expression for the cross product, we have:

[math]\displaystyle{ \mathbf{a\times b} = -\mathbf{b\times a} }[/math].

If c = (c¹, c², c³) is a third vector, then the triple scalar product equals

[math]\displaystyle{ \mathbf{a}\cdot(\mathbf{b\times c}) = \varepsilon_{ijk} a^i b^j c^k. }[/math]

From this expression, it can be seen that the triple scalar product is antisymmetric when exchanging any pair of arguments. For example,

[math]\displaystyle{ \mathbf{a}\cdot(\mathbf{b\times c})= -\mathbf{b}\cdot(\mathbf{a\times c}) }[/math].

Curl (one vector field)

If F = (F¹, F², F³) is a vector field defined on some open set of [math]\displaystyle{ \mathbb{R}^3 }[/math] as a function of position x = (x¹, x², x³) (using Cartesian coordinates). Then the ith component of the curl of F equals

[math]\displaystyle{ (\nabla \times \mathbf{F})^i(\mathbf{x}) = \varepsilon_{ijk}\frac{\partial}{\partial x^j} F^k(\mathbf{x}), }[/math]

which follows from the cross product expression above, substituting components of the gradient vector operator (nabla).

Tensor density

In any arbitrary curvilinear coordinate system and even in the absence of a metric on the manifold, the Levi-Civita symbol as defined above may be considered to be a tensor density field in two different ways. It may be regarded as a contravariant tensor density of weight +1 or as a covariant tensor density of weight −1. In n dimensions using the generalized Kronecker delta,

[math]\displaystyle{ \begin{align} \varepsilon^{\mu_1 \dots \mu_n} &= \delta^{\mu_1 \dots \mu_n}_{\,1 \,\dots \,n} \, \\ \varepsilon_{\nu_1 \dots \nu_n} &= \delta^{\,1 \,\dots \,n}_{\nu_1 \dots \nu_n} \,. \end{align} }[/math]

Notice that these are numerically identical. In particular, the sign is the same.

Levi-Civita tensors

On a pseudo-Riemannian manifold, one may define a coordinate-invariant covariant tensor field whose coordinate representation agrees with the Levi-Civita symbol wherever the coordinate system is such that the basis of the tangent space is orthonormal with respect to the metric and matches a selected orientation. This tensor should not be confused with the tensor density field mentioned above. The presentation in this section closely follows Carroll 2004.

The covariant Levi-Civita tensor (also known as the Riemannian volume form) in any coordinate system that matches the selected orientation is

[math]\displaystyle{ E_{a_1\dots a_n} = \sqrt{\left|\det [g_{ab}]\right|}\, \varepsilon_{a_1\dots a_n} \,, }[/math]

where g_ab is the representation of the metric in that coordinate system. We can similarly consider a contravariant Levi-Civita tensor by raising the indices with the metric as usual,

[math]\displaystyle{ E^{a_1\dots a_n} = E_{b_1\dots b_n} \prod_{i=1}^n g^{a_i b_i} \,, }[/math]

but notice that if the metric signature contains an odd number of negative eigenvalues q, then the sign of the components of this tensor differ from the standard Levi-Civita symbol:

[math]\displaystyle{ E^{a_1\dots a_n} = \frac{ \sgn \left( \det [g_{ab}] \right) }{ \sqrt{ \left| \det [g_{ab}] \right| }} \, \varepsilon^{a_1\dots a_n} , }[/math]

where sgn(det[g_ab]) = (−1)^q, and [math]\displaystyle{ \varepsilon^{a_1\dots a_n} }[/math] is the usual Levi-Civita symbol discussed in the rest of this article. More explicitly, when the tensor and basis orientation are chosen such that [math]\displaystyle{ E_{01\dots n} = +\sqrt{\left|\det [g_{ab}]\right|} }[/math], we have that [math]\displaystyle{ E^{01\dots n} = \frac{\sgn(\det [g_{ab}])}{\sqrt{\left|\det [g_{ab}]\right|}} }[/math].

From this we can infer the identity,

[math]\displaystyle{ E^{\mu_1\dots\mu_p\alpha_1\dots\alpha_{n-p}}E_{\mu_1\dots\mu_p\beta_1\dots\beta_{n-p}} = (-1)^q p!\delta^{\alpha_1\dots\alpha_{n-p}}_{\beta_1\dots\beta_{n-p}} \,, }[/math]

where

[math]\displaystyle{ \delta^{\alpha_1 \dots \alpha_{n-p}}_{\beta_1 \dots \beta_{n-p}} = (n-p)! \delta^{\lbrack \alpha_1}_{\beta_1} \dots \delta^{\alpha_{n-p} \rbrack}_{\beta_{n-p}} }[/math]

is the generalized Kronecker delta.

Example: Minkowski space

In Minkowski space (the four-dimensional spacetime of special relativity), the covariant Levi-Civita tensor is

[math]\displaystyle{ E_{\alpha \beta \gamma \delta} = \pm \sqrt{ \left| \det [g_{\mu \nu}] \right| } \, \varepsilon_{\alpha \beta \gamma \delta} \,, }[/math]

where the sign depends on the orientation of the basis. The contravariant Levi-Civita tensor is

[math]\displaystyle{ E^{\alpha \beta \gamma \delta} = g^{\alpha \zeta} g^{\beta \eta} g^{\gamma \theta} g^{\delta \iota} E_{\zeta \eta \theta \iota} \,. }[/math]

The following are examples of the general identity above specialized to Minkowski space (with the negative sign arising from the odd number of negatives in the signature of the metric tensor in either sign convention):

[math]\displaystyle{ \begin{align} E_{\alpha \beta \gamma \delta} E_{\rho \sigma \mu \nu} & = -g_{\alpha \zeta} g_{\beta \eta} g_{\gamma \theta} g_{\delta \iota} \delta^{\zeta \eta \theta \iota}_{\rho \sigma \mu \nu} \\ E^{\alpha \beta \gamma \delta} E^{\rho \sigma \mu \nu} & = -g^{\alpha \zeta} g^{\beta \eta} g^{\gamma \theta} g^{\delta \iota} \delta^{\rho \sigma \mu \nu}_{\zeta \eta \theta \iota} \\ E^{\alpha \beta \gamma \delta} E_{\alpha \beta \gamma \delta} & = - 24 \\ E^{\alpha \beta \gamma \delta} E_{\rho \beta \gamma \delta} & = - 6 \delta^{\alpha}_{\rho} \\ E^{\alpha \beta \gamma \delta} E_{\rho \sigma \gamma \delta} & = - 2 \delta^{\alpha \beta}_{\rho \sigma} \\ E^{\alpha \beta \gamma \delta} E_{\rho \sigma \theta \delta} & = - \delta^{\alpha \beta \gamma}_{\rho \sigma \theta} \,. \end{align} }[/math]

See also

• List of permutation topics
• Symmetric tensor

References

• Misner, C.; Thorne, K. S.; Wheeler, J. A. (1973). Gravitation. W. H. Freeman & Co. pp. 85–86, §3.5. ISBN 0-7167-0344-0. 
• Neuenschwander, D. E. (2015). Tensor Calculus for Physics. Johns Hopkins University Press. pp. 11, 29, 95. ISBN 978-1-4214-1565-9. 
• Carroll, Sean M. (2004), Spacetime and Geometry, Addison-Wesley, ISBN 0-8053-8732-3, https://www.preposterousuniverse.com/spacetimeandgeometry/ 

External links

• Weisstein, Eric W.. "Permutation Tensor". http://mathworld.wolfram.com/PermutationTensor.html.