# 8.  A New Way to Write the Action Integral

Michael Fowler

## Introduction

Following Landau, we'll first find how the action integral responds to incremental changes in the endpoint coordinates and times, then use the result to write the action integral itself in a new, more intuitive way.  This new formulation shows very directly the link to quantum mechanics, and variation of the action in this form gives Hamilton's equations immediately.

## Function of Endpoint Position

We’ll now think of varying the action in a slightly different way.  (Note:  We’re using Landau's notation.)  Previously, we considered the integral of the Lagrangian over all possible different paths from the initial place and time ${q}^{\left(1\right)},\text{ }{t}_{1}$ to the final place and time ${q}^{\left(2\right)},\text{ }{t}_{2}$ and found the path of minimum action.  Now, though, we’ll start with that path, the actual physical path, and investigate the corresponding action as a function of the final endpoint variables, given a fixed beginning place and time.

Taking one degree of freedom (the generalization is straightforward), for a small path variation the incremental change in action

$\delta S={\left[\frac{\partial L}{\partial \stackrel{˙}{q}}\delta q\right]}_{{t}_{1}}^{{t}_{2}}+\underset{{t}_{1}}{\overset{{t}_{2}}{\int }}\left(\frac{\partial L}{\partial q}-\frac{d}{dt}\frac{\partial L}{\partial \stackrel{˙}{q}}\right)\delta qdt$

(Recall that first term comes from the calculus of variations when we allow the end point to vary -- it’s exactly the same point we previously discussed in the brachistochrone problem of fastest time for a given horizontal distance, allowing the vertical position of the endpoint to be a free parameter.)

With the incremental variation, we’ve gone from the physical path $P$ (followed by the system in configuration space from ${q}^{\left(1\right)},{t}_{1}$ to ${q}^{\left(2\right)},{t}_{2}$ ) to a second path ${P}^{\prime }$ beginning at the same place and time, and ending at the same time ${t}_{2}$ as $P,$ but at a slightly different place ${q}^{\left(2\right)}+\delta q\left({t}_{2}\right)$.

Both paths $P,{P}^{\prime }$ are fully determined by their initial and final positions and times, so $P,{P}^{\prime }$ must correspond to slightly different initial velocities.  The important point is that since both paths describe the natural dynamical development of the system from the initial conditions, the system obeys the equations of motion at all times along both paths, and therefore the integral term in the above equation is identically zero

Writing $\delta q\left({t}_{2}\right)=\delta q,\text{ }{p}^{\left(2\right)}=p$ the action, regarded as a function of the final position variable, with the final time fixed at ${t}_{2}$, has the differential

$\delta S\left({q}^{\left(2\right)},{t}_{2}\right)={\left[\frac{\partial L}{\partial \stackrel{˙}{q}}\delta q\right]}_{{t}_{1}}^{{t}_{2}}={p}^{\left(2\right)}\delta {q}^{\left(2\right)}=p\delta q.$

For the multidimensional case, the incremental change in the action on varying the final position variable is given by (dropping the superscript)

$\partial S/\partial {q}_{i}={p}_{i},$

## Function of Endpoint Time

What about the action as a function of the final point arrival time

Since $S=\underset{{t}_{1}}{\overset{{t}_{2}}{\int }}Ldt,$ the total time derivative $dS/d{t}_{2}=L\left({q}^{\left(2\right)},{t}_{2}\right),$ the value of the Lagrangian at the endpoint.  Remember we are defining the action at a point as that from integrating along the true path up to that point

Landau denotes ${t}_{2}$ by just $t$, so he writes $dS/dt=L$, and we’ll be doing this, but it’s crucial to keep in mind that the endpoint position and time are the variables here!

If we now allow an incremental time increase, ${t}_{2}\to {t}_{2}+dt$, with the final coordinate position as a free parameter, the  dynamical path will now continue on, to an incrementally different finishing point.

This will give (with $t$ understood from now on to mean ${t}_{2}$, and ${q}_{i}$ means ${q}_{i}^{\left(2\right)}$ )

$\frac{dS\left({q}_{i},t\right)}{dt}=\frac{\partial S}{\partial t}+\sum _{i}\frac{\partial S}{\partial {q}_{i}}{\stackrel{˙}{q}}_{i}=\frac{\partial S}{\partial t}+\sum _{i}{p}_{i}{\stackrel{˙}{q}}_{i}.$

Putting this together with $dS/dt=L$ gives immediately the partial time derivative

$\partial S/\partial t=L-\sum _{i}{p}_{i}{\stackrel{˙}{q}}_{i}=-H$

and therefore, combining this with the result $\partial S/\partial {q}_{i}={p}_{i}$ from the previous section,

$dS\left({q}_{i},t\right)=\sum _{i}{p}_{i}d{q}_{i}-Hdt.$

This, then, is the total differential of the action as a function of the spatial and time coordinates of the end of the path.

## Varying Both Ends

The argument given above for the incremental change in action from varying the endpoint is clearly equally valid for varying the beginning point of the integral (there will be a sign change, of course), so

$dS\left({q}_{i}^{\left(2\right)},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{t}_{2},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{q}_{i}^{\left(1\right)},\text{ }\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{t}_{1}\right)=\sum _{i}{p}_{i}^{\left(2\right)}d{q}_{i}^{\left(2\right)}-{H}^{\left(2\right)}d{t}_{2}-\sum _{i}{p}_{i}^{\left(1\right)}d{q}_{i}^{\left(1\right)}+{H}^{\left(1\right)}d{t}_{1}.$

The initial and final coordinates and times specify the action and the time development of the system uniquely.

(Note: We’ll find this equation again in the section on canonical transformations -- the action will be seen there to be the generating function of the time-development canonical transformation, this will become clear when we get to it.)

## Another Way of Writing the Action Integral

Up to this point, we’ve always written the action as an integral of the Lagrangian with respect to time along the path,

$S\left({q}_{i}^{\left(2\right)},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{t}_{2},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{q}_{i}^{\left(1\right)},\text{ }\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{t}_{1}\right)=\underset{{q}^{\left(1\right)},\text{\hspace{0.17em}}{t}_{1}}{\overset{{q}^{\left(2\right)},\text{\hspace{0.17em}}{t}_{2}}{\int }}Ldt$.

However, the expression derived in the last section for the increment of action generated by an incremental change in the path endpoint is clearly equally valid for the contribution to the action from some interior increment of the path, say from $\left(q,t\right)$ to $\left(q+dq,t+dt\right)$, so we can write the total action integral as the sum of these increments:

$S\left({q}_{i},t\right)=\int dS=\int \left(\sum _{i}{p}_{i}d{q}_{i}-Hdt\right).$

In this integral, of course, the $d{q}_{i}$ add up to cover the whole path.

(In writing $\left({q}_{i},t\right)$ we’re following Landau’s default practice of taking the action as a function of the final endpoint coordinates and time, assuming the beginning point to be fixed. This is almost always fine$—$we’ll make clear when it isn’t.)

## How this Classical Action Relates to Phase in Quantum Mechanics

The link between classical and quantum mechanics is particularly evident in the expression for the action integral given above.  In the so-called semi-classical regime of quantum mechanics, the de Broglie waves oscillate with wavelengths much smaller than typical sizes in the system.  This means that locally it’s an adequate approximation to treat the Schrödinger wave function as a plane wave,

$\psi \left(x,t\right)=A\left(x,t\right){e}^{i\left(kx-\omega t\right)}=A\left(x,t\right){e}^{\left(i/\hslash \right)\left(px-Et\right)}$

where the amplitude function $A\left(x,t\right)$ only varies over distances much greater than the wavelength, and times far longer than the oscillation period.  This expression is valid in almost all the classically accessible regions, invalid in the neighborhood of turning points, but the size of those neighborhoods goes to zero in the classical limit.

As we’ve discussed earlier, in the Dirac-Feynman formulation of quantum mechanics, to find the probability amplitude of a particle propagating from one point to another, we add contributions from all possible paths between the two points, each path contributing a term with phase equal to $i/\hslash$ times the action integral along the path.

From the semi-classical Schrödinger wave function above, it’s clear that the change in phase from a small change in the endpoint is $\left(i/\hslash \right)\left(pdx-Edt\right)$, coinciding exactly with the incremental contribution to the action in  $S=\int dS=\int \left(\sum _{i}{p}_{i}d{q}_{i}-Hdt\right).$

So again we see, here very directly, how the action along a classical path is a multiple of the quantum mechanical phase change along the path.

## Hamilton’s Equations from Action Minimization

For arbitrary small path variations $\delta q,\delta p$ in phase space, the minimum action condition using the form of action given above generates Hamilton’s equations.

(Note for nitpickers: This may seem a bit surprising, since we generated this form of the action using the equations along the actual dynamical path, how can we vary it and still use them? Bear with me, you’ll see.)

We’ll prove this for a one dimensional system, it’s trivial to go to many variables, but it clutters up the equations.

For a small path deviation $\delta q,\delta p$, the change in the action $S=\int \left(pdq-Hdt\right)$ is

$\delta S=\int \left[\delta pdq+pd\left(\delta q\right)-\left(\partial H/\partial q\right)\delta qdt-\left(\partial H/\partial p\right)\delta pdt\right]=0$

and integrating $pd\left(\delta q\right)$ by parts, with $\delta p=\delta q=0$ at the endpoints,

$\delta S=\int \delta p\left\{dq-\left(\partial H/\partial p\right)dt\right\}+\left[p\delta q\right]-\int \delta q\left\{dp+\left(\partial H/\partial q\right)dt\right\}=0.$

The path variations $\delta p,\delta q$ are independent and arbitrary, so must have identically zero coefficients$—$Hamilton’s equations follow immediately, $\stackrel{˙}{q}=\partial H/\partial p,$ $\stackrel{˙}{p}=-\partial H/\partial q.$

Again, it’s worth emphasizing the close parallel with quantum mechanics: Hamilton’s equations written using Poisson brackets are:

$\stackrel{˙}{q}=\left[H,q\right],\text{ }\stackrel{˙}{p}=\left[H,p\right].$

In quantum mechanics, the corresponding Heisenberg equations of motion for position and momentum operators in terms of commutators are

$\stackrel{˙}{q}=\left(1/i\hslash \right)\left[H,q\right],\text{ }\stackrel{˙}{p}=\left(1/i\hslash \right)\left[H,p\right].$

## How Can p, q Really Be Independent Variables?

It may seem a little odd at first that varying $p,q$ as independent variables leads to the same equations as the Lagrangian minimization, where we only varied $q,$ and that variation “locked in” the variation of $\stackrel{˙}{q}.$ And, isn’t $p$ defined in terms of $q,\stackrel{˙}{q}$ by $p=\partial L/\partial \stackrel{˙}{q},$ which is some function of $q,\stackrel{˙}{q}$?  So wouldn’t varying $q$ automatically determine the variation of $p$?

The answer is, no, $p$ is not defined as $p=\partial L/\partial \stackrel{˙}{q}$ from the start in Hamilton’s formulation.  In this Hamiltonian approach, $p,q$ really are taken as independent variables, then varying them to find the minimum path gives the equations of motion, including the relation between $p$ and $q,\stackrel{˙}{q}$

This comes about as follows: Along the minimum action path, we just established that

$dH\left(p,q\right)=\stackrel{˙}{q}dp-\stackrel{˙}{p}dq.$

We also have that $L=p\stackrel{˙}{q}-H,$ so (Legendre transformation!)

$dL\left(q,\stackrel{˙}{q}\right)=pd\stackrel{˙}{q}+\stackrel{˙}{p}dq,$

from which, along the physical path,   So this identity, previously written as the definition of $p,$ now arises as a consequence of the action minimization in phase space.