解析力学から様々な理論分野を組み立てる 第1部 第12章 §1.12.2 ★★★ { Q i , P j } = δ i j \{Q_i,P_j\}=\delta_{ij} { Q i , P j } = δ ij などを保つ変換を正準変換と定義し、それが Hamilton 方程式の形を保つことを示す
定義と行列の形
定義 1.12.1(正準変換)
f f f 自由度の相空間の変数の取り替え ( q , p ) ↦ ( Q , P ) (q,p)\mapsto(Q,P) ( q , p ) ↦ ( Q , P ) (Q i = Q i ( q , p , t ) Q_i=Q_i(q,p,t) Q i = Q i ( q , p , t ) , P i = P i ( q , p , t ) P_i=P_i(q,p,t) P i = P i ( q , p , t ) 、滑らかで逆写像を持つ)が、各時刻で
{ Q i , Q j } = 0 , { P i , P j } = 0 , { Q i , P j } = δ i j ( i , j = 1 , … , f ) (1.12.4) \{Q_i,Q_j\}=0,\qquad\{P_i,P_j\}=0,\qquad\{Q_i,P_j\}=\delta_{ij}\qquad(i,j=1,\dots,f)\tag{1.12.4} { Q i , Q j } = 0 , { P i , P j } = 0 , { Q i , P j } = δ ij ( i , j = 1 , … , f ) ( 1.12.4 )
を満たすとき、正準変換 と呼ぶ。括弧は古い変数 ( q , p ) (q,p) ( q , p ) で計算する。
式 (1.12.4) は、古い変数が満たす { q i , q j } = 0 \{q_i,q_j\}=0 { q i , q j } = 0 , { p i , p j } = 0 \{p_i,p_j\}=0 { p i , p j } = 0 , { q i , p j } = δ i j \{q_i,p_j\}=\delta_{ij} { q i , p j } = δ ij と同じ形の関係を、新しい変数にも要求している。これを行列で書いておく。2 f 2f 2 f 個の変数を1列に並べて z = ( q 1 , … , q f , p 1 , … , p f ) z=(q_1,\dots,q_f,p_1,\dots,p_f) z = ( q 1 , … , q f , p 1 , … , p f ) とし、2 f × 2 f 2f\times2f 2 f × 2 f 行列
J = ( 0 1 f − 1 f 0 ) (1.12.5) \mathsf J=\begin{pmatrix}0&1_f\\-1_f&0\end{pmatrix}\tag{1.12.5} J = ( 0 − 1 f 1 f 0 ) ( 1.12.5 )
を使うと(1 f 1_f 1 f は f × f f\times f f × f の単位行列)、Poisson 括弧は勾配の列ベクトル ∇ z A = ( ∂ A / ∂ z a ) \nabla_zA=(\partial A/\partial z_a) ∇ z A = ( ∂ A / ∂ z a ) を使って
{ A , B } = ∑ i = 1 f ( ∂ A ∂ q i ∂ B ∂ p i − ∂ A ∂ p i ∂ B ∂ q i ) = ( ∇ z A ) T J ∇ z B (1.12.6) \{A,B\}=\sum_{i=1}^f\Bigl(\frac{\partial A}{\partial q_i}\frac{\partial B}{\partial p_i}-\frac{\partial A}{\partial p_i}\frac{\partial B}{\partial q_i}\Bigr)=(\nabla_zA)^{\mathsf T}\,\mathsf J\,\nabla_zB\tag{1.12.6} { A , B } = i = 1 ∑ f ( ∂ q i ∂ A ∂ p i ∂ B − ∂ p i ∂ A ∂ q i ∂ B ) = ( ∇ z A ) T J ∇ z B ( 1.12.6 )
と書くことができる。新しい変数も Z = ( Q 1 , … , Q f , P 1 , … , P f ) Z=(Q_1,\dots,Q_f,P_1,\dots,P_f) Z = ( Q 1 , … , Q f , P 1 , … , P f ) と並べ、Jacobi 行列を M a b = ∂ Z a / ∂ z b M_{ab}=\partial Z_a/\partial z_b M ab = ∂ Z a / ∂ z b とすると、∇ z Z a \nabla_zZ_a ∇ z Z a は M M M の第 a a a 行なので、新しい変数どうしの括弧を並べた行列は ( { Z a , Z b } ) = M J M T \bigl(\{Z_a,Z_b\}\bigr)=M\mathsf JM^{\mathsf T} ( { Z a , Z b } ) = M J M T になる。式 (1.12.4) は ( { Z a , Z b } ) = J \bigl(\{Z_a,Z_b\}\bigr)=\mathsf J ( { Z a , Z b } ) = J のことなので、
正準変換 ⟺ M J M T = J ⟺ M T J M = J . (1.12.7) \text{正準変換}\iff M\mathsf JM^{\mathsf T}=\mathsf J\iff M^{\mathsf T}\mathsf JM=\mathsf J .\tag{1.12.7} 正準変換 ⟺ M J M T = J ⟺ M T J M = J . ( 1.12.7 )
2つ目の同値は、J − 1 = − J \mathsf J^{-1}=-\mathsf J J − 1 = − J を使って M J M T = J M\mathsf JM^{\mathsf T}=\mathsf J M J M T = J の両辺の逆行列をとると ( M T ) − 1 J M − 1 = J (M^{\mathsf T})^{-1}\mathsf J M^{-1}=\mathsf J ( M T ) − 1 J M − 1 = J となり、左から M T M^{\mathsf T} M T 、右から M M M を掛ければ出る。M T J M = J M^{\mathsf T}\mathsf JM=\mathsf J M T J M = J を満たす行列をシンプレクティック行列 と呼ぶ。
1自由度では 2 × 2 2\times2 2 × 2 行列について M J M T = ( det M ) J M\mathsf JM^{\mathsf T}=(\det M)\,\mathsf J M J M T = ( det M ) J が恒等的に成り立つので、条件は det M = 1 \det M=1 det M = 1 だけになり、1.12.1 節の式 (1.12.3) に戻る。1自由度の正準変換とは、向きを変えずに面積を保つ変換のことである。
括弧はどちらの変数で計算しても同じ
定理 1.12.2(Poisson 括弧の不変性 ★★★)
( q , p ) ↦ ( Q , P ) (q,p)\mapsto(Q,P) ( q , p ) ↦ ( Q , P ) が正準変換なら、相空間の任意の関数 A A A , B B B について、括弧を古い変数で計算しても新しい変数で計算しても同じ値になる。
{ A , B } q , p = { A , B } Q , P . \{A,B\}_{q,p}=\{A,B\}_{Q,P}. { A , B } q , p = { A , B } Q , P .
証明 A A A を新しい変数の関数と見ると、連鎖律で ∂ A / ∂ z b = ∑ a ( ∂ A / ∂ Z a ) ( ∂ Z a / ∂ z b ) \partial A/\partial z_b=\sum_a(\partial A/\partial Z_a)(\partial Z_a/\partial z_b) ∂ A / ∂ z b = ∑ a ( ∂ A / ∂ Z a ) ( ∂ Z a / ∂ z b ) 、すなわち ∇ z A = M T ∇ Z A \nabla_zA=M^{\mathsf T}\nabla_ZA ∇ z A = M T ∇ Z A 。これを式 (1.12.6) に入れて
{ A , B } q , p = ( ∇ Z A ) T M J M T ∇ Z B = ( ∇ Z A ) T J ∇ Z B = { A , B } Q , P . \{A,B\}_{q,p}=(\nabla_ZA)^{\mathsf T}M\mathsf JM^{\mathsf T}\nabla_ZB=(\nabla_ZA)^{\mathsf T}\mathsf J\,\nabla_ZB=\{A,B\}_{Q,P}. { A , B } q , p = ( ∇ Z A ) T M J M T ∇ Z B = ( ∇ Z A ) T J ∇ Z B = { A , B } Q , P .
2つ目の等号で正準変換の条件 (1.12.7) を使った。
この定理から、正準変換が Hamilton 方程式の形を保つことが出る。ここでは変換が時刻を陽に含まない場合を扱い、時刻を含む場合は1.12.5 節で扱う。
定理 1.12.3(Hamilton 方程式の形が保たれる ★★★)
時刻を陽に含まない正準変換 ( q , p ) ↦ ( Q , P ) (q,p)\mapsto(Q,P) ( q , p ) ↦ ( Q , P ) と Hamiltonian H ( q , p , t ) H(q,p,t) H ( q , p , t ) について、H H H を新しい変数で書き直したものを K ( Q , P , t ) = H ( q ( Q , P ) , p ( Q , P ) , t ) K(Q,P,t)=H\bigl(q(Q,P),p(Q,P),t\bigr) K ( Q , P , t ) = H ( q ( Q , P ) , p ( Q , P ) , t ) とする。( q , p ) (q,p) ( q , p ) が H H H の Hamilton 方程式に従うなら、( Q , P ) (Q,P) ( Q , P ) は K K K の Hamilton 方程式
Q ˙ i = ∂ K ∂ P i , P ˙ i = − ∂ K ∂ Q i \dot Q_i=\frac{\partial K}{\partial P_i},\qquad\dot P_i=-\frac{\partial K}{\partial Q_i} Q ˙ i = ∂ P i ∂ K , P ˙ i = − ∂ Q i ∂ K
に従う。
証明 任意の関数 A ( q , p ) A(q,p) A ( q , p ) の時間変化は d A d t = { A , H } q , p \frac{\mathrm{d} A}{\mathrm{d} t}=\{A,H\}_{q,p} d t d A = { A , H } q , p である(第1部第11章 )。A = Q i A=Q_i A = Q i とおき、定理 1.12.2 で括弧を新しい変数で計算すると
Q ˙ i = { Q i , H } q , p = { Q i , K } Q , P = ∑ j ( ∂ Q i ∂ Q j ∂ K ∂ P j − ∂ Q i ∂ P j ∂ K ∂ Q j ) = ∂ K ∂ P i . \dot Q_i=\{Q_i,H\}_{q,p}=\{Q_i,K\}_{Q,P}=\sum_j\Bigl(\frac{\partial Q_i}{\partial Q_j}\frac{\partial K}{\partial P_j}-\frac{\partial Q_i}{\partial P_j}\frac{\partial K}{\partial Q_j}\Bigr)=\frac{\partial K}{\partial P_i}. Q ˙ i = { Q i , H } q , p = { Q i , K } Q , P = j ∑ ( ∂ Q j ∂ Q i ∂ P j ∂ K − ∂ P j ∂ Q i ∂ Q j ∂ K ) = ∂ P i ∂ K .
A = P i A=P_i A = P i でも同じ計算で P ˙ i = { P i , K } Q , P = − ∂ K / ∂ Q i \dot P_i=\{P_i,K\}_{Q,P}=-\partial K/\partial Q_i P ˙ i = { P i , K } Q , P = − ∂ K / ∂ Q i を得る。
正準変換の例
点変換 :座標だけを Q i = Q i ( q ) Q_i=Q_i(q) Q i = Q i ( q ) と取り替え、運動量を P i = ∑ j ( ∂ q j / ∂ Q i ) p j P_i=\sum_j(\partial q_j/\partial Q_i)\,p_j P i = ∑ j ( ∂ q j / ∂ Q i ) p j とおく変換は正準変換である。1自由度なら Q = g ( q ) Q=g(q) Q = g ( q ) , P = p / g ′ ( q ) P=p/g'(q) P = p / g ′ ( q ) で、{ Q , P } = g ′ ( q ) ⋅ 1 g ′ ( q ) − 0 = 1 \{Q,P\}=g'(q)\cdot\frac{1}{g'(q)}-0=1 { Q , P } = g ′ ( q ) ⋅ g ′ ( q ) 1 − 0 = 1 となる。Lagrangian の側で一般化座標を取り替える操作(第1部第2章 )は、相空間ではこの点変換にあたる。運動量の変換の式は、p j = ∂ L / ∂ q ˙ j p_j=\partial L/\partial\dot q_j p j = ∂ L / ∂ q ˙ j と q ˙ j = ∑ i ( ∂ q j / ∂ Q i ) Q ˙ i \dot q_j=\sum_i(\partial q_j/\partial Q_i)\dot Q_i q ˙ j = ∑ i ( ∂ q j / ∂ Q i ) Q ˙ i から P i = ∂ L / ∂ Q ˙ i P_i=\partial L/\partial\dot Q_i P i = ∂ L / ∂ Q ˙ i を計算すると出てくる。
座標と運動量の入れ替え :Q = p Q=p Q = p , P = − q P=-q P = − q は { Q , P } = 0 ⋅ 0 − 1 ⋅ ( − 1 ) = 1 \{Q,P\}=0\cdot0-1\cdot(-1)=1 { Q , P } = 0 ⋅ 0 − 1 ⋅ ( − 1 ) = 1 で正準変換である。前の節の例題で、調和振動子の4分の1周期の流れがほぼこの形になっていた。点変換は座標と運動量を混ぜないが、正準変換は混ぜてよい。Lagrangian の側の座標変換より、正準変換のほうがずっと広い。
反例 :Q = 2 q Q=2q Q = 2 q , P = 2 p P=2p P = 2 p は一見ただの目盛りの取り直しだが、{ Q , P } = 4 \{Q,P\}=4 { Q , P } = 4 なので正準変換でない。面積が4倍になるためである。H H H を書き直しただけの K K K では Hamilton 方程式の形が崩れ、K K K を 4 4 4 倍すると形が戻る。この種の定数倍の変換は1.12.6 節で扱う。
2自由度の回転
平面の回転 Q 1 = q 1 cos α − q 2 sin α Q_1=q_1\cos\alpha-q_2\sin\alpha Q 1 = q 1 cos α − q 2 sin α , Q 2 = q 1 sin α + q 2 cos α Q_2=q_1\sin\alpha+q_2\cos\alpha Q 2 = q 1 sin α + q 2 cos α に、運動量にも同じ回転 P 1 = p 1 cos α − p 2 sin α P_1=p_1\cos\alpha-p_2\sin\alpha P 1 = p 1 cos α − p 2 sin α , P 2 = p 1 sin α + p 2 cos α P_2=p_1\sin\alpha+p_2\cos\alpha P 2 = p 1 sin α + p 2 cos α を組み合わせると、M = d i a g ( R , R ) M=\mathrm{diag}(R,R) M = diag ( R , R ) (R R R は回転行列)で、R T R = 1 R^{\mathsf T}R=1 R T R = 1 から M T J M = J M^{\mathsf T}\mathsf JM=\mathsf J M T J M = J が成り立つ。これは点変換の運動量の式からも出る。座標だけ回して運動量を回さない M = d i a g ( R , 1 ) M=\mathrm{diag}(R,1) M = diag ( R , 1 ) では M T J M ≠ J M^{\mathsf T}\mathsf JM\neq\mathsf J M T J M = J となり、正準変換でない。
例題 1.12.2 4つの候補を判定して、Hamilton 方程式を書き直す ★★
オリジナル 1自由度で、次の4つの変換が正準変換かどうかを判定せよ((a) では q > 0 q>0 q > 0 とする)。
(a) Q = q 2 Q=q^2 Q = q 2 , P = p / ( 2 q ) P=p/(2q) P = p / ( 2 q )
(b) Q = p Q=p Q = p , P = − q P=-q P = − q
(c) Q = q Q=q Q = q , P = p + q 2 P=p+q^2 P = p + q 2
(d) Q = 2 q Q=2q Q = 2 q , P = 2 p P=2p P = 2 p
さらに (a) で H = 1 2 p 2 + 1 2 q 4 H=\frac12p^2+\frac12q^4 H = 2 1 p 2 + 2 1 q 4 を新しい変数で書いた K ( Q , P ) K(Q,P) K ( Q , P ) を求め、古い変数の Hamilton 方程式から計算した Q ˙ \dot Q Q ˙ , P ˙ \dot P P ˙ が K K K の Hamilton 方程式と一致することを確かめよ。
解答 1自由度なので { Q , P } = 1 \{Q,P\}=1 { Q , P } = 1 かどうかを調べればよい。(a) { Q , P } = 2 q ⋅ 1 2 q − 0 ⋅ ∂ P ∂ q = 1 \{Q,P\}=2q\cdot\frac{1}{2q}-0\cdot\frac{\partial P}{\partial q}=1 { Q , P } = 2 q ⋅ 2 q 1 − 0 ⋅ ∂ q ∂ P = 1 。正準変換(g ( q ) = q 2 g(q)=q^2 g ( q ) = q 2 の点変換)。(b) { Q , P } = 0 ⋅ 0 − 1 ⋅ ( − 1 ) = 1 \{Q,P\}=0\cdot0-1\cdot(-1)=1 { Q , P } = 0 ⋅ 0 − 1 ⋅ ( − 1 ) = 1 。正準変換。(c) { Q , P } = 1 ⋅ 1 − 0 ⋅ 2 q = 1 \{Q,P\}=1\cdot1-0\cdot2q=1 { Q , P } = 1 ⋅ 1 − 0 ⋅ 2 q = 1 。正準変換。運動量だけを位置の関数でずらす変換で、点変換の形には入らない。(d) { Q , P } = 2 ⋅ 2 − 0 = 4 \{Q,P\}=2\cdot2-0=4 { Q , P } = 2 ⋅ 2 − 0 = 4 。正準変換でない。
(a) の逆変換は q = Q q=\sqrt Q q = Q , p = 2 q P = 2 Q P p=2qP=2\sqrt Q\,P p = 2 q P = 2 Q P なので、K = 1 2 ( 2 Q P ) 2 + 1 2 Q 2 = 2 Q P 2 + 1 2 Q 2 K=\frac12(2\sqrt Q\,P)^2+\frac12Q^2=2QP^2+\frac12Q^2 K = 2 1 ( 2 Q P ) 2 + 2 1 Q 2 = 2 Q P 2 + 2 1 Q 2 。K K K の Hamilton 方程式は Q ˙ = ∂ K / ∂ P = 4 Q P \dot Q=\partial K/\partial P=4QP Q ˙ = ∂ K / ∂ P = 4 QP , P ˙ = − ∂ K / ∂ Q = − 2 P 2 − Q \dot P=-\partial K/\partial Q=-2P^2-Q P ˙ = − ∂ K / ∂ Q = − 2 P 2 − Q 。一方、古い変数では q ˙ = p \dot q=p q ˙ = p , p ˙ = − 2 q 3 \dot p=-2q^3 p ˙ = − 2 q 3 なので
Q ˙ = 2 q q ˙ = 2 q p = 2 Q ⋅ 2 Q P = 4 Q P , P ˙ = p ˙ 2 q − p q ˙ 2 q 2 = − q 2 − p 2 2 q 2 = − Q − 2 P 2 . \dot Q=2q\dot q=2qp=2\sqrt Q\cdot2\sqrt Q\,P=4QP,\qquad
\dot P=\frac{\dot p}{2q}-\frac{p\,\dot q}{2q^2}=-q^2-\frac{p^2}{2q^2}=-Q-2P^2 . Q ˙ = 2 q q ˙ = 2 q p = 2 Q ⋅ 2 Q P = 4 QP , P ˙ = 2 q p ˙ − 2 q 2 p q ˙ = − q 2 − 2 q 2 p 2 = − Q − 2 P 2 .
どちらも K K K の Hamilton 方程式と一致する。新しい変数で書いた K K K には Q P 2 QP^2 Q P 2 という座標と運動量の積が現れ、「運動エネルギー+ポテンシャル」の形ではなくなるが、Hamilton 方程式の形は保たれている。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com