14Manifold_Differential_Forms.ipynb

Kernel: SageMath 9.6

14. Differential $k$-forms

This notebook is part of the Introduction to manifolds in SageMath by Andrzej Chrzeszczyk (Jan Kochanowski University of Kielce, Poland).

version()

'SageMath version 9.6, Release Date: 2022-05-15'

Warning: in this notebook there are many repetitions with respect to notebook 10. Although mathematically alternating forms on modules from notebook 10 are generalizations of differential forms, the SageMath code in the present notebook differs significantly from that in notebook 10.

Let $M$ be a smooth manifold.

A differential $k$-form or simply $k$-form on a manifold $M$ is a smooth map $ω$ that, to each point $p ∈ M$, assigns $k$-linear antisymmetric form $ω(p)=ω_p$ on the tangent space $T_pM$.

Reminder. $k$-linear antisymmetric forms on arbitrary modules were presented in notebook 10, and smoothness of covariant tensor fields in notebook 13.

One can use an equivalent definition:

A differential form of degree $k$, or differential $k$-form, $ω$ on a manifold $M$ is a smooth tensor field of type $(0,k)$ on $M$, with the antisymmetry property:

ParseError: KaTeX parse error: Undefined control sequence: \label at position 122: …. . . , X_k ), \̲l̲a̲b̲e̲l̲{}\tag{14.1} \e…

for $1\leq i < j \leq k$, $X_1 , . . . , X_k ∈ \mathfrak{X}(M)$.

Note that from (14.1) it follows that if $X_i=X_j$ then $ω(X_1 , . . . , X_i , . . . , X_j , . . . , X_k ) = 0$.
The relation (14.1) can be reformulated as

$$\omega(X_{\sigma(1)},\ldots,X_{\sigma(k)})=\mathrm{sign}(\sigma) \omega(X_1,\ldots,X_k),$$

for $\sigma\in S_k$ (permutations and their signs were defined in notebook 10).

A 0-form is a smooth real-valued function on $M$.
The set of the k-forms on a manifold $M$, is denoted by $\Omega^k(M)$. It is a submodule of the module of covariant tensor fields $T^{(0,k)} M$ on the manifold $M$.

Example 14.1

Let us define a differential 2-form on a 3-dimensional manifold.

M = Manifold(3, 'M')             # manifold M
X.<x,y,z> = M.chart()
a = M.diff_form(2, name='a') ;   # define 2-form
print(a)                         # show information on a

2-form a on the 3-dimensional differentiable manifold M

# parent?

# mathematical object of which "a" is an element
show(a.parent())

$\displaystyle \Omega^{2}\left(M\right)$

The obtained module can be defined independently:

Om2 = M.diff_form_module(2)  # define module of 2-forms
show(Om2)

$\displaystyle \Omega^{2}\left(M\right)$

Antisymmetrization operation Alt

If $t$ is a covariant tensor field from $T^{(0,k)}M$, then we can define a differential $k$-form called Alt($t$) in the following way

ParseError: KaTeX parse error: Undefined control sequence: \label at position 144: …_{\sigma(k)}), \̲l̲a̲b̲e̲l̲{}\tag{14.2} \e…

for $X_i\in\mathfrak{X}(M)$.

Alt is a linear map on the module $T^{(0,k)}M$ (over $C^\infty(M)$) into $\Omega^k(M)$.

To prove that for $t\in T^{(0,k)}M$ we have $\text{Alt}(t)\in \Omega^k(M)$ note that for arbitrary permutation $\tau\in S_k$

$$\text{Alt}(t)(v_{\tau(1)},\ldots,v_{\tau(k)}) =\frac{1}{k!}\sum_{\sigma\in S_k}\text{sign}\sigma\, t(v_{\sigma(\tau(1))},\ldots,v_{\sigma(\tau(k))})\\ =\frac{1}{k!}\sum_{\sigma\in S_k}\text{sign}\sigma\, \text{sign}\sigma\,t(v_{\tau(1)},\ldots,v_{\tau(k)})=\text{sign}\tau\,t(v_{1},\ldots,v_{k}).$$

We can check also that if $t\in \Omega^k(M)$, then Alt($t$)=$t$.
In fact, assume that $t(v_{\sigma(1)},\ldots,v_{\sigma(k)})=\text{sign}\sigma\, t(v_1,\ldots,v_k)$. Then $$\text{Alt}(t)(v_1,\ldots,v_k)= \frac{1}{k!}\sum_{\sigma\in S_k}(\text{sign} \sigma)t(v_{\sigma(1)},\ldots,v_{\sigma(k)})\\ =\frac{1}{k!}\sum_{\sigma\in S_k}(\text{sign} \sigma)(\text{sign} \sigma)\, t(v_{1},\ldots,v_{k})=t(v_{1},\ldots,v_{k}).$$

The last two observations lead to the relation $$\text{Alt}(\text{Alt}(t))=\text{Alt}(t), \quad \text{or equivalently}\quad \text{Alt}^2=\text{Alt}.$$

Wedge product

If $ω$ is a $k$-differential form and $η$ is an $l$-differential form on $M$, the exterior, or wedge, product $ω ∧ η\in\Omega^{(k+l)}(M)$, is defined by ParseError: KaTeX parse error: Undefined control sequence: \label at position 87: …a\otimes\eta). \̲l̲a̲b̲e̲l̲{}\tag{14.3} \e…

Remark. Some authors use different coefficients in this formula. Our version is in accordance with that in SageMath Manifolds.

Example 14.2

If $α$ and $β$ are 1-forms, applying (14.3), (14.2) and using the fact that there are only two permutations of two elements with opposite signs we have $$(α ∧ β)(X_1 , X_2 ) = \frac{(1+1)!}{1!1!}\mathrm{Alt} (α ⊗ β)(X_1 , X_2 )\\ = \frac{2!}{1!1!}\frac{1}{2!} [ (α ⊗ β)(X_1 , X_2 ) − (α ⊗ β)(X_2 , X_1 )]\\ = α(X_1 )β(X_2 ) − α(X_2 )β(X_1 ) = (α ⊗ β − β ⊗ α)(X_1 , X_2 ),$$

for $X_1 , X_2 ∈ \mathfrak{X}(M).$ Consequently $$α ∧ β = (α ⊗ β − β ⊗ α) = −β ∧ α,\quad \mbox{for }\quad α, β ∈ \Omega^1 (M).$$

Basic algebraic properties of wedge product

From the multilinearity of the tensor product and linearity of antisymmetrization it follows.

For $\alpha,\alpha_1,\alpha_2\in \Omega^k(M),\beta,\beta_1,\beta_2\in \Omega^m(M)$ and $a\in C^\infty(M)$

$$\begin{equation} \begin{matrix} (\alpha_1+\alpha_2)\wedge\beta=\alpha_1\wedge\beta+\alpha_2\wedge\beta,\\ \alpha\wedge(\beta_1+\beta_2)=\alpha\wedge\beta_1+\alpha\wedge\beta_2,\\ (a\alpha)\wedge\beta=a(\alpha\wedge\beta)=\alpha\wedge(a\beta). \end{matrix} \tag{14.4} \end{equation}$$

Associativity of the wedge product

To check the associativity of the wedge product we need some properties of the antisymmetrization operation.

First observe that if $t\in T^{(0,k)}M$ and $s\in T^{(0,m)}M$, then

$$\begin{equation} \text{Alt}(t)=0\quad \Longrightarrow\quad \text{Alt}(t\otimes s)=\text{Alt}(s\otimes t)=0. \tag{Alt1} \end{equation}$$

To prove this, let us note that

$$(k + m)!\, \text{Alt}(t ⊗ s)(v_1 , . . . , v_{k+m} ) = \sum_{\sigma\in S_{k+m}} (\text{sign}\, σ) t (v_{σ(1)} , . . . , v_{σ(k)} )s(v_{σ(k+1)} , . . . , v_{σ(k+m)} ).$$

First define the subgroup of $S_{k+m}$,

$$G=\{\sigma\in S_{k+m}: (\sigma(k+1),\ldots,\sigma(k+m))=(k+1,\ldots,k+m)\},$$

and compute the sum from the right hand side of the previous equality restricted to this subgroup.

$$\sum_{\sigma\in G}(\text{sign}\, σ)\,t (v_{σ(1)} , . . . , v_{σ(k)} )s(v_{σ(k+1)} , . . . , v_{σ(k+m)} )\\ =\Big[\sum_{\sigma\in G}(\text{sign}\, σ)\,t (v_{σ(1)} , . . . , v_{σ(k)} )\Big] s(v_{k+1} , . . . , v_{k+m} )\\= k!(\text{Alt}(t))(v_{1} , . . . , v_{k} ) s(v_{k+1} , . . . , v_{k+m} )= k!(\text{Alt}(t)\otimes s)(v_1,\ldots,v_{k+m}),$$

so the entire sum vanishes if the sum in brackets vanishes, i.e. Alt($t$) vanishes.

Now for fixed $\tilde\sigma\in S_{k+m}$ let us sum over the coset $G\tilde\sigma=\{\sigma\tilde\sigma: \sigma \in G$}.

$$\sum_{\sigma'\in G\tilde\sigma}(\text{sign}\, σ')(t ⊗ s)(v_{σ'(1)} , . . . , v_{σ'(k+m)} )\\ =\sum_{\sigma\in G}(\text{sign}\, σ\tilde\sigma)(t ⊗ s)(v_{\sigma(\tilde\sigma(1))} , . . . , v_{\sigma(\tilde \sigma(k+m))} )\\ =\text{sign}\tilde\sigma \Big[\sum_{\sigma\in G}\text{sign}\sigma \, t(v_{\sigma(\tilde\sigma(1))} , . . . , v_{\sigma(\tilde \sigma(k))} )\Big] s(v_{\sigma(\tilde\sigma(k+1))} , . . . , v_{\sigma(\tilde \sigma(k+m))} )\\ =\text{sign}\tilde\sigma\, k!\,(\text{Alt}(t)) (v_{\tilde\sigma(1)} , . . . , v_{\tilde \sigma(k)} ) s(v_{\tilde\sigma(k+1)} , . . . , v_{\tilde \sigma(k+m)} )\\ \text{sign}\tilde\sigma\, k!\,\big[\text{Alt}(t)\otimes s\big] (v_{\tilde\sigma(1)} , . . . , v_{\tilde \sigma(k+m)} ) .$$

Thus if $\text{Alt}(t)$ vanishes then the sums over all cosets $G\tilde\sigma$ vanish. Since the group $S_{k+m}$ is a sum of cosets $G\tilde\sigma$, the implication (Alt1) holds true.

The second ingredient of our associativity proof is the following

$$\begin{equation} \text{Alt}(\text{Alt}(t ⊗ s) ⊗ r) = \text{Alt}(t ⊗ s ⊗ r) = \text{Alt}(t ⊗ \text{Alt}(s ⊗ r)). \tag{Alt2} \end{equation}$$

In fact, the linearity of Alt and the relation $\text{Alt}^2=\text{Alt}$ implies

$$\text{Alt}(\text{Alt}(s ⊗ r) − s ⊗ r)=\text{Alt}(\text{Alt}(s ⊗ r))-\text{Alt}(s\otimes r)=0.$$

Applying the implication (Alt1) to $\ \text{Alt}(s ⊗ r) − s ⊗ r\$ we obtain

$$\text{Alt}(t\otimes\text{Alt}(s\otimes r))-\text{Alt}(t\otimes s\otimes r)= \text{Alt}(t\otimes[\text{Alt}(s\otimes r)-s\otimes r]).$$

Since we have checked that the value of Alt on the expression in brackets vanishes, then the implication (Alt1) gives us

$$\text{Alt}(t\otimes\text{Alt}(s\otimes r)-\text{Alt}(t\otimes s\otimes r)) =\text{Alt}(t\otimes[\text{Alt}(s\otimes r)-s\otimes r])=0,$$

i.e., the second equality from (Alt2) holds true

$$\text{Alt}(t ⊗ \text{Alt}(s ⊗ r)) = \text{Alt}(t ⊗ s ⊗ r).$$

The first equality in (Alt2) can be proved analogously.

From (Alt2) it follows for $\alpha\in\Omega^k(M),\ \beta\in \Omega^l(M), \ \gamma\in \Omega^m(M)$

$$(\alpha\wedge \beta)\wedge \gamma= \frac{(k+m+l)!}{(k+m)!l!}\text{Alt}((\alpha\wedge\beta)\otimes\gamma)\\ =\frac{(k+m+l)!}{(k+m)!l!}\text{Alt}(\frac{(k+m)!}{k!m!}(\text{Alt}(\alpha\otimes\beta)\otimes\gamma)\\= \frac{(k+m+l)!}{k!m!l!}\text{Alt}(\alpha\otimes\beta\otimes\gamma).$$

Analogously for $α ∧ (β ∧ γ )$. We have proved associativity of the wedge product:

$$(α ∧ β) ∧ γ=α ∧ (β ∧ γ ) = \frac{(k+l+m)!}{k!l!m!} \mathrm{Alt} (α ⊗ β ⊗ γ ),$$

for $\alpha\in\Omega^k(M),\ \beta\in \Omega^l(M), \ \gamma\in \Omega^m(M)$, and more generally

ParseError: KaTeX parse error: Undefined control sequence: \label at position 136: …s \alpha_r ), \̲l̲a̲b̲e̲l̲{}\tag{14.6} \e…

for $\alpha_i\in\Omega^{k_i}(M).$

For example if $x^1,\ldots,x^n$ are local coordinates on $M$, then since $dx^i$ are 1-forms, then $dx^{i_1}\wedge\ldots\wedge dx^{i_k}$ are $k$-forms on $M$.

Differential forms in local coordinates

Let $(x^1 , . . . , x^n )$ be a local coordinate system on $M$. According to (13.7) a $k$-differential form possesses the local representation

ParseError: KaTeX parse error: Undefined control sequence: \label at position 66: … · ⊗ dx^{i_k}, \̲l̲a̲b̲e̲l̲{}\tag{14.7} \e…

where

ParseError: KaTeX parse error: Undefined control sequence: \label at position 122: …x^{i_k}}\Big). \̲l̲a̲b̲e̲l̲{}\tag{14.8} \e…

Since $\omega$ is antisymmetric, $ω_{i_1 ...i_k}$ is also antisymmetric.
Since Alt is linear and Alt($\omega$)=$\omega\$ we have

$$\begin{matrix} \omega=\mathrm{Alt}(\omega)=\sum_{i_1,\ldots,i_k=1}^n ω_{i_1 ...i_k}\mathrm{Alt}(dx^{i_1}\otimes\ldots\otimes dx^{i_k}).\\ \end{matrix}$$

From (14.6) it follows that (note that $dx^{i_j}$ are 1-forms) $$\mathrm{Alt}(dx^{i_1}\otimes\ldots\otimes dx^{i_k})=\frac{1}{k!}dx^{i_1}\wedge\ldots\wedge dx^{i_k}.$$ Thus we have ParseError: KaTeX parse error: Undefined control sequence: \label at position 228: …. \end{matrix} \̲l̲a̲b̲e̲l̲{}\tag{14.9} \e…

Note that every element of the sum with increasing indices has $k!$ counterparts in the sum with arbitrarily ordered indices. For arbitrary permutation $\sigma\in S_k$ (no summation here)

$$\omega_{\sigma(i_1)\ldots\sigma(i_k)}dx^{\sigma(i_1)}\wedge\ldots\wedge dx^{\sigma(i_k)}=\mathrm{sign}(\sigma)^2ω_{i_1 ...i_k}dx^{i_1}\wedge\ldots\wedge dx^{i_k}\\=ω_{i_1 ...i_k}dx^{i_1}\wedge\ldots\wedge dx^{i_k}.$$

Elements of the set

$$\{dx^{i_1}\wedge\dots\wedge dx^{i_k}: 1\leq i_1<\ldots<i_k\leq n\}$$

are linearly independent, since if

$$0=\omega=\sum_{i_1<...<i_k}\omega_{i_1\ldots i_k}dx^{i_1}\wedge\dots\wedge dx^{i_k},$$

then taking the basis $\{\frac{\partial}{\partial x^{1}},\ldots,\frac{\partial}{\partial x^{n}}\}$ and an increasing index sequence $(j_1,\ldots,j_k)$, we obtain

$$0=\omega(\frac{\partial}{\partial x^{j_1}},\ldots,\frac{\partial}{\partial x^{j_k}})=k!\sum_{i_1<...<i_k}\omega_{i_1\ldots i_k} \text{Alt}(\{dx^{i_1}\otimes\dots\otimes dx^{i_k})\big( \frac{\partial}{\partial x^{j_1}},\ldots,\frac{\partial}{\partial x^{j_k}}\big)\\ =\sum_{i_1<...<i_k}\omega_{i_1\ldots i_k} \sum_{\sigma\in S_k}(\text{sign}\, \sigma) dx^{i_1}(\frac{\partial}{\partial x^{j_{\sigma(1)}}})\ldots dx^{i_k}(\frac{\partial}{\partial x^{j_{\sigma(k)}}}).$$

Since $(i_1 , . . . , i_k )$ and $( j_1 , . . . , j_k )$ are both increasing, there is only one choice of $i_1,..., i_k\$ which gives a nonzero value, namely $\ i_1 = j_1,\ldots,i_k =j_k.$ Since

$$dx^{j_1}(\frac{\partial}{\partial x^{j_{1}}})\ldots dx^{j_k}(\frac{\partial}{\partial x^{j_{k}}})=1$$

we have $$0=\omega(\frac{\partial}{\partial x^{j_1}},\ldots,\frac{\partial}{\partial x^{j_k}})=\omega_{j_1\ldots j_k}.$$

Anticommutativity

If $\omega ∈ \Omega^k(M)$ and $\eta ∈ \Omega^l(M)$, then

$$\omega ∧ \eta = (−1)^{kl} \eta ∧ \omega.$$

This is consequence of (we use the unordered version of (14.9))

$$(k+l)!ω ∧ η = ω_{i_1 ...i_k} η_{ j_1 ... j_l} dx^{i_1} ∧ · · · ∧ dx^{i_k} ∧ dx^{j_1} ∧ · · · ∧ dx^{j_l}\\ = (−1)^{kl} ω_{i_1 ...i_k} η_{j_1 ... j_l} dx^{j_1} ∧ · · · ∧ dx^{j_l} ∧ dx^{i_1} ∧ · · · ∧ dx^{i_k} = (k+l)!(−1)^{kl} η ∧ ω.$$

Abbreviated notations for $k$-forms

For arbitrary differential forms in $\Omega^k(M)$ sometimes we would not want to actually write out all of the elements of the local frames of $\Omega^k(M)$, so instead we write ParseError: KaTeX parse error: Undefined control sequence: \label at position 43: …m_I a_I dx^I . \̲l̲a̲b̲e̲l̲{}\tag{14.10} \… Here the $I$ stands for the sequence of $k$ increasing indices $i_1 i_2 \ldots i_k$: $1 ≤ i_1 < i_2 < \ldots < i_k ≤ n$. That is, we sum over $I ∈ J_{k,n} = \{(i_1 i_2 \ldots i_k ) : 1 ≤ i_1 < i_2 < \ldots < i_k ≤ n\}.$

For example, for $k = 3$ and $n = 4$ we have $I\in\{123,124,134,234\}.$

If $I$ and $J$ are disjoint, then we have $dx^I ∧ dx^J = ±dx^K$ where $K = I ∪ J$, but is reordered to be in increasing order. Elements with repeated indices are dropped. Using this notation we can compute the wedge product as follows ParseError: KaTeX parse error: Undefined control sequence: \label at position 152: …Ib_J\Big)dx^K. \̲l̲a̲b̲e̲l̲{}\tag{14.11} \…

Example 14.3

Let us define a differential 2-form on a 3-dimensional manifold.

%display latex
N = 3                 # dimension of manifold M 
M = Manifold(N, 'M')  # manifold M
X = M.chart(' '.join(['x'+str(i)+':x^{'+str(i)+'}' for i in range(N)]))  # chart on M
x0, x1, x2 = X[:]     # coordinates x^0, x^1, x^2 of chart X as the Python variables x0, x1, x2
a = M.diff_form(2, name='a')  # differential 2-form on M
a01 = M.scalar_field(function('a01')(x0,x1,x2), name='a01')  # component a_{01}
a02 = M.scalar_field(function('a02')(x0,x1,x2), name='a02')  # component a_{02}
a12 = M.scalar_field(function('a12')(x0,x1,x2), name='a12')  # component a_{12}
a[0,1] = a01; a[0,2] = a02        # by antisymmetry, 6 nonzero components defined
a[1,2] = a12  
a.disp()

$\displaystyle a = a_{01}\left({x^{0}}, {x^{1}}, {x^{2}}\right) \mathrm{d} {x^{0}}\wedge \mathrm{d} {x^{1}} + a_{02}\left({x^{0}}, {x^{1}}, {x^{2}}\right) \mathrm{d} {x^{0}}\wedge \mathrm{d} {x^{2}} + a_{12}\left({x^{0}}, {x^{1}}, {x^{2}}\right) \mathrm{d} {x^{1}}\wedge \mathrm{d} {x^{2}}$

The arguments of components can be omitted

Manifold.options.omit_function_arguments=True
a.disp()

$\displaystyle a = a_{01} \mathrm{d} {x^{0}}\wedge \mathrm{d} {x^{1}} + a_{02} \mathrm{d} {x^{0}}\wedge \mathrm{d} {x^{2}} + a_{12} \mathrm{d} {x^{1}}\wedge \mathrm{d} {x^{2}}$

# mathematical object  of which "a" is an element
a.parent()

$\displaystyle \Omega^{2}\left(M\right)$

The matrix of 2-form must be antisymmetric.

a[:]

$\displaystyle \left(\begin{array}{rrr} 0 & a_{01} & a_{02} \\ -a_{01} & 0 & a_{12} \\ -a_{02} & -a_{12} & 0 \end{array}\right)$

We can check the formula $\ a_{ij}=a(\frac{\partial}{\partial x^0},\frac{\partial}{\partial x^1})$.

# check the matrix a(d/dxi,d/dxj)
fr = X.frame();
matrix([[a(fr[i0],fr[i1]).expr() for i1 in range(3)] for i0 in range(3)])  

$\displaystyle \left(\begin{array}{rrr} 0 & a_{01}\left({x^{0}}, {x^{1}}, {x^{2}}\right) & a_{02}\left({x^{0}}, {x^{1}}, {x^{2}}\right) \\ -a_{01}\left({x^{0}}, {x^{1}}, {x^{2}}\right) & 0 & a_{12}\left({x^{0}}, {x^{1}}, {x^{2}}\right) \\ -a_{02}\left({x^{0}}, {x^{1}}, {x^{2}}\right) & -a_{12}\left({x^{0}}, {x^{1}}, {x^{2}}\right) & 0 \end{array}\right)$

Example 14.4

Let us show how the antisymmetrization operation works for tensor fields $t\in T^{(0,2)}M$.

N = 3                         # dimension of manifold M 
M = Manifold(N, 'M')          # manifold M
X = M.chart(' '.join(['x'+str(i)+':x^{'+str(i)+'}' for i in range(N)]))  # chart on M
x0, x1, x2 = X[:]     # coordinates x^0, x^1, x^2 of chart X as the Python variables x0, x1, x2
t = M.tensor_field(0,2, name='t')  # tensor field of type (0,2)
f2 = [['f'+str(i)+str(j) for j in range(N)] for i in range(N)]   # table of comp. names
f = [[M.scalar_field(function(f2[j][k])(x0,x1,x2), name=f2[j][k]) # table of comp. funcs
      for k in range(N)] for j in range(N)]
t[:] = f                        # define all components

t.disp()                      # show tensor t

$\displaystyle t = f_{00} \mathrm{d} {x^{0}}\otimes \mathrm{d} {x^{0}} + f_{01} \mathrm{d} {x^{0}}\otimes \mathrm{d} {x^{1}} + f_{02} \mathrm{d} {x^{0}}\otimes \mathrm{d} {x^{2}} + f_{10} \mathrm{d} {x^{1}}\otimes \mathrm{d} {x^{0}} + f_{11} \mathrm{d} {x^{1}}\otimes \mathrm{d} {x^{1}} + f_{12} \mathrm{d} {x^{1}}\otimes \mathrm{d} {x^{2}} + f_{20} \mathrm{d} {x^{2}}\otimes \mathrm{d} {x^{0}} + f_{21} \mathrm{d} {x^{2}}\otimes \mathrm{d} {x^{1}} + f_{22} \mathrm{d} {x^{2}}\otimes \mathrm{d} {x^{2}}$

Now we apply the antisymmetrization operation, the result is 2-form at.

at = t.antisymmetrize()       # antisymmetrization of t
print(at)                     # show info on at

2-form on the 3-dimensional differentiable manifold M

Manifold.options.omit_function_arguments=True
at.disp()                     # show at

$\displaystyle \left( \frac{1}{2} \, f_{01} - \frac{1}{2} \, f_{10} \right) \mathrm{d} {x^{0}}\wedge \mathrm{d} {x^{1}} + \left( \frac{1}{2} \, f_{02} - \frac{1}{2} \, f_{20} \right) \mathrm{d} {x^{0}}\wedge \mathrm{d} {x^{2}} + \left( \frac{1}{2} \, f_{12} - \frac{1}{2} \, f_{21} \right) \mathrm{d} {x^{1}}\wedge \mathrm{d} {x^{2}}$

The matrix of components of $t$ is a full 3-3 matrix:

t[:]                          # components of t

$\displaystyle \left(\begin{array}{rrr} f_{00} & f_{01} & f_{02} \\ f_{10} & f_{11} & f_{12} \\ f_{20} & f_{21} & f_{22} \end{array}\right)$

and the antisymmetrized tensor has the antisymmetric matrix of components:

at[:]                        # components of antisymmetrization

$\displaystyle \left(\begin{array}{rrr} 0 & \frac{1}{2} \, f_{01} - \frac{1}{2} \, f_{10} & \frac{1}{2} \, f_{02} - \frac{1}{2} \, f_{20} \\ -\frac{1}{2} \, f_{01} + \frac{1}{2} \, f_{10} & 0 & \frac{1}{2} \, f_{12} - \frac{1}{2} \, f_{21} \\ -\frac{1}{2} \, f_{02} + \frac{1}{2} \, f_{20} & -\frac{1}{2} \, f_{12} + \frac{1}{2} \, f_{21} & 0 \end{array}\right)$

Example 14.5

Consider a tensor field of type $(0,3)$ on a 4-dimensional manifold and its antisymmetrization.

%display latex
N = 4                           # dimension of the manifold M
M = Manifold(N, 'M')          # manifold M
X = M.chart(' '.join(['x'+str(i)+':x^{'+str(i)+'}' for i in range(N)]))  # chart on M
t3 = M.tensor_field(0,3, name='t3') # tensor field of type (0,3)
print(t3)                     # show info on t3
                              # table of component names:
f3 = [[['a'+str(i)+str(j)+str(k) for k in range(N)] 
      for j in range(N)] for i in range(N)]
                              # table  of component functions:
f = [[[M.scalar_field(function(f3[i][j][k])(*X), name=f3[i][j][k]) 
       for k in range(N)] for j in range(N)] for i in range(N)]
t3[:] = f                     # define all components

Tensor field t3 of type (0,3) on the 4-dimensional differentiable manifold M

Let $e$ be the frame:

e = X.frame(); e                # frame of M

$\displaystyle \left(M, \left(\frac{\partial}{\partial {x^{0}} },\frac{\partial}{\partial {x^{1}} },\frac{\partial}{\partial {x^{2}} },\frac{\partial}{\partial {x^{3}} }\right)\right)$

and let's show how the definition of the antisymmetrization works if we take $(v_1,v_2,v_3)=(\frac{\partial}{\partial x^0}, \frac{\partial}{\partial x^1},\frac{\partial}{\partial x^2})$ in (14.2):

S3 = Permutations(3).list()     # S_3 group
s = 1/factorial(3)*sum([sign(p)*t3(e[p[0]-1],e[p[1]-1],e[p[2]-1])
                        for p in S3])     # special case of (14.2)

Manifold.options.omit_function_arguments=True
s.disp()                      # show the result
#s.expr()

$\displaystyle \begin{array}{llcl} & M & \longrightarrow & \mathbb{R} \\ & \left({x^{0}}, {x^{1}}, {x^{2}}, {x^{3}}\right) & \longmapsto & \frac{1}{6} \, a_{012} - \frac{1}{6} \, a_{021} - \frac{1}{6} \, a_{102} + \frac{1}{6} \, a_{120} + \frac{1}{6} \, a_{201} - \frac{1}{6} \, a_{210} \end{array}$

As we can see the sum contains 3!=6 elements with signs equal to permutations signs.

For comparison one can use the antisymmetrize method:

                                # built in
t3a = t3.antisymmetrize()       # antisymmetrization

t3a(e[0],e[1],e[2]).disp()

$\displaystyle \begin{array}{llcl} & M & \longrightarrow & \mathbb{R} \\ & \left({x^{0}}, {x^{1}}, {x^{2}}, {x^{3}}\right) & \longmapsto & \frac{1}{6} \, a_{012} - \frac{1}{6} \, a_{021} - \frac{1}{6} \, a_{102} + \frac{1}{6} \, a_{120} + \frac{1}{6} \, a_{201} - \frac{1}{6} \, a_{210} \end{array}$

Example 14.6

Define two 1-forms with symbolic components and compute their wedge product.

N = 3                               # dimension of manifold M
M = Manifold(N, 'M')                # manifold M
X = M.chart(' '.join(['x'+str(i)+':x^{'+str(i)+'}' for i in range(N)]))  # chart on M
a = M.diff_form(1,name='a')           # 1-form a
b = M.diff_form(1,name='b')           # 1-form b
ast = ['a'+str(j) for j in range(N)]  # list of component names 
bst = ['b'+str(j) for j in range(N)]  # list of component names
af = [M.scalar_field(function(ast[j])(*X),name=ast[j]) 
      for j in range(N)]             # list of component functions
bf = [M.scalar_field(function(bst[j])(*X),name=bst[j]) 
      for j in range(N)]             # list of component functions
a[:] = af                             # define all components of a
b[:] = bf                             # define all components of b

Manifold.options.omit_function_arguments=True
%display latex
a.disp()                            # show a

$\displaystyle a = a_{0} \mathrm{d} {x^{0}} + a_{1} \mathrm{d} {x^{1}} + a_{2} \mathrm{d} {x^{2}}$

b.disp()                            # show b

$\displaystyle b = b_{0} \mathrm{d} {x^{0}} + b_{1} \mathrm{d} {x^{1}} + b_{2} \mathrm{d} {x^{2}}$

# mathematical object of which "b" is an element.
b.parent()

$\displaystyle \Omega^{1}\left(M\right)$

Wedge product $a\wedge b$:

ab = a.wedge(b)                     # wedge product
ab.disp()                           # in SageMath

$\displaystyle a\wedge b = \left( -a_{1} b_{0} + a_{0} b_{1} \right) \mathrm{d} {x^{0}}\wedge \mathrm{d} {x^{1}} + \left( -a_{2} b_{0} + a_{0} b_{2} \right) \mathrm{d} {x^{0}}\wedge \mathrm{d} {x^{2}} + \left( -a_{2} b_{1} + a_{1} b_{2} \right) \mathrm{d} {x^{1}}\wedge \mathrm{d} {x^{2}}$

Compare with the result obtained from definition (14.3):

ab1 = factorial(1+1)/factorial(1)/factorial(1)*(a*b).antisymmetrize()
ab1.disp()                         # use formula (14.3)

$\displaystyle \left( -a_{1} b_{0} + a_{0} b_{1} \right) \mathrm{d} {x^{0}}\wedge \mathrm{d} {x^{1}} + \left( -a_{2} b_{0} + a_{0} b_{2} \right) \mathrm{d} {x^{0}}\wedge \mathrm{d} {x^{2}} + \left( -a_{2} b_{1} + a_{1} b_{2} \right) \mathrm{d} {x^{1}}\wedge \mathrm{d} {x^{2}}$

ab1 == ab                            # check if both methods agree

$\displaystyle \mathrm{True}$

Define the matrix containing the components of 1-forms in rows:

ma = matrix([a[:],b[:]])
ma

$\displaystyle \left(\begin{array}{rrr} a_{0} & a_{1} & a_{2} \\ b_{0} & b_{1} & b_{2} \end{array}\right)$

The components of the wedge product of 1-forms are minors of the matrix of components of these 1-forms.

ma.minors(2)

$\displaystyle \left[-a_{1} b_{0} + a_{0} b_{1}, -a_{2} b_{0} + a_{0} b_{2}, -a_{2} b_{1} + a_{1} b_{2}\right]$

Example 14.7

Consider the wedge product of three 1-forms (in 3-dimensional manifold).

# continuation
Manifold.options.omit_function_arguments=True
cst = ['c'+str(j) for j in range(N)]  # list of component names
cf = [M.scalar_field(function(cst[j])(*X),name=cst[j]) 
      for j in range(N)]              # list of component functions
c = M.diff_form(1,name='c')           # 1-form c
c[:] = cf                             # define all coefficients
abc = (a.wedge(b)).wedge(c)           # wedge prod. a /\ b /\ c
abc.disp()

$\displaystyle a\wedge b\wedge c = \left( -{\left(a_{2} b_{1} - a_{1} b_{2}\right)} c_{0} + {\left(a_{2} b_{0} - a_{0} b_{2}\right)} c_{1} - {\left(a_{1} b_{0} - a_{0} b_{1}\right)} c_{2} \right) \mathrm{d} {x^{0}}\wedge \mathrm{d} {x^{1}}\wedge \mathrm{d} {x^{2}}$

We can recognize in the result the Laplace expansion of the determinant $\$ det $\left(\begin{matrix} a_0,a_1,a_2\\ b_0,b_1,b_2\\ c_0,c_1,c2 \end{matrix}\right),\ \$ so we have
$\quad a\wedge b\wedge c=$ det $\left(\begin{matrix} a_0,a_1,a_2\\ b_0,b_1,b_2\\ c_0,c_1,c2 \end{matrix}\right) dx^0\wedge dx^1\wedge dx^2.$

We can check the observation using SageMath:

d = X.coframe()
d

$\displaystyle \left(M, \left(\mathrm{d} {x^{0}},\mathrm{d} {x^{1}},\mathrm{d} {x^{2}}\right)\right)$

# compare abc with the 3-form with the determinant as coefficient
abc == det(matrix([a[:],b[:],c[:]]))*(d[0].wedge(d[1])).wedge(d[2])

$\displaystyle \mathrm{True}$

The unique component of the wedge product $abc=a\wedge b\wedge c$ must be equal to $\ \ abc(\frac{\partial}{\partial x^0},\frac{\partial}{\partial x^1},\frac{\partial}{\partial x^2})$.

Let us check it with SageMath Manifolds.

X.frame()                         # default frame

$\displaystyle \left(M, \left(\frac{\partial}{\partial {x^{0}} },\frac{\partial}{\partial {x^{1}} },\frac{\partial}{\partial {x^{2}} }\right)\right)$

# check if the value of abc on the frame 
# is equal to det([a[:],b[:],c[:]])
abc(*X.frame()) == det(matrix([a[:],b[:],c[:]]))

$\displaystyle \mathrm{True}$

# the same in an other way
abc[0,1,2] == det(matrix([a[:],b[:],c[:]]))

$\displaystyle \mathrm{True}$

Example 14.8

Consider two 2-forms on a 4-dimensional manifold.

N = 4                         # dimension of manifold M
M = Manifold(N, 'M')        # manifold M
X = M.chart(' '.join(['x'+str(i)+':x^{'+str(i)+'}' for i in range(N)]))  # chart on M
a=M.diff_form(2,name='a')   # 2-form a
b=M.diff_form(2,name='b')   # 2-form b

We define first the names of functions,

def ast(i,j): return 'a'+str(i)+str(j) # names of comp. for a
    
def bst(i,j): return 'b'+str(i)+str(j) # names of comp. for b  

and the upper triangles of component matrices.

for i0 in range(N):         # component functions for a
    for i1 in range(N):      
        if i1>i0:
            a[i0,i1] = M.scalar_field(function(ast(i0,i1))(*X),
                                      name=ast(i0,i1))
for i0 in range(N):         # component functions for b
    for i1 in range(N):      
        if i1>i0:
            b[i0,i1] = M.scalar_field(function(bst(i0,i1))(*X),
                                      name=bst(i0,i1))

Manifold.options.omit_function_arguments=True
%display latex
a[:]                       # show components of a

$\displaystyle \left(\begin{array}{rrrr} 0 & a_{01} & a_{02} & a_{03} \\ -a_{01} & 0 & a_{12} & a_{13} \\ -a_{02} & -a_{12} & 0 & a_{23} \\ -a_{03} & -a_{13} & -a_{23} & 0 \end{array}\right)$

b[:]                       # show components of b

$\displaystyle \left(\begin{array}{rrrr} 0 & b_{01} & b_{02} & b_{03} \\ -b_{01} & 0 & b_{12} & b_{13} \\ -b_{02} & -b_{12} & 0 & b_{23} \\ -b_{03} & -b_{13} & -b_{23} & 0 \end{array}\right)$

Recall that according to the general rule, to obtain the wedge product we use all possible strictly increasing and disjoint sequences $I=(i_1,i_2),\ J=(j_1,j_2)$ and define $dx^K=e^1\wedge e^2\wedge e^3\wedge e^4,$ where $K=(1,2,3,4)$ is reordered disjoint union $I\cup J$. We have

$$\Big(\sum_I a_Idx^I\Big)\wedge\Big(\sum_J b_Jdx^J\Big)= \sum_K\Big(\sum_{\substack{I\cup J=K\\I\cap J=\emptyset}}\pm a_Ib_J\Big)dx^K.$$

To be less formal, to compute the wedge product we

a) take all possible components of the first 2-form with increasing permutations $I$ of two indices from (0,1,2,3);

b) multiply them by components of the second 2-form with increasing permutations $J$ of the remaining 2 indices;

c) the products are taken with + or - depending on the sign of the joined permutation $(I,J)$ of (0,1,2,3).

w = a.wedge(b)                    # w=a/\b
w.display()                       # show w

$\displaystyle a\wedge b = \left( a_{23} b_{01} - a_{13} b_{02} + a_{12} b_{03} + a_{03} b_{12} - a_{02} b_{13} + a_{01} b_{23} \right) \mathrm{d} {x^{0}}\wedge \mathrm{d} {x^{1}}\wedge \mathrm{d} {x^{2}}\wedge \mathrm{d} {x^{3}}$

What's next?

Take a look at the notebook Pullback of tensor fields.