tecktom.com

この章のまとめ

2026年10月17日12分で読めますAI生成教科書『虎の巻』第2部 第3章 Heisenberg の運動方程式と古典の正準方程式の対応

AI 生成・木下が未勉強

量子力学の様々な構成方法 第2部 第3章 §2.3.7 ★★★ 要請の2つの形・保存するもの・線形の場合の共通例題・非線形の場合のずれ・橋を1つにまとめ、総合問題で磁場中の荷電粒子まで通す

行列力学の時間発展の要請は、正準方程式の形 (E1) と交換子の形 (E2) の2通りに書け、交換子が形式的な微分を与えるという補題 2.3.3 によって同値である(2.3.1 節)。正準変数の方程式は、Hamiltonian の順序の選び方によらず、古典の Hamilton 方程式に帽子を付けた形で厳密に成り立つ。

時間発展は Jacobi 恒等式によって交換関係を保ち、古典の正準変換の量子版になっている(2.3.2 節)。H^\hat H と交換する量が保存し、エネルギーの基底では行列要素が Bohr の振動数 ωnm\omega_{nm} で振動する。これは第1章で Heisenberg が置いた要請を、この章の運動方程式から導き直したものである。互いに交換しない保存量があれば準位は縮退する。

Hamiltonian が2次以下なら、運動方程式は古典と同じ係数の線形方程式で、古典の解の初期値を演算子にしたものが厳密解になる(2.3.3・2.3.4 節)。共通例題③の波束の広がり σ(t)2=σ02+(ℏt/2mσ0)2\sigma(t)^2=\sigma_0^2+(\hbar t/2m\sigma_0)^2、①の行列要素 qn+1,n(t)∝eiωtq_{n+1,n}(t)\propto\mathrm{e}^{\mathrm{i}\omega t}、②の Rabi 振動 Pe=(Ω/Ω′)2sin⁡2(Ω′t/2)P_e=(\Omega/\Omega')^2\sin^2(\Omega't/2) は、いずれも状態を動かさずに演算子の解から出た。3次以上のポテンシャルでは期待値の方程式が閉じず、量子の ⟨q^(t)⟩\langle \hat q(t) \rangle は1本の古典軌道でなく古典の集団の平均に沿って崩れ、準位の離散性のために後で復活する(2.3.5 節)。

橋の節(2.3.6 節)では、この運動方程式が基準の Heisenberg 描像、古典の Hamilton 方程式、位相空間の Moyal 方程式、密度行列の Liouville 型の方程式とそれぞれ式で対応することを見た。次の章(第2部第4章)では、時間発展でなく交換関係の代数だけから、調和振動子と角運動量のスペクトルを出す。

総合問題 磁場と電場の中の荷電粒子を行列力学で解く ★★★

オリジナル平面(xyxy 面)内を動く質量 mm、電荷 qe>0q_{\mathrm e}>0 の粒子が、一様な磁場 B=(0,0,B)\bm{B}=(0,0,B) と一様な電場 E=(E,0,0)\bm{E}=(E,0,0) の中にある。ベクトルポテンシャルを A=(0,Bx,0)\bm{A}=(0,Bx,0) にとり、Hamiltonian を

H^=12m(p^x2+(p^y−qeBx^)2)−qeEx^\hat H=\frac{1}{2m}\Bigl(\hat p_x^2+(\hat p_y-q_{\mathrm e}B\hat x)^2\Bigr)-q_{\mathrm e}E\hat x

とする。速度の演算子を mv^x=p^xm\hat v_x=\hat p_x、mv^y=p^y−qeBx^m\hat v_y=\hat p_y-q_{\mathrm e}B\hat x、サイクロトロン角振動数を ωc=qeB/m\omega_{\mathrm c}=q_{\mathrm e}B/m と書く。

  1. (1) [v^x, v^y]=iℏωc/m\left[\hat v_x,\,\hat v_y\right]=\mathrm{i}\hbar\omega_{\mathrm c}/m を示し、x^\hat x、y^\hat y、v^x\hat v_x、v^y\hat v_y の Heisenberg 方程式を書け。古典の Lorentz 力の運動方程式と比べよ。(使う道具:2.3.1 節・第2部第2章の処方)
  2. (2) (1) の方程式は線形である。u^=v^y+E/B\hat u=\hat v_y+E/B と置いて v^x(t)\hat v_x(t)、v^y(t)\hat v_y(t) を求め、速度の1周期の時間平均が (0,−E/B)(0,-E/B) であることを示せ。(使う道具:2.3.3 節)
  3. (3) 旋回中心 X^=x^+v^y/ωc\hat X=\hat x+\hat v_y/\omega_{\mathrm c}、Y^=y^−v^x/ωc\hat Y=\hat y-\hat v_x/\omega_{\mathrm c} について、X^˙=0\dot{\hat X}=0、Y^˙=−E/B\dot{\hat Y}=-E/B、[X^, Y^]=−iℏ/(qeB)\left[\hat X,\,\hat Y\right]=-\mathrm{i}\hbar/(q_{\mathrm e}B) を示せ。E=0E=0 のとき命題 2.3.5 から何が言えるか。(使う道具:2.3.2 節)
  4. (4) H^=m2(v^x2+u^2)−qeEX^−mE22B2\hat H=\frac m2(\hat v_x^2+\hat u^2)-q_{\mathrm e}E\hat X-\frac{mE^2}{2B^2} と書けること、X^\hat X が v^x\hat v_x、u^\hat u と交換することを示し、a^=m/(2ℏωc) (v^x+iu^)\hat a=\sqrt{m/(2\hbar\omega_{\mathrm c})}\,(\hat v_x+\mathrm{i}\hat u) を使ってエネルギーを求めよ。E=0E=0 のとき、1つの準位に属する状態の数が面積 SS あたり qeBS/(2πℏ)q_{\mathrm e}BS/(2\pi\hbar) であることを、[X^, Y^]\left[\hat X,\,\hat Y\right] を正準交換関係と見なして説明せよ。(使う道具:2.3.2 節・第2部第4章の調和振動子・第2部第1章)
  5. (5) GaAs の伝導電子(∣qe∣=e\lvert q_{\mathrm e}\rvert =e、m=0.067 mem=0.067\,m_{\mathrm e})で B=1 TB=1\,\mathrm T、E=103 V m−1E=10^3\,\mathrm{V\,m^{-1}} のとき、ℏωc\hbar\omega_{\mathrm c}、磁気長 ℓB=ℏ/(eB)\ell_B=\sqrt{\hbar/(eB)}、1準位あたりの面密度、ドリフトの速さ、mE2/(2B2)mE^2/(2B^2) を数で求めよ。(使う道具:(1)〜(4))
解答

(1) [p^x, x^]=−iℏ\left[\hat p_x,\,\hat x\right]=-\mathrm{i}\hbar で、p^y\hat p_y は x^\hat x、p^x\hat p_x と交換するから[mv^x, mv^y]=[p^x, −qeBx^]=iℏqeB\left[m\hat v_x,\,m\hat v_y\right]=\left[\hat p_x,\,-q_{\mathrm e}B\hat x\right]=\mathrm{i}\hbar q_{\mathrm e}B、すなわち [v^x, v^y]=iℏωc/m\left[\hat v_x,\,\hat v_y\right]=\mathrm{i}\hbar\omega_{\mathrm c}/m。H^=m2(v^x2+v^y2)−qeEx^\hat H=\frac m2(\hat v_x^2+\hat v_y^2)-q_{\mathrm e}E\hat x と書けるので、[x^, H^]=iℏv^x\left[\hat x,\,\hat H\right]=\mathrm{i}\hbar\hat v_x、[y^, H^]=iℏv^y\left[\hat y,\,\hat H\right]=\mathrm{i}\hbar\hat v_y([y^, p^y]=iℏ\left[\hat y,\,\hat p_y\right]=\mathrm{i}\hbar)。Leibniz 則から [v^x, v^y2]=2iℏωcv^y/m\left[\hat v_x,\,\hat v_y^2\right]=2\mathrm{i}\hbar\omega_{\mathrm c}\hat v_y/m、[v^x, −qeEx^]=iℏqeE/m\left[\hat v_x,\,-q_{\mathrm e}E\hat x\right]=\mathrm{i}\hbar q_{\mathrm e}E/m、[v^y, v^x2]=−2iℏωcv^x/m\left[\hat v_y,\,\hat v_x^2\right]=-2\mathrm{i}\hbar\omega_{\mathrm c}\hat v_x/m、[v^y, x^]=0\left[\hat v_y,\,\hat x\right]=0 なので

x^˙=v^x,y^˙=v^y,v^˙x=ωcv^y+qeEm,v^˙y=−ωcv^x.\dot{\hat x}=\hat v_x,\quad \dot{\hat y}=\hat v_y,\quad \dot{\hat v}_x=\omega_{\mathrm c}\hat v_y+\frac{q_{\mathrm e}E}{m},\quad \dot{\hat v}_y=-\omega_{\mathrm c}\hat v_x .

古典の mv˙=qe(E+v×B)m\dot{\bm{v}}=q_{\mathrm e}(\bm{E}+\bm{v}\times\bm{B}) は成分で mv˙x=qeE+qeBvym\dot v_x=q_{\mathrm e}E+q_{\mathrm e}Bv_y、mv˙y=−qeBvxm\dot v_y=-q_{\mathrm e}Bv_x で、帽子を外せば同じ式である。p^y\hat p_y と x^\hat x が交換するので、(p^y−qeBx^)2(\hat p_y-q_{\mathrm e}B\hat x)^2 に順序の曖昧さは無い(前の章の処方がそのまま当たる)。

(2) qeE/m=ωcE/Bq_{\mathrm e}E/m=\omega_{\mathrm c}E/B だから v^˙x=ωcu^\dot{\hat v}_x=\omega_{\mathrm c}\hat u、u^˙=−ωcv^x\dot{\hat u}=-\omega_{\mathrm c}\hat v_x で、(v^x,u^)(\hat v_x,\hat u) は角速度 ωc\omega_{\mathrm c} で回る。

v^x(t)=v^xcos⁡ωct+u^sin⁡ωct,v^y(t)=−EB+u^cos⁡ωct−v^xsin⁡ωct.\hat v_x(t)=\hat v_x\cos\omega_{\mathrm c}t+\hat u\sin\omega_{\mathrm c}t,\qquad \hat v_y(t)=-\frac EB+\hat u\cos\omega_{\mathrm c}t-\hat v_x\sin\omega_{\mathrm c}t .

1周期 2π/ωc2\pi/\omega_{\mathrm c} で平均すると三角関数の項が消え、速度の平均は (0,−E/B)(0,-E/B)、すなわち E×B/B2\bm{E}\times\bm{B}/B^2 の向きのドリフトになる。古典の解と同じ式で、状態によらず演算子として成り立つ。

(3) X^˙=v^x+v^˙y/ωc=v^x−v^x=0\dot{\hat X}=\hat v_x+\dot{\hat v}_y/\omega_{\mathrm c}=\hat v_x-\hat v_x=0、Y^˙=v^y−v^˙x/ωc=v^y−v^y−E/B=−E/B\dot{\hat Y}=\hat v_y-\dot{\hat v}_x/\omega_{\mathrm c}=\hat v_y-\hat v_y-E/B=-E/B。交換子は [X^, Y^]=[x^, −v^x]/ωc+[v^y, y^]/ωc−[v^y, v^x]/ωc2\left[\hat X,\,\hat Y\right]=\left[\hat x,\,-\hat v_x\right]/\omega_{\mathrm c}+\left[\hat v_y,\,\hat y\right]/\omega_{\mathrm c}-\left[\hat v_y,\,\hat v_x\right]/\omega_{\mathrm c}^2 で、[x^, v^x]=iℏ/m\left[\hat x,\,\hat v_x\right]=\mathrm{i}\hbar/m、[v^y, y^]=−iℏ/m\left[\hat v_y,\,\hat y\right]=-\mathrm{i}\hbar/m、[v^y, v^x]=−iℏωc/m\left[\hat v_y,\,\hat v_x\right]=-\mathrm{i}\hbar\omega_{\mathrm c}/m を入れると1ωc(−iℏm−iℏm+iℏm)=−iℏmωc=−iℏqeB\frac{1}{\omega_{\mathrm c}}\bigl(-\frac{\mathrm{i}\hbar}{m}-\frac{\mathrm{i}\hbar}{m}+\frac{\mathrm{i}\hbar}{m}\bigr)=-\frac{\mathrm{i}\hbar}{m\omega_{\mathrm c}}=-\frac{\mathrm{i}\hbar}{q_{\mathrm e}B}。E=0E=0 なら X^\hat X と Y^\hat Y はどちらも保存し、互いに交換しないので、命題 2.3.5 により準位は縮退する(面積が有限の試料で考えればスペクトルは離散的)。古典では、旋回中心がどこにあっても同じエネルギーの円運動ができることに当たる。

(4) x^=X^−v^y/ωc\hat x=\hat X-\hat v_y/\omega_{\mathrm c} を代入すると −qeEx^=−qeEX^+m(E/B)v^y-q_{\mathrm e}E\hat x=-q_{\mathrm e}E\hat X+m(E/B)\hat v_y で、m2v^y2+mEBv^y=m2u^2−mE22B2\frac m2\hat v_y^2+m\frac EB\hat v_y=\frac m2\hat u^2-\frac{mE^2}{2B^2} だから与えられた形になる。[X^, v^x]=[x^, v^x]+[v^y, v^x]/ωc=iℏ/m−iℏ/m=0\left[\hat X,\,\hat v_x\right]=\left[\hat x,\,\hat v_x\right]+\left[\hat v_y,\,\hat v_x\right]/\omega_{\mathrm c}=\mathrm{i}\hbar/m-\mathrm{i}\hbar/m=0、[X^, v^y]=0\left[\hat X,\,\hat v_y\right]=0(x^\hat x と v^y\hat v_y は交換)で、X^\hat X は u^\hat u とも交換する。[v^x+iu^, v^x−iu^]=−2i[v^x, u^]=2ℏωc/m\left[\hat v_x+\mathrm{i}\hat u,\,\hat v_x-\mathrm{i}\hat u\right]=-2\mathrm{i}\left[\hat v_x,\,\hat u\right]=2\hbar\omega_{\mathrm c}/m から [a^, a^†]=1\left[\hat a,\,\hat a^\dagger\right]=1、a^†a^=m2ℏωc(v^x2+u^2)−12\hat a^\dagger\hat a=\frac{m}{2\hbar\omega_{\mathrm c}}(\hat v_x^2+\hat u^2)-\frac12 だから m2(v^x2+u^2)=ℏωc(a^†a^+12)\frac m2(\hat v_x^2+\hat u^2)=\hbar\omega_{\mathrm c}(\hat a^\dagger\hat a+\frac12)。a^†a^\hat a^\dagger\hat a の固有値は n=0,1,2,…n=0,1,2,\dots(第2部第4章)、X^\hat X はそれと交換して固有値 XX を持つので

En,X=ℏωc(n+12)−qeEX−mE22B2.E_{n,X}=\hbar\omega_{\mathrm c}\Bigl(n+\frac12\Bigr)-q_{\mathrm e}EX-\frac{mE^2}{2B^2}.

E=0E=0 では XX によらない Landau 準位 ℏωc(n+12)\hbar\omega_{\mathrm c}(n+\frac12) になる。[X^, Y^]=−iℓ2\left[\hat X,\,\hat Y\right]=-\mathrm{i}\ell^2(ℓ2=ℏ/(qeB)\ell^2=\hbar/(q_{\mathrm e}B))は、ℏ\hbar を ℓ2\ell^2 に置き換えた正準交換関係である(符号は Y^\hat Y を座標、X^\hat X を運動量と読めばそろう)。1自由度の相空間で状態1つが面積 2πℏ2\pi\hbar を占めるのと同じく(第2部第1章の Bohr–Sommerfeld の数え方)、(X,Y)(X,Y) の平面では状態1つが面積 2πℓ22\pi\ell^2 を占めるので、面積 SS あたりの状態の数は S/(2πℓ2)=qeBS/(2πℏ)S/(2\pi\ell^2)=q_{\mathrm e}BS/(2\pi\hbar) である。

(5) ℏωc=ℏeB/m=1.73 meV\hbar\omega_{\mathrm c}=\hbar eB/m=1.73\,\mathrm{meV}、ℓB=25.7 nm\ell_B=25.7\,\mathrm{nm}、面密度 eB/(2πℏ)=eB/h=2.42×1014 m−2eB/(2\pi\hbar)=eB/h=2.42\times10^{14}\,\mathrm{m^{-2}}(スピン1つあたり)、ドリフトの速さ E/B=103 m s−1E/B=10^3\,\mathrm{m\,s^{-1}}、mE2/(2B2)=1.9×10−7 eVmE^2/(2B^2)=1.9\times10^{-7}\,\mathrm{eV}(ℏωc\hbar\omega_{\mathrm c} の約 10−410^{-4})。電子は qe=−e<0q_{\mathrm e}=-e<0 なので回る向きは逆になるが、ドリフトの向き E×B/B2\bm{E}\times\bm{B}/B^2 と、上の大きさは変わらない。

木下の勉強ステータス

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

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

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