物理のかぎしっぽ 記事ソース/共役な物理量とポアソン括弧 のバックアップ差分(No.4)

#rst2hooktail_source
 ============================================================
 共役な物理量とポアソン括弧
 ============================================================
 
 この記事では、通常
 
 <tex>
 \left[ A , B \right] = \dfrac{\partial A}{\partial q} \dfrac{\partial B}{\partial p} - \dfrac{\partial A}{\partial p} \dfrac{\partial B}{\partial q}
 \tag{##}
 </tex>
 
 で定義されるポアソン括弧について考えることで、共役な変数とは何か?を考えます。
 古典力学です。解析力学では一般的な共役な物理量と言うと、 $pq$ 等の様に、積の次元がエネルギー×時間の次元を持つと言う定義です。
 
 復習(q,pに対するポアソン括弧)
 ===============================
 
 まずは、 $q,p$ を基にしたポアソン括弧を区別する為、
 
 <tex>
 \left[ A , B \right]_{q,p} \equiv \dfrac{\partial A}{\partial q} \dfrac{\partial B}{\partial p} - \dfrac{\partial A}{\partial p} \dfrac{\partial B}{\partial q}
 \tag{##}
 </tex>
 
 と書くことにします。ある物理量 $X$ と変数 $q$ or $p$ のどちらか単独でポアソン括弧に入れると、それは正負の符号を無視すれば、共役な変数での微分をした事になるのでした。
 
 <tex>
 \left[ q , X \right]_{q,p} 
 &= \dfrac{\partial q}{\partial q} \dfrac{\partial X}{\partial p} - \dfrac{\partial q}{\partial p} \dfrac{\partial X}{\partial q} = \dfrac{\partial X}{\partial p}
 \tag{##}
 </tex>
 
 <tex>
 \left[ p , X \right]_{q,p} 
 &= \dfrac{\partial p}{\partial q} \dfrac{\partial X}{\partial p} - \dfrac{\partial p}{\partial p} \dfrac{\partial X}{\partial q} = - \dfrac{\partial X}{\partial q}
 \tag{##}
 </tex>
 
 実験(t,Eに対するポアソン括弧)
 ===============================
 
 ここで、共役な変数として有名な時間 $t$ とエネルギー $E$ を $q$ と $p$ の代わりに使ってみたらどうなるでしょう?自由粒子の時 $E=\dfrac{p^2}{2m}$ として、 $\left[ q , E \right]_{t,E}$ を計算すると、 $q=vt$ , $p=mv$ ですから、
 
 <tex>
 \left[ q , E \right]_{t,E} 
 &= \dfrac{\partial q}{\partial t} \dfrac{\partial E}{\partial E} - \dfrac{\partial q}{\partial E} \dfrac{\partial E}{\partial t} \\
 &= \dfrac{\partial (vt)}{\partial t} \dfrac{\partial E}{\partial E} - \dfrac{\partial (vt)}{\partial E} \dfrac{\partial E}{\partial t} \\
 &= v \cdot 1 - 0 \\
 &= \dfrac{p}{m} \\
 &= \dfrac{\partial \dfrac{p^2}{2m}}{\partial p} \\
 &= \left[ q , E \right]_{q,p}
 \tag{##}
 </tex>
 
 おや、これは偶然でしょうか?
 
 調和振動子の時( $A=\dfrac{p^2}{2m}$ )もやってみましょう。もう一つを $B=pq$ にしてみます。この
 解の一つは、 $\omega = \sqrt{\dfrac{k}{m}}$ として、 $q = C \cos (\omega t)$ 、 $p = m\dot{q} = -m \omega C \sin (\omega t)$ です。
 
 <tex>
 \left[ \dfrac{p^2}{2m} + \dfrac{k}{2}q^2 , pq \right]_{q,p} 
 &= \dfrac{\partial (p^2/2m + kq^2/2)}{\partial q} \dfrac{\partial (pq)}{\partial p} - \dfrac{\partial (p^2/2m + kq^2/2)}{\partial p} \dfrac{\partial (pq)}{\partial q} \\
 &= kq \cdot q - \dfrac{p}{m} \cdot p \\
 &= -\dfrac{p^2}{m} + kq^2 \\
 &= -m \omega^2 C^2 \sin^2 (\omega t) + k C^2 \cos (\omega t) \\
 &= k C^2 \cos (2 \omega t)
 \tag{##}
 </tex>
 
 です。一方、 $t,E$ のポアソン括弧は、
 
 <tex>
 \left[ \dfrac{p^2}{2m} + \dfrac{k}{2}q^2 , pq \right]_{t,E} 
 &= \dfrac{\partial E}{\partial t} \dfrac{\partial (pq)}{\partial E} - \dfrac{\partial E}{\partial E} \dfrac{\partial (pq)}{\partial t} \\
 &= 0 \cdot ? - 1 \cdot \dfrac{\partial (-m \omega C^2 \sin \omega t \cos \omega t)}{\partial t} \\
 &= \dfrac{\partial (\dfrac{m}{2} \omega C^2 \sin 2 \omega t)}{\partial t} \\
 &= m \omega^2 C^2 \cos (2 \omega t) \\
 &= k C^2 \cos (2 \omega t) \\
 &= \left[ \dfrac{p^2}{2m} + \dfrac{k}{2}q^2 , pq \right]_{q,p}
 \tag{##}
 </tex>
 
 なるほど、一致しましたね。これは偶然ではなさそうです。
 
 q,pの括弧とt,Eの括弧が同じ値を取ることの証明
 ================================================
 
 個々の例に囚われていてはいけません。ある程度一般化して考えましょう。
 
 <tex>
 t &= A(q,p) \\
 E &= B(q,p)
 \tag{##}
 </tex>
 
 と置いて計算してみましょう!(ちなみに $E$ は陽に $t$ に依らないものとします。実はこれは「正準変換」に含まれます。)
 
 <tex>
 \dfrac{\partial}{\partial t} &= \dfrac{\partial q}{\partial t} \dfrac{\partial }{\partial q} + \dfrac{\partial p}{\partial t} \dfrac{\partial }{\partial p} \\
 \dfrac{\partial}{\partial E} &= \dfrac{\partial q}{\partial E} \dfrac{\partial }{\partial q} + \dfrac{\partial p}{\partial E} \dfrac{\partial }{\partial p} \tag{##}
 </tex>
 
 ポアソン括弧 $t,E$ を計算すると、
 
 <tex>
 \left[ A, B \right]_{t,E} 
 &=  \dfrac{\partial A}{\partial t} \dfrac{\partial q}{\partial E} + \dfrac{\partial A}{\partial E} \dfrac{\partial B}{\partial t} \\
 &=  \left( \dfrac{\partial A}{\partial q} \dfrac{\partial q}{\partial t} + \dfrac{\partial A}{\partial p} \dfrac{\partial p}{\partial t} \right) \cdot \left( \dfrac{\partial B}{\partial q} \dfrac{\partial q}{\partial E} + \dfrac{\partial B}{\partial p} \dfrac{\partial p}{\partial E} \right) \\
 &-  \left( \dfrac{\partial A}{\partial q} \dfrac{\partial q}{\partial E} + \dfrac{\partial A}{\partial p} \dfrac{\partial p}{\partial E} \right) \cdot \left( \dfrac{\partial B}{\partial q} \dfrac{\partial q}{\partial t} + \dfrac{\partial B}{\partial p} \dfrac{\partial p}{\partial t} \right) \\
 &= \dfrac{\partial A}{\partial q} \dfrac{\partial B}{\partial q} \left(\dfrac{\partial q}{\partial t} \dfrac{\partial q}{\partial E} - \dfrac{\partial q}{\partial E} \dfrac{\partial q}{\partial t} \right) \\
 &+ \dfrac{\partial A}{\partial q} \dfrac{\partial B}{\partial p} \left(\dfrac{\partial q}{\partial t} \dfrac{\partial p}{\partial E} - \dfrac{\partial q}{\partial E} \dfrac{\partial p}{\partial t} \right) \\
 &- \dfrac{\partial A}{\partial p} \dfrac{\partial B}{\partial q} \left(\dfrac{\partial q}{\partial t} \dfrac{\partial p}{\partial E} - \dfrac{\partial q}{\partial E} \dfrac{\partial p}{\partial t} \right) \\
 &+ \dfrac{\partial A}{\partial q} \dfrac{\partial B}{\partial p} \left(\dfrac{\partial p}{\partial t} \dfrac{\partial p}{\partial E} - \dfrac{\partial p}{\partial E} \dfrac{\partial p}{\partial t} \right) \\
 &= \left( \dfrac{\partial A}{\partial q} \dfrac{\partial B}{\partial p} - \dfrac{\partial A}{\partial p} \dfrac{\partial B}{\partial q} \right) \cdot \left(\dfrac{\partial q}{\partial t} \dfrac{\partial p}{\partial E} - \dfrac{\partial q}{\partial E} \dfrac{\partial p}{\partial t} \right) \\
 &= \left[ A, B \right]_{q,p} \cdot \left[ q, p \right]_{t,E} \tag{##}
 </tex>
 
 よって、
 
 <tex>
 \left[ A, B \right]_{t,E} = \left[ A, B \right]_{q,p} \tag{##}
 </tex>
 
 を満たすには、
 
 <tex>
 \left[ q, p \right]_{t,E} = 1 \tag{##}
 </tex>
 
 が成立していなければなりません。
 
 最後の仕上げ
 ==============
 
 <tex>
 \left[ q, p \right]_{t,E} = 1 \tag{##}
 </tex>
 
 を示します。今、ハミルトニアンは $H=E$ とします。記号こそ違うものの同じエネルギーを表しています。
 
 <tex>
 \dfrac{dq}{dt} \ \ \ &= \dfrac{\partial H }{\partial p} = \ \ \ \dfrac{\partial q }{\partial t} + \dfrac{\partial q}{\partial E} \dfrac{\partial E}{\partial t} \\
 \dfrac{dp}{dt} \ \ \ &= -\dfrac{\partial H }{\partial q} = \ \ \ - \dfrac{\partial p }{\partial t} - \dfrac{\partial p}{\partial E} \dfrac{\partial E}{\partial t} \\
 \tag{##}
 </tex>
 
 より、
 
 <tex>
 \dfrac{\partial q }{\partial t} &= \dfrac{dq}{dt} - \dfrac{\partial q}{\partial E} \dfrac{\partial E}{\partial t} \\
 \dfrac{\partial p }{\partial t} &= \dfrac{dp}{dt} - \dfrac{\partial p}{\partial E} \dfrac{\partial E}{\partial t} 
 \tag{##}
 </tex>
 
 ですから、
 
 <tex>
 \left[ q, p \right]_{t,E}
 &= \dfrac{\partial q}{\partial t} \dfrac{\partial p}{\partial E} - \dfrac{\partial q}{\partial E} \dfrac{\partial p}{\partial t} \\
 &= \left( \dfrac{dq}{dt} - \dfrac{\partial q}{\partial E} \dfrac{\partial E}{\partial t} \right) \dfrac{\partial p}{\partial E} - \left( \dfrac{dp}{dt} - \dfrac{\partial p}{\partial E} \dfrac{\partial E}{\partial t} \right) \dfrac{\partial q}{\partial E} \\
 &= \dfrac{dq}{dt} \dfrac{\partial p}{\partial E} - \dfrac{dp}{dt} \dfrac{\partial q}{\partial E} \\
 &= \dfrac{\partial H }{\partial p} \dfrac{\partial p}{\partial E} + \dfrac{\partial H }{\partial q} \dfrac{\partial q}{\partial E} \\
 &= \dfrac{\partial H }{\partial E} \\
 &= 1 \ \ \ (\because H=E)
 \tag{##}
 </tex>
 
 ふう、証明できました。ここで、共役な変数とは、最初に一対 $p,q$ を共役と決めたら、確かめたい変数 $A,B$ として、そのポアソン括弧 $\left[ A,B \right]_{q,p}$ 、もしくは、 $\left[ q,p \right]_{A,B}$ が1に等しくなるような変数の組と言えそうです。逆に1に等しい時、それらの変数は共役であるかどうかは分かりませんでした。
 
 ちなみに、僕が気になっていることを話しておくと、電荷 $Q[C]$ と磁荷 $Q_m[Wb]$ の積は、 $\mathrm{rot} \bm{E} = - \dfrac{\partial \bm{B}}{\partial t}$ から確かめると、なんと $[C][Wb]=[J \cdot s]$ で(少なくとも)次元については共役かもしれないということが言えます。実際に、 $[Q,Q_m]_{q,p}=1$ になるかどうかは分かりません。そもそもこんな計算は現実に対応する物理的な量ではないかもしれませんが、興味の一つに持っておこうと思います。
 
 今日はここまで、お疲れさまでした!
 
 @@reference: 大貫義郎,物理テキストシリーズ2『解析力学』,岩波書店,1987,第9章,4000077422@@
 @@reference: 大貫義郎,物理テキストシリーズ2 解析力学,岩波書店,1987,第9章,4000077422@@
 
 @@author:クロメル@@
 @@accept:2019-09-15@@
 @@category:解析力学@@
 @@id:poissonBracketOfEt@@
トップ   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS
Modified by 物理のかぎプロジェクト PukiWiki 1.4.6 Copyright © 2001-2005 PukiWiki Developers Team. License is GPL.
Based on "PukiWiki" 1.3 by yu-ji Powered by PHP 5.3.29 HTML convert time to 0.007 sec.