解析力学から様々な理論分野を組み立てる 第4部 第10章 §4.10.4 ★★ 一様な電場の下で k k k は一定の速さで BZ を回り、電子は位置が振動する。周期 T B = 2 π ℏ / ( e E a ) T_{\mathrm B}=2\pi\hbar/(eEa) T B = 2 π ℏ/ ( e E a ) 、振れ幅は帯の幅/ ( e E ) /(eE) / ( e E )
k k k は一定の速さで進む
1次元の結晶に + x +x + x 向きの一様な電場 E E E をかける。ϕ = − E x \phi=-Ex ϕ = − E x で、式 (4.10.6) は ℏ k ˙ = − e E \hbar\dot k=-eE ℏ k ˙ = − e E だから
k ( t ) = k 0 − e E ℏ t (4.10.14) k(t)=k_0-\frac{eE}{\hbar}t\tag{4.10.14} k ( t ) = k 0 − ℏ e E t ( 4.10.14 )
である。自由な電子なら速度も一定の割合で増え続ける。しかしバンドの中では ε ( k ) \varepsilon (k) ε ( k ) と v ( k ) v(k) v ( k ) が k k k の周期関数なので、k k k が 2 π / a 2\pi/a 2 π / a 進むたびに速度は元に戻る。k k k が BZ を1周する時間は
T B = 2 π / a e E / ℏ = 2 π ℏ e E a = h e E a , ω B = 2 π T B = e E a ℏ (4.10.15) T_{\mathrm B}=\frac{2\pi/a}{eE/\hbar}=\frac{2\pi\hbar}{eEa}=\frac{h}{eEa},\qquad \omega_{\mathrm B}=\frac{2\pi}{T_{\mathrm B}}=\frac{eEa}{\hbar}\tag{4.10.15} T B = e E /ℏ 2 π / a = e E a 2 π ℏ = e E a h , ω B = T B 2 π = ℏ e E a ( 4.10.15 )
で、これを Bloch 振動の周期と呼ぶ。
位置はエネルギー保存から直接求まる
位置 x ( t ) x(t) x ( t ) は x ˙ = v ( k ( t ) ) \dot x=v(k(t)) x ˙ = v ( k ( t )) を積分すれば求まるが、エネルギー保存 (4.10.7) を使うと積分しなくてよい。− e ϕ = e E x -e\phi=eEx − e ϕ = e E x なのでε ( k ) + e E x \varepsilon (k)+eEx ε ( k ) + e E x が一定で、
x ( t ) = x 0 − ε ( k ( t ) ) − ε ( k 0 ) e E (4.10.16) x(t)=x_0-\frac{\varepsilon (k(t))-\varepsilon (k_0)}{eE}\tag{4.10.16} x ( t ) = x 0 − e E ε ( k ( t )) − ε ( k 0 ) ( 4.10.16 )
となる。ε \varepsilon ε は幅 4 t h 4t_{\mathrm h} 4 t h の帯に収まるので、x x x も幅 4 t h / ( e E ) 4t_{\mathrm h}/(eE) 4 t h / ( e E ) の範囲に閉じ込められる。電子は電場の中で一方向に流れず、同じ所を往復する。tight-binding で k 0 = 0 k_0=0 k 0 = 0 から出発すると ε ( k ( t ) ) − ε ( 0 ) = 2 t h ( 1 − cos ω B t ) \varepsilon (k(t))-\varepsilon (0)=2t_{\mathrm h}\bigl(1-\cos\omega_{\mathrm B}t\bigr) ε ( k ( t )) − ε ( 0 ) = 2 t h ( 1 − cos ω B t ) で、
x ( t ) = x 0 − 2 t h e E ( 1 − cos ω B t ) x(t)=x_0-\frac{2t_{\mathrm h}}{eE}\bigl(1-\cos\omega_{\mathrm B}t\bigr) x ( t ) = x 0 − e E 2 t h ( 1 − cos ω B t )
である(図 4.10.2 )。電子はまず力の向き(− x -x − x )へ進み、k k k が BZ の端を越えたあたりで向きを変えて戻ってくる。力の向きへ進む間はポテンシャルエネルギー e E x eEx e E x を失ってバンドのエネルギー ε \varepsilon ε を得る。ε \varepsilon ε がバンドの上の端に達すると、それ以上エネルギーを受け取る先が無いので、電子は戻るしかない。
図 4.10.2 一様な電場の中の Bloch 振動(1次元の tight-binding、k 0 = 0 k_0=0 k 0 = 0 )。(a) k k k は一定の速さで減り、BZ の端 − π / a -\pi/a − π / a に達すると + π / a +\pi/a + π / a に戻る(同じ状態の別名)。(b) 位置は周期 T B T_{\mathrm B} T B で振動し、振れ幅は 4 t h / ( e E ) 4t_{\mathrm h}/(eE) 4 t h / ( e E ) 。式 (4.10.16) と、Hamilton 方程式を数値積分した軌道(verify_ch4-10.py)は 10 − 6 10^{-6} 1 0 − 6 の精度で一致する
Wannier–Stark の階段
Bloch 振動は周期 T B T_{\mathrm B} T B の周期運動なので、第1部第14章の言葉で言えば、作用角変数を当てて角振動数 ω B \omega_{\mathrm B} ω B を持つ運動である。量子力学では、周期運動の準位は間隔 ℏ ω = ℏ ⋅ 2 π / T \hbar\omega=\hbar\cdot2\pi/T ℏ ω = ℏ ⋅ 2 π / T で並ぶ(Bohr–Sommerfeld の量子化で d E / d J = ω \mathrm{d} E/\mathrm{d} J=\omega d E / d J = ω を J J J の間隔 ℏ \hbar ℏ に掛けたもの)。Bloch 振動では ℏ ω B = e E a \hbar\omega_{\mathrm B}=eEa ℏ ω B = e E a で、準位は
E ν = ν e E a + 一定 , ν = 0 , ± 1 , ± 2 , … E_\nu=\nu\,eEa+\text{一定},\qquad \nu=0,\pm1,\pm2,\dots E ν = ν e E a + 一定 , ν = 0 , ± 1 , ± 2 , …
と等間隔に並ぶ。これを Wannier–Stark の階段と呼ぶ。隣の準位の間隔が、電子が格子1つ分 a a a ずれたときのポテンシャルエネルギーの差 e E a eEa e E a に等しいのは、電場の中で並進 x → x + a x\to x+a x → x + a を行うとエネルギーが e E a eEa e E a ずれることの表れである。
なぜ普通の金属では見えないか
普通の結晶では a a a が 0.3 n m 0.3\,\mathrm{nm} 0.3 nm 程度なので、T B T_{\mathrm B} T B はかなり長い。その間に電子は不純物やフォノンに散乱され、k k k が1周する前に元の分布へ戻されてしまう。Bloch 振動を見るには、a a a を人工的に大きくした半導体の超格子を使う。Feldmann らは GaAs/AlGaAs の超格子で、電場に比例する周期の振動を観測した1 1 J. Feldmann, K. Leo, J. Shah, D. A. B. Miller, J. E. Cunningham, T. Meier, G. von Plessen, A. Schulze, P. Thomas, S. Schmitt-Rink (1992), Phys. Rev. B 46 , 7252–7255. 。
例題 4.10.4 超格子と普通の結晶の Bloch 周期 ★★
オリジナル 電場の大きさを E = 1.0 × 10 5 V m − 1 E=1.0\times10^5\,\mathrm{V\,m^{-1}} E = 1.0 × 1 0 5 V m − 1 とする。
(1) 周期 a = 10 n m a=10\,\mathrm{nm} a = 10 nm 、ミニバンドの幅 4 t h = 20 m e V 4t_{\mathrm h}=20\,\mathrm{meV} 4 t h = 20 meV の超格子で、Bloch 振動の周期 T B T_{\mathrm B} T B 、振動数 1 / T B 1/T_{\mathrm B} 1/ T B 、位置の振れ幅を求めよ。
(2) 同じ電場で a = 0.30 n m a=0.30\,\mathrm{nm} a = 0.30 nm の普通の結晶では T B T_{\mathrm B} T B はいくらか。散乱の時間を τ = 1 × 10 − 14 s \tau=1\times10^{-14}\,\mathrm s τ = 1 × 1 0 − 14 s として T B / τ T_{\mathrm B}/\tau T B / τ を求め、振動が見えるかを判断せよ。
解答 (1) 式 (4.10.15) から
T B = h e E a = 6.626 × 10 − 34 1.602 × 10 − 19 × 1.0 × 10 5 × 10 × 10 − 9 = 4.14 × 10 − 12 s , 1 T B = 242 G H z . T_{\mathrm B}=\frac{h}{eEa}=\frac{6.626\times10^{-34}}{1.602\times10^{-19}\times1.0\times10^5\times10\times10^{-9}}=4.14\times10^{-12}\,\mathrm s,\qquad \frac1{T_{\mathrm B}}=242\,\mathrm{GHz}. T B = e E a h = 1.602 × 1 0 − 19 × 1.0 × 1 0 5 × 10 × 1 0 − 9 6.626 × 1 0 − 34 = 4.14 × 1 0 − 12 s , T B 1 = 242 GHz .
振れ幅は式 (4.10.16) から帯の幅を e E eE e E で割ったもので、20 × 10 − 3 e V / ( e × 1.0 × 10 5 V m − 1 ) = 2.0 × 10 − 7 m = 200 n m 20\times10^{-3}\,\mathrm{eV}/(e\times1.0\times10^5\,\mathrm{V\,m^{-1}})=2.0\times10^{-7}\,\mathrm m=200\,\mathrm{nm} 20 × 1 0 − 3 eV / ( e × 1.0 × 1 0 5 V m − 1 ) = 2.0 × 1 0 − 7 m = 200 nm 。超格子の周期 20 20 20 個分を往復する。
(2) T B = h / ( e E a ) = 1.38 × 10 − 10 s T_{\mathrm B}=h/(eEa)=1.38\times10^{-10}\,\mathrm s T B = h / ( e E a ) = 1.38 × 1 0 − 10 s で、T B / τ = 1.38 × 10 4 T_{\mathrm B}/\tau=1.38\times10^4 T B / τ = 1.38 × 1 0 4 。k k k が BZ を1周するまでに 10 4 10^4 1 0 4 回散乱されるので、振動は見えない。T B T_{\mathrm B} T B は a a a に反比例するので、a a a を30倍以上にした超格子では周期が同じ割合で短くなり、散乱される前に k k k が BZ を1周する条件を満たしやすくなる。
1 J. Feldmann, K. Leo, J. Shah, D. A. B. Miller, J. E. Cunningham, T. Meier, G. von Plessen, A. Schulze, P. Thomas, S. Schmitt-Rink (1992), Phys. Rev. B 46 , 7252–7255.
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com