tecktom.com

出てくるもの (1):正準運動量と Hamiltonian

2026年10月15日9分で読めますAI生成教科書『虎の巻』第2部 第2章 電磁場中の荷電粒子
目次
  1. 正準運動量
  2. Hamiltonian
  3. エネルギーと時間による変化
  4. 速度の成分どうしの Poisson 括弧
  5. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第2部 第2章 §2.2.3 ★★★ p=mr˙+qA\bm{p}=m\dot{\bm{r}}+q\bm{A}、H=(p−qA)2/(2m)+qϕH=(\bm{p}-q\bm{A})^2/(2m)+q\phi。HH の値に A\bm{A} は入らず、速度の成分どうしの Poisson 括弧は磁場になる

正準運動量

正準運動量は式 (2.2.5) そのもので、

p=mr˙+qA(r,t)(2.2.10)\bm{p}=m\dot{\bm{r}}+q\bm{A}(\bm{r},t)\tag{2.2.10}

である。p\bm{p} は mr˙m\dot{\bm{r}} と qAq\bm{A} の和で、速度に質量を掛けた量(運動学的な運動量)とは一致しない。第1部第7章の例題(第1部第7章の「速度の1次の項は運動量を変えるがエネルギーを変えない」)で、速度の1次の項 b(q)q˙b(q)\dot q が運動量を bb だけずらすことを見た。ここではその bb が qAq\bm{A} である。

Hamiltonian

LL の速度についての2階微分は mm 倍の単位行列で正則なので、式 (2.2.10) は r˙=(p−qA)/m\dot{\bm{r}}=(\bm{p}-q\bm{A})/m と解ける。これを H=p⋅r˙−LH=\bm{p}\cdot\dot{\bm{r}}-L に入れる(第1部第8章)。p⋅r˙−qr˙⋅A=(p−qA)⋅r˙=(p−qA)2/m\bm{p}\cdot\dot{\bm{r}}-q\dot{\bm{r}}\cdot\bm{A}=(\bm{p}-q\bm{A})\cdot\dot{\bm{r}}=(\bm{p}-q\bm{A})^2/m を使うと

H(r,p,t)=(p−qA(r,t))22m+qϕ(r,t)(2.2.11)H(\bm{r},\bm{p},t)=\frac{\bigl(\bm{p}-q\bm{A}(\bm{r},t)\bigr)^2}{2m}+q\phi(\bm{r},t)\tag{2.2.11}

を得る。p−qA=mr˙\bm{p}-q\bm{A}=m\dot{\bm{r}} を入れると、HH の値は

H=12mr˙2+qϕ(2.2.12)H=\frac12m\dot{\bm{r}}^2+q\phi\tag{2.2.12}

で、A\bm{A} を含まない。磁場は HH の関数の形(p\bm{p} を p−qA\bm{p}-q\bm{A} にずらすこと)には入るが、エネルギーの値には入らない。自由粒子の Hamiltonian p2/(2m)\bm{p}^2/(2m) の p\bm{p} を p−qA\bm{p}-q\bm{A} に置き換え、qϕq\phi を足したものになっている。この置き換えを最小結合と呼ぶ。

Hamilton 方程式を書くと、

r˙=∂H∂p=p−qAm,p˙i=−∂H∂xi=q∑jpj−qAjm∂Aj∂xi−q∂ϕ∂xi=q∑jx˙j∂Aj∂xi−q∂ϕ∂xi.\dot{\bm{r}}=\frac{\partial H}{\partial \bm{p}}=\frac{\bm{p}-q\bm{A}}{m},\qquad \dot p_i=-\frac{\partial H}{\partial x_i}=q\sum_j\frac{p_j-qA_j}{m}\frac{\partial A_j}{\partial x_i}-q\frac{\partial \phi}{\partial x_i}=q\sum_j\dot x_j\frac{\partial A_j}{\partial x_i}-q\frac{\partial \phi}{\partial x_i}.

第1式は式 (2.2.10) の書き直しである。第2式の右辺は ∂L/∂xi\partial L/\partial x_i(式 (2.2.7))に等しく、p˙i=ddt(mx˙i+qAi)\dot p_i=\frac{\mathrm{d} }{\mathrm{d} t}(m\dot x_i+qA_i) と組み合わせると式 (2.2.9) が再び出る。

エネルギーと時間による変化

第1部第8章の dH/dt=∂H/∂t\mathrm{d} H/\mathrm{d} t=\partial H/\partial t を当てると、

dHdt=∂H∂t=q∂ϕ∂t−q r˙⋅∂A∂t.(2.2.13)\frac{\mathrm{d} H}{\mathrm{d} t}=\frac{\partial H}{\partial t}=q\frac{\partial \phi}{\partial t}-q\,\dot{\bm{r}}\cdot\frac{\partial \bm{A}}{\partial t}.\tag{2.2.13}

ϕ\phi と A\bm{A} が時刻によらなければ、H=12mr˙2+qϕH=\frac12m\dot{\bm{r}}^2+q\phi が保存する。磁場がどれだけ強くても、静的な場の中の粒子のエネルギーは運動エネルギーと静電的なエネルギーの和で、磁場はそれを変えない。A\bm{A} が時刻とともに変わると、−qr˙⋅∂A/∂t-q\dot{\bm{r}}\cdot\partial\bm{A}/\partial t の分だけエネルギーが変わる。これは電場のうち −∂A/∂t-\partial\bm{A}/\partial t の部分(電磁誘導の電場)が粒子にする仕事率である(2.2.6節)。

速度の成分どうしの Poisson 括弧

相空間の関数としての運動学的な運動量 mx˙i=pi−qAi(r,t)m\dot x_i=p_i-qA_i(\bm{r},t) の Poisson 括弧(第1部第11章)を計算する。{pi,Aj(r)}=−∂Aj/∂xi\{p_i,A_j(\bm{r})\}=-\partial A_j/\partial x_i を使うと

{px−qAx, py−qAy}=−q{px,Ay}−q{Ax,py}=q(∂Ay∂x−∂Ax∂y)=qBz(2.2.14)\{p_x-qA_x,\ p_y-qA_y\}=-q\{p_x,A_y\}-q\{A_x,p_y\}=q\Bigl(\frac{\partial A_y}{\partial x}-\frac{\partial A_x}{\partial y}\Bigr)=qB_z\tag{2.2.14}

となり、ほかの組も同じ計算で {mx˙i,mx˙j}=q∑kϵijkBk\{m\dot x_i,m\dot x_j\}=q\sum_k\epsilon_{ijk}B_k になる。正準運動量の成分どうしは {px,py}=0\{p_x,p_y\}=0 だが、速度の成分どうしは磁場があると Poisson 括弧が 00 にならない。一方、位置と速度の組は {xi,mx˙j}={xi,pj}=δij\{x_i,m\dot x_j\}=\{x_i,p_j\}=\delta_{ij} のまま変わらない。

式 (2.2.14) の右辺は磁場だけで書けている。A\bm{A} を取り替えても変わらない量どうしの括弧が、同じく変わらない量になった。量子力学では Poisson 括弧が交換子 [⋅,⋅]/(iℏ)[\cdot,\cdot]/(\mathrm{i}\hbar) に置き換わるので、磁場の中では速度の xx 成分と yy 成分を同時に確定させることができない。この事実が一様な磁場の中のエネルギー準位(Landau 準位)の出発点になる(第4部第6章)。

例題 2.2.3 円運動の中心の座標の Poisson 括弧 ★★★

オリジナル一様な磁場 B=Bez\bm{B}=B\bm{e}_z(q>0q>0, B>0B>0)の中の xyxy 平面内の運動を考える。A\bm{A} は ∇×A=Bez\nabla\times\bm{A}=B\bm{e}_z を満たす任意のものとする。円運動の中心の座標を、相空間の関数として

xg≔x+py−qAyqB,yg≔y−px−qAxqBx_{\mathrm g}\coloneqq x+\frac{p_y-qA_y}{qB},\qquad y_{\mathrm g}\coloneqq y-\frac{p_x-qA_x}{qB}

と定める。

  1. (1) 例題 2.2.2 の円運動で、(xg,yg)(x_{\mathrm g},y_{\mathrm g}) が円の中心であることを確かめよ。
  2. (2) {xg, px−qAx}\{x_{\mathrm g},\,p_x-qA_x\} と {xg, py−qAy}\{x_{\mathrm g},\,p_y-qA_y\} を求め、xgx_{\mathrm g} が保存することを示せ(H=[(px−qAx)2+(py−qAy)2]/(2m)H=[(p_x-qA_x)^2+(p_y-qA_y)^2]/(2m))。
  3. (3) {xg,yg}\{x_{\mathrm g},y_{\mathrm g}\} を求めよ。
解答

(1) py−qAy=my˙p_y-qA_y=m\dot y, px−qAx=mx˙p_x-qA_x=m\dot x なので xg=x+y˙/ωcx_{\mathrm g}=x+\dot y/\omega_{\mathrm c}, yg=y−x˙/ωcy_{\mathrm g}=y-\dot x/\omega_{\mathrm c}。例題 2.2.2 の解 x+iy=(xg+iyg)+iuωcx+\mathrm{i} y=(x_{\mathrm g}+\mathrm{i} y_{\mathrm g})+\frac{\mathrm{i} u}{\omega_{\mathrm c}}(u=x˙+iy˙u=\dot x+\mathrm{i}\dot y)から、(x+iy)−i(x˙+iy˙)ωc=(x+y˙ωc)+i(y−x˙ωc)(x+\mathrm{i} y)-\frac{\mathrm{i}(\dot x+\mathrm{i}\dot y)}{\omega_{\mathrm c}}=\bigl(x+\frac{\dot y}{\omega_{\mathrm c}}\bigr)+\mathrm{i}\bigl(y-\frac{\dot x}{\omega_{\mathrm c}}\bigr) が中心に等しい。

(2) {x,px−qAx}=1\{x,p_x-qA_x\}=1, {x,py−qAy}=0\{x,p_y-qA_y\}=0 と式 (2.2.14) の {py−qAy,px−qAx}=−qB\{p_y-qA_y,p_x-qA_x\}=-qB を使うと

{xg,px−qAx}=1+−qBqB=0,{xg,py−qAy}=0+{py−qAy,py−qAy}qB=0.\{x_{\mathrm g},p_x-qA_x\}=1+\frac{-qB}{qB}=0,\qquad\{x_{\mathrm g},p_y-qA_y\}=0+\frac{\{p_y-qA_y,p_y-qA_y\}}{qB}=0.

HH は2つの速度の成分だけの関数なので {xg,H}=0\{x_{\mathrm g},H\}=0 で、xgx_{\mathrm g} は保存する。ygy_{\mathrm g} も同様である。

(3) {x,y}=0\{x,y\}=0, {x,px−qAx}=1\{x,p_x-qA_x\}=1, {py−qAy,y}=−1\{p_y-qA_y,y\}=-1 と式 (2.2.14) から

{xg,yg}=−1qB{x,px−qAx}+1qB{py−qAy,y}−1(qB)2{py−qAy,px−qAx}=−1qB−1qB+1qB=−1qB.\begin{aligned}\{x_{\mathrm g},y_{\mathrm g}\}&=-\frac{1}{qB}\{x,p_x-qA_x\}+\frac{1}{qB}\{p_y-qA_y,y\}-\frac{1}{(qB)^2}\{p_y-qA_y,p_x-qA_x\}\\ &=-\frac{1}{qB}-\frac{1}{qB}+\frac{1}{qB}=-\frac{1}{qB}.\end{aligned}

円の中心の2つの座標は、互いに共役な正準変数の組(−1/(qB)-1/(qB) 倍を除いて)になっている。量子力学では [x^g,y^g]=−iℏ/(qB)[\hat x_{\mathrm g},\hat y_{\mathrm g}]=-\mathrm{i}\hbar/(qB) で、中心の位置は面積 2πℏ/(qB)2\pi\hbar/(qB) より細かく決められない。面積 SS の試料に入る中心の数はおよそ qBS/(2πℏ)qBS/(2\pi\hbar) で、これが Landau 準位1つあたりの状態の数になる(第4部第6章)。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ