The Principle of Least Action

Jason Gross, December 7, 2010
Last Updated September 23, 2023

Introduction

Recall that we defined the Lagrangian to be the kinetic energy less potential energy, L=K-U, at a point.  The action is then defined to be the integral of the Lagrangian along the path,

$S{\int }_{t_0}^{t_1}Ldt{\int }_{t_0}^{t_1}K-Udt$

It is (remarkably!) true that, in any physical system, the path an object actually takes minimizes the action.  It can be shown that the extrema of action occur at

$\frac{\partial L}{\partial q}-\frac{d}{dt}\frac{\partial L}{\partial \stackrel{.}{q}}0$

This is called the Euler equation, or the Euler-Lagrange Equation.

Derivation

Least action: F=m a

Suppose we have the Newtonian kinetic energy, $K=\frac{1}{2}m{v}^2$, and a potential that depends only on position, $U=U\left(\stackrel{\rightharpoonup }{r}\right)$.  Then the Euler-Lagrange equations tell us the following:

$-U^{\prime }\left[r\left[t\right]\right]-m{r}^{\mathrm{\prime \prime }}\left[t\right]⩵0$

Rearrangement gives

Least action with no potential

Suppose we have no potential, U=0.  Then L=K, so the Euler-Lagrange equations become

$\frac{\partial K}{\partial q}-\frac{d}{dt}\frac{\partial K}{\partial \stackrel{.}{q}}0$

For Newtonian kinetic energy, $K=\frac{1}{2}m{\stackrel{.}{x}}^2$, this is just

This is a straight line, as expected.

Least action with gravitational potential

Suppose we have gravitational potential close to the surface of the earth, U=m g y, and Newtonian kinetic energy, $K=\frac{1}{2}m{\stackrel{.}{y}}^2$.  Then the Euler-Lagrange equations become

This is a parabola, as expected.

Constants of motion: Momenta

We may rearrange the Euler-Lagrange equations to obtain

$\frac{\partial L}{\partial q}=\frac{d}{dt}\frac{\partial L}{\partial \stackrel{.}{q}}$

If it happens that $\frac{\partial L}{\partial q}=0$, then $\frac{d}{dt}\frac{\partial L}{\partial \stackrel{.}{q}}$ is also zero.  This means that $\frac{\partial L}{\partial \stackrel{.}{q}}$ is a constant (with respect to time).  We call $\frac{\partial L}{\partial \stackrel{.}{q}}$ a (conserved) momentum of the system.

Linear Momentum

By noting that Newtonian kinetic energy, $K=\frac{1}{2}m{v}^2$, is independent of the time derivatives of position, if potential energy depends only on position, we can infer that $\frac{\partial L}{\partial \stackrel{.}{x}}$ (and, similarly, $\frac{\partial L}{\partial \stackrel{.}{y}}$ and $\frac{\partial L}{\partial \stackrel{.}{z}}$) are constant.  Then $\frac{\partial L}{\partial \stackrel{.}{x}}=\frac{\partial }{\partial \stackrel{.}{x}}\left(\frac{1}{2}m{\stackrel{.}{x}}^2\right)=m\stackrel{.}{x}$.  This is just standard linear momentum, m v.

Angular Momentum

Let us change to polar coordinates.

$\frac{1}{2}m{r^{\prime }\left(t\right)}^2+\frac{1}{2}m{r\left(t\right)}^2{{\theta }^{\prime }\left(t\right)}^2$

Using dot notation, this is

$K/.\mathrm{r\_}\text{'}\left[t\right]\to \mathrm{OverDot}\left[r\right]/.\mathrm{r\_}\left[t\right]\to r//\mathrm{TraditionalForm}$

$\frac{1}{2}{\stackrel{.}{\theta }}^{2}m{r}^2+\frac{m{\stackrel{.}{r}}^2}{2}$

Note that θ does not appear in this expression.  If potential energy is not a function of θ (is only a function of r), then $\frac{\partial L}{\partial \stackrel{.}{\theta }}=m{r}^2\stackrel{.}{\theta }$ is constant.  This is standard angular momentum, $m{r}^2\omega =rmr\omega =r\times mv$.

Classic Problem: Brachistochrone (“shortest time”)

Problem

A bead starts at x=0, y=0, and slides down a wire without friction, reaching a lower point $\left(x_f,y_f\right)$.  What shape should the wire be in order to have the bead reach $\left(x_f,y_f\right)$ in as little time as possible.

Solution

Idea

Use the Euler equation to minimize the time it takes to get from $\left(x_i,y_i\right)$ to $\left(x_f,y_f\right)$.

Implementation

Letting ds be the infinitesimal distance element and v be the travel speed,

Now we apply the Euler equation to $f=\sqrt{\frac{1+{\left(x\text{'}\right)}^2}{2gy}}$ and change t→y, $\stackrel{.}{x}\to x\text{'}$.

Squaring both sides and making a special choice for the constant gives

To solve this, change variables:

$y=A\left(1-\cos \left(\theta \right)\right),dy=A\sin \left(\theta \right)d\theta$

$\mathrm{FullSimplify}\left[2Ay-y^2/.y\to A\left(1-\mathrm{Cos}\left[\theta \right]\right)\right]$

$A^2{\mathrm{Sin}\left[\theta \right]}^2$

Full solution: The brachistochrone is described by

$x=A\left(\theta -\sin \left(\theta \right)\right)\text{}y=A\left(1-\cos \left(\theta \right)\right)$

There’s no analytic solution, but we can compute them.

Classic Problem: Catenary

Problem

Suppose we have a rope of length l and linear mass density λ.  Suppose we fix its ends at points $\left(x_0,y_0\right)$ and $\left(x_f,y_f\right)$.  What shape does the rope make, hanging under the influence of gravity?

Solution

Idea

Calculate the potential energy of the rope as a function of the curve, y(x), and minimize this quantity using the Euler-Lagrange equations.

Implementation

Suppose we have curve parameterized by t, (x(t), y(t)).  The potential energy associated with this curve is

Note that if we choose to factor ds the other way (for y'), we get a mess.

Now we apply the Euler-Lagrange equation to $f=\lambda gy\sqrt{1+{\left(x\text{'}\right)}^2}$ and change t→y, $\stackrel{.}{x}\to x\text{'}$.

Since $\frac{\partial f}{\partial x}=0$, $\frac{\partial f}{\partial x\text{'}}$ is constant, say $a=\frac{1}{\lambda g}\frac{\partial f}{\partial x\text{'}}=\frac{yx\text{'}}{\sqrt{1+{\left(x\text{'}\right)}^2}}$.  Then

$x\text{'}=\frac{dx}{dy}=\pm \frac{a}{\sqrt{y^2-a^2}}$

Using the fact that

$\int \frac{dy}{\sqrt{y^2-a^2}}{\cosh }^{-1}\left(\frac{y}{a}\right)+b,$

integration of x' gives

$x\left(y\right)=\pm a{\cosh }^{-1}\left(\frac{y}{a}\right)+b$

where b is a constant of integration.

Plotting this for a=1, b=0 gives:

$\mathrm{Clear}\left[y\right];\mathrm{Manipulate}\left[\mathrm{ParametricPlot}\left[\left\{\left\{-a\mathrm{ArcCosh}\left[\frac{t}{a}\right]+b,t\right\},\left\{a\mathrm{ArcCosh}\left[\frac{t}{a}\right]+b,t\right\}\right\},\left\{t,\mathrm{ymin},\mathrm{ymax}\right\},\mathrm{PlotStyle}\to \mathrm{Black},\mathrm{AspectRatio}\to \mathrm{Automatic}\right],\left\{\left\{a,1\right\},-5,5\right\},\left\{\left\{b,0\right\},-5,5\right\},\left\{\left\{\mathrm{ymin},0,y_{\min }\right\},-5,5\right\},\left\{\left\{\mathrm{ymax},2,y_{\max }\right\},-5,5\right\}\right]$

Problem: Bead on a Ring

From 8.033 Quiz #2

Problem

A bead of mass m slides without friction on a circular hoop of radius R. The angle θ is defined so that when the bead is at the bottom of the hoop, θ=0. The hoop is spun about its vertical axis with angular velocity ω. Gravity acts downward with acceleration g.
Find an equation describing how θ evolves with time.
Find the minimum value of ω for the bead to be in equilibrium at some value of θ other than zero.
(“equilibrium” means that $\stackrel{.}{\theta }$ and $\ddot{\theta }$ are both zero.) How large must ω be in order to make θ=π/2?

Solution

The general Lagrangian for the object in Cartesian coordinates is

$\mathrm{Clear}\left[x,y,z,t\right];\left(L=\frac{1}{2}m\left({x\text{'}\left[t\right]}^2+{y\text{'}\left[t\right]}^2+{z\text{'}\left[t\right]}^2\right)-mgz\left[t\right]\right)//\mathrm{TraditionalForm}$

$\frac{1}{2}m\left({x^{\prime }\left(t\right)}^2+{y^{\prime }\left(t\right)}^2+{z^{\prime }\left(t\right)}^2\right)-gmz\left(t\right)$

Converting to polar coordinates, and using the constraints that φ=ω t and r=R, using the conversion

gives

$LgmR\cos \left(\theta \right)-gmR-\frac{1}{4}m{R}^2{\omega }^2\cos \left(2\theta \right)+\frac{1}{2}{\stackrel{.}{\theta }}^{2}m{R}^2+\frac{1}{4}m{R}^2{\omega }^2$

$0\frac{\partial L}{\partial \theta }-\frac{d}{dt}\frac{\partial L}{\partial \stackrel{.}{\theta }}-gmR\sin \left(\theta \right)+m{R}^2{\omega }^2\sin \left(\theta \right)\cos \left(\theta \right)-m{R}^2{\theta }^{\mathrm{\prime \prime }}\left(t\right)$

${\theta }^{\mathrm{\prime \prime }}\left(t\right)\frac{R{\omega }^2\sin \left(\theta \left(t\right)\right)\cos \left(\theta \left(t\right)\right)-g\sin \left(\theta \left(t\right)\right)}{R}$

Finding the minimum value of ω for the bead to be in equilibrium gives

$\frac{R{\omega }^2\sin \left(\theta \left(t\right)\right)\cos \left(\theta \left(t\right)\right)-g\sin \left(\theta \left(t\right)\right)}{R}0$

$\cos \left(\theta \right)\frac{g}{R{\omega }^2}$

In order for this to have a solution, we must have

$\omega \ge \sqrt{\frac{g}{R}}$

If θ=π/2, then cos(θ)=0, so ω=∞.

Problem 11.8: K & K 8.12

Problem

A pendulum is rigidly fixed to an axle held by two supports so that it can only swing in a plane perpendicular to the axle. The pendulum consists of a mass m attached to a massless rod of length l. The supports are mounted on a platform which rotates with constant angular velocity Ω. Find the pendulum’s frequency assuming the amplitude is small.

Solution by torque

(From the problem set solutions)

The torque about the pivot point is

${\stackrel{\rightharpoonup }{\tau }}_p\stackrel{\rightharpoonup }{\alpha }{I}_p$

$\hat{k}:-g\ell m\sin \left(\theta \right)+\ell F_{\mathrm{cent}}\cos \left(\theta \right)\ddot{\theta }{I}_p$

The centrifugal effective force is

$F_{\mathrm{cent}}=m\left(\ell \sin \left(\theta \right)\right){\Omega }^2$

For small angles, sin(θ)=θ, cos(θ)=1.  Then equation (1) becomes

If ${\Omega }^2>\frac{g}{\ell }$, the motion is no longer harmonic.

Solution by least action

The general Lagrangian for the object in Cartesian coordinates is

$\mathrm{Clear}\left[x,y,z,t\right];\left(L=\frac{1}{2}m\left({x\text{'}\left[t\right]}^2+{y\text{'}\left[t\right]}^2+{z\text{'}\left[t\right]}^2\right)-mgz\left[t\right]\right)//\mathrm{TraditionalForm}$

$\frac{1}{2}m\left({x^{\prime }\left(t\right)}^2+{y^{\prime }\left(t\right)}^2+{z^{\prime }\left(t\right)}^2\right)-gmz\left(t\right)$

Converting to polar coordinates, and using the constraints that φ=Ω t and r=ℓ, using the conversion

gives

$Lgm\ell \cos \left(\theta \right)-gm\ell -\frac{1}{4}m{\Omega }^2{\ell }^2\cos \left(2\theta \right)+\frac{1}{2}{\stackrel{.}{\theta }}^{2}m{\ell }^2+\frac{1}{4}m{\Omega }^2{\ell }^2$

$0\frac{\partial L}{\partial \theta }-\frac{d}{dt}\frac{\partial L}{\partial \stackrel{.}{\theta }}-gm\ell \sin \left(\theta \right)-m{\ell }^2{\theta }^{\mathrm{\prime \prime }}\left(t\right)+m{\Omega }^2{\ell }^2\sin \left(\theta \right)\cos \left(\theta \right)$

${\theta }^{\mathrm{\prime \prime }}\left(t\right)\frac{{\Omega }^2\ell \sin \left(\theta \left(t\right)\right)\cos \left(\theta \left(t\right)\right)-g\sin \left(\theta \left(t\right)\right)}{\ell }$

Note that this is, after minor changes of variable, the exact same equation that we found in the previous problem.  We should(’ve) expect(ed) this.
Making the first order approximation that θ≈0 (Taylor expanding around θ=0 to the first order), we get

${\theta }^{\mathrm{\prime \prime }}\left(t\right)-\left(\frac{g}{\ell }-{\Omega }^2\right)\theta \left(t\right)$

This is the differential equation for a harmonic oscillator, with

$\omega =\sqrt{\frac{g}{\ell }-{\Omega }^2}$

If ${\Omega }^2>\frac{g}{\ell }$, the motion is no longer harmonic.