tecktom.com

この章のまとめ

2026年10月15日13分で読めますAI生成教科書『虎の巻』第4部 第10章 バンドの中の電子の半古典運動

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第4部 第10章 §4.10.7 ★★★ H=ε(p+eA)−eϕH=\varepsilon (\bm{p}+e\bm{A})-e\phi の Hamilton 方程式が半古典運動のすべて。有効質量は2階微分、Bloch 振動は帯の有界性、Onsager の量子化は k\bm{k} 空間の面積=作用変数

結晶の中の電子は、周期ポテンシャルの効果をバンド ε(k)\varepsilon (\bm{k}) にまとめると、Hamiltonian H=ε((p+eA)/ℏ)−eϕH=\varepsilon \bigl((\bm{p}+e\bm{A})/\hbar\bigr)-e\phi の1粒子になる。Hamilton 方程式から、速度は群速度 r˙=1ℏ∇kε\dot{\bm{r}}=\frac1\hbar\nabla_{\bm{k}}\varepsilon、運動方程式は ℏk˙=−e(E+r˙×B)\hbar\dot{\bm{k}}=-e(\bm{E}+\dot{\bm{r}}\times\bm{B}) になり、ゲージポテンシャルは消える。ℏk\hbar\bm{k} は運動量でなく、離散並進の演算子 T^R\hat T_{\bm{R}} の固有値 eik⋅R\mathrm{e}^{\mathrm{i}\bm{k}\cdot\bm{R}} から作った結晶運動量で、ℏG\hbar\bm{G} の分だけ不定である。古典力学では離散的な対称性から保存量は出ないが、量子力学ではユニタリな T^R\hat T_{\bm{R}} の固有値が保存する。一様な力の下では ℏk˙=F\hbar\dot{\bm{k}}=\bm{F} が厳密に成り立つ。

有効質量は H=ε(k)H=\varepsilon (\bm{k}) の2階微分の行列で、バンドの上の端では負になる。満ちたバンドは ε\varepsilon の周期性から電流を運ばず、上の端の近くの空いた状態は電荷 +e+e、正の質量の正孔として動く。一様な電場の下では kk が一定の速さで BZ を回り、電子は周期 TB=2πℏ/(eEa)T_{\mathrm B}=2\pi\hbar/(eEa)、振れ幅 ((帯の幅)/(eE))/(eE) で振動する。磁場の下では ℏkx\hbar k_x と ℏky\hbar k_y が共役な組になり、k\bm{k} 空間の閉軌道の面積 A\mathcal A が作用変数 J=ℏ2A/(2πeB)J=\hbar^2\mathcal A/(2\pi eB) になる。1周の時間からサイクロトロン質量 mc=ℏ22π∂A∂εm_{\mathrm c}=\frac{\hbar^2}{2\pi}\frac{\partial \mathcal A}{\partial \varepsilon }、Bohr–Sommerfeld の条件から Onsager の関係 A=2πeB(ν+γ)/ℏ\mathcal A=2\pi eB(\nu+\gamma)/\hbar と量子振動の周期が出る。

この章では Berry 曲率が 00 のバンドを考えた。Berry 曲率があると Poisson 括弧 (4.10.8) の {ri,rj}\{r_i,r_j\} が 00 でなくなり、速度に異常速度の項が加わる。それを第4部第11章(第4部第11章)で扱う。バンドの中の電子の集団の輸送は、この章の運動方程式を Boltzmann 方程式(第3部第8章)に入れて、固体の巻 第4部(ハミルトニアンの選び方で固体を考える 第4部第2章)で計算する。

この章の結果

  • 選ぶ Hamiltonian:H(r,p)=ε(k)−eϕ(r)H(\bm{r},\bm{p})=\varepsilon \bigl(\bm{k}\bigr)-e\phi(\bm{r})、ℏk=p+eA(r)\hbar\bm{k}=\bm{p}+e\bm{A}(\bm{r})(Peierls 置換)。
  • 半古典運動方程式:r˙=1ℏ∇kε\dot{\bm{r}}=\frac1\hbar\nabla_{\bm{k}}\varepsilon、ℏk˙=−e(E+r˙×B)\hbar\dot{\bm{k}}=-e(\bm{E}+\dot{\bm{r}}\times\bm{B})。{ℏki,ℏkj}=−e∑lϵijlBl\{\hbar k_i,\hbar k_j\}=-e\sum_l\epsilon_{ijl}B_l。
  • 結晶運動量:T^Rψnk=eik⋅Rψnk\hat T_{\bm{R}}\psi_{n\bm{k}}=\mathrm{e}^{\mathrm{i}\bm{k}\cdot\bm{R}}\psi_{n\bm{k}} の位相。G\bm{G} の分だけ不定。一様な力では ℏk˙=F\hbar\dot k=F が厳密。
  • 有効質量 (m∗−1)ij=1ℏ2∂ki∂kjε(m^{*-1})_{ij}=\frac1{\hbar^2}\partial_{k_i}\partial_{k_j}\varepsilon。tight-binding で m∗(0)=ℏ2/(2tha2)m^*(0)=\hbar^2/(2t_{\mathrm h}a^2)。満ちたバンドの電流は 00、正孔は電荷 +e+e。
  • Bloch 振動:TB=2πℏ/(eEa)T_{\mathrm B}=2\pi\hbar/(eEa)、x(t)=x0−[ε(k(t))−ε(k0)]/(eE)x(t)=x_0-[\varepsilon (k(t))-\varepsilon (k_0)]/(eE)、Wannier–Stark の間隔 eEaeEa。
  • 磁場:J=ℏ2A/(2πeB)J=\hbar^2\mathcal A/(2\pi eB)、mc=ℏ22π∂A∂εm_{\mathrm c}=\frac{\hbar^2}{2\pi}\frac{\partial \mathcal A}{\partial \varepsilon }、Aν=2πeB(ν+γ)/ℏ\mathcal A_\nu=2\pi eB(\nu+\gamma)/\hbar、Δ(1/B)=2πe/(ℏAext)\Delta(1/B)=2\pi e/(\hbar\mathcal A_{\mathrm{ext}})。

総合問題 2次元の正方格子の電子:有効質量・Bloch 振動・量子振動 ★★★

オリジナル格子間隔 aa の2次元の正方格子の tight-binding のバンド

ε(k)=−2th(cos⁡kxa+cos⁡kya),−πa<kx,ky≤πa\varepsilon (\bm{k})=-2t_{\mathrm h}\bigl(\cos k_xa+\cos k_ya\bigr),\qquad -\frac\pi a<k_x,k_y\le\frac\pi a

を考える(スピンは1つだけ数える)。BZ の中心を Γ=(0,0)\Gamma=(0,0)、辺の中点を X=(π/a,0)\mathrm X=(\pi/a,0)、角を M=(π/a,π/a)\mathrm M=(\pi/a,\pi/a) と書く。

  1. (1) (使う道具:4.10.2節・4.10.3節、第1部第8章の Hamilton 方程式)群速度 v(k)\bm{v}(\bm{k}) を求めよ。Γ\Gamma・X・M での有効質量テンソルの逆数 (m∗−1)ij(m^{*-1})_{ij} を求め、それぞれの点がバンドの底・鞍点・上の端のどれかを答えよ。
  2. (2) (使う道具:4.10.4節)Γ\Gamma にいた電子に、xx 軸に沿った一様な電場 EE をかけたときと、対角線 (1,1)/2(1,1)/\sqrt2 に沿った電場 EE をかけたときで、Bloch 振動の周期と、場に沿った方向の位置の振れ幅を比べよ。
  3. (3) (使う道具:4.10.5節)B=Bez\bm{B}=B\bm{e}_z をかける。バンドの底の近くで ε≈−4th+tha2∣k∣2\varepsilon \approx-4t_{\mathrm h}+t_{\mathrm h}a^2\lvert \bm{k}\rvert ^2 を示し、等エネルギー線の面積 A(ε)\mathcal A(\varepsilon ) とサイクロトロン質量 mcm_{\mathrm c} を求めよ。mcm_{\mathrm c} を (1) の Γ\Gamma の有効質量と比べよ。
  4. (4) (使う道具:4.10.5節、第1部第14章の作用変数)ε=0\varepsilon =0 の等エネルギー線が ∣kx∣+∣ky∣=π/a\lvert k_x\rvert +\lvert k_y\rvert =\pi/a の正方形になることを示し、その面積と、a=0.40 nma=0.40\,\mathrm{nm} のときの量子振動の BdHB_{\mathrm{dH}} を求めよ。また ε=+th\varepsilon =+t_{\mathrm h} の軌道はどの点を囲み、どちら向きに回るか。
  5. (5) (使う道具:(4)・4.10.5節、第4部第4章の総合問題の縮退度)εF=0\varepsilon _{\mathrm F}=0 まで電子が詰まっているとき(バンドのちょうど半分)、B=10 TB=10\,\mathrm T で εF\varepsilon _{\mathrm F} より下にある Landau 準位の数 BdH/BB_{\mathrm{dH}}/B を求め、それに1つの準位あたりの面積あたりの状態数 eB/heB/h を掛けると電子の面密度 1/(2a2)1/(2a^2) に一致することを確かめよ。
解答

(1) Hamilton 方程式 (4.10.5) から

v(k)=2thaℏ(sin⁡kxa, sin⁡kya).\bm{v}(\bm{k})=\frac{2t_{\mathrm h}a}{\hbar}\bigl(\sin k_xa,\ \sin k_ya\bigr).

2階微分は対角で ∂2ε∂kx2=2tha2cos⁡kxa\frac{\partial^2\varepsilon }{\partial k_x^2}=2t_{\mathrm h}a^2\cos k_xa、∂2ε∂ky2=2tha2cos⁡kya\frac{\partial^2\varepsilon }{\partial k_y^2}=2t_{\mathrm h}a^2\cos k_ya、交差項は 00 なので

Γ: 2tha2ℏ2(1001),X: 2tha2ℏ2(−1001),M: 2tha2ℏ2(−100−1).\Gamma:\ \frac{2t_{\mathrm h}a^2}{\hbar^2}\begin{pmatrix}1&0\\0&1\end{pmatrix},\qquad \mathrm X:\ \frac{2t_{\mathrm h}a^2}{\hbar^2}\begin{pmatrix}-1&0\\0&1\end{pmatrix},\qquad \mathrm M:\ \frac{2t_{\mathrm h}a^2}{\hbar^2}\begin{pmatrix}-1&0\\0&-1\end{pmatrix}.

Γ\Gamma は両方向に正の質量 ℏ2/(2tha2)\hbar^2/(2t_{\mathrm h}a^2) の底(ε=−4th\varepsilon =-4t_{\mathrm h})、M は両方向に負の質量の上の端(ε=+4th\varepsilon =+4t_{\mathrm h})、X は xx 方向に負・yy 方向に正の鞍点(ε=0\varepsilon =0)である。どの点でも v=0\bm{v}=0 で、ε\varepsilon の停留点になっている。

(2) xx 方向の電場では ℏk˙x=−eE\hbar\dot k_x=-eE、ky=0k_y=0 のままで、運動は1次元の場合と同じになる。周期は TB=2πℏ/(eEa)T_{\mathrm B}=2\pi\hbar/(eEa)、ε(kx,0)=−2thcos⁡kxa−2th\varepsilon (k_x,0)=-2t_{\mathrm h}\cos k_xa-2t_{\mathrm h} は幅 4th4t_{\mathrm h} を動くので、振れ幅は 4th/(eE)4t_{\mathrm h}/(eE) である。

対角の電場では ℏk˙x=ℏk˙y=−eE/2\hbar\dot k_x=\hbar\dot k_y=-eE/\sqrt2 で、kx=ky=−eEt/(2ℏ)k_x=k_y=-eEt/(\sqrt2\hbar)。ε\varepsilon が元に戻るのは kxk_x が 2π/a2\pi/a 進んだときで、周期は2π2ℏ/(eEa)=2 TB2\pi\sqrt2\hbar/(eEa)=\sqrt2\,T_{\mathrm B} である。対角線に沿った格子の周期は 2a\sqrt2a なので、式 (4.10.15) の aa を 2a\sqrt2a に替えた形と言うこともできる(ただし k\bm{k} が Γ\Gamma から M までの全部を通る点が違う)。ε(k,k)=−4thcos⁡ka\varepsilon (k,k)=-4t_{\mathrm h}\cos ka は幅 8th8t_{\mathrm h} を動くので、場に沿った方向の振れ幅はエネルギー保存 ε+eE s=\varepsilon +eE\,s= 一定(ss は場に沿った座標)から 8th/(eE)8t_{\mathrm h}/(eE) で、xx 方向の場合の2倍になる。

(3) cos⁡u≈1−u2/2\cos u\approx1-u^2/2 から ε≈−2th(2−a2kx2/2−a2ky2/2)=−4th+tha2∣k∣2\varepsilon \approx-2t_{\mathrm h}(2-a^2k_x^2/2-a^2k_y^2/2)=-4t_{\mathrm h}+t_{\mathrm h}a^2\lvert \bm{k}\rvert ^2。等エネルギー線は半径 k=(ε+4th)/(tha2)k=\sqrt{(\varepsilon +4t_{\mathrm h})/(t_{\mathrm h}a^2)} の円で

A(ε)=πε+4ththa2,mc=ℏ22π∂A∂ε=ℏ22tha2.\mathcal A(\varepsilon )=\pi\frac{\varepsilon +4t_{\mathrm h}}{t_{\mathrm h}a^2},\qquad m_{\mathrm c}=\frac{\hbar^2}{2\pi}\frac{\partial \mathcal A}{\partial \varepsilon }=\frac{\hbar^2}{2t_{\mathrm h}a^2}.

底の近くでは、サイクロトロン質量は (1) の Γ\Gamma の有効質量に一致する。等方的な放物線のバンドでは、軌道全体で決まる mcm_{\mathrm c} と1点で決まる m∗m^* が同じになるからである。エネルギーが上がると放物近似は外れる。たとえば ε=−3th\varepsilon =-3t_{\mathrm h} の厳密な面積は数値積分で 3.36/a23.36/a^2 で、放物近似の π/a2\pi/a^2 より 7% 大きい(図 4.10.3)。

(4) ε=0\varepsilon =0 は cos⁡kxa=−cos⁡kya=cos⁡(π−kya)\cos k_xa=-\cos k_ya=\cos(\pi-k_ya) で、BZ の中では kxa=±(π−∣ky∣a)k_xa=\pm(\pi-\lvert k_y\rvert a)、つまり ∣kx∣+∣ky∣=π/a\lvert k_x\rvert +\lvert k_y\rvert =\pi/a の正方形である。対角線の長さ 2π/a2\pi/a の正方形なので面積は A=12(2π/a)2=2π2/a2\mathcal A=\frac12(2\pi/a)^2=2\pi^2/a^2 で、BZ の面積 4π2/a24\pi^2/a^2 のちょうど半分である。量子振動の周期は式 (4.10.23) から

BdH=ℏ2πe⋅2π2a2=πℏea2=3.1416×1.0546×10−341.602×10−19×(0.40×10−9)2=1.29×104 T.B_{\mathrm{dH}}=\frac{\hbar}{2\pi e}\cdot\frac{2\pi^2}{a^2}=\frac{\pi\hbar}{ea^2}=\frac{3.1416\times1.0546\times10^{-34}}{1.602\times10^{-19}\times(0.40\times10^{-9})^2}=1.29\times10^4\,\mathrm T.

ε=+th\varepsilon =+t_{\mathrm h}(cos⁡kxa+cos⁡kya=−12\cos k_xa+\cos k_ya=-\frac12)の軌道は M を囲む(BZ の4隅にある曲線が周期性 (4.10.1) で1つの閉曲線につながる)。k→k+(π/a,π/a)\bm{k}\to\bm{k}+(\pi/a,\pi/a) で ε→−ε\varepsilon \to-\varepsilon となるので、M のまわりの ε=+th\varepsilon =+t_{\mathrm h} の軌道の面積は、Γ\Gamma のまわりの ε=−th\varepsilon =-t_{\mathrm h} の軌道の面積に等しい。回る向きは、M の近くで m∗<0m^*<0 だから速度が k−kM\bm{k}-\bm{k}_{\mathrm M} と逆向きで、ℏk˙=−ev×B\hbar\dot{\bm{k}}=-e\bm{v}\times\bm{B} の向きも逆になり、+z+z 側から見て時計回りである(正孔的な軌道)。

(5) BdH/B=1.29×104/10=1.29×103B_{\mathrm{dH}}/B=1.29\times10^4/10=1.29\times10^3 本の準位が εF\varepsilon _{\mathrm F} より下にある。1つの準位あたりの面積あたりの状態数はeB/h=1.602×10−19×10/6.626×10−34=2.42×1015 m−2eB/h=1.602\times10^{-19}\times10/6.626\times10^{-34}=2.42\times10^{15}\,\mathrm m^{-2} で、掛けると

BdHB⋅eBh=πℏea2⋅e2πℏ=12a2=3.13×1018 m−2\frac{B_{\mathrm{dH}}}{B}\cdot\frac{eB}{h}=\frac{\pi\hbar}{ea^2}\cdot\frac{e}{2\pi\hbar}=\frac{1}{2a^2}=3.13\times10^{18}\,\mathrm m^{-2}

となる。BZ に 1/a21/a^2(面積あたり)の状態があり、その半分が詰まっているので、電子の面密度は確かに 1/(2a2)1/(2a^2) である。BB が式から消えるのは、磁場が強いほど1つの準位に入る状態が増え、そのぶん εF\varepsilon _{\mathrm F} より下の準位の数が減るからで、Onsager の関係が状態の数を正しく数えていることの確かめになっている。

木下の勉強ステータス

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

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

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