• ベストアンサー

フェルミー準位

フェルミー準位ってなんですか?式において、 n(x)=n(0)exp{-(E(x)-Ef)/kT} のEfってことはわかるのですが、それが何を意味するかわかりません。 言葉で教えてください。

質問者が選んだベストアンサー

  • ベストアンサー
  • Umada
  • ベストアンサー率83% (1169/1405)
回答No.2

まず ・固体中の電子の有するエネルギーはどんな値でも自由にとれるわけでなく、いくつかのエネルギー準位に限られる。 ・電子は上記の式に従ってエネルギー準位の低い方から入っていく。ただし下の準位から全部詰めて入っていくわけでなく、上記の式に従って分布する。(従って低い準位でも若干の空きがあり、一方で高い準位にも少しは入っている) ということを思い出して下さい。これが分からないとFermi準位については分かりませんから、統計力学か電子論、固体物性の教科書をもう一度復習して下さい。 Fermi準位の理解には「大学の大部屋の講義」を思い浮かべて頂くと分かりやすいと思います。 まず座席は、前の列から後ろの列までたくさんあるとします。とりあえず30列くらいあるとしましょうか。 そこへ学生が入ってきます。学生はあまり前の方に座りたがりませんから、席は大体後ろから詰まっていきます。もちろん後ろの1列が完全に埋まってから次の1列が埋まり始めるわけでなく、後ろから2列目や3列目も早い段階から埋まり始めることでしょう。人数は少ないながらも、前の方にも何人かは熱心な学生が座っていることでしょう。図に描くと下のようになります。 座席の埋まっている数 ↑ │          │         ■ │         ■   │       ■■■ │  ■ ■ ■■■■ └─────────→教授との距離  前        * 後 もう見当が付いたことと思いますが、 座席の場所(教授からの距離):エネルギー準位 横一列の中で、埋まっている座席の割合:占有確率 に対応するわけです。「前の席ほどエネルギー準位が高い」ということですね。 Fermi準位とは、「だいたいどの辺の準位まで埋まっているか」「(電子は)平均としてどの辺の準位にいるか」を表す指標です。(厳密に言うならFermi-Dirac統計で「占有確率が1/2になるエネルギー準位」) 上の図では部屋の総収容数に比べて学生が少ないですから、後ろの方に学生が少しだけ座っておしまいです。ですから「Fermi準位は低い」ということになります。だいたい*を付けた付近がFermi準位でしょう。 さて期末試験が近付いて、学生の数が増えたとしましょう。後ろは埋まってしまって、あふれた学生は必然的に前の方にも座ることになります。従って分布は下の図のようになるでしょう。 座席の埋まっている数 ↑       ■■■ │      ■■■■ │     ■■■■■ │    ■■■■■■   │   ■■■■■■■ │ ■■■■■■■■■ └─────────→教授との距離  前    *    後 この場合は学生はやむなく前に座るわけで、Fermi準位も上がります。大体*の辺りでしょうか。「学生の平均の位置」あるいは「座席はどの辺まで埋まっているか」は、前寄りになるわけですね。 最初の n(x)=n(0)exp{-(E(x)-Ef)/kT} に立ち返れば、これは n(0)exp(Ef/kT)exp{-E(x)/kT} と書き直せますから、Fermiエネルギーは一種の規格化因子とも言えます。 エネルギー準位の数(正確には状態密度)×占有確率を前エネルギーに亘って積分すれば、それは系の全電子数に等しくなければなりません。そこで規格化因子exp(Ef/kT)が入るわけです。 私自身、多少いいかげんな理解をしている部分がありますので、上記の回答の信頼性は適当に割引いてお読みください。

tera242
質問者

お礼

凄くわかりやすい説明で感動しました。 これで、テストでフェルミー準位について説明するとききっちりと書けます。 本当にどうもありがとうございました。

その他の回答 (4)

  • siegmund
  • ベストアンサー率64% (701/1090)
回答No.5

siegmund です. Umada さんの座席のご説明はとてもおもしろかったです. 今度,フェルミ分布の話をするときに使わせていただこうかと思っています. > 半導体ですとバンドギャップ部分だけ座席が取り払われている教室を > イメージすればよろしいのでしょうか・ なるほど,good idea ですね. 大講義室ですと,真ん中当たりに黒板と平行に通路がある部屋もありますね. そういう部屋の話にしますか. 通路の前と後ろとでは座席占有確率が急に変化しそうで, そこらへんも価電子帯と伝導帯の性質をよく反映しそうです. > あとうまい具合にFermi準位がそこに来ているような図にしないといけないですね そこがちょっと難しそうです.

tera242
質問者

お礼

ぜひとも、座席の説明をしてみてください。 日常的なことを授業に入れれば、学生の反応もいいはずです。 また、何かお世話になることがあるかもしれませんが、 そのときはよろしくお願いします。

  • Umada
  • ベストアンサー率83% (1169/1405)
回答No.4

あちゃー、やっぱりボロがありましたか。siegmundさん、いつもご指摘ありがとうございます。 座席の数と状態密度の関係について「座席定員×占有確率=座っている人数」になるのは全くご指摘の通りです。座席の数はどの列でも一定として図を描きましたが、実際には列によって座席の数が異なりますよね。 半導体-金属の違いについてのご指摘もありがとうございました、勉強不足が露呈してしまいました。(半導体ですとバンドギャップ部分だけ座席が取り払われている教室をイメージすればよろしいのでしょうか・・・あとうまい具合にFermi準位がそこに来ているような図にしないといけないですね) 学生が前に座らないのは不変の真理(?)ですね。私も学生の時分は「前の方はポテンシャルが高いから、私の波動関数はその辺にはあまりないのだ」などとうそぶいて、いつも後ろに座っておりました(笑)。

tera242
質問者

お礼

僕も授業はどちらかというと、後ろの方に座るタイプです。黒板の文字が小さい時は前列にのり出しますが・・・。 学生の中でも前に座る人と後ろに座る人でかなりわかれてるような気がします。活発な電子とそうでない電子がいるんですね。僕はどちらかというとそうでないほうかもしれないですけど・・・。

  • siegmund
  • ベストアンサー率64% (701/1090)
回答No.3

(1)  n(x)=n(0)exp{-(E(x)-Ef)/kT} の式からすると,半導体の話のようですね. Umada さんの解説はなかなかわかりやすく,面白い解説でした. 私の講義も後ろの方から学生が埋まります(^^;). さて,Umada さんのご説明は,本質的に金属に対する説明です. 教授との距離がエネルギーに対応し (後ろの方がエネルギーが低い,逆であって欲しいんですがね), 座席が横にいくつあるか(Umada さんの図では6つ)が同じエネルギーにある準位の数 (状態密度)に対応しています. 本当は状態密度はエネルギーの関数なのですが, Umada さんの図ではわかりやすく6に固定しています. Umada さんの説明に適合する分布関数はフェルミ-ディラックの分布関数 (2)  f(x) = 1/{exp[(Ef-E(x))/kT] + 1} です. (a) (2)で E(x) = Ef とおくと,f(x) = 1/2 になること (b) (2)で T→0 とすると,E(x)<Ef なら f(x) = 1, E(x)>Ef なら f(x) = 0.   すなわち,フェルミ準位以下は完全に占有されていて,フェルミ準位以上は空. であることを確認してください. なお,(2)の分母の +1 を -1にすると, ボーズ・アインシュタイン分布関数になります. (2)の分母の +1 を無視すると(1)のタイプになりますが, 無視できるためには (3)  exp[(Ef-E(x))/kT] >> 1 が必要で,通常の金属ではこの条件は満たされません (縮退している,といっています). ところが,半導体で Ef が禁制帯内にあり,しかも上下のバンド端から十分離れていると (3)の条件が満たされます. このとき(2)を (4)  f(x) = exp[-(E(x)-Ef)/kT] の形に書くことができて,これが(1)の正体です. 繰り返しになりますが,(4)が使えるのは(3)の条件の下(つまり(4)の f <<1 )の ときのみです. そこら辺を考えずに(4)を使うと,フェルミ粒子なのに占有確率が1を越える というおかしなことになります.

tera242
質問者

お礼

Umadaさんの話がある特殊な場合に成り立つことがよくわかりました。 補足の説明大変ありがとうございます。 回答の文を見るに当たり、大学で授業をなさってるようですね。 あなたの授業のフェルミー準位が上がることを祈っています。(笑)

  • may-may-jp
  • ベストアンサー率26% (324/1203)
回答No.1

「半導体」のところで習ったような・・・。

関連するQ&A

  • フェルミ準位

    エネルギーEとE+dEの間にある電子の数n(E)dEは n(E)dE=Z(E)F(E)dE Z(E):単位体積、単位エネルギーあたりの状態密度 F(E):フェルミ・ディラックの分布関数 F(E)=1/(1+exp[(E-E_f)/kT]) T:絶対温度 E_f:フェルミ準位 電子の状態密度は Z(E)dE=(4π/h^3)*(2m)^(3/2)e^(1/2)dE m:固体中での電子の有効質量 h:プランク定数 T=0Kでは n=(4π/h^3)*(2m)^(3/2)∫[0~E_f0]e^(1/2)dE E_f0:T=0KのときのE_f 変形するとE_f0=(h^2/2m)(3n/8π)^(2/3) T>0Kのときは n(E)=∫[0~∞]Z(E)dE/(1+exp[(E-E_f)/kT]) ここでE_f>>kTとすると E_f≒E_f0[1-(π^2/12)(kT/E_f0)^2] となるのですが、この最終式を導くにはどうしたらいいのでしょうか。とくに「12」という分母がどこから出てくるのかが気になります。

  • フェルミ準位について

    フェルミ・ディラックの分布関数は f(E)=1/{exp(E-EF/kT)+1}である。 1)Ef/k=5×10^4、T=5×10^2K とするならば -∂f/∂E と E/k の関係をグラフに書け。 2)E=Ef+δ とすると   f(δ)=1-f(-δ) である。これを証明せよ。 1)は f(E) を微分して ∂f/∂E=-1/4kT まではできたんですけど、そこからどうやったらいいのかわかりません。 固体物理学の本を見ても証明とかは省いてあってできないんです。 どちらかひとつだけでもいいので教えてください。お願いします。

  • 2準位系の粒子数

    2準位系の粒子数 2準位形のカノニカル分布を考えています。 全粒子数はN、状態1がE=0、状態2がE=εとします。 このとき分配関数Zは(1)、(2)式のようにかけます。 ここで、なぜ以下のように書けるのか分かりません。 「よって、A,Bの状態をとる粒子数をそれぞれn_A,n_Bとすると、 (3)、(4)式のようにかける」

  • ゾンマーフェルト展開

    エネルギーEとE+dEの間にある電子の数n(E)dEは n(E)dE=Z(E)F(E)dE Z(E):単位体積、単位エネルギーあたりの状態密度 F(E):フェルミ・ディラックの分布関数 F(E)=1/(1+exp[(E-E_f)/kT]) T:絶対温度 E_f:フェルミ準位 電子の状態密度は Z(E)dE=(4π/h^3)*(2m)^(3/2)e^(1/2)dE m:固体中での電子の有効質量 h:プランク定数 T=0Kでは n=(4π/h^3)*(2m)^(3/2)∫[0~E_f0]e^(1/2)dE E_f0:T=0KのときのE_f 変形するとE_f0=(h^2/2m)(3n/8π)^(2/3) T>0Kのときは n(E)=∫[0~∞]Z(E)dE/(1+exp[(E-E_f)/kT]) ここでE_f>>kTとすると E_f≒E_f0[1-(π^2/12)(kT/E_f0)^2] この式を導こうとしていたところです。 先日、回答者の方からのお力をいただきまして、 以下のように計算してみました。 フェルミ分布関数fはT=0でステップ関数なので、df/dE はδ関数。ところが有限温度だとステップがぼやけるため、df/dE はガウス関数で近似できる。 n(E)=∫[0~∞]Z(E)dE/(1+exp[(E-E_f)/kT]) dF(E)/dE≒-(1/sqrt(2π))exp(-((E-E_f)^2/2σ^2)) F(E)=1/(1+exp[(E-E_f)/kT])の導関数にE=E_fを代入した式=-1/4kT≒(-1/sqrt(2π)) σ=2sqrt(2)kT/sqrt(π) 部分積分を行う。 n(E)=4π/h^3(2m)^(3/2){[F(E)・(2/3)E^(3/2)]_0^∞-∫(0→∞)(dF/dE)(2/3)E^(3/2)dE} =(2/(3sqrt(2π)σ))∫(0→∞)exp(-(E-E_f)^2/2σ^2)E^(3/2)dE ここでいきづまっています。3/2乗に2乗の指数関数が出てきていて、どう積分したものやらと思っております。ゾンマーフェルト展開についてのっている本だけでも紹介していただけないでしょうか。少しでも助言をお願いします。

  • 半導体工学の問題

    半導体電子中の電子電流密度Jnとすると Jn = e * n *μn * E + e * Dn * (∂n/∂x) (1) で与えられ、ここでμnは電子の移動度、Dnは電子の拡散係数である。 電子密度nを n = ni * exp{e*(Vf-V) / k*T } (2) と書き、Vfはフェルミ準位に対応する電位で、場所は一定である。磁界の強さEは E = - ∂V / ∂x (3) で与えられる。 (2)式をxで微分して、(3)でを用い、平衡状態でJn=0と置いて、アインシュタインの関係式を導きなさい。 どう頑張ってもアインシュタインの関係式(D=μkT)にならないのでお力をお貸しください。

  • 分配関数と状態密度

    参考書によって分配関数Z=Σw(E)exp(-βE) (1) と書いてあったり、Z=Σexp(-βE) (2) とあったりします。またZ=1/{h^(3N)N!}∫exp(-βH)d^NΓ (3) と表している時もあります。(ここではw(E)は微視的状態密度、Hはハミルトニアン、β=1/(kT)とする。) ・微視的状態密度と物性物理で使う状態密度とは違いがありますか? ・(1)式と(2)式では何が違うのですか? ・当然、問題によって異なるでしょうが(1)、(2)、(3)式はどういう時にどれを使えば良いのか、問題を解くときのコツをお願いします。

  • 積分の問題

    d/dx∫(x^2→2x+1)e^(-kt^2)dt この問題は、先に∫(x^2→2x+1)e^(-kt^2)dtを計算して、その答えをxで微分するという意味でしょうか?

  • <mv^2/2>を求める

    こんばんは。熱力学での質問です。 ∫exp(-ax^2)dx=(π/a)^1/2・・・・・・・・(ⅰ) ∫x^2exp(-ax^2)dx=a^-3/2(π/4)^1/2・・・(ⅱ) ρ(0)=π^-3/2(m/2kT)^3/2・・・・・・・・(ⅲ) が与えられた状態で、 <mv^2/2>=∫∫∫(mv^2/2)ρ(0)exp(-mv^2/2kT)d^3v・・・・(ⅳ) を求めよ、という問題で、(ⅳ)を変形して (m/2)ρ(0)∫∫∫v^2 exp(-mv^2/2kT)d^3v とし、 (ⅱ)でv=x , a=m/2kT として、(ⅱ)、(ⅲ)を用いて(ⅳ)を計算したところ、答えがどうしても (kT)^3/2m^2 となってしまい、公式の<mv^2/2>=3kT/2 と一致しません。どこで間違えたか分からず困っています。何度か計算してみたのですが、ただの計算ミスでしょうか…? どなたか回答よろしくお願いいたします。

  • 積分の問題を教えていただけませんか。

    質問(1)<V>= 4π(m/2πkT)^(3/2)∫[0,∞]V^3exp(-mV^2/2kT)dV =√(8kT/mπ) こうなるようなのですが、これはなぜでしょうか?途中式など教えてください。

  • 積分の問題

    d/dx∫(x^2→2x+1)e^(-kt^2)dt の途中式を教えてください。