Topics in Geometric Mechanics: Week 8

$ %% not provided in align* environment \def\newsavebox#1{} \def\intertext#1{\text{#1}} \def\let#1#2{} \def\makeatletter{} \def\makeatother{} \newenvironment{theorem}{\textbf{Theorem.}\rm}{} %%% Local Variables: %%% mode: latex %%% TeX-master: t %%% End: %% sets \def\R{\mathbf{R}} \def\Z{\mathbf{Z}} \def\C{\mathbf{C}} \def\Q{\mathbf{Q}} \def\sphere#1{\mathbf{S}^{#1}} \def\rp#1{\R P^{#1}} \def\set#1{\left\{ #1 \right\}} \def\st{\,\mathrm{s.t.}\,} \def\extalg#1#2{\Lambda^{#1}(#2)} \def\symalg#1#2{\mathsf{S}^{#1}(#2)} \def\tenalg#1#2{\mathsf{T}^{#1}(#2)} \def\hom#1#2{\mathrm{Hom}(#1;#2)} \def\extproduct{\wedge} \def\symmetricproduct{\cdot} \def\tensorproduct{\otimes} \def\image#1{\mathrm{Im}\,#1} \def\kernel#1{\mathrm{Ker}\,#1} \def\Trace[#1]#2{\mathrm{Tr}\ifx#1\empty\else{}_{#1}\fi\ifx#2\empty\else\,(#2)\fi} \def\orbit#1#2{{#1}\cdot{#2}} \def\stab#1#2{\mathrm{stab}_{#1}(#2)} \def\Span{\mathrm{span}} %% Lie groups \def\GL#1{\mathrm{GL}(#1)} \def\SL#1{\mathrm{SL}(#1)} \def\Symp#1{\mathrm{Sp}(#1)} \def\Orth#1{\mathrm{O}(#1)} \def\SOrth#1{\mathrm{SO}(#1)} \def\Unitary#1{\mathrm{U}(#1)} \def\SUnitary#1{\mathrm{SU}(#1)} \def\Torus#1{\mathbf{T}^{#1}} \def\H#1{\mathbf{H}^{#1}} \def\Diagonal#1{\mathbf{Diag}(#1)} \def\diag#1{\mathbf{diag}\left(#1\right)} \def\Liegp#1{\mathrm{#1}} %% Lie algebras \def\Mat#1#2{\mathrm{Mat}_{{#1} \times {#1}}(#2)} \def\Matrect#1#2#3{\mathrm{Mat}_{{#1} \times {#2}}(#3)} \def\gl#1{\mathfrak{gl}(#1)} \def\spl#1{\mathfrak{sl}(#1)} \def\symp#1{\mathfrak{sp}(#1)} \def\orth#1{\mathfrak{o}(#1)} \def\sorth#1{\mathfrak{so}(#1)} \def\unitary#1{\mathfrak{u}(#1)} \def\sunitary#1{\mathfrak{su}(#1)} \def\torus#1{\mathfrak{t}^{#1}} \def\liealg#1{\mathfrak{#1}} \def\liebracket#1#2{[#1,#2]} \def\Ad#1#2{\mathrm{Ad}_{#1}{#2}} \def\ad#1#2{\mathrm{ad}_{#1}{#2}} \def\coAd#1#2{\mathrm{Ad}^*_{#1}{#2}} \def\coad#1#2{\mathrm{ad}^*_{#1}{#2}]} %% operators \def\D#1{\,\mathrm{d}{#1}\,} \def\Dat#1#2{\,\mathrm{d}_{#1}{#2}\,} \def\lieder#1#2{{\sf{L}}_{#1}{#2}} \def\covder#1#2{\,\frac{D #2}{d{#1}}} \def\liebracket#1#2{\left[#1,#2\right]} \def\crossproduct{{\boldsymbol{\times}}} \def\ip#1#2{\langle #1, #2 \rangle} \def\interiorproduc#1#2{\iota_{#1}{#2}} \def\norm#1{|#1|} \def\lt{<} \def\gt{>} \def\ddt#1#2{\frac{d{#2}}{d{#1}}} \def\dndtn#1#2#3{\frac{d^{#3}{#2}}{d{#1}^{#3}}} \def\wint#1#2{\int\limits_{#1}^{#2}} \def\didi#1#2{\frac{\partial{#1}}{\partial{#2}}} \def\dindin#1#2#3{\frac{\partial^{#3}{#1}}{\partial{#2}^{#3}}} \def\eulerlagrange#1#2#3#4{\ddt{#1}{}\didi{#4}{#3} - \didi{#4}{#2}} %% labels, etc. \def\labelenumi{{\bf\Alph{enumi}.}} \def\solutioncolor{blue} \def\solutiontopsep{3mm} \def\solutionbotsep{0mm} \newenvironment{solution}{\vspace{\solutiontopsep}\noindent\textbf{Solution}. \textcolor{\solutioncolor}\bgroup}{\egroup\vspace{\solutionbotsep}} %% redo quote environment \def\quote{} \def\endquote{} %% redo enumerate environment \makeatletter \let\oldenumerate\enumerate \let\oldendenumerate\endenumerate \def\enumerate{\bgroup\oldenumerate\setlength{\itemsep}{3mm}\setlength{\topsep}{3mm}} \def\endenumerate{\oldendenumerate\egroup} \makeatother \def\museincludegraphicsoptions{% width=0.25\textwidth% } \def\museincludegraphics{% \begingroup \catcode`\|=0 \catcode`\\=12 \catcode`\#=12 \expandafter\includegraphics\expandafter[\museincludegraphicsoptions] } \def\musefigurewidth{!} \def\musefigureheight{!} \newsavebox{\mypicboc}% \def\includefigure#1#2#3#4{% \edef\fileending{#4}% \def\svgend{svg}% \def\pdftend{pdf_t}% \def\pdftexend{pdf_tex}% \def\bang{!}% \begin{lrbox}{\mypicboc}% \ifx\fileending\pdftend% \input{#3.#4}% \else% \ifx\fileending\pdftexend \input{#3.#4}% \else% \ifx\fileending\svgend% \input{#3.#4}% \fi\fi% \includegraphics{#3.#4}% \fi% \end{lrbox}% \edef\thiswidth{#1}\edef\thisheight{#2}% \ifx\thiswidth\bang% \ifx\thisheight\bang% \resizebox{!}{!}{\usebox{\mypicboc}}% \else% \resizebox{!}{\thisheight}{\usebox{\mypicboc}}% \fi% \else% \ifx\thisheight\bang% \resizebox{\thiswidth}{!}{\usebox{\mypicboc}}% \else% \resizebox{\thiswidth}{\thisheight}{\usebox{\mypicboc}}% \fi\fi} $
Navigation
MTH-696a
Course outline
Assignments & Solutions
Schedule
Lecture Notes
Sources
Home Page of Leo Butler

Week 8

This week we will look at

The Euler-Lagrange Equations

Recall

$L : TM \to \R$ is a smooth function of the position $x$ and velocity $\dot{x}$. The Euler-Lagrange equations are \begin{align*} \eulerlagrange{t}{x}{\dot{x}}{L} &= 0 \end{align*}

The Central Force Problem

Conservative, central force \begin{align*} F(x) &= -\lambda(r)x &&\implies&& F(x) = -\nabla V(r), &&&& r=\norm{x} \end{align*}

The central force problem.
The central force problem.

Conservative, central force \begin{align*} L &= \frac{1}{2} \norm{\dot{x}}^2 - V(r)\\ \eulerlagrange{t}{x}{\dot{x}}{L} &= 0 &&\implies&& \ddot{x} = - \lambda(r)x. \end{align*}

The central force problem.
The central force problem.

Quadratic Potential

Suppose that

\begin{align*} V(x) &= \frac{1}{2} \norm{x}^2 = \frac{1}{2} r^2, && \ddot{x} = -x. \end{align*}

Newtonian Potential

Suppose that

\begin{align*} V(x) &= -\dfrac{1}{r}, && \ddot{x} = - \dfrac{x}{r^3}. \end{align*} Then \begin{align*} L &= \frac{1}{2} (\dot{r}^2 + r^2 \dot{\theta}^2) + \dfrac{1}{r}. \end{align*} So \begin{align*} \ddt{t}{} r^2 \dot{\theta} &= 0 && \ddot{r} = r \dot{\theta}^2 - \dfrac{1}{r^2}. \end{align*}

Solutions

Substitute \begin{align*} \ddt{t}{\theta} &= \dfrac{c}{r^2}, && r = u^{-1} && \ddt{t}{r} = c u^2 \ddt{\theta}{r} = -c \ddt{\theta}{u}. \end{align*} This gives \begin{align*} %%-c^2 u^2 \dndtn{\theta}{u}{2} &= c^2 u^3 - u^2,\\ \dndtn{\theta}{u}{2} &= c^{-2} - u. \end{align*}

The solutions are conic sections

\[ r = \dfrac{1}{c^{-2} + e \cos(\theta-\theta_0)} \]

An elliptical, parabolic and hyperbolic orbit.
An elliptical, parabolic and hyperbolic orbit.

Tensoriality

The Euler-Lagrange equations \[ \eulerlagrange{t}{x}{\dot{x}}{L} = 0 \qquad (*) \] are tensorial. That is, if $y = \phi(x)$ is a change of coordinates, then (*) takes the same form in the $y$ system of coordinates.

Proof 1

Boring.

Let $v=\dot{x}$ and $w=\dot{y}$, so $w=\D{\phi}\cdot v$. \begin{align*} \didi{L}{v_i} &= \didi{L}{w_{\alpha}} \didi{\phi_{\alpha}}{x_i} && \didi{L}{x_i} = \didi{L}{y_{\alpha}} \didi{\phi_{\alpha}}{x_i} + \didi{L}{w_{\alpha}} \didi{w_{\alpha}}{x_i} \end{align*} Therefore \begin{align*} \underbrace{% \eulerlagrange{t}{x_i}{v_i}{L}% }_{\text{LHS}} &= \underbrace{% \left[ \eulerlagrange{t}{y_{\alpha}}{w_{\alpha}}{L} \right]% }_{\text{RHS}}% \didi{\phi_{\alpha}}{x_i} + \didi{L}{w_{\alpha}}% \underbrace{% \left[ \ddt{t}{}\didi{\phi_{\alpha}}{x_i} - \didi{w_{\alpha}}{x_i} \right]% }_{\text{NIL}} \end{align*} But $y_{\alpha} = \phi_{\alpha}$ and $w_{\alpha} = \displaystyle\ddt{t}{y_{\alpha}}$, so NIL vanishes. Since $\phi$ is a diffeomorphism, the LHS vanishes for all $i$ iff RHS vanishes for all $\alpha$.

Proof 2

Let $a \lt b$ and let $p, q \in M$. Define \[ \Omega_{p, q}^{a, b} = \set{ c : [a, b] \to M \st c(a)=p, c(b)=q, c \textrm{ is piecewise AC} }. \] Define the action functional $A : \Omega_{p, q}^{a, b} \to \R$ by \[ A(c) := \int_a^b L(c(t), \dot{c}(t)) \D{t}. \]

If $c$ is a critical point of $A$, then $c$ solves the Euler-Lagrange equations.

This implies that the Euler-Lagrange equations are tensorial.

Variations of a curve.
Variations of a curve.

Proof 2 - cont'd

Let $x=(x_i)$ be a a coordinate system on $M$ and let $c_s \in \Omega_{p, q}^{a, b}$ be a smooth path of curves. Then \begin{align*} \ddt{s}{} A(c_s) &= \int_a^b \left[ \didi{L}{x_{\alpha}} \didi{}{s} c_s^{\alpha} + \overbrace{% \didi{L}{\dot{x}_{\alpha}}% }^{u} \underbrace{% \ddt{t}{} \didi{}{s} c^{\alpha}_s% }_{\D{v}}% \right] \D{t} \\ &= \int_a^b \underbrace{% \left[ \didi{L}{x_{\alpha}} - \ddt{t}{} \didi{L}{\dot{x}_{\alpha}} \right]% }_{E_{\alpha}} \underbrace{% \didi{c^{\alpha}_s}{s}% }_{v^{\alpha}} \D{t} + \underbrace{% \left[ \didi{L}{\dot{x}_{\alpha}} \didi{c^{\alpha}_s}{s} \right]_{t=a}^{t=b} }_{\text{NIL}} \end{align*}

Fundamental Theorem

Let $f : [a, b] \to \R$ be a continuous function. If, for all continuous functions $g$ such that $g(a)=g(b)=0$, \[ \int_a^b f(t)g(t) \D{t} = 0 \] then $f \equiv 0$.

If there is an $x \in [a, b]$ such that $f(x) \gt 0$, then by continuity, there is an $x_0 \in (a, b)$ such that $f(x_0) \gt 0$. By continuity, there is a $\delta \gt 0$ such that $f(x)^2 \gt \frac{1}{2} f(x_0)^2$ for all $x \in [x_0 - 2\delta, x_0 + 2\delta] \subset (a, b)$. Let $g : [a, b] \to \R$ be a non-negative continuous function such that $g = f$ on $[x_0 - \delta, x_0 + \delta]$ and $g(a)=g(b)=0$. Therefore \[ \int_a^b f(t)g(t) \D{t} \geq \delta f(x_0)^2 \gt 0. \]

Proof 2 - finished

The variation vector field $v = (v_{\alpha})$ can be an arbitrary continuous function that vanishes at $a$ and $b$. Now apply the fundamental theorem.

Legendre Transform

Let $L : TM \to \R$ be $C^2$. Define the Legendre transform $F : TM \to T^*M$ by \[ F(x, v) = (x, \didi{L}{v}). \]

If $L$ is $C^2$ and convex, then $F$ is a $C^1$ diffeomorphism. Moreover, the function \[ H(x, p) = \max_{v \in T_xM} \ip{p}{v} - L(x, v) \] is $C^2$ and $c$ solves the Euler-Lagrange equations for $L$ iff $\gamma = F \circ c$ satisfies Hamilton's equations for $H$.

Proof 3

If $L$ is convex, i.e. \[ \dindin{L}{v}{2} \gt 0 \] then there is a unique $v \in T_xM$ maximizing $\ip{p}{v} - L$. It solves \[ p = \didi{L}{v} \] for $v$. This implies $H$ is $C^2$. Also, $L$ satisfies \[ L(x, v) = \max_{p \in T^*_xM} \ip{p}{v} - H(x, p). \]

Proof 3 - cont'd

In sum \begin{align*} p &= \didi{L}{v} && v = \didi{H}{p} && \didi{H}{x} = -\didi{L}{x} \end{align*} Which implies \begin{align*} \dot{p} &= \underbrace{% \ddt{t}{} \didi{L}{v} = \didi{L}{x}% }_{\text{Euler-Lagrange}}% = -\didi{H}{x} && \dot{x} = v = \didi{H}{p}. \end{align*}

Let \begin{align*} \Omega &= \sum_i \D{x}_i \wedge \D{p}_i && X_H = \sum_i \didi{H}{p_i} \didi{}{x_i} - \didi{H}{x_i} \didi{}{p_i} \end{align*} We see that $\interiorproduc{X_H}{\Omega} = \D{H}$.

It suffices to prove that $\Omega$ is a tensor field.

The Symplectic Form

The $2$-form $\Omega = \sum_i \D{x}_i \wedge \D{p}_i$ is a tensor field on $T^*M$.

Let $y = \phi(x)$ be a change of coordinates, so that $p = (\D{\phi})^* q$ where $\D{\phi}^* : T^*_y M \to T^*_x M$. Then \begin{align*} \sum_{\alpha} \D{y_\alpha} \wedge \D{q_{\alpha}} &= \sum_{\alpha, \beta} \didi{\phi_{\alpha}}{x_{\beta}} \D{x_{\beta}} \wedge \D{q_{\alpha}} \\ &= \sum_{\beta} \D{x_{\beta}} \wedge \left[ \sum_{\alpha} \didi{\phi_{\alpha}}{x_{\beta}} \D{q_{\alpha}} \right] \\ &= \sum_{\beta} \D{x_{\beta}} \wedge \D{p_{\beta}}. \end{align*}