ハミルトニアンの選び方で固体を考える 第5部 第5章 §5.5.4 ★★★ 質量 0 0 0 の Dirac コーンを1周すると Berry 位相は π \pi π で、Onsager の量子化の位相のずれが消え、ε = 0 \varepsilon =0 ε = 0 に Landau 準位が1つ残る。この準位を電子と正孔で半分ずつ分けるので、Hall 伝導度は 4 ( N + 1 2 ) e 2 / h 4(N+\frac12)e^2/h 4 ( N + 2 1 ) e 2 / h の半整数の並びになる
Dirac 点のまわりの Berry 位相
Haldane 模型で Δ = 0 \Delta=0 Δ = 0 、t 2 = 0 t_2=0 t 2 = 0 とすると、質量 0 0 0 のグラフェンに戻る。このとき Berry 曲率は Dirac 点に集中した δ \delta δ 関数になる。例題 5.5.1 の質量 m m m を 0 0 0 に近づけると、Ω − \Omega_- Ω − は幅 ∣ m ∣ / ( ℏ v ) \lvert m\rvert /(\hbar v) ∣ m ∣ / ( ℏ v ) の領域に集まり、流束は ± π \pm\pi ± π のまま変わらない。したがって K のまわりの閉じた等エネルギー線を1周する Berry 位相は ± π \pm\pi ± π である。π \pi π と − π -\pi − π は位相としては同じなので、質量を正に取っても負に取っても、円を1周すると固有ベクトルの符号が反転する。これは擬スピンが q \bm{q} q の向きに沿って x y xy x y 平面の中で1回転し、スピン 1 2 \frac12 2 1 を 2 π 2\pi 2 π 回すと − 1 -1 − 1 倍になるのと同じ事情である。
Landau 準位
解析力学から様々な理論分野を組み立てる 第4部第11章 で、Berry 位相 Γ n \Gamma_n Γ n を持つ軌道の Onsager の量子化がA ν = 2 π e B ℏ ( ν + 1 2 − Γ n 2 π ) \mathcal A_\nu=\frac{2\pi eB}{\hbar}\bigl(\nu+\frac12-\frac{\Gamma_n}{2\pi}\bigr) A ν = ℏ 2 π e B ( ν + 2 1 − 2 π Γ n ) になることを見た。A ν \mathcal A_\nu A ν は k \bm{k} k 空間の軌道が囲む面積で、ν = 0 , 1 , 2 , … \nu=0,1,2,\dots ν = 0 , 1 , 2 , … である。質量 0 0 0 の Dirac コーンでは Γ n = ± π \Gamma_n=\pm\pi Γ n = ± π なので、位相のずれ 1 2 − Γ n 2 π \frac12-\frac{\Gamma_n}{2\pi} 2 1 − 2 π Γ n は 0 0 0 (または 1 1 1 。ν \nu ν の付け替えで同じ)になる。半径 q q q の円の面積 π q 2 \pi q^2 π q 2 と ε = ℏ v q \varepsilon =\hbar vq ε = ℏ v q を入れると
ε ν = sgn ( ν ) 2 e ℏ v 2 B ∣ ν ∣ , ν = 0 , ± 1 , ± 2 , … (5.5.12) \varepsilon _\nu=\operatorname{sgn}(\nu)\sqrt{2e\hbar v^2B\lvert \nu\rvert },\qquad\nu=0,\pm1,\pm2,\dots\tag{5.5.12} ε ν = sgn ( ν ) 2 e ℏ v 2 B ∣ ν ∣ , ν = 0 , ± 1 , ± 2 , … ( 5.5.12 )
を得る(ν < 0 \nu<0 ν < 0 は下のコーンの準位)。解析力学から様々な理論分野を組み立てる 第4部第11章 は上のコーンだけを書いたが、下のコーンでも同じ量子化が成り立つ。式 (5.5.12) は Dirac 方程式の Landau 準位を厳密に解いた結果と一致する(verify_ch5-05.py で調和振動子の基底の行列を対角化して確かめた)。普通の放物線のバンドの Landau 準位 ℏ ω c ( ν + 1 2 ) \hbar\omega_{\mathrm c}(\nu+\frac12) ℏ ω c ( ν + 2 1 ) との違いは2つある。1つは間隔が等しくなく ∣ ν ∣ B \sqrt{\lvert \nu\rvert B} ∣ ν ∣ B に比例すること、もう1つは ε = 0 \varepsilon =0 ε = 0 にちょうど1つの準位(ν = 0 \nu=0 ν = 0 )があることである。
1つの準位に入る状態の数は、面積あたり e B / h eB/h e B / h にスピンの2と谷の2を掛けた 4 e B / h 4eB/h 4 e B / h である。ν = 0 \nu=0 ν = 0 の準位は上のコーンと下のコーンの境目にあり、電子と正孔が入れ替わる対称性(ε → − ε \varepsilon \to-\varepsilon ε → − ε )から、電荷を中性にしたとき半分だけ満たされている。
半整数の Hall 伝導度
電荷が中性の点から測った電子の面密度を n e n_{\mathrm e} n e とする。Fermi 準位が ν = N \nu=N ν = N と N + 1 N+1 N + 1 の準位の間(N ≥ 0 N\ge0 N ≥ 0 )にあるとき、ν = 0 \nu=0 ν = 0 の準位の半分と、ν = 1 , … , N \nu=1,\dots,N ν = 1 , … , N の準位が中性の点から足された分なので
n e = 4 e B h ( N + 1 2 ) n_{\mathrm e}=\frac{4eB}{h}\Bigl(N+\frac12\Bigr) n e = h 4 e B ( N + 2 1 )
である。解析力学から様々な理論分野を組み立てる 第4部第11章 の Středa の式 σ x y = − e ( ∂ n e / ∂ B ) μ \sigma_{xy}=-e(\partial n_{\mathrm e}/\partial B)_\mu σ x y = − e ( ∂ n e / ∂ B ) μ を当てると
σ x y = − 4 e 2 h ( N + 1 2 ) = − e 2 h × 2 , 6 , 10 , … (5.5.13) \sigma_{xy}=-\frac{4e^2}{h}\Bigl(N+\frac12\Bigr)=-\frac{e^2}{h}\times2,\ 6,\ 10,\ \dots\tag{5.5.13} σ x y = − h 4 e 2 ( N + 2 1 ) = − h e 2 × 2 , 6 , 10 , … ( 5.5.13 )
となる。正孔の側(n e < 0 n_{\mathrm e}<0 n e < 0 )では符号が逆の + 2 , + 6 , … +2,+6,\dots + 2 , + 6 , … になる。2次元の電子気体の量子 Hall 効果が ν e 2 / h \nu e^2/h ν e 2 / h (ν \nu ν は整数)であるのに対し、グラフェンでは e 2 / h e^2/h e 2 / h の単位で 4 4 4 ずつ、半分ずれた並びになる(図 5.5.3 )。
この「半分」は谷ごとの Chern 数で読むことができる。磁場の中の Landau 準位はそれぞれ Chern 数 1 1 1 のバンドで、普通の電子気体では満ちた準位の数がそのまま Chern 数の和になる。グラフェンの1つの Dirac コーンは、例題 5.5.1 で見たように、それ自体が ± 1 2 \pm\frac12 ± 2 1 の Chern 数を持つ(質量を 0 0 0 に近づけた極限で)。1つの谷、1つのスピンあたりの Hall 伝導度は、ν = 0 \nu=0 ν = 0 の準位が空のとき + 1 2 e 2 / h +\frac12\,e^2/h + 2 1 e 2 / h 、満ちたとき − 1 2 e 2 / h -\frac12\,e^2/h − 2 1 e 2 / h で、そこから上の準位を1つ満たすごとに − e 2 / h -e^2/h − e 2 / h ずつ変わる。スピンと谷の4倍で式 (5.5.13) になる。Haldane 模型では、磁場の代わりに t 2 sin ϕ t_2\sin\phi t 2 sin ϕ によって2つの谷の質量の符号が逆になり、谷ごとの 1 2 \frac12 2 1 が足し合って磁場なしで 1 1 1 になった。どちらも、Dirac コーン1つが持つ半整数が見えている。
図 5.5.3 グラフェンの Hall 伝導度の模式図(温度 0 0 0 、乱れによって準位の中央付近の状態だけが電流を運ぶとした理想化)。横軸は電荷が中性の点から測った電子の面密度を1つの準位の状態の数 4 e B / h 4eB/h 4 e B / h で割ったもの。点線の台地は ± 2 , ± 6 , ± 10 , … \pm2,\pm6,\pm10,\dots ± 2 , ± 6 , ± 10 , … に並び、各台地の中央は n e / ( 4 e B / h ) = N + 1 2 n_{\mathrm e}/(4eB/h)=N+\frac12 n e / ( 4 e B / h ) = N + 2 1 にあって、n e = 0 n_{\mathrm e}=0 n e = 0 をまたいで − 2 -2 − 2 から + 2 +2 + 2 へ跳ぶ。実線は古典の Hall 効果 σ x y = − n e e / B \sigma_{xy}=-n_{\mathrm e}e/B σ x y = − n e e / B で、台地の中央を通る。
半整数の並びは 2005 年に2つのグループが独立に測った(K. S. Novoselov et al. (2005), Nature 438 , 197;Y. Zhang, Y.-W. Tan, H. L. Stormer and P. Kim (2005), Nature 438 , 201)。ν = 1 \nu=1 ν = 1 の準位と ν = 0 \nu=0 ν = 0 の準位の間隔は普通の半導体の Landau 準位の間隔より大きく、グラフェンでは室温でも量子 Hall 効果が見えることが報告されている(K. S. Novoselov et al. (2007), Science 315 , 1379)。
例題 5.5.4 グラフェンの Landau 準位の大きさ ★★
オリジナル v = 1.0 × 10 6 m s − 1 v=1.0\times10^6\,\mathrm{m\,s^{-1}} v = 1.0 × 1 0 6 m s − 1 、B = 10 T B=10\,\mathrm T B = 10 T とする。
(1) ε 1 \varepsilon _1 ε 1 を eV と温度(ε 1 / k B \varepsilon _1/ k_{\mathrm B} ε 1 / k B )で求め、ε 2 − ε 1 \varepsilon _2-\varepsilon _1 ε 2 − ε 1 と比べよ。
(2) GaAs の2次元電子(有効質量 m ∗ = 0.067 m e m^*=0.067m_{\mathrm e} m ∗ = 0.067 m e )の Landau 準位の間隔 ℏ ω c = ℏ e B / m ∗ \hbar\omega_{\mathrm c}=\hbar eB/m^* ℏ ω c = ℏ e B / m ∗ と比べよ。
(3) 最初の台地 σ x y = − 2 e 2 / h \sigma_{xy}=-2e^2/h σ x y = − 2 e 2 / h の中央での電子の面密度と、そのときの ∣ ρ y x ∣ \lvert \rho_{yx}\rvert ∣ ρ y x ∣ を求めよ。
解答 (1) ε 1 = v 2 e ℏ B = 1.0 × 10 6 × 2 × 1.602 × 10 − 19 × 1.0546 × 10 − 34 × 10 = 1.84 × 10 − 20 J = 0.115 e V \varepsilon _1=v\sqrt{2e\hbar B}=1.0\times10^6\times\sqrt{2\times1.602\times10^{-19}\times1.0546\times10^{-34}\times10}=1.84\times10^{-20}\,\mathrm J=0.115\,\mathrm{eV} ε 1 = v 2 e ℏ B = 1.0 × 1 0 6 × 2 × 1.602 × 1 0 − 19 × 1.0546 × 1 0 − 34 × 10 = 1.84 × 1 0 − 20 J = 0.115 eV 、ε 1 / k B = 1.33 × 10 3 K \varepsilon _1/ k_{\mathrm B}=1.33\times10^3\,\mathrm K ε 1 / k B = 1.33 × 1 0 3 K である。ε 2 = 2 ε 1 \varepsilon _2=\sqrt2\,\varepsilon _1 ε 2 = 2 ε 1 なので ε 2 − ε 1 = 0.414 ε 1 = 0.0475 e V \varepsilon _2-\varepsilon _1=0.414\,\varepsilon _1=0.0475\,\mathrm{eV} ε 2 − ε 1 = 0.414 ε 1 = 0.0475 eV で、準位の間隔はエネルギーが上がるほど狭くなる。
(2) ℏ ω c = 1.0546 × 10 − 34 × 1.602 × 10 − 19 × 10 / ( 0.067 × 9.109 × 10 − 31 ) = 2.77 × 10 − 21 J = 17.3 m e V \hbar\omega_{\mathrm c}=1.0546\times10^{-34}\times1.602\times10^{-19}\times10/(0.067\times9.109\times10^{-31})=2.77\times10^{-21}\,\mathrm J=17.3\,\mathrm{meV} ℏ ω c = 1.0546 × 1 0 − 34 × 1.602 × 1 0 − 19 × 10/ ( 0.067 × 9.109 × 1 0 − 31 ) = 2.77 × 1 0 − 21 J = 17.3 meV で、ε 1 \varepsilon _1 ε 1 はその 6.6 6.6 6.6 倍である。Dirac コーンでは ε = 0 \varepsilon =0 ε = 0 の近くの状態の密度が小さい(第2部第9章 の D ( ε ) ∝ ∣ ε ∣ D(\varepsilon )\propto\lvert \varepsilon \rvert D ( ε ) ∝ ∣ ε ∣ )ので、同じ数の状態を集めて1つの準位にするために広いエネルギーの幅が要る。そのため最初の間隔が大きい。
(3) 台地 N = 0 N=0 N = 0 の中央は ν = 0 \nu=0 ν = 0 と ν = 1 \nu=1 ν = 1 の準位の間の中央で、面密度は n e = 4 e B h ⋅ 1 2 = 2 e B h = 4.84 × 10 15 m − 2 n_{\mathrm e}=\frac{4eB}{h}\cdot\frac12=\frac{2eB}{h}=4.84\times10^{15}\,\mathrm{m^{-2}} n e = h 4 e B ⋅ 2 1 = h 2 e B = 4.84 × 1 0 15 m − 2 である。σ x x = 0 \sigma_{xx}=0 σ xx = 0 なので ∣ ρ y x ∣ = 1 / ∣ σ x y ∣ = h / ( 2 e 2 ) = 12.9 k Ω \lvert \rho_{yx}\rvert =1/\lvert \sigma_{xy}\rvert =h/(2e^2)=12.9\,\mathrm{k}\Omega ∣ ρ y x ∣ = 1/ ∣ σ x y ∣ = h / ( 2 e 2 ) = 12.9 k Ω である。(検算 scripts/verify_ch5-05.py)
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com