This article is part of a set of four; the common factor is Calculus of Variations. In classical physics Calculus of Variations is applied in three areas: Optics, Statics, and Dynamics. Each article in the set is written as a standalone article, resulting in some degree of overlap.
The other three articles:

Foundation:	Calculus of Variations
Optics:	Fermat's stationary time
Statics:	The catenary

The Energy-Position equation

1. From F = ma to the Work-Energy theorem

2. From Work-Energy to stationary action

3. Trial trajectories: response to variation
    3.1 Potential increases linear with height
    3.2 Potential increases quadratic with displacement
    3.3 Potential increases proportional to cube of displacement

Appendix I Verification

Appendix II The meaning of stationary

Appendix III Jacob's Lemma

This exposition is about Hamilton's stationary action. The usual presentation is to posit Hamilton's stationary action, and to proceed with showing that F = ma can be recovered from it. However, it is also possible to proceed in the other direction. You can start with F = ma, and proceed to Hamilton's stationary action in all forward steps. The development is in two stages:
- Derivation of the Work-Energy theorem from F = ma
- Demonstration that in cases where the Work-Energy theorem holds good Hamilton's stationary action will hold good also.

The Work-Energy theorem and Hamilton's stationary action have the following in common: the physics taking place is described in terms of kinetic energy and potential energy. As an overarching name for both the Work-Energy theorem and Hamilton's stationary action I will use the expression 'energy mechanics'.

Interactive diagrams
The intention of this article is to let the interactive diagrams tell the story. You can choose to jump ahead to the interactive diagrams However, I strongly recommend you return to these paragraphs at some point; for overall understanding they are necessary.

Requirement for well defined potential energy

Given a particular force: in order to have a well defined expression for potential energy the force must have the property that the work done in moving from some point A to some point B is independent of the path taken between the two points. The validity of energy mechanics is limited to forces with that property.

(It may be that in terms of a deeper theory such as quantum physics all interactions actually have that potential-is-independent-of-the-path property, but we should not blindly assume that.)

As we know: the principle of conservation of energy is asserted without any restriction. That's why the concept is referred as a 'principle', it is a blanket statement.

The discussion in this article avoids the blanket statement. The discussion in this article is limited to the classes of cases where by way of experimental corroboration it is evident that the work done in moving from some point A to some point B is independent of the path taken between the two points.

1. From F = ma to the Work-Energy theorem

The derivation starts with F=ma, and evaluating the integral from position coordinate s₀ to position coordinate s:

$\int_{s_0}^s F \ ds = \int_{s_0}^s ma \ ds$

(1.1)

In the steps (1.2) to (1.8) the right hand side of (1.1) is evaluated.

For the time being I will not write the factor m, it is a multiplicative factor that is just carried from step to step. In the final expression I will include the factor m again.

As we know, acceleration and position are related according to (1.2) and (1.3).

$ds = v \ dt$

(1.2)

$dv = a \ dt$

(1.3)

The integral for acceleration from a starting point s₀ to a final point s

$\int_{s_0}^s a \ ds$

(1.4)

Intermediary step: change of the differential according to (1.2), with corresponding change of limits.

$\int_{t_0}^t a \ v \ dt$

(1.5)

Change the order:

$\int_{t_0}^t v \ a \ dt$

(1.6)

Change of the differential according to (1.3), with corresponding change of limits.

$\int_{v_0}^v v \ dv$

(1.7)

So we have:

$\int_{s_0}^s a \ ds = \tfrac{1}{2}v^2 - \tfrac{1}{2}v_0^2$

(1.8)

The elements that give rise to the relation (1.8) are: (1.2) and (1.3), and the mathematics of differentiation/integration. That is: no physics is referred to, hence (1.8) is independent from any physics.

After multiplying with m we use the right hand side of (1.8) for the right hand side of (1.1), thus arriving at the Work-Energy theorem:

The Work-Energy theorem

To obtain the work-energy theorem: set up the integral of force with respect to the position coordinate, substitute the force F with ma, and then restate according to Torricelli's formula:

$\int_{s_0}^s F \ ds = \int_{s_0}^s ma \ ds = \tfrac{1}{2}mv^2 - \tfrac{1}{2}mv_0^2$

$\int_{s_0}^s F \ ds = \tfrac{1}{2}mv^2 - \tfrac{1}{2}mv_0^2$

(1.9)

(About the quantity ½mv² on the right hand side of (1.9): we have that earlier in the history of theory of motion a quantity of the form mv² was recognized in another context: collisions. In elastic collision a quantity proportional to mv² is conserved. This quantity-proportional-to-mv² was referred to as 'the living force'. As we know: today a single concept is used, kinetic energy, with the value supplied by the Work-Energy theorem: ½mv².)

To emphasize the structure of the Work-Energy theorem I give the elements arranged in a table. In each row the statement follows mathematically from the statements in the row(s) above it.

	$v = \frac{ds}{dt}$	$a = \frac{dv}{dt}$
	$\int_{s_0}^s a \ ds = \tfrac{1}{2}v^2 - \tfrac{1}{2}v_0^2$
$\int_{s_0}^s F \ ds = \tfrac{1}{2}mv^2 - \tfrac{1}{2}mv_0^2$

Potential energy

Historically the concept of potential energy arose before the work-energy theorem reached today's standard form. With the benefit of hindsight we recognize that the concept of potential energy capitalizes on the work-energy theorem.

Potential energy is defined as the negative of the left hand side of (1.9); the integral of force over distance.

$\Delta E_p = - \int_{s_0}^s F \ ds$

(1.10)

In order for the potential energy to be well defined a specific condition must be met: the value of the integral from point s₀ to point s must be independent of how the test object moves from point s₀ to point s.

When the potential energy is well defined:

$\Delta (E_k) = - \Delta (E_p)$

(1.11)

That is: with potential energy defined according to (1.10): when an object is moving down a potential gradient the kinetic energy will increase by the amount that the potential energy decreases. (With everything reversed when moving in the opposite direction, of course.)

(1.12) combines three statements. I present these three statements as a unit to emphasize the interconnection between them. While the statements are different mathematically, the physics content of these three is the same.

$\begin{array}{rcl} F & = & ma \\[+10pt] \int_{s_0}^s F \ ds & = & \tfrac{1}{2}mv^2 - \tfrac{1}{2}mv_0^2 \\[+10pt] -\Delta E_p & = & \Delta E_k \end{array}$

(1.12)

2. From Work-Energy to stationary action

Visual demonstration

graphlet 2.1 offers a visual representation of the connection between the work-energy theorem and Hamilton's stationary action

In the section following this section the connection will be given in equation form.

Picture 2.1 Graphlet

The upturned parabola in the upper left sub-panel represents height as a function of time, of an object that is thrown vertically upward. The object ascends, decelerated by a uniform force, and is then accelerated back down.

The slider at the bottom of the diagram represents sweeping out variation. When the variation is at the value zero the trial trajectory coincides with the true trajectory.

Note especially that variation of a trial trajectory is exclusively variation of position coordinate. As a consequence: differentiation with respect to variation can always be re-expressed as differentiation with respect to position coordinate.

The following three relations are plotted:
upper-left subpanel: height of the trial trajectory as a function of time
upper-right and lower-left: Kinetic energy (red) and Work done (green) as a function of time
lower-right: the value of the integral-with-respect-to-time as a function of variational parameter p_v.

Use the slider to sweep out variation.
As you move the slider:
· Upper-left: The height of the trial trajectory changes linear with the slider position.
· Upper-right and lower-left: the shape of the curve for each energy changes; for the respective energies the rate of change is different.
· Lower-right:
·Green curve: the value changes linear with the variation.
·Red curve: the value changes quadratically.

Now focus on the lower-right subpanel, which displays the value of the integral as a function of the variational parameter p_v. The moving dot represents the integral value that corresponds to the current p_v value. Examine the point in variation space where the trial trajectory coincides with the true trajectory: at that point the rate of change of the red integral value matches the rate of change of the green integral value.

The lower-right subpanel has a checkbox labeled 'sum'. With that checkbox checked the sum of the two integrals is displayed, in the color blue. At the point in variation space where the trial trajectory coincides with the true trajectory: at that point the derivative of the blue curve is zero.

The usual notation for Hamilton's action is to represent it with the letter S. In the graphlet the red integral has the label S_k for 'kinetic energy component of the action S'. Potential energy is defined as the negative of work done, hence the green integral has the label '-S_k' for 'potential energy component of the action S'. The blue curve represents Hamilton's action:

(2.1)

Derivative with respect to position.

We pick things up from equation (1.12) and we build onto it result by making use of the fundamental theorem of Calculus: differentiation and integration are each other's inverse.

In section 1. the derivation of the work-energy theorem from F=ma was by way of integration with respect to the position coordinate. The following equation specifies the inverse: differentiation with respect to the position coordinate

$\frac{\int_{s_0}^s F \ ds}{ds} = \frac{\tfrac{1}{2}mv^2}{ds}$

(2.2)

Given that differentiation is the inverse of integration we know in advance that processing this differentiation will recover F=ma.

In this section the property that is demonstrated visually in graphlet 2.1 is demonstrated in equation form.

The following property is key:
Differentiation and integration are linear operations.
As we know: that means we can move additions and subtractions in and out.

$\frac{d \ f(x)}{dx} + \frac{d \ g(x)}{dx} = \frac{d \big( \ f(x) + g(x) \ \big)}{dx}$

$\int f(x)dx + \int g(x)dx = \int \big( \ f(x) + g(x) \ \big) dx$

About the two equations below: (2.3) is a repeat of (2.2)

Given that integration and differentiation are linear operations: the property expressed in (2.3) propagates to the form expressed in (2.4).

$\frac{\int_{s_0}^s F \ ds}{ds} = \frac{\tfrac{1}{2}mv^2}{ds}$

(2.3)

$\frac{d \Big( \int_{t_1}^{t_2} \ \big(\int F \ ds \big) \ dt \Big)}{ds} = \frac{d \Big( \int_{t_1}^{t_2} \ \big(\tfrac{1}{2}mv^2 \big) \ dt \Big)}{ds}$

(2.4)

In the above the forms ∫F ds and ½mv² have been retained to emphasize the continuity with the work-energy theorem.

About the right hand side of (2.3), the differentiation of kinetic energy with respect to position coordinate. That differentiation can be restated in a way that makes it easier to process. The following two differentiation operations are equivalent:

$\frac{\tfrac{1}{2}mv^2}{ds} \ \Leftrightarrow \ \frac{d}{dt}\frac{d (\tfrac{1}{2}mv^2)}{d v}$

(2.5)

In the following expression is a restatement of (2.3), with the right hand side restated according to (2.5)

$\frac{d(-E_p)}{ds} = \frac{d}{dt}\frac{d (E_k)}{d v}$

(2.6)

Rearrange to have both terms on the same side:

$\frac{d(-E_p)}{ds} - \frac{d}{dt}\frac{d (E_k)}{d v} = 0$

(2.7)

The following notation says the same thing as (2.7) does, since ∂E_k/∂s=0, and ∂E_k/∂v=0.

$\frac{\partial(E_k - E_p)}{\partial s} - \frac{d}{dt}\frac{\partial (E_k - E_p)}{\partial v} = 0$

(2.8)

Of course: with: (2.8) is the equation that we get when the Lagrangian (E_k-E_p) is inserted in the Euler-Lagrange operator

3. Trial trajectories: response to variation

The following three graphlets are three instances of the same graphlet, for three successive classes of cases.

As you manipulate the sliders: the process to home in on the true trajectory is to manipulate the trial worldline such that over the entire trajectory the slope of the kinetic energy curve (red) matches the slope of the potential energy curve (green). When the trial trajectory is at thet the point in variation space such that those two slopes match the derivative of Hamilton's action is zero.

3.1 Potential increases linear with height

Picture 3.1. Graphlet
Potential energy increases linear with height; the trajectory is a parabola

Graphlet 3.1 is for the case of a uniform force, causing an acceleration of 2 m/s².
As we know: with a uniform force the curve that represents the height as a function of time is a parabola.

In the 'energy' sub-panel the green curve represents the minus potential energy. In effect the potential energy curve has been flipped upside down. That way you can see directly that when the trial trajectory coincides with the true trajectory the red and green curve are parallel to each other everywhere.

(Any form of calculating the true trajectory uses the derivative of the potential energy. That is, in calculation the value of the potential energy itself is not used. Because of that the choice of zero point of potential energy is arbitrary.)

In the lower right subpanel:
The blue dot represents Hamilton's action.
The label of the horizontal axis is p_v, which stands for 'variational parameter'. The positioning of the dots corresponds to the value of the main slider.

In the case of linear potential: when the trial trajectory coincides with the true trajectory the value of Hamilton's action is minimal/minimized. In the case implemented in this graphlet: any change of the trial trajectory results in raising the value of Hamilton's action.

3.2 Potential increases quadratic with displacement

Picture 3.2. Graphlet
Potential energy increases quadratic with displacement

Graphlet 3.2 is for the case where the force increases in linear proportion to the displacement. This case is commonly referred to as Hooke's law.
As we know: With Hooke's law the resulting motion is harmonic oscillation and the curve that represents the displacement as a function of time is the sine function.

With Hooke's law the potential energy increases with the square of the displacement. That is: with Hooke's law we have that the rate of change of potential energy is on par with that of the kinetic energy: with Hooke's law the expressions for the energies are both quadratic expressions.

As we know: with Hooke's law the amplitude and period of the resulting oscillation are independent from each other. Hamilton's stationary action corroborates that: when the trial trajectory is the harmonic oscillation function then if you change the amplitude the value of the action remains the same.

(In the diagram the actual curve is the cosine function; the point is that it is the harmonic oscillation function.)

3.3 Potential increases proportional to cube of displacement

Picture 3.3. Graphlet
Potential energy proportional to cube of displacement.

Graphlet 3.3 is for the case where the force increases proportional to the square of the displacement, hence the potential energy increases with the cube.

Click the button 'Show numerical' to show how the curve displayed in graphlet 6 was independently verified

Comparing the trajectories

In graphlet 3.1 the potential energy as a function of position is linear. Whenever the potential energy as a function of position is of lower order than the expression for the kinetic energy the action reaches a minimum when the derivative of the action is zero. In graphlet 3.2 the potential energy function and the kinetic energy function are both quadratic expressions, so the action is neither a minimum nor a maximum. Here in graphlet 3.3 the potential energy as a function of position is a higher order expression than the one for the kinetic energy (cubic versus quadratic), and consequently when the derivative of the action is zero the action reaches a maximum.

I recommend enlarging the graphlet to the full width of the browser window. Visibility of the navigation column of this page can be toggled. When the navigation column has been hidden: use the button 'Larger' and zoom in on the page as a whole to make the graphlet fill the entire width of the browser window.

How to operate the demonstration:

The graphlet set 3.1, 3.2, and 3.3 is designed to be discoverable, but for the sake of completeness I provide a description.

The main slider, located at the bottom of the graphlet, executes a global variation sweep. The value displayed in the "knob" of the main slider is a value to implement the variation, from here on I will refer to it as the 'variational parameter p_v' In the 'integral' panel the variational parameter p_v is along the horizontal axis.

Names for the 4 sub-panels:
- upper left: Control panel
- lower left: Height panel
- upper right: Energies panel
- lower right: Integrals panel

Note: The 3 sub-panels with a coordinate system are named after their vertical axis name: Height, Energy, Integral.

In the starting configuration the trajectory points (height panel) have been placed such that they coincide (to a very good approximation) with the true trajectory of the object.

When the object is moving along the true trajectory the rate of change of kinetic energy matches the rate of change of potential energy (work-energy theorem). This match of rate of change has been emphasized by turning the graph of the potential energy (green) upside down. When the object is moving along the true trajectory the red and green curve have the same slope along the entire trajectory.

The energies panel is where it happens. The energies panel represents how the equation that you are using solves for the true trajectory.

Hamilton's Action

$S = \int (E_k - E_p) dt$

Hamilton's Action is represented by the blue dot in the Integrals panel. Use the main slider to sweep out variation. The blue dot follows a curve; it is a curve in variation space. The slope of that curve represents the derivative of Hamilton's Action with respect to variation.

At the point where over the entire trajectory the derivative of the kinetic energy is equal to the derivative of the minus potential energy the derivative of the blue dot is zero.

In the control panel the two sliders on the far left are adjusters for the trajectory. The upper adjuster morphs the trajectory towards a curve that is more blunt than the true trajectory; the lower adjuster morphs the trajectory towards a triangle shape.

These adjusters allow morphing of the trial trajectory while maintaining that during the entire ascent the velocity never reverses from decreasing to increasing again, which would be unphysical.

The row of ten sliders is for adjustment of individual nodes. The three radio buttons toggle between three sets of node adjusting sliders.
The '× 1' button: for a ratio of 1-to-1 of moving the slider and movement of the corresponding node of the trajectory.
The '× 0.1' button: ratio of 10 to 1
The '× 0.01' button: ratio of 100 to 1, for fine adjustment.

The button 'Consolidate':
Clicking the button 'Consolidate' does the following: the current value of the '× 0.1' slider is transferred to the '× 1' slider. That is: after clicking the button 'Consolidate' the position of the '× 0.1' slider is reset to its zero position, and the '× 1' slider has been incremented accordingly.

In the energies panel: the grey point is draggable, and it drags the entire curve of the kinetic energy (red) with it. By shifting the kinetic energy curve over to the potential energy curve the user can verify that the red and green curve are parallel to each other along the entire trajectory. (The evaluation of the integral uses the unshifted position of the kinetic energy curve.)

In the integrals panel:
· Red/Green point: value of the integral of the red/green curve of the energies panel
· Blue point: the sum of the values of the red point and the green point

Finding the true trajectory by adjusting the sliders

As a seed move the '× 1' sliders so that the nodes are positioned along a straight line, say at 45 degrees.
In the 'Energies' subpanel, move the graph of the kinetic energy (red dots) down, on average a bit below the potential energy graph (green dots)
Use the '× 0.1' sliders to move the nodes upward bringing the red dots and green dots into alignment.
The true trajectory has the property that the rate of change of kinetic energy matches the rate of change of potential energy, so you want to make those two curves parallel to each other.

Picture I.1. Graphlet
Infinitesimal variation δs

Appendix I: Verification

The Work-Energy theorem is constructed by stating the integral with respect to position of F=ma

Hence taking the derivative with respect to position will recover F=ma.

Simplify to initial position coordinate zero and initial velocity zero:

$\int F \ ds = \tfrac{1}{2}mv^2$

(I.12)

the derivative with respect to position:

$\frac{\int F \ ds}{ds} = \frac{\tfrac{1}{2}mv^2}{ds} \qquad \Longleftrightarrow \qquad F = ma$

(I.13)

In the Euler-Lagrange equation, (I.10) and (I.11), it looks as if the way the kinetic energy term is processed is different from how the kinetic energy term is processed in the Energy-Position equation, but that is not actually the case.

In classical mechanics the Euler-Lagrange equation takes the derivative of the energy with respect to the position coordinate. Proof: (I.13) proceeds with the differentiation of kinetic energy specified by the Euler-Lagrange equation:

$\frac{d}{dt} \frac{dE_k}{dv} = \frac{d}{dt}mv = ma$

(I.14)

(I.14) evaluates to ma, showing that in classical mechanics the Energy-Position equation (1.14), and the Euler-Lagrange equation (I.10)/(I.11) are the same equation.

Appendix II: The meaning of 'stationary'

Picture II.1. Graphlet
At the x-coordinate where the red and green curve have equal slope the blue curve is at zero slope.

Red and green are both ascending functions. We want to identify the coordinate where the rate of change of the red value matches the rate of change of the green value. The shortest way to get there is by setting up an equation with the derivative of the red curve on one side and the derivative of the green curve on the other side.

$\frac{d\left( \ f(x) \ \right)}{dx} = \frac{d \left( \ g(x) \ \right)}{dx}$

(II.1)

$\frac{d(\tfrac{1}{2}x^2)}{dx} = \frac{d(ln(x)-\tfrac{1}{2})}{dx}$

(II.2)

$x = \frac{1}{x}$

(II.3)

In the case of Hamilton's stationary action the same comparison is performed, but in a form that takes more steps. We see the following: a third function is defined, which is called the Lagrangian L: L=(E_k-E_p), and the point to be identified is called the point of 'stationary action'. 'Stationary action' is another way of saying: identify the point where the derivative of the Lagrangian is zero.

Summerizing:
Hamilton's stationary action works by imposing the constraint that was established with (1.11): the rate of change of kinetic energy must match the rate of change of potential energy.

Back to the paragraph 'stationary' in the article.

Appendix III: Jacob's Lemma

When Johann Bernoulli had presented the Brachistochrone problem to the mathematicians of the time Jacob Bernoulli was among the few who was able to find the solution independently. The treatment by Jacob Bernoulli is in the Acta Eruditorum, May 1697, pp. 211-217

Jacob opens his treatment with an observation concerning the fact that the curve that is sought is a minimum curve.

Picture III.1. Image
Jacob's Lemma

Image source:
Acta Eruditorum 1697 p. 205

Lemma. Let ACEDB be the desired curve along which a heavy point falls from A to B in the shortest time, and let C and D be two points on it as close together as we like. Then the segment of arc CED is among all segments of arc with C and D as end points the segment that a heavy point falling from A traverses in the shortest time. Indeed, if another segment of arc CFD were traversed in a shorter time, then the point would move along ACFDB in a shorter time than along ACEDB, which is contrary to our supposition.

Jacob's lemma generalizes to all cases where the curve that you want to find is an extremum; either a maximum or a minimum. If the evaluation is an extremum for the entire curve, then it is also an extremum for any sub-section of the curve, down to infinitesimally short subsections.

The interactive diagrams on this page are created with the Javascript Library JSXGraph. JSXGraph is developed at the Lehrstuhl für Mathematik und ihre Didaktik, University of Bayreuth, Germany.

This work is licensed under a Creative Commons Attribution-ShareAlike 3.0 Unported License.

Last time this page was modified: September 14 2024