Physlib

Physlib.ClassicalMechanics.HamiltonsEquations

Hamilton's equations

In this module, given a Hamiltonian function `H : Time → X → X → ℝ`, we define the operator `hamiltonEqOp` which when equals zero implies hamilton's equations.

We show that the variational derivative of the action functional `∫ ⟪p, dq/dt⟫ - H(t, p, q) dt` is equal to the `hamiltonEqOp` applied to `(p, q)`.

References

- G. J. Sussman and J. Wisdom, "Structure and Interpretation of Classical Mechanics", Section 3.1.2. <https://groups.csail.mit.edu/mac/users/gjs/6946/sicm-html/book-Z-H-36.html#%_sec_3.1.2>

4 declarations

definition

Hamilton equation operator hamiltonEqOp(H,p,q)\text{hamiltonEqOp}(H, p, q)

Given a Hamiltonian function H:Time→X→X→RH: \text{Time} \to X \to X \to \mathbb{R}, a momentum trajectory p:Time→Xp: \text{Time} \to X, and a position trajectory q:Time→Xq: \text{Time} \to X, the Hamilton equation operator is defined as the function mapping each time tt to the pair: (dqdt(t)−∇pH(t,p(t),q(t)), −dpdt(t)−∇qH(t,p(t),q(t))) \left( \frac{dq}{dt}(t) - \nabla_p H(t, p(t), q(t)), \ -\frac{dp}{dt}(t) - \nabla_q H(t, p(t), q(t)) \right) In this expression, dqdt\frac{dq}{dt} and dpdt\frac{dp}{dt} denote the time derivatives of the trajectories, ∇pH\nabla_p H denotes the gradient of the Hamiltonian with respect to the momentum (the first XX argument), and ∇qH\nabla_q H denotes the gradient with respect to the position (the second XX argument). The vanishing of this operator is equivalent to the classical Hamilton's equations q˙=∂H∂p\dot{q} = \frac{\partial H}{\partial p} and p˙=−∂H∂q\dot{p} = -\frac{\partial H}{\partial q}.

theorem

Explicit formula for the Hamilton equation operator hamiltonEqOp\text{hamiltonEqOp}

Let H:Time→X→X→RH: \text{Time} \to X \to X \to \mathbb{R} be a Hamiltonian function, p:Time→Xp: \text{Time} \to X be a momentum trajectory, and q:Time→Xq: \text{Time} \to X be a position trajectory. The Hamilton equation operator hamiltonEqOp(H,p,q)\text{hamiltonEqOp}(H, p, q) is the function mapping each time tt to the pair: (dqdt(t)−∇pH(t,p(t),q(t)), −dpdt(t)−∇qH(t,p(t),q(t))) \left( \frac{dq}{dt}(t) - \nabla_p H(t, p(t), q(t)), \ -\frac{dp}{dt}(t) - \nabla_q H(t, p(t), q(t)) \right) where dqdt\frac{dq}{dt} and dpdt\frac{dp}{dt} denote the time derivatives of the trajectories, and ∇pH\nabla_p H and ∇qH\nabla_q H denote the gradients of the Hamiltonian with respect to the momentum (the first XX argument) and the position (the second XX argument), respectively.

theorem

hamiltonEqOp(H,p,q)=0\text{hamiltonEqOp}(H, p, q) = 0 is equivalent to Hamilton's equations

Let H:Time→X→X→RH: \text{Time} \to X \to X \to \mathbb{R} be a Hamiltonian function, where the first XX argument corresponds to momentum and the second to position. For a momentum trajectory p:Time→Xp: \text{Time} \to X and a position trajectory q:Time→Xq: \text{Time} \to X, the Hamilton equation operator hamiltonEqOp(H,p,q)\text{hamiltonEqOp}(H, p, q) vanishes if and only if the trajectories satisfy Hamilton's canonical equations for all tt: dqdt(t)=∇pH(t,p(t),q(t))anddpdt(t)=−∇qH(t,p(t),q(t)) \frac{dq}{dt}(t) = \nabla_p H(t, p(t), q(t)) \quad \text{and} \quad \frac{dp}{dt}(t) = -\nabla_q H(t, p(t), q(t)) where ∇pH\nabla_p H and ∇qH\nabla_q H denote the gradients of the Hamiltonian with respect to its momentum and position arguments, respectively.

theorem

The variational gradient of the action ∫(⟨p,q˙⟩−H)dt\int (\langle p, \dot{q} \rangle - H) dt is the Hamilton equation operator

Let H:Time→X→X→RH : \text{Time} \to X \to X \to \mathbb{R} be a C∞C^\infty smooth Hamiltonian function, where the two XX arguments represent momentum and position respectively. Let pq:Time→X×Xpq : \text{Time} \to X \times X be a C∞C^\infty smooth trajectory in phase space, denoted as pq(t)=(p(t),q(t))pq(t) = (p(t), q(t)). The variational gradient (functional derivative) of the action functional S[p,q]=∫(⟨p(t),q˙(t)⟩−H(t,p(t),q(t)))dt S[p, q] = \int \left( \langle p(t), \dot{q}(t) \rangle - H(t, p(t), q(t)) \right) dt evaluated at the trajectory pqpq is equal to the Hamilton equation operator hamiltonEqOp(H,p,q)\text{hamiltonEqOp}(H, p, q), which maps each time tt to the pair: (q˙(t)−∇pH(t,p(t),q(t)), −p˙(t)−∇qH(t,p(t),q(t))) \left( \dot{q}(t) - \nabla_p H(t, p(t), q(t)), \ -\dot{p}(t) - \nabla_q H(t, p(t), q(t)) \right) where q˙\dot{q} and p˙\dot{p} are time derivatives, and ∇pH\nabla_p H and ∇qH\nabla_q H are the gradients of the Hamiltonian with respect to momentum and position respectively.