summaryrefslogtreecommitdiff
path: root/source/know/concept/hamiltonian-mechanics
diff options
context:
space:
mode:
authorPrefetch2022-10-14 23:25:28 +0200
committerPrefetch2022-10-14 23:25:28 +0200
commit6ce0bb9a8f9fd7d169cbb414a9537d68c5290aae (patch)
treea0abb6b22f77c0e84ed38277d14662412ce14f39 /source/know/concept/hamiltonian-mechanics
Initial commit after migration from Hugo
Diffstat (limited to 'source/know/concept/hamiltonian-mechanics')
-rw-r--r--source/know/concept/hamiltonian-mechanics/index.md308
1 files changed, 308 insertions, 0 deletions
diff --git a/source/know/concept/hamiltonian-mechanics/index.md b/source/know/concept/hamiltonian-mechanics/index.md
new file mode 100644
index 0000000..610f8bd
--- /dev/null
+++ b/source/know/concept/hamiltonian-mechanics/index.md
@@ -0,0 +1,308 @@
+---
+title: "Hamiltonian mechanics"
+date: 2021-07-03
+categories:
+- Physics
+- Classical mechanics
+layout: "concept"
+---
+
+**Hamiltonian mechanics** is an alternative formulation of classical mechanics,
+which equivalent to Newton's laws,
+but often mathematically advantageous.
+It is built on the shoulders of [Lagrangian mechanics](/know/concept/lagrangian-mechanics/),
+which is in turn built on [variational calculus](/know/concept/calculus-of-variations/).
+
+
+## Definitions
+
+In Lagrangian mechanics, use a Lagrangian $L$,
+which depends on position $q(t)$ and velocity $\dot{q}(t)$,
+to define the momentum $p(t)$ as a derived quantity.
+Hamiltonian mechanics switches the roles of $\dot{q}$ and $p$:
+the **Hamiltonian** $H$ is a function of $q$ and $p$,
+and the velocity $\dot{q}$ is derived from it:
+
+$$\begin{aligned}
+ \pdv{L(q, \dot{q})}{\dot{q}} = p
+ \qquad \quad
+ \pdv{H(q, p)}{p} \equiv \dot{q}
+\end{aligned}$$
+
+Conveniently, this switch turns out to be
+[Legendre transformation](/know/concept/legendre-transform/):
+$H$ is the Legendre transform of $L$,
+with $p = \partial L / \partial \dot{q}$ taken as
+the coordinate to replace $\dot{q}$.
+Therefore:
+
+$$\begin{aligned}
+ \boxed{
+ H(q, p) \equiv \dot{q} \: p - L(q, \dot{q})
+ }
+\end{aligned}$$
+
+This almost always works,
+because $L$ is usually a second-order polynomial of $\dot{q}$,
+and thus convex as required for Legendre transformation.
+In the above expression,
+$\dot{q}$ must be rewritten in terms of $p$ and $q$,
+which is trivial, since $p$ is proportional to $\dot{q}$ by definition.
+
+The Hamiltonian $H$ also has a direct physical meaning:
+for a mass $m$, and for $L = T - V$,
+it is straightforward to show that $H$ represents the total energy $T + V$:
+
+$$\begin{aligned}
+ H
+ = \dot{q} \: p - L
+ = m \dot{q}^2 - L
+ = 2 T - (T - V)
+ = T + V
+\end{aligned}$$
+
+Just as Lagrangian mechanics,
+Hamiltonian mechanics scales well for large systems.
+Its definition is generalized as follows to $N$ objects,
+where $p$ is shorthand for $p_1, ..., p_N$:
+
+$$\begin{aligned}
+ \boxed{
+ H(q, p)
+ \equiv \bigg( \sum_{n = 1}^N \dot{q}_n \: p_n \bigg) - L(q, \dot{q})
+ }
+\end{aligned}$$
+
+The positions and momenta $(q, p)$ form a phase space,
+i.e. they fully describe the state.
+
+An extremely useful concept in Hamiltonian mechanics
+is the **Poisson bracket** (PB),
+which is a binary operation on two quantities $A(q, p)$ and $B(q, p)$,
+denoted by $\{A, B\}$:
+
+$$\begin{aligned}
+ \boxed{
+ \{ A, B \}
+ \equiv \sum_{n = 1}^N \Big( \pdv{A}{q_n} \pdv{B}{p_n} - \pdv{A}{p_n} \pdv{B}{q_n} \Big)
+ }
+\end{aligned}$$
+
+
+## Canonical equations
+
+Lagrangian mechanics has a single Euler-Lagrange equation per object,
+yielding $N$ second-order equations of motion in total.
+In contrast, Hamiltonian mechanics has $2 N$ first-order equations of motion,
+known as **Hamilton's canonical equations**:
+
+$$\begin{aligned}
+ \boxed{
+ - \pdv{H}{q_n} = \dot{p}_n
+ \qquad
+ \pdv{H}{p_n} = \dot{q}_n
+ }
+\end{aligned}$$
+
+<div class="accordion">
+<input type="checkbox" id="proof-canoneq"/>
+<label for="proof-canoneq">Proof</label>
+<div class="hidden">
+<label for="proof-canoneq">Proof.</label>
+For the first equation,
+we differentiate $H$ with respect to $q_n$,
+and use the chain rule:
+
+$$\begin{aligned}
+ \pdv{H}{q_n}
+ &= \pdv{}{q_n}\Big( \sum_{j} \dot{q}_j \: p_j - L \Big)
+ \\
+ &= \sum_{j} \bigg( \Big( \dot{q}_j \pdv{p_j}{q_n} + p_j \pdv{\dot{q}_j}{q_n} \Big)
+ - \Big( \pdv{L}{q_n} + \pdv{L}{\dot{q}_j} \pdv{\dot{q}_j}{q_n} \Big) \bigg)
+ \\
+ &= \sum_{j} \Big( p_j \pdv{\dot{q}_j}{q_n} - \pdv{L}{q_n} - p_j \pdv{\dot{q}_j}{q_n} \Big)
+ = - \pdv{L}{q_n}
+\end{aligned}$$
+
+We use the Euler-Lagrange equation here,
+leading to the desired equation:
+
+$$\begin{aligned}
+ - \pdv{L}{q_n} = - \dv{}{t}\Big( \pdv{L}{\dot{q}_n} \Big) = - \dv{p_n}{t} = - \dot{p}_n
+\end{aligned}$$
+
+The second equation is somewhat trivial,
+since $H$ is defined to satisfy it in the first place.
+Nevertheless, we can prove it by brute force,
+using the same approach as above:
+
+$$\begin{aligned}
+ \pdv{H}{p_n}
+ &= \pdv{}{p_n}\Big( \sum_{j} \dot{q}_j \: p_j - L \Big)
+ \\
+ &= \sum_{j} \bigg( \Big( \dot{q}_j \pdv{p_j}{p_n} + p_j \pdv{\dot{q}_j}{p_n} \Big)
+ - \Big( \pdv{L}{q_j} \pdv{q_j}{p_n} + \pdv{L}{\dot{q}_j} \pdv{\dot{q}_j}{p_n} \Big) \bigg)
+ \\
+ &= \dot{q}_n + \sum_{j} \Big( p_j \pdv{\dot{q}_j}{p_n}
+ - 0 \pdv{L}{q_j} - p_j \pdv{\dot{q}_j}{p_n} \Big)
+ = \dot{q}_n
+\end{aligned}$$
+</div>
+</div>
+
+Just like in Lagrangian mechanics, if $H$ does not explicitly contain $q_n$,
+then $q_n$ is called a **cyclic coordinate**, and leads to the conservation of $p_n$:
+
+$$\begin{aligned}
+ \dot{p}_n = - \pdv{H}{q_n} = 0
+ \quad \implies \quad
+ p_n = \mathrm{conserved}
+\end{aligned}$$
+
+Of course, there may be other conserved quantities.
+Generally speaking, the $t$-derivative of an arbitrary quantity $A(q, p, t)$ is as follows,
+where $\ipdv{}{t}$ is a "soft" derivative
+(only affects explicit occurrences of $t$),
+and $\idv{}{t}$ is a "hard" derivative
+(also affects implicit $t$ inside $q$ and $p$):
+
+$$\begin{aligned}
+ \boxed{
+ \dv{A}{t}
+ = \{ A, H \} + \pdv{A}{t}
+ }
+\end{aligned}$$
+
+<div class="accordion">
+<input type="checkbox" id="proof-diff-t"/>
+<label for="proof-diff-t">Proof</label>
+<div class="hidden">
+<label for="proof-diff-t">Proof.</label>
+We differentiate via the multivariate chain rule,
+insert the canonical equations,
+and eventually recognize the PB definition:
+
+$$\begin{aligned}
+ \dv{A}{t}
+ &= \sum_{n} \Big( \pdv{A}{q_n} \pdv{q_n}{t} + \pdv{A}{p_n} \pdv{p_n}{t} \Big) + \pdv{A}{t}
+ \\
+ &= \sum_{n} \Big( \pdv{A}{q_n} \dot{q}_n + \pdv{A}{p_n} \dot{p}_n \Big) + \pdv{A}{t}
+ \\
+ &= \sum_{n} \Big( \pdv{A}{q_n} \pdv{H}{p_n} - \pdv{A}{p_n} \pdv{H}{q_n} \Big) + \pdv{A}{t}
+\end{aligned}$$
+</div>
+</div>
+
+Assuming that $H$ does not explicitly depend on $t$,
+the above property naturally leads us to an alternative
+way of writing Hamilton's canonical equations:
+
+$$\begin{aligned}
+ \dot{q}_n = \{ q_n, H \}
+ \qquad \quad
+ \dot{p}_n = \{ p_n, H \}
+\end{aligned}$$
+
+
+
+## Canonical coordinates
+
+So far, we have assumed that the phase space coordinates $(q, p)$
+are the *positions* and *canonical momenta*, respectively,
+and that led us to Hamilton's canonical equations.
+
+In theory, we could make a transformation of the following general form:
+
+$$\begin{aligned}
+ q \to Q(q, p)
+ \qquad \quad
+ p \to P(q, p)
+\end{aligned}$$
+
+However, most choices of $(Q, P)$ would not preserve Hamilton's equations.
+Any $(Q, P)$ that do keep this form
+are known as **canonical coordinates**,
+and the corresponding transformation is a **canonical transformation**.
+That is, any $(Q, P)$ that satisfy:
+
+$$\begin{aligned}
+ - \pdv{H}{Q_n} = \dot{P}_n
+ \qquad \quad
+ \pdv{H}{P_n} = \dot{Q}_n
+\end{aligned}$$
+
+Then we might as well write $H(q, p)$ as $H(Q, P)$.
+So, which $(Q, P)$ fulfill this?
+It turns out that the following must be satisfied for all $n, j$,
+where $\delta_{nj}$ is the Kronecker delta:
+
+$$\begin{aligned}
+ \boxed{
+ \{ Q_n, Q_j \} = \{ P_n, P_j \} = 0
+ \qquad
+ \{ Q_n, P_j \} = \delta_{nj}
+ }
+\end{aligned}$$
+
+<div class="accordion">
+<input type="checkbox" id="proof-cantrans"/>
+<label for="proof-cantrans">Proof</label>
+<div class="hidden" markdown="1">
+<label for="proof-cantrans">Proof.</label>
+Assuming that $Q_n$, $P_n$ and $H$ do not explicitly depend on $t$,
+we use our expression for the $t$-derivative of an arbitrary quantity,
+and apply the multivariate chain rule to it:
+
+$$\begin{aligned}
+ \dot{Q}_n
+ &= \{Q_n, H\}
+ = \sum_{n} \bigg( \pdv{Q_n}{q_n} \pdv{H}{p_n} - \pdv{Q_n}{p_n} \pdv{H}{q_n} \bigg)
+ \\
+ &= \sum_{n, j} \bigg( \pdv{Q_n}{q_n} \Big( \pdv{H}{Q_j} \pdv{Q_j}{p_n} + \pdv{H}{P_j} \pdv{P_j}{p_n} \Big)
+ - \pdv{Q_n}{p_n} \Big( \pdv{H}{Q_j} \pdv{Q_j}{q_n} + \pdv{H}{P_j} \pdv{P_j}{q_n} \Big) \bigg)
+ \\
+ &= \sum_{n, j} \bigg( \pdv{H}{Q_j} \Big( \pdv{Q_n}{q_n} \pdv{Q_j}{p_n} - \pdv{Q_n}{p_n} \pdv{Q_j}{q_n} \Big)
+ + \pdv{H}{P_j} \Big( \pdv{Q_n}{q_n} \pdv{P_j}{p_n} - \pdv{Q_n}{p_n} \pdv{P_j}{q_n} \Big) \bigg)
+ \\
+ &= \sum_{j} \bigg( \pdv{H}{Q_j} \{Q_n, Q_j\} + \pdv{H}{P_j} \{Q_n, P_j\} \bigg)
+\end{aligned}$$
+
+This is equivalent to Hamilton's equation $\dot{Q}_n = \ipdv{H}{P_n}$
+if and only if $\{Q_n, Q_j\} = 0$ for all $n$ and $j$,
+and if $\{Q_n, P_j\} = \delta_{nj}$.
+
+Next, we do the exact same thing with $P_n$ instead of $Q_n$,
+giving an analogous result:
+
+$$\begin{aligned}
+ \dot{P}_n
+ &= \{P_n, H\}
+ = \sum_{n} \bigg( \pdv{P_n}{q_n} \pdv{H}{p_n} - \pdv{P_n}{p_n} \pdv{H}{q_n} \bigg)
+ \\
+ &= \sum_{n, j} \bigg( \pdv{P_n}{q_n} \Big( \pdv{H}{Q_j} \pdv{Q_j}{p_n} + \pdv{H}{P_j} \pdv{P_j}{p_n} \Big)
+ - \pdv{P_n}{p_n} \Big( \pdv{H}{Q_j} \pdv{Q_j}{q_n} + \pdv{H}{P_j} \pdv{P_j}{q_n} \Big) \bigg)
+ \\
+ &= \sum_{n, j} \bigg( \pdv{H}{Q_j} \Big( \pdv{P_n}{q_n} \pdv{Q_j}{p_n} - \pdv{P_n}{p_n} \pdv{Q_j}{q_n} \Big)
+ + \pdv{H}{P_j} \Big( \pdv{P_n}{q_n} \pdv{P_j}{p_n} - \pdv{P_n}{p_n} \pdv{P_j}{q_n} \Big) \bigg)
+ \\
+ &= \sum_{j} \bigg( \pdv{H}{Q_j} \{P_n, Q_j\} + \pdv{H}{P_j} \{P_n, P_j\} \bigg)
+\end{aligned}$$
+
+Which is equivalent to Hamilton's equation $\dot{P}_n = -\ipdv{H}{Q_n}$
+if and only if $\{P_n, P_j\} = 0$,
+and $\{Q_n, P_j\} = - \delta_{nj}$.
+The PB is anticommutative,
+i.e. $\{A, B\} = - \{B, A\}$.
+</div>
+</div>
+
+If you have experience with quantum mechanics,
+the latter equation should look suspiciously similar
+to the *canonical commutation relation* $[\hat{Q}, \hat{P}] = i \hbar$.
+
+
+
+## References
+1. R. Shankar,
+ *Principles of quantum mechanics*, 2nd edition,
+ Springer.