Geodesics of left invariant metrics on matrix Lie groups – Part 2 Conservation laws

In the last post, Geodesics of left invariant metrics on matrix Lie groups – Part 1,we have derived Arnold’s equation – that is a half of the problem of finding geodesics on a Lie group endowed with left-invariant metric.

Suppose G is a Lie group, and g(\xi,\eta) is a scalar product (i.e. a nondegenerate bilinear form) on its Lie algebra Lie(G). Then, using left translations g defines a left invariant (Riemannian or pseudo-Riemannian) metric on the whole group G. If t\mapsto a(t)\in G is a path in G, we use left translations to define the image \omega(t) of the tangent vector \dot{a}(t) in Lie(G)

(1)   \begin{equation*}\omega(t)=a(t)^{-1}\dot{a}(t).\end{equation*}

Usually it is written more carefully, using L_{a^{-1}} instead of a^{-1}, but I am using a simplified notation, well adapted to dealing with problems for matrix groups.

If a(t) is a geodesic for the metric g, then \omega(t) satisfies the Arnold’s equation

(2)   \begin{equation*}\dot{\omega}=B(\omega,\omega),\end{equation*}

where

    \[B:Lie(G)\times Lie(G)\rightarrow Lie(G)\]

is defined as

(3)   \begin{equation*}B^k_{ij}=g^{kl}B_{ij,l}=g^{kl}C_{jl,i}=g^{kl}g_{im}C_{jl}^m,\end{equation*}

and C_{jl}^m are the structure constants of G, that is

(4)   \begin{equation*}[\xi_i,\xi_j]=C_{ij}^k\, \xi_k\end{equation*}

where \xi_i form a basis for Lie(G) and

(5)   \begin{equation*}B_{ij,l}=g_{lk}B_{ij}^k,\,C_{ij,l}=g_{lk}C_{ij}^k.\end{equation*}

Equation (2) is a system of nonlinear ordinary differential equations with constant coefficients. In this form it can be found in Arnold’s 1966 paper “Sur la g\’eom\’etrie differentielle des groupes de Lie de dimension infinie et ses applications \`a l’hydrodynamique des fluides parfaits”, Ann. Inst.
Fourier (Grenoble) 16 (1966), 319-361.

In his blog post “The Euler-Arnold equation” Terrence Tao mentions that some people call this equation the “Euler-Arnold” equation, while some other prefer to skip Arnold’s name completely and call it “Euler-Poincare equation”. Go-figure!

Solving equations (2) is just one half of the whole problem of finding a geodesic. Once \omega^{i}(t) are known, we need to solve the linear differential equation with variable coefficients that results from the definition of \omega (??):

(6)   \begin{equation*} \dot{a}(t)=a(t)\omega(t),\end{equation*}

For the case of the rotation group and free rigid body we did it in Taming the T-handle continued.

B. Kolev in his paper “Lie groups and mechanics. An introduction” shows that for the rotation group O(n) the equations of motion for the free rigid body are completely integrable. We will not need this result, but is useful to know that in general case there are always two quadratic constants of motion corresponding to “kinetic energy” and “square of the angular momentum”. We were discussing these constants of motion for the group O(3) in
“Asymmetric Spinning Top – The Hardest Concept To Grasp In Physics” – they are used in so-called Poinsot construction. We will discuss a version of it for the group O(2,1) in the future.

The first observation is that the “kinetic energy” T(t)=g_{ij}\omega(t)^i\omega(t)^j is, in fact, a constant, independent of t. To see that this is the case, we differentiate:

(7)   \begin{equation*}\frac{dT(t)}{dt}=g_{kl}\dot{\omega}^{k}\omega^l+g_{kl}\omega^{k}\dot{\omega^l}=2g_{kl}\omega^k\dot{\omega}^l.\end{equation*}

In the previous post we wrote Eq. (2) as

(8)   \begin{equation*}g_{kl}\,\dot{\omega}^l=C_{jk,i}\,\omega^{i}\omega^{j}.\end{equation*}

Substituting into Eq. (7) we get

(9)   \begin{equation*}\frac{dT(t)}{dt}=2C_{jk,i}\omega^{i}\omega^j \omega^{k}=0.\end{equation*}

The result is zero because C_{jk,i} is antisymmetric in j,k while \omega^j\omega^k is symmetric. That is an often used property: if A^{ij}=-A^{ji} and B_{ij}=B_{ji}, then the contraction A^{ij}B_{ij}=0. Indeed A^{ij}B_{ij}=-A^{ji}B_{ij}=-A^{ji}B_{ji}=-A^{ij}B_{ij}, where in the last equality we have exchanged dummy indices names i\leftrightarrow j. If a number is equal to its negative, it must be zero.

To get the formula for the second quadratic invariant we need to return to the Ad-invariant scalar product that we have denoted \mathring{g} in Killing vectors, geodesics, and “Noether’s theorem”:

(10)   \begin{equation*}\mathring{g}(\eta_1,\eta_2)=\mbox{const}\, \frac{1}{2}Re(\mbox{Tr}(\eta_1\eta_2))).\end{equation*}

The fact that \mathring{g} is Ad-invariant implies an important relation between the matrix \mathring{g}_{ij} and the structure constants C_{ij}^k
that we are going to use. Ad-invariance means that:

(11)   \begin{equation*}\mathring{g}(e^{t\eta}\xi_1 e^{-t\eta},e^{t\eta}\xi_2 e^{-t\eta})=\mathring{g}(\xi_1,\xi_2)\end{equation*}

for all t\in\mathbf{R} and \xi_1,\xi_2,\eta\in Lie(G).
Differentiating at t=0 we get

(12)   \begin{equation*}\mathring{g}([\eta,\xi_1],\xi_2)+\mathring{g}(\xi_1,[\eta,\xi_2])=0.\end{equation*}

Setting \xi_1=\xi_i,\xi_2=\xi_j,\eta=\xi_k we get

(13)   \begin{equation*}C_{ki}^l\mathring{g}_{lj}+C_{kj}^l\mathring{g}_{li}=0.\end{equation*}

Multiplying both sides by g^{im}g^{jn} we obtain

(14)   \begin{equation*} C_{ki}^n\,g^{im}+C_{kj}^m\,g^{jn}=0\end{equation*}

We can now derive the second conservation law. The angular momentum m_i is defined as

(15)   \begin{equation*}m_i=g_{ij}\omega^{j}.\end{equation*}

Notice that m is a covector, a one-form on Lie(G), it is in the dual Lie(G)^* of Lie(G). It is the metric that connects the space to its dual. While vectors in Lie(G) play an active role, they generate transformations, elements in the dual, one-forms from Lie(G)^*, are “passive”, they evaluate vectors to numbers. It is the metric that is the third element here, that allows the active principle to connect to the passive principle. The metric depends on the mass distribution. In application to rigid bodies the inertia tensor is encoded in the metric on the rotation group.

The second conservation law states that the square of the angular momentum evaluated with the Ad-invariant metric \mathring{g} is constant:

(16)   \begin{equation*} m_0^2(t)=\mathring{g}^{ij}m_i(t)m_j(t)=\mbox{const}.\end{equation*}

To verify we differentiate and use Eq. (8) rewritten as

(17)   \begin{equation*}\dot{m}_k=C_{jk}^{i}\,m_{i}\omega^{j}.\end{equation*}

(18)   \begin{equation*}\frac{dm_0^2(t)}{dt}=2\mathring{g}^{kl}\dot{m}_k m_l=2\mathring{g}^{kl}C_{jk}^{i}m_im_l\omega^j.\end{equation*}

Now, according to Eq. (14) \mathring{g}^{kl}C_{jk}^{i} is antisymmetric in (i,l), while the product m_im_l is symmetric, therefore we get zero:

(19)   \begin{equation*}\frac{dm_0^2(t)}{dt}=0.\end{equation*}

In the following posts we will first return to the case of the rotation group in three dimensions and the rigid body, and then try to apply a similar reasoning to the case of the Lorentz group O(2,1) in 2+1 dimensions.