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:TimeXXRH: \text{Time} \to X \to X \to \mathbb{R}, a momentum trajectory p:TimeXp: \text{Time} \to X, and a position trajectory q:TimeXq: \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˙=Hp\dot{q} = \frac{\partial H}{\partial p} and p˙=Hq\dot{p} = -\frac{\partial H}{\partial q}.

theorem

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

Let H:TimeXXRH: \text{Time} \to X \to X \to \mathbb{R} be a Hamiltonian function, p:TimeXp: \text{Time} \to X be a momentum trajectory, and q:TimeXq: \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:TimeXXRH: \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:TimeXp: \text{Time} \to X and a position trajectory q:TimeXq: \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:TimeXXRH : \text{Time} \to X \to X \to \mathbb{R} be a CC^\infty smooth Hamiltonian function, where the two XX arguments represent momentum and position respectively. Let pq:TimeX×Xpq : \text{Time} \to X \times X be a CC^\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.