| この記事は検証可能な参考文献や出典が全く示されていないか、不十分です。出典を追加して記事の信頼性向上にご協力ください。(このテンプレートの使い方) 出典検索?: "正準変数" – ニュース · 書籍 · スカラー · CiNii · J-STAGE · NDL · dlib.jp · ジャパンサーチ · TWL(2013年5月) |
正準変数(せいじゅんへんすう、英: canonical variable)とは、ハミルトン形式の解析力学において、物体の運動を記述する基本変数として用いられる一般化座標と一般化運動量の組をいう[1][2][3]。しばしば一般化座標は文字 q 、一般化運動量は p で表される。正準(カノニカル、英: canonical)という語は標準的、慣例的という意味を表す[4]。ウィリアム・ローワン・ハミルトンによって導入された正準変数による形式に正準(仏: canonique)という語を充てたのは、カール・グスタフ・ヤコブ・ヤコビである[5][6]。
ニュートン力学やラグランジュ力学においては、基本変数が一般化座標 q とその時間微分である一般化速度 ·q であったが、ハミルトン力学においては、一般化座標と一般化運動量が用いられる。ラグランジアンL=L(q,·q,t) は一般化座標、一般化速度、時間の関数である。ここで L にルジャンドル変換
![{\displaystyle H=p{\dot {q}}-L}](https://wikimedia.org/api/rest_v1/media/math/render/svg/de8a51af06dabd30181d7b82ab40102dbfb26f40)
を施すことで一般化座標、一般化運動量、時間を変数とする関数ハミルトニアンH=H(q,p,t)が得られ、正準方程式
![{\displaystyle {\dot {q}}={\frac {\partial H}{\partial p}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/72268721e438f08bd3f79505999e16672b3e31dc)
![{\displaystyle {\dot {p}}=-{\frac {\partial H}{\partial q}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1ca7d48dc4a7a4bd8b9b6b0f4a4aa819b672d569)
が得られる。
一般化座標q=(q1,.., qn)と正準共役な一般化運動量 p=(p1,.., pn)の組による2n個の変数
(q, p)=(q1,.., qn, p1,.., pn)を系の状態を指定する独立な変数と見なしたときに、(q, p)を正準変数という[1][2][3]。このとき、一般化座標 q を正準座標、一般化運動量 p を正準運動量とも呼ぶ[2]。正準変数(q, p)を座標とする2n次元の空間を相空間という[1][2][3]。系の状態は相空間上の1点で指定される。ハミルトニアンをH=H(q,p,t)とするとき、物体の運動を記述する運動方程式は
![{\displaystyle {\dot {q}}^{i}(t)={\frac {\partial H}{\partial p_{i}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/765d9cac24a073b76b4749d23d74567855267ce3)
![{\displaystyle {\dot {p}}_{i}(t)=-{\frac {\partial H}{\partial q^{i}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4dc530f6adc1a9caaf0e56065c43e01a078a9068)
で与えられる。但し、ドット記号は時間微分を表す。この方程式をハミルトンの正準方程式という。この正準方程式で時間発展が定まる力学系を自由度nのハミルトン力学系、またはハミルトン系という。ハミルトン力学系での系の時間発展は相空間上の軌道(q(t), p(t))で与えられる。
2n個の変数 z=(z1,.., zn, zn+1,.., z2n)を
![{\displaystyle (z^{1},\cdots ,z^{n},z^{n+1},\cdots ,z^{2n}):=(q^{1},\cdots ,q^{n},p_{1},\cdots ,p_{n})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f343f99e5ab262fe9580739d43f052feb687dc9d)
で定義すると正準変数をまとめて、z=(q, p)で表記することができる。列ベクトルでの表記を z=(q1,.., qn, p1,.., pn)Tとすると、正準方程式は
![{\displaystyle {\dot {\boldsymbol {z}}}(t)=\Omega \nabla H({\boldsymbol {z}}(t))}](https://wikimedia.org/api/rest_v1/media/math/render/svg/577efe7ced5455a9a5b08d0ac9d2516372822e1c)
となる。ここでは∇はナブラ演算子である。Ω=(ωi,j)は
![{\displaystyle \Omega ={\begin{pmatrix}0_{n}&I_{n}\\-I_{n}&0_{n}\end{pmatrix}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2da4b0d7e12754fe9736a2be9fbdaa2fdf213248)
で定義される2n × 2n行列である。Ω内の0nはn 次の零行列、Inはn 次の単位行列である。
相空間上の関数f=f(q, p, t)、g=g(q, p, t)に対し
![{\displaystyle \{f,g\}:=\sum _{i=1}^{n}{\Big (}{\frac {\partial f}{\partial q_{i}}}{\frac {\partial g}{\partial p_{i}}}-{\frac {\partial f}{\partial p_{i}}}{\frac {\partial g}{\partial q_{i}}}{\Big )}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3b5db2f88366b600a4a465c504746958556ecaea)
で定義される{f, g}をポアソン括弧という。正準方程式による時間発展(q(t), p(t))に対し、f=f(q(t), p(t), t)の時間変化は
![{\displaystyle {\frac {d}{dt}}f(q(t),p(t),t)={\frac {\partial }{\partial t}}f(q(t),p(t),t)+\{f,H\}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/874b2897368e7f555383b905bf9eb8913b397b25)
とハミルトニアン H とのポアソン括弧で表される。特に正準変数の時間発展を記述する正準方程式は
![{\displaystyle {\dot {q}}_{i}(t)=\{q_{i},H\}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6d43b26786b1f9a99316da7897056ff4fb936546)
![{\displaystyle {\dot {p}}_{i}(t)=\{p_{i},H\}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d999fee36c01422560ab42bdb93b8c2cabfd62fa)
となる。
正準変数をまとめた表記z=(q, p)と行列Ω=(ωij)を用いるとポアソン括弧は
![{\displaystyle \{f,g\}=\sum _{i,j}\omega _{ij}{\frac {\partial f}{\partial z_{i}}}{\frac {\partial g}{\partial z_{j}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fc16104984b03e076fc724da18906d071eb70941)
となる。
電荷量eをとする質量mの荷電粒子の電磁場中における運動を考える[4]。3次元空間での粒子の位置座標x =(x, y, z )を一般化座標にとる。スカラーポテンシャルをφ(x,t)、ベクトルポテンシャルをA(x,t)とすると、荷電粒子のラグランジアンは、
![{\displaystyle L({\boldsymbol {x}},{\boldsymbol {\dot {x}}},t)={\frac {1}{2}}m{\boldsymbol {\dot {x}}}^{2}-e\phi ({\boldsymbol {x}},t)+e{\boldsymbol {\dot {x}}}\cdot {\boldsymbol {A}}({\boldsymbol {x}},t)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/40e1e6098f600f01a6710f8fdc7dba7933c2cbbb)
で与えられる。ここでx =(x, y, z )に正準共役な運動量p =(px, py, pz )は
![{\displaystyle p_{x}:={\frac {\partial L}{\partial {\dot {x}}}}=m{\dot {x}}+eA_{1}({\boldsymbol {x}},t)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b53c4da21fe88212e5d4511ac3b725f6c51dd59a)
![{\displaystyle p_{y}:={\frac {\partial L}{\partial {\dot {y}}}}=m{\dot {y}}+eA_{2}({\boldsymbol {x}},t)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2db1583444ba5e57367b93f6bac59c2223aea235)
![{\displaystyle p_{z}:={\frac {\partial L}{\partial {\dot {z}}}}=m{\dot {z}}+eA_{3}({\boldsymbol {x}},t)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a386e1b42b976d121e8b2190ecb78b87097f8a36)
である。これをベクトルで表記すると
![{\displaystyle {\boldsymbol {p}}={\frac {\partial L}{\partial {\boldsymbol {\dot {x}}}}}=m{\boldsymbol {\dot {x}}}+e{\boldsymbol {A}}({\boldsymbol {x}},t)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/813f5b4ff1e71281a7bfe8c5d74f0f7ea01fa13b)
となる。ハミルトニアンは
![{\displaystyle {\begin{aligned}H({\boldsymbol {x}},{\boldsymbol {p}},t)&={\boldsymbol {p}}\cdot {\boldsymbol {\dot {x}}}-L\\&={\frac {1}{2m}}{\bigl (}{\boldsymbol {p}}-e{\boldsymbol {A}}({\boldsymbol {x}},t){\bigr )}^{2}+e\phi ({\boldsymbol {x}},t)\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/673b3801fbab975d1d980893269dafc8945313a8)
である。
距離r=√x2+y2+z2のみに依存する中心力ポテンシャルV=V(r)の下での質量mの粒子の運動を考える[4]。3次元空間での粒子の位置の極座標表示を(x,y,z)=(rsinθ cosφ, rsinθ sinφ, rsinφ)とし、極座標(r,θ,φ)を一般化座標にとる。このとき、粒子のラグランジアンは
![{\displaystyle L(r,\theta ,\phi ,{\dot {r}},{\dot {\theta }},{\dot {\phi }})={\frac {1}{2}}m{\bigl (}{\dot {r}}^{2}+r^{2}{\dot {\theta }}^{2}+r^{2}\sin ^{2}{\theta }\,{\dot {\phi }}^{2}{\bigr )}-V(r)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/59c6a3f2a7af0de2b30513d236f3648eddf68502)
で与えられる。(r,θ,φ)に正準共役な運動量(pr,pθ,pφ)は
![{\displaystyle p_{r}:={\frac {\partial L}{\partial {\dot {r}}}}=m{\dot {r}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0a35ba56ac8a7501b0e9b3b8861fdeadc3d2021d)
![{\displaystyle p_{\theta }:={\frac {\partial L}{\partial {\dot {\theta }}}}=mr^{2}{\dot {\theta }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/77ee4c08c1872a4a020f6b0f708a2d3b48054da5)
![{\displaystyle p_{\phi }:={\frac {\partial L}{\partial {\dot {\phi }}}}=mr^{2}\sin ^{2}{\theta }\,{\dot {\phi }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bb16cc8554a51ab2f7b2e2d6a42d140ae801e53b)
である。ハミルトニアンは
![{\displaystyle {\begin{aligned}H(r,\theta ,\phi ,p_{r},p_{\theta },p_{\phi })&=p_{r}{\dot {r}}+p_{\theta }{\dot {\theta }}+p_{\phi }{\dot {\phi }}-L\\&={\frac {p_{r}^{\,2}}{2m}}+{\frac {p_{\theta }^{\,2}}{2mr^{2}}}+{\frac {p_{\phi }^{\,2}}{2mr^{2}\sin ^{2}{\theta }}}+V(r)\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2c53dbc1c4248790367a6cfaaf8f98abbf09c765)
である。
Wrecke 4: 129-126.