記事ソース/平行移動、回転、時間発展とポアソン括弧
メニュー現在 9 名がオンラインです。 最新の25件2023-12-12
2023-11-11
2023-11-06
2023-07-15
2022-09-14
2022-07-01
2022-06-12
2022-04-13
2021-12-03
2021-10-07
2021-08-12
2021-07-26
2021-06-30
2021-06-06
2021-05-02
2021-04-17
2021-03-20
2021-03-19
2021-03-11
|
記事ソース/平行移動、回転、時間発展とポアソン括弧†これはrst2hooktailの記事ソース保存・変換用です（詳細）． コンバート公開・更新メニュー ▼▲記事ソースの内容============================================================
平行移動、回転、時間発展とポアソン括弧
============================================================
ハミルトニアンは時間発展、運動量は平行移動、角運動量は回転と関係していると聞いたことはありませんか？
僕が知っていたのは、ネーターの定理あたりの話からだったのですが、
ポアソン括弧を使ってもっと直接的関係が得られることに気づきました。
量子力学では有名な話です。ただ、私はその古典力学での対応を知りませんでした。
この記事は証明と言うよりは、検証・確認です。
ポアソン括弧
==============================
ポアソン括弧とは、運動量 $p$ 、位置 $q$ 、時間 $t$ と正準変数の関数 $A \equiv A(q,p),B \equiv B(q,p)$ を用いた、
<tex>
[A,B] = \dfrac{\partial A}{\partial q}\dfrac{\partial B}{\partial p} -\dfrac{\partial A}{\partial p}\dfrac{\partial B}{\partial q} \tag{##}
</tex>
と言う計算の事です。
時間発展とハミルトニアン
==============================
これは有名だと思います。関数 $F(q,p,t)$ の時間発展は、ハミルトニアン $H$ を用いて、
<tex>
\dfrac{d}{dt}F(q,p,t) = [F,H] + \dfrac{\partial F}{\partial t} \tag{##}
</tex>
という方程式になります。特に $F$ が時間を陽に含まない時、
<tex>
\dfrac{d}{dt}F(q,p) = [F,H] \tag{##}
</tex>
となります。試しに
<tex>
F=q,H=\dfrac{p^2}{2m}+\dfrac{kq^2}{2} \tag{##}
</tex>
として式を計算すると、
<tex>
\dfrac{dq}{dt} &= [q,\dfrac{p^2}{2m}+\dfrac{kq^2}{2}] \\
&= \dfrac{\partial q}{\partial q}\dfrac{\partial \left( \dfrac{p^2}{2m}+\dfrac{kq^2}{2} \right)}{\partial p} -\dfrac{\partial q}{\partial p}\dfrac{\partial \left( \dfrac{p^2}{2m}+\dfrac{kq^2}{2} \right)}{\partial q} \\
&= p/m
\tag{##}
</tex>
となり、確かに正しい運動方程式です。 $H=\dfrac{p^2}{2m}+\dfrac{kq^2}{2}$ としたハミルトン方程式
<tex>
\dfrac{dq}{dt} &= \dfrac{\partial H}{\partial p} \\
&= p/m
</tex>
と等価です。
平行移動と運動量
=================================
式 $(3)$ に対応させて以下の式を考えてみました。
<tex>
\dfrac{d}{dq}F(q,p) = [F(q,p),p] \tag{##}
</tex>
右辺を計算すると、
<tex>
[F(q,p),p] = \dfrac{\partial F}{\partial q}\dfrac{\partial p}{\partial p} -\dfrac{\partial F}{\partial p}\dfrac{\partial p}{\partial q}
\tag{##}
</tex>
特に $\dfrac{k}{2}q^2$ としてみると、左辺が、
<tex>
\dfrac{d}{dq} \left( \dfrac{k}{2}q^2 \right) = kq \tag{##}
</tex>
であり、右辺が
<tex>
[\dfrac{k}{2}q^2,p] &= \dfrac{\partial \left( \dfrac{k}{2}q^2 \right) }{\partial q}\dfrac{\partial p}{\partial p} -\dfrac{\partial \left( \dfrac{k}{2}q^2 \right)}{\partial p}\dfrac{\partial p}{\partial q} \\
&= kq
</tex>
となり、一致していますね。
回転と角運動量
=============================
最後に二次元平面での回転を見ましょう。
<tex>
\bm{L} = \bm{q} \times \bm{p} \tag{##}
</tex>
ですから、二次元に制限すると、 $L_3=q_1p_2-q_2p_1$ を考えれば良く、
<tex>
\dfrac{d}{d \theta}F(q_1,q_2,p_1,p_2) &= \sum_{i} [F,L_3] \\
&= \sum_{i} \left( \dfrac{\partial F}{\partial q_i}\dfrac{\partial L_3}{\partial p_i} - \dfrac{\partial F}{\partial p_i}\dfrac{\partial L_3}{\partial q_i} \right)
\tag{##}
</tex>
です。 $\theta$ は回転角です。例えば、
<tex>
[F(q_1,q_2,p_1,p_2),L_3] &= \dfrac{\partial F}{\partial q_1}\dfrac{\partial L_3}{\partial p_1} -\dfrac{\partial F}{\partial p_1}\dfrac{\partial L_3}{\partial q_1} + \dfrac{\partial F}{\partial q_2}\dfrac{\partial L_3}{\partial p_2} -\dfrac{\partial F}{\partial p_2}\dfrac{\partial L_3}{\partial q_2} \\
&= \dfrac{\partial F}{\partial q_1} \cdot (-q_2) - \dfrac{\partial F}{\partial p_1} \cdot p_2 + \dfrac{\partial F}{\partial q_2} \cdot q_1 - \dfrac{\partial F}{\partial p_2} \cdot (-p_1)
\tag{##}
</tex>
となります。ここで、 $F=q_1 \mathrm{or} F=q_2$ とすると、
<tex>
\dfrac{d q_1}{d \theta} &= [q_1,L_3] = -q_2 \\
\dfrac{d q_2}{d \theta} &= [q_2,L_3] = q_1 \tag{##}
</tex>
はて、これはどうしたらいいでしょう？
行列を使って書き直してみましょう。
<tex>
\dfrac{d}{d \theta} \begin{pmatrix} q_1 \\ q_2 \end{pmatrix} = \begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix} \begin{pmatrix} q_1 \\ q_2 \end{pmatrix} \tag{##}
</tex>
となり、行列の指数関数を使えば良さそうです。 $J = \begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix}$ として、
解くと、
<tex>
\begin{pmatrix} q_1(\theta) \\ q_2(\theta) \end{pmatrix} &= \exp (J \theta) \begin{pmatrix} q_1(0) \\ q_2(0) \end{pmatrix} \\
&= \begin{pmatrix} \cos \theta & - \sin \theta \\ \sin \theta & \cos \theta \end{pmatrix} \begin{pmatrix} q_1(0) \\ q_2(0) \end{pmatrix}
\tag{##}
</tex>
となりました。これは正に（能動的）回転の行列ですね。
それでは、今日はこの辺で。お疲れ様でした。
@@author:クロメル@@
@@accept:2018-04-01@@
@@category:解析力学@@
@@id:poissonAndOperators@@
|
