tecktom.com

つなぎ:解析力学から様々な理論分野を組み立てる 第4部第11章から借りるものと、2バンドの公式

2026年10月15日12分で読めますAI生成教科書『虎の巻』第5部 第5章 つなぎ:Berry 曲率と異常 Hall 効果
目次
  1. 移したものと、ここで使う結論
  2. 当てる道具:2バンドの模型の Berry 曲率
  3. 木下の層

AI 生成・木下が未勉強

ハミルトニアンの選び方で固体を考える 第5部 第5章 §5.5.1 ★★ 借りるのは Berry 曲率の定義・異常速度・Hall 伝導度の式・Chern 数の整数性・対称性の制約の5つ。2バンドの模型では Berry 曲率が単位ベクトル n(k)\bm{n}(\bm{k}) の張る立体角で書け、Chern 数は n\bm{n} が球面を覆う回数になる

強磁性の金属に電流を流すと、外から磁場をかけなくても、電流と垂直な向きに電圧が現れる。これを異常 Hall 効果と呼ぶ。授業で習う半古典運動の式では、磁場が無ければ電子の速度は電場に沿った成分しか持たないので、この横向きの電流は出てこない。バンドの電子の波束が k\bm{k} 空間を動くときに Berry 位相を得て、その回転である Berry 曲率 Ωn(k)\Omega_n(\bm{k}) が速度に横向きの補正(異常速度)を加える、というのが答えである。この答えを1階の Lagrangian から導く部分は、解析力学の道具がそのまま当たるので、解析力学から様々な理論分野を組み立てる 第4部第11章へ移した。この章の最初の節は、そこから何を借りるかを短くまとめる「つなぎ」である(5.5.1節)。

残りの節で、この巻の主題である「どの H^\hat H を選び、何を削ったか」に沿って Berry 曲率の現れる現象を読む。選ぶ H^\hat H は、第2部第9章のグラフェンの tight-binding 模型に、副格子の準位の差と、時間反転を破る第2近接の複素数の飛び移りを足した Haldane 模型である(5.5.2節)。この模型では、磁場の正味の流束が 00 のまま Chern 数 11 の絶縁体ができ、Hall 伝導度が e2/he^2/h に量子化される。出てくるものの第一は、Chern 数が 00 でない絶縁体の端に、ギャップを横切る一方向にだけ進む状態が必ず現れることである(5.5.3節)。第二は、質量 00 のグラフェンでは Dirac 点のまわりの Berry 位相が π\pi になり、磁場の中の量子 Hall 効果が 4(N+12)e2/h4(N+\frac12)e^2/h という半整数の並びになることである(5.5.4節)。そのあと、削ったものを並べて近似の裏付けの札を貼り、強磁性金属の異常 Hall 効果と量子異常 Hall 効果の測定で確かめる(5.5.5節)。最後に授業の Hall 効果の式と並べ(5.5.6節)、まとめと総合問題で終える。スピンと軌道の結合で時間反転を保ったまま現れるトポロジカル絶縁体は、第5部第6章で扱う。

この章では、電子の電荷を −e-e(e>0e>0)とする。Hall 伝導度は jx=σxyEyj_x=\sigma_{xy}E_y で、Hall 抵抗率は電流を xx 向きに流して yy 向きの電場を測る ρyx=Ey/jx\rho_{yx}=E_y/j_x で定義する(この本の流儀)。解析力学から様々な理論分野を組み立てる 第4部第11章は同じ量を ρxy\rho_{xy} と書く(添字の順が逆なだけで量は同じ)。時間は時間微分の点 r˙\dot{\bm{r}}、k˙\dot{\bm{k}} としてしか現れないので、ホッピングは記法の正本どおり tt と書く。スカラーポテンシャルは出てこないので、ϕ\phi は Haldane 模型の位相に使う。

移したものと、ここで使う結論

バンド nn の Bloch 関数を ψnk(r)=eik⋅runk(r)\psi_{n\bm{k}}(\bm{r})=\mathrm{e}^{\mathrm{i}\bm{k}\cdot\bm{r}}u_{n\bm{k}}(\bm{r}) と書く。解析力学から様々な理論分野を組み立てる 第4部第11章で導いた結論のうち、この章で使うのは次の5つである。

1つめは定義で、Berry 接続と Berry 曲率を

An(k)=i⟨unk∣∇kunk⟩,Ωn(k)=∂An,y∂kx−∂An,x∂ky(5.5.1)\bm{A}_n(\bm{k})=\mathrm{i}\langle u_{n\bm{k}} | \nabla_{\bm{k}}u_{n\bm{k}} \rangle ,\qquad \Omega_n(\bm{k})=\frac{\partial A_{n,y}}{\partial k_x}-\frac{\partial A_{n,x}}{\partial k_y}\tag{5.5.1}

とする(2次元。内積は単位胞の上でとる)。An\bm{A}_n は unku_{n\bm{k}} の位相の選び方で変わるが、Ωn\Omega_n は変わらない。2つめは半古典運動で、波束の中心 r\bm{r} と波数 k\bm{k} は

ℏk˙=−e(E+r˙×B),r˙=1ℏ∇kεn(k)−k˙×Ωn(k)(5.5.2)\hbar\dot{\bm{k}}=-e\bigl(\bm{E}+\dot{\bm{r}}\times\bm{B}\bigr),\qquad \dot{\bm{r}}=\frac1\hbar\nabla_{\bm{k}}\varepsilon _n(\bm{k})-\dot{\bm{k}}\times\bm{\Omega}_n(\bm{k})\tag{5.5.2}

に従う(Ωn=Ωnez\bm{\Omega}_n=\Omega_n\bm{e}_z)。第2式の最後の項が異常速度で、電場の中では eℏE×Ωn\frac e\hbar\bm{E}\times\bm{\Omega}_n となって電場に垂直に向く。3つめは Hall 伝導度で、磁場の無い2次元の結晶の平衡状態について

σxy=−e2ℏ∑n∫BZd2k(2π)2 f0(εn(k)) Ωn(k)(5.5.3)\sigma_{xy}=-\frac{e^2}{\hbar}\sum_n\int_{\mathrm{BZ}}\frac{\mathrm{d}^2k}{(2\pi)^2}\,f_0(\varepsilon _n(\bm{k}))\,\Omega_n(\bm{k})\tag{5.5.3}

である。散乱の緩和時間を含まないこの部分を、異常 Hall 伝導度の内因的な部分と呼ぶ。4つめは量子化で、バンド nn がほかのバンドとギャップで隔たっているとき、Chern 数

Cn=12π∫BZΩn(k) d2k(5.5.4)C_n=\frac{1}{2\pi}\int_{\mathrm{BZ}}\Omega_n(\bm{k})\,\mathrm{d}^2k\tag{5.5.4}

は整数になり、Fermi 準位がギャップの中にある絶縁体では σxy=−(e2/h)∑n 詰まったCn\sigma_{xy}=-(e^2/h)\sum_{n\ \text{詰まった}}C_n となる。5つめは対称性の制約で、時間反転対称性があれば Ωn(−k)=−Ωn(k)\Omega_n(-\bm{k})=-\Omega_n(\bm{k})、空間反転対称性があれば Ωn(−k)=Ωn(k)\Omega_n(-\bm{k})=\Omega_n(\bm{k})、両方あれば Ωn≡0\Omega_n\equiv0 である。時間反転だけでも σxy=0\sigma_{xy}=0 になる。

この巻ではこの5つを結果として使い、導き直さない。式 (5.5.2) の第2式は、第2部第9章で使った半古典運動の式に異常速度を足したものである。第2部第9章では、グラフェンの Dirac 点のように2本のバンドが接する所では半古典運動の前提が崩れると書き、扱いをこの章へ送った。この章では、Dirac 点に質量を入れてギャップを開いた模型を主に扱い、ギャップが閉じる点は「Chern 数が変わる境目」として現れる。

当てる道具:2バンドの模型の Berry 曲率

この章の模型はどれも、波数ごとに 2×22\times2 の Bloch 行列

H(k)=g0(k) 1+g(k)⋅σ,ε±(k)=g0±∣g∣,n(k)≔g(k)∣g(k)∣(5.5.5)H(\bm{k})=g_0(\bm{k})\,\mathbb 1+\bm{g}(\bm{k})\cdot\bm{\sigma},\qquad\varepsilon _\pm(\bm{k})=g_0\pm\lvert \bm{g}\rvert ,\qquad\bm{n}(\bm{k})\coloneqq\frac{\bm{g}(\bm{k})}{\lvert \bm{g}(\bm{k})\rvert }\tag{5.5.5}

の形をしている。Pauli 行列 σ\bm{\sigma} は第2部第9章と同じく A・B の副格子の2成分(擬スピン)に作用する。g0g_0 は2つのバンドを一緒に上下させるだけで固有ベクトルを変えないので、Berry 曲率は単位ベクトル n(k)\bm{n}(\bm{k}) だけで決まる。

解析力学から様々な理論分野を組み立てる 第4部第9章で、スピン 12\frac12 の H^=R⋅σ\hat H=\bm{R}\cdot\bm{\sigma} の固有状態の Berry 曲率を、パラメータ R\bm{R} の空間で求めた。下の状態(固有値 −∣R∣-\lvert \bm{R}\rvert)では Ω−(R)=+R/(2∣R∣3)\bm{\Omega}_-(\bm{R})=+\bm{R}/(2\lvert \bm{R}\rvert ^3) で、原点に置いた単極子の場になる。R\bm{R} を g(k)\bm{g}(\bm{k}) とし、k\bm{k} 平面の微小な四辺形が g\bm{g} の空間に写る面要素 ∂kxg×∂kyg dkxdky\partial_{k_x}\bm{g}\times\partial_{k_y}\bm{g}\,\mathrm{d} k_x\mathrm{d} k_y を貫く流束を数えると、

Ω−(k)=g⋅(∂kxg×∂kyg)2∣g∣3=12 n⋅(∂kxn×∂kyn),Ω+(k)=−Ω−(k)(5.5.6)\Omega_-(\bm{k})=\frac{\bm{g}\cdot(\partial_{k_x}\bm{g}\times\partial_{k_y}\bm{g})}{2\lvert \bm{g}\rvert ^3}=\frac12\,\bm{n}\cdot\bigl(\partial_{k_x}\bm{n}\times\partial_{k_y}\bm{n}\bigr),\qquad\Omega_+(\bm{k})=-\Omega_-(\bm{k})\tag{5.5.6}

を得る。2つめの等号は、n⋅(∂n×∂n)\bm{n}\cdot(\partial\bm{n}\times\partial\bm{n}) を g=∣g∣n\bm{g}=\lvert \bm{g}\rvert \bm{n} で書き直すと、∣g∣\lvert \bm{g}\rvert の微分を含む項が n⋅(n×⋅)=0\bm{n}\cdot(\bm{n}\times\cdot)=0 で消えることから出る。n⋅(∂kxn×∂kyn) dkxdky\bm{n}\cdot(\partial_{k_x}\bm{n}\times\partial_{k_y}\bm{n})\,\mathrm{d} k_x\mathrm{d} k_y は、k\bm{k} の小さな四辺形を n\bm{n} で単位球面へ写した像の立体角である。したがって下のバンドの Chern 数は

C−=12π∫BZΩ− d2k=14π∫BZn⋅(∂kxn×∂kyn) d2k(5.5.7)C_-=\frac{1}{2\pi}\int_{\mathrm{BZ}}\Omega_-\,\mathrm{d}^2k=\frac{1}{4\pi}\int_{\mathrm{BZ}}\bm{n}\cdot\bigl(\partial_{k_x}\bm{n}\times\partial_{k_y}\bm{n}\bigr)\,\mathrm{d}^2k\tag{5.5.7}

で、Brillouin 域の全体を n\bm{n} で写したとき、単位球面(立体角 4π4\pi)を向きを込めて何回覆うかを数える。

第2部第8章の SSH 模型では、g\bm{g} が xyxy 平面の中にあり、kk が Brillouin 域を1周するあいだに g\bm{g} が原点のまわりを何回まわるか(巻き数)が端の状態の数を決めた。2次元では g\bm{g} が3次元の向きを持ち、巻き数の代わりに球面を覆う回数が現れる。g\bm{g} が k\bm{k} のどこでも 00 にならない(ギャップが開いている)限り n\bm{n} は滑らかに決まり、覆う回数は整数のまま連続には変えられない。C−C_- は、パラメータを動かす途中で g(k)=0\bm{g}(\bm{k})=0 となる点、つまりギャップが閉じる点が現れたときにだけ変わる。

例題 5.5.1 立体角の公式で質量のある Dirac 模型を解く ★★★

オリジナルg(k)=(ℏvkx, ℏvky, m)\bm{g}(\bm{k})=(\hbar vk_x,\ \hbar vk_y,\ m)(v>0v>0、mm は実数で 00 でない)とする。

  1. (1) 式 (5.5.6) から Ω−(k)\Omega_-(\bm{k}) を求め、解析力学から様々な理論分野を組み立てる 第4部第11章の例題で固有ベクトルから求めた ℏ2v2m/(2g3)\hbar^2v^2m/(2g^3)(g=∣g∣g=\lvert \bm{g}\rvert)と一致することを確かめよ。
  2. (2) k\bm{k} が平面全体を動くとき、n\bm{n} は単位球面のどの部分を覆うか。そこから ∫Ω− d2k\int\Omega_-\,\mathrm{d}^2k を積分せずに求めよ。
  3. (3) g=(−ℏvkx, ℏvky, m)\bm{g}=(-\hbar vk_x,\ \hbar vk_y,\ m) に変えると、(2) の答えはどう変わるか。
解答

(1) ∂kxg=(ℏv,0,0)\partial_{k_x}\bm{g}=(\hbar v,0,0)、∂kyg=(0,ℏv,0)\partial_{k_y}\bm{g}=(0,\hbar v,0) なので ∂kxg×∂kyg=(0,0,ℏ2v2)\partial_{k_x}\bm{g}\times\partial_{k_y}\bm{g}=(0,0,\hbar^2v^2)、g⋅(⋯ )=ℏ2v2m\bm{g}\cdot(\cdots)=\hbar^2v^2m である。式 (5.5.6) の最初の形に入れると Ω−=ℏ2v2m/(2g3)\Omega_-=\hbar^2v^2m/(2g^3) で、固有ベクトルから求めた値と一致する。固有ベクトルの位相の取り方を考えずに済むのが、この公式の利点である。

(2) m>0m>0 なら nz=m/g>0n_z=m/g>0 で、∣k∣→∞\lvert \bm{k}\rvert \to\infty で nz→0n_z\to0 になる。n\bm{n} の xyxy 成分は k\bm{k} と同じ向きを向くので、k\bm{k} の平面全体は北半球にちょうど1回、向きを保って写る。北半球の立体角は 2π2\pi なので、∫Ω− d2k=12×2π=π\int\Omega_-\,\mathrm{d}^2k=\frac12\times2\pi=\pi、12π∫Ω−d2k=12\frac1{2\pi}\int\Omega_-\mathrm{d}^2k=\frac12 である。m<0m<0 なら南半球を同じ向きに覆う。南半球を外側から見た向きで数えると、k\bm{k} の向きを保つ写像は向きを逆にした覆い方になり、−π-\pi、つまり −12-\frac12 になる。(1) の式で mm の符号を変えても同じ結論になる。

(3) kxk_x の符号を反転すると、k\bm{k} 平面を鏡に映してから写すことになり、覆う部分は同じ半球だが向きが逆になる。よって 12π∫Ω−d2k=−12sgn⁡m\frac1{2\pi}\int\Omega_-\mathrm{d}^2k=-\frac12\operatorname{sgn}m である。1つの Dirac 点は質量の符号と回り方(kxk_x の符号)の組で ±12\pm\frac12 を持つ。次の節では、この半分ずつを2つの谷で足して Chern 数を作る。(検算 scripts/verify_ch5-05.py)

ほかの巻では

木下の勉強ステータス

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

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

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