共役な物理量とポアソン括弧

この記事では,通常

\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{1}

で定義されるポアソン括弧について考えることで,共役な変数とは何か?を考えます. 古典力学です.解析力学では一般的な共役な物理量と言うと, pq 等の様に,積の次元がエネルギー×時間の次元を持つと言う定義です.

復習(q,pに対するポアソン括弧)

まずは, q,p を基にしたポアソン括弧を区別する為,

\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{2}

と書くことにします.ある物理量 X と変数 q or p のどちらか単独でポアソン括弧に入れると,それは正負の符号を無視すれば,共役な変数での微分をした事になるのでした.

\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{3} \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{4}

実験(t,Eに対するポアソン括弧)

ここで,共役な変数として有名な時間 t とエネルギー Eqp の代わりに使ってみたらどうなるでしょう?自由粒子の時 E=\dfrac{p^2}{2m} として, \left[ q , E \right]_{t,E} を計算すると, q=vt , p=mv ですから,

\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{5}

おや,これは偶然でしょうか?

調和振動子の時( A=\dfrac{p^2}{2m}+\dfrac{k}{2}q^2 )もやってみましょう.もう一つを B=pq にしてみます.この 解の一つは, \omega = \sqrt{\dfrac{k}{m}} として, q = C \cos (\omega t)p = m\dot{q} = -m \omega C \sin (\omega t) です.

\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{6}

です.一方, t,E のポアソン括弧は( E は運動の定数なので \dfrac{\partial E}{\partial t}=0 ),

\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{7}

なるほど,一致しましたね.これは偶然ではなさそうです.

q,pの括弧とt,Eの括弧が同じ値を取ることの証明

個々の例に囚われていてはいけません.ある程度一般化して考えましょう.

A &= A(q,p) = A(t,E) \\B &= B(q,p) = B(t,E)\tag{8}

と置いて計算してみましょう!

\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{9}

ポアソン括弧 t,E を計算すると,

\left[ A, B \right]_{t,E} &=  \dfrac{\partial A}{\partial t} \dfrac{\partial B}{\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{10}

よって,

\left[ A, B \right]_{t,E} = \left[ A, B \right]_{q,p} \tag{11}

を満たすには,

\left[ q, p \right]_{t,E} = 1 \tag{12}

が成立していなければなりません.

最後の仕上げ

\left[ q, p \right]_{t,E} = 1 \tag{13}

を示します.今,ハミルトニアンは H=E とします.記号こそ違うものの同じエネルギーを表しています. 一応, Et 依存性を考慮しておくと,

\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{14}

より,

\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{15}

ですから,

\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{16}

ふう,証明できました.ここで,共役な変数とは,最初に一対 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 になるかどうかは分かりません.そもそもこんな計算は現実に対応する物理的な量ではないかもしれませんが,興味の一つに持っておこうと思います.

今日はここまで,お疲れさまでした!