ハミルトニアンの選び方で固体を考える 第4部 第9章 §4.9.1 ★★★ 有効質量のバンドに静電ポテンシャルを足した1電子の H ^ \hat H H ^ を選び、その静電ポテンシャルを電荷から Poisson 方程式で決める。平衡では化学ポテンシャルが一定なので、電子と正孔の数は局所の電位の Boltzmann 因子になる
第4部第5章 では、p 型と n 型の半導体をつなぐと境目に空乏層ができ、バンドが曲がることを見た。そこでバンドを曲げたのは、2つの半導体の化学ポテンシャルの違いである。この章では、バンドを外から曲げる。半導体の表面に薄い絶縁膜(酸化膜)を付け、その上に金属の電極(ゲート)を置く。この金属(Metal)・酸化膜(Oxide)・半導体(Semiconductor)の重なりを MOS 構造と呼ぶ。ゲートに電圧をかけると、酸化膜はコンデンサーの誘電体として働き、半導体の表面に電荷が誘起される。その電荷は Poisson 方程式に従って半導体の電位を曲げ、曲がったバンドの中で電子と正孔の数は Boltzmann 因子に従って変わる。この電荷と電位の釣り合いを解くことが、この章の計算の中心である。
選ぶハミルトニアンは新しくない。伝導帯と価電子帯の端の近くを有効質量で書いた1電子のバンド(第2部第7章 )に、ゆっくり変わる静電ポテンシャル − e ϕ ( r ) -e\phi(\bm{r}) − e ϕ ( r ) を足すだけである。新しく加わるのは、その ϕ \phi ϕ を電子と正孔と不純物のイオンの電荷そのものから Poisson 方程式で決める、という自己無撞着の条件である(4.9.1節)。
結論を先に書いておく。p 型の半導体の表面では、ゲートの電圧を上げていくと、まず正孔が押し出されて負のアクセプターのイオンだけが残る空乏層ができる。さらに上げて表面の電位が 2 ϕ B 2\phi_{\mathrm B} 2 ϕ B (ϕ B \phi_{\mathrm B} ϕ B は半導体の中の化学ポテンシャルが禁制帯の中央からどれだけ下にあるかを電圧で表したもの)を越えると、表面には正孔の代わりに電子が集まる。これを反転層と呼ぶ。反転層ができるゲートの電圧がしきい値電圧 V t h V_{\mathrm{th}} V th で、それを越えた分の電圧はほとんどそのまま反転層の電子の電荷 C o x ( V G − V t h ) C_{\mathrm{ox}}(V_{\mathrm G}-V_{\mathrm{th}}) C ox ( V G − V th ) になる(4.9.2 節)。反転層の両端に n 型の電極(ソースとドレイン)を付けると、ゲートの電圧で電子の数を、したがって電流を変えることができる。これが MOS 型の電界効果トランジスタ(MOSFET)で、ドレインの電圧を上げると電流はやがて増えなくなる(飽和)。しきい値より下では、電流は熱で励起された電子の拡散で流れ、ゲートの電圧に対して指数関数的に変わる。その傾きには、Fermi 分布の裾で決まる下限がある(4.9.3 節)。反転層の厚さは数 nm で、電子の de Broglie 波長と同じ程度である。そこで表面に垂直な運動は量子化され、電子は面内だけを自由に動く2次元電子系になる。2次元電子系に強い磁場をかけると Landau 準位ができ、ゲートの電圧で準位の詰まり方を変えることができる。量子 Hall 効果が最初に見つかったのは、この Si の MOSFET の反転層である(4.9.4 節)。そのあと、削ったものを並べて近似の裏付けの札を貼り(4.9.5 節)、容量と電流の測定で確かめ、授業の式と並べる(4.9.6 節)。トランジスタを組み合わせた回路は扱わない。
この章では、電子の電荷を − e -e − e (e > 0 e>0 e > 0 )と書き、静電ポテンシャルを ϕ \phi ϕ 、化学ポテンシャルを μ \mu μ 、電子の移動度を μ e \mu_{\mathrm e} μ e と書く。Si と酸化膜(SiO2 _2 2 )の比誘電率を ε r s \varepsilon_{\mathrm{rs}} ε rs ・ε r o \varepsilon_{\mathrm{ro}} ε ro と書き、数値例ではそれぞれ 11.7 11.7 11.7 ・3.9 3.9 3.9 とする。電荷と容量は、断らない限り単位面積あたりの量である。しきい値電圧は温度 T T T と紛れないよう V t h V_{\mathrm{th}} V th と書く。
選ぶ H ^ \hat H H ^
MOS 構造の半導体の部分を、表面を x = 0 x=0 x = 0 、半導体の中へ向かう向きを x > 0 x>0 x > 0 とする座標で考える。ゲートは表面全体を一様に覆っているとし、すべての量は x x x だけの関数とする。伝導帯の電子に対して選ぶハミルトニアンは
H ^ = ε c + ∑ i = x , y , z p ^ i 2 2 m i − e ϕ ( x ) (4.9.1) \hat H=\varepsilon _{\mathrm c}+\sum_{i=x,y,z}\frac{\hat p_i^2}{2m_i}-e\phi(x)\tag{4.9.1} H ^ = ε c + i = x , y , z ∑ 2 m i p ^ i 2 − e ϕ ( x ) ( 4.9.1 )
である。ε c \varepsilon _{\mathrm c} ε c は半導体の奥(バルク)の伝導帯の底、m i m_i m i は伝導帯の谷の有効質量で、Si では回転楕円体の谷の長軸方向に m l m_{\mathrm l} m l 、短軸方向に m t m_{\mathrm t} m t になる。正孔に対しても、価電子帯の端 ε v \varepsilon _{\mathrm v} ε v と正孔の有効質量で同じ形のものを選ぶ。式 (4.9.1) は、結晶の周期ポテンシャルの中の電子を、バンドの端の近くで有効質量を持つ自由な粒子に置き換え、そこへ格子間隔よりずっとゆっくり変わる外のポテンシャルを足したものである(有効質量近似。第2部第7章 ・第4部第5章 )。
静電ポテンシャル ϕ ( x ) \phi(x) ϕ ( x ) は外から与えるものでなく、電荷の分布から決まる。半導体の中の電荷密度は、正孔の数密度 p ( x ) p(x) p ( x ) ・電子の数密度 n ( x ) n(x) n ( x ) ・アクセプターの数密度 N A N_{\mathrm A} N A (すべてイオン化しているとする)でe ( p ( x ) − n ( x ) − N A ) e\bigl(p(x)-n(x)-N_{\mathrm A}\bigr) e ( p ( x ) − n ( x ) − N A ) と書けるので、Poisson 方程式は
d 2 ϕ d x 2 = − e ε 0 ε r s ( p ( x ) − n ( x ) − N A ) (4.9.2) \frac{\mathrm{d}^2\phi}{\mathrm{d} x^2}=-\frac{e}{\varepsilon_0\varepsilon_{\mathrm{rs}}}\bigl(p(x)-n(x)-N_{\mathrm A}\bigr)\tag{4.9.2} d x 2 d 2 ϕ = − ε 0 ε rs e ( p ( x ) − n ( x ) − N A ) ( 4.9.2 )
になる。ここで p ( x ) p(x) p ( x ) と n ( x ) n(x) n ( x ) は、式 (4.9.1) の固有状態に電子を詰めた結果として決まる。つまり、電子が感じるポテンシャルを、その電子自身を含めた電荷の平均から決めている。これは電子どうしの Coulomb 相互作用を平均の場で置き換える Hartree の近似で、この章の計算はすべてこの自己無撞着の枠の中で行う。
当てる道具:平衡では化学ポテンシャルが一定
ゲートに電流が流れない(酸化膜は絶縁体)ので、ゲートの電圧をかけたあとの半導体は熱平衡にある。熱平衡では化学ポテンシャル μ \mu μ が場所によらない。一方、バンドの端は ε c − e ϕ ( x ) \varepsilon _{\mathrm c}-e\phi(x) ε c − e ϕ ( x ) のように場所によって曲がる。ϕ \phi ϕ が格子間隔や電子の波長に比べてゆっくり変わるなら、各点の近くでは平らなバンドがその点の高さにあると考えてよく、第2部第7章 の非縮退の式の ε c \varepsilon _{\mathrm c} ε c を ε c − e ϕ ( x ) \varepsilon _{\mathrm c}-e\phi(x) ε c − e ϕ ( x ) に置き換えたものがその点の電子の数になる。バルクで ϕ = 0 \phi=0 ϕ = 0 とすると
n ( x ) = n 0 e e ϕ ( x ) / ( k B T ) , p ( x ) = p 0 e − e ϕ ( x ) / ( k B T ) , n 0 p 0 = n i 2 (4.9.3) n(x)=n_0\,\mathrm{e}^{e\phi(x)/( k_{\mathrm B} T)},\qquad p(x)=p_0\,\mathrm{e}^{-e\phi(x)/( k_{\mathrm B} T)},\qquad n_0p_0=n_{\mathrm i}^2\tag{4.9.3} n ( x ) = n 0 e e ϕ ( x ) / ( k B T ) , p ( x ) = p 0 e − e ϕ ( x ) / ( k B T ) , n 0 p 0 = n i 2 ( 4.9.3 )
である。n 0 n_0 n 0 と p 0 p_0 p 0 はバルクの電子と正孔の数密度で、n i n_{\mathrm i} n i は真性キャリア密度である。p 型では p 0 ≃ N A p_0\simeq N_{\mathrm A} p 0 ≃ N A 、n 0 = n i 2 / N A n_0=n_{\mathrm i}^2/N_{\mathrm A} n 0 = n i 2 / N A になる。ϕ > 0 \phi>0 ϕ > 0 の所ではバンドが下へ曲がり、電子が増えて正孔が減る。式 (4.9.3) を式 (4.9.2) に入れたものを Poisson–Boltzmann 方程式と呼ぶ。
バルクの化学ポテンシャルの位置を表す量として、
e ϕ B ≔ ε i − μ = k B T ln N A n i (4.9.4) e\phi_{\mathrm B}\coloneqq\varepsilon _{\mathrm i}-\mu= k_{\mathrm B} T\ln\frac{N_{\mathrm A}}{n_{\mathrm i}}\tag{4.9.4} e ϕ B : = ε i − μ = k B T ln n i N A ( 4.9.4 )
を定める。ε i \varepsilon _{\mathrm i} ε i は電子と正孔の数が等しくなるときの化学ポテンシャル(禁制帯のほぼ中央)で、p 型では μ \mu μ が ε i \varepsilon _{\mathrm i} ε i より e ϕ B e\phi_{\mathrm B} e ϕ B だけ下にある。2つ目の等号は、p 0 = n i e ( ε i − μ ) / ( k B T ) p_0=n_{\mathrm i}\mathrm{e}^{(\varepsilon _{\mathrm i}-\mu)/( k_{\mathrm B} T)} p 0 = n i e ( ε i − μ ) / ( k B T ) (第2部第7章 )と p 0 = N A p_0=N_{\mathrm A} p 0 = N A から出る。
式 (4.9.2) の右辺に現れる長さの目安が、Debye の遮蔽長(第4部第5章 で pn 接合の n 側について定義したものと同じ量で、ここではアクセプターの密度で書く)
λ D = ε 0 ε r s k B T e 2 N A (4.9.5) \lambda_{\mathrm D}=\sqrt{\frac{\varepsilon_0\varepsilon_{\mathrm{rs}} k_{\mathrm B} T}{e^2N_{\mathrm A}}}\tag{4.9.5} λ D = e 2 N A ε 0 ε rs k B T ( 4.9.5 )
である。表面の電位が k B T / e k_{\mathrm B} T/e k B T / e より小さいとき、式 (4.9.3) を1次まで展開すると式 (4.9.2) は ϕ ′ ′ = ϕ / λ D 2 \phi''=\phi/\lambda_{\mathrm D}^2 ϕ ′′ = ϕ / λ D 2 になり、表面の電位は e − x / λ D \mathrm{e}^{-x/\lambda_{\mathrm D}} e − x / λ D で減る。正孔が少し動いて、表面の電荷を λ D \lambda_{\mathrm D} λ D ほどの厚さで打ち消すからである。
ゲートの電圧と表面の電位
ゲートの金属と半導体の化学ポテンシャルの差を、ゲートの電圧 V G V_{\mathrm G} V G と呼ぶ(金属の化学ポテンシャルが半導体より e V G eV_{\mathrm G} e V G だけ下にあるとき V G > 0 V_{\mathrm G}>0 V G > 0 )。V G V_{\mathrm G} V G は、酸化膜にかかる電圧 V o x V_{\mathrm{ox}} V ox と、半導体の表面とバルクの電位差 ϕ s ≔ ϕ ( 0 ) \phi_{\mathrm s}\coloneqq\phi(0) ϕ s : = ϕ ( 0 ) に分かれる。ただし、金属と半導体の仕事関数が違ったり酸化膜に固定された電荷があったりすると、V G = 0 V_{\mathrm G}=0 V G = 0 でもバンドは曲がっている。バンドが平らになる電圧を平坦帯電圧 V F B V_{\mathrm{FB}} V FB と呼び、その分をずらして
V G − V F B = ϕ s + V o x , V o x = − Q s C o x , C o x = ε 0 ε r o d o x (4.9.6) V_{\mathrm G}-V_{\mathrm{FB}}=\phi_{\mathrm s}+V_{\mathrm{ox}},\qquad V_{\mathrm{ox}}=-\frac{Q_{\mathrm s}}{C_{\mathrm{ox}}},\qquad C_{\mathrm{ox}}=\frac{\varepsilon_0\varepsilon_{\mathrm{ro}}}{d_{\mathrm{ox}}}\tag{4.9.6} V G − V FB = ϕ s + V ox , V ox = − C ox Q s , C ox = d ox ε 0 ε ro ( 4.9.6 )
と書く。Q s Q_{\mathrm s} Q s は半導体の中の全電荷(単位面積あたり)、d o x d_{\mathrm{ox}} d ox は酸化膜の厚さである。2つ目の式は、ゲートの金属に − Q s -Q_{\mathrm s} − Q s の電荷が載り、酸化膜を平行板コンデンサーとみなしたものである。Q s Q_{\mathrm s} Q s は ϕ s \phi_{\mathrm s} ϕ s の関数として Poisson–Boltzmann 方程式から決まるので、式 (4.9.6) を ϕ s \phi_{\mathrm s} ϕ s について解けば、ゲートの電圧から表面の電位が決まる。
例題 4.9.1 p 型 Si の基本の長さと電圧 ★★
オリジナル 室温(T = 300 K T=300\,\mathrm K T = 300 K 、k B T / e = 0.02585 V k_{\mathrm B} T/e=0.02585\,\mathrm V k B T / e = 0.02585 V )の p 型 Si で、N A = 1.0 × 10 23 m − 3 N_{\mathrm A}=1.0\times10^{23}\,\mathrm{m^{-3}} N A = 1.0 × 1 0 23 m − 3 、n i = 1.0 × 10 16 m − 3 n_{\mathrm i}=1.0\times10^{16}\,\mathrm{m^{-3}} n i = 1.0 × 1 0 16 m − 3 、ε r s = 11.7 \varepsilon_{\mathrm{rs}}=11.7 ε rs = 11.7 とする。
(1) バルクの電子の数密度 n 0 n_0 n 0 を求めよ。
(2) ϕ B \phi_{\mathrm B} ϕ B を求めよ。
(3) Debye の遮蔽長 λ D \lambda_{\mathrm D} λ D を求めよ。
(4) アクセプターの平均の間隔 N A − 1 / 3 N_{\mathrm A}^{-1/3} N A − 1/3 を λ D \lambda_{\mathrm D} λ D と比べ、式 (4.9.2) でアクセプターの電荷を一様な密度 e N A eN_{\mathrm A} e N A に均したことがどの長さで正しいかを述べよ。
解答 (1) n 0 = n i 2 / N A = ( 1.0 × 10 16 ) 2 / ( 1.0 × 10 23 ) = 1.0 × 10 9 m − 3 n_0=n_{\mathrm i}^2/N_{\mathrm A}=(1.0\times10^{16})^2/(1.0\times10^{23})=1.0\times10^{9}\,\mathrm{m^{-3}} n 0 = n i 2 / N A = ( 1.0 × 1 0 16 ) 2 / ( 1.0 × 1 0 23 ) = 1.0 × 1 0 9 m − 3 である。正孔の 10 − 14 10^{-14} 1 0 − 14 倍しかない。
(2) ln ( N A / n i ) = ln 10 7 = 16.12 \ln(N_{\mathrm A}/n_{\mathrm i})=\ln10^7=16.12 ln ( N A / n i ) = ln 1 0 7 = 16.12 なので ϕ B = 0.02585 × 16.12 = 0.417 V \phi_{\mathrm B}=0.02585\times16.12=0.417\,\mathrm V ϕ B = 0.02585 × 16.12 = 0.417 V である。
(3) ε 0 ε r s = 11.7 × 8.854 × 10 − 12 = 1.036 × 10 − 10 F m − 1 \varepsilon_0\varepsilon_{\mathrm{rs}}=11.7\times8.854\times10^{-12}=1.036\times10^{-10}\,\mathrm{F\,m^{-1}} ε 0 ε rs = 11.7 × 8.854 × 1 0 − 12 = 1.036 × 1 0 − 10 F m − 1 、e N A = 1.602 × 10 4 C m − 3 eN_{\mathrm A}=1.602\times10^{4}\,\mathrm{C\,m^{-3}} e N A = 1.602 × 1 0 4 C m − 3 から
λ D = 1.036 × 10 − 10 × 0.02585 1.602 × 10 4 = 1.672 × 10 − 16 m = 12.9 n m \lambda_{\mathrm D}=\sqrt{\frac{1.036\times10^{-10}\times0.02585}{1.602\times10^{4}}}=\sqrt{1.672\times10^{-16}}\,\mathrm m=12.9\,\mathrm{nm} λ D = 1.602 × 1 0 4 1.036 × 1 0 − 10 × 0.02585 = 1.672 × 1 0 − 16 m = 12.9 nm
である。
(4) N A − 1 / 3 = ( 1.0 × 10 23 ) − 1 / 3 m = 21.5 n m N_{\mathrm A}^{-1/3}=(1.0\times10^{23})^{-1/3}\,\mathrm m=21.5\,\mathrm{nm} N A − 1/3 = ( 1.0 × 1 0 23 ) − 1/3 m = 21.5 nm で、λ D \lambda_{\mathrm D} λ D より長い。λ D \lambda_{\mathrm D} λ D の大きさの立方体に入るアクセプターは平均 N A λ D 3 = 0.22 N_{\mathrm A}\lambda_{\mathrm D}^3=0.22 N A λ D 3 = 0.22 個しかない。したがって、λ D \lambda_{\mathrm D} λ D の長さで電位が変わる様子を一様な電荷で計算した結果は、個々のイオンの位置で揺らぐ実際の電位を平均したものにすぎない。一様に均してよいのは、4.9.2 節で出てくる空乏層の厚さ(0.1 μ m 0.1\,\mathrm{\mu m} 0.1 μ m ほど)のように、アクセプターの間隔よりずっと長い長さで電荷を数える量である。この点は4.9.5 節でもう一度取り上げる。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com