Simple Example: Particle on the Surface of a Cylinder

Consider a particle constrained to move on the surface of a cylinder of radius \bgroup\color{black}$ R$\egroup, with a force toward the origin \bgroup\color{black}$ \vec{F}=-k\vec{r}$\egroup. We will work the problem in cylindrical coordinates, \bgroup\color{black}$ (R,\theta,z)$\egroup with \bgroup\color{black}$ R$\egroup fixed. It is rather simple to work this problem from Newton's laws for from the Lagrangian so the point of this exercise is only to see that the Hamiltonian formalism has a lot in common with the methods we know.

Of course the potential energy is

\bgroup\color{black}$\displaystyle U={1\over 2}kr^2={1\over 2}k(R^2+z^2) $\egroup

and the kinetic energy is

\bgroup\color{black}$\displaystyle T={1\over 2}m(\dot{z}^2+R^2\dot{\theta}^2) . $\egroup

Since the potential does not depend on the velocities,

\bgroup\color{black}$\displaystyle H=T+U={1\over 2}m(\dot{z}^2+R^2\dot{\theta}^2)+{1\over 2}k(R^2+z^2) . $\egroup

We need to write this as a function of the coordinates and momenta. So first find the momenta conjugate to \bgroup\color{black}$ z$\egroup and \bgroup\color{black}$ \theta$\egroup.

\bgroup\color{black}$\displaystyle p_i={\partial L\over \partial q_i} = {\partial T\over \partial q_i} $\egroup

since the potential doesn't depend on the velocities.

$\displaystyle p_z$ $\displaystyle = {\partial T\over \partial \dot{z}}=m\dot{z}$    
$\displaystyle p_\theta$ $\displaystyle = {\partial T\over \partial \dot{\theta}}=mR^2\dot{\theta}$    

(Eventually, you should know the momenta conjugate to the standard coordinates.)

Next we write the Hamiltonian in terms of the momenta.

\bgroup\color{black}$\displaystyle H={1\over 2m}\left(p_z^2+{p_\theta^2\over R^2}\right)+{1\over 2}k(R^2+z^2) . $\egroup

Now we look at Hamilton's equations. The equations for \bgroup\color{black}$ \dot{q}_i$\egroup give us the momentum in terms of \bgroup\color{black}$ \dot{q}$\egroup again and so there is nothing new there.

$\displaystyle \dot{z}$ $\displaystyle ={\partial H\over\partial p_z}={p_z\over m}$    
$\displaystyle \dot{\theta}$ $\displaystyle ={\partial H\over\partial p_\theta}={p_\theta\over mR^2}$    

The equations of motion can be found in the equations for \bgroup\color{black}$ \dot{p}_i$\egroup

$\displaystyle \dot{p}_z$ $\displaystyle =-{\partial H\over\partial z } = -kz$    
$\displaystyle \dot{p}_\theta$ $\displaystyle =-{\partial H\over\partial \theta } = 0$    

The \bgroup\color{black}$ \theta$\egroup momentum (angular momentum) is conserved while \bgroup\color{black}$ z$\egroup executes harmonic motion.

$\displaystyle \dot{p}_z$ $\displaystyle = -kz$    
$\displaystyle m\ddot{z}$ $\displaystyle = -kz$    

Jim Branson 2012-10-21