フェルミオンの挙動の統計的記述
フェルミ=ディラック統計は、 パウリの排他原理 に従う 、相互作用しない多数の 同一粒子からなる 系 の 物理学 に適用される 量子統計 の一種です。その結果、 エネルギー状態 上の粒子のフェルミ=ディラック分布が得られます 。これは、1926年に エンリコ・フェルミ と ポール・ディラック がそれぞれ独立にこの分布を導出したことにちなんで名付けられました。 [1] [2]フェルミ=ディラック統計は 統計 力学の分野の一部であり、 量子力学 の原理を用いています 。
フェルミ・ディラック統計は、 熱力学的平衡状態 にある、 半整数 スピン (1/2、3/2など)を持つ、同一で区別のつかない粒子( フェルミオン )に適用されます。粒子間の相互作用が無視できる場合、システムは単一粒子の エネルギー状態 によって記述できます。その結果、これらの状態にわたる粒子のフェルミ・ディラック分布が生成され、2つの粒子が同じ状態を占めることができなくなり、システムの特性に大きな影響を与えます。フェルミ・ディラック統計は、 スピン1/2 を持つフェルミオンの一種である 電子 に最も一般的に適用されます。
フェルミ=ディラック統計に対応するのが ボーズ=アインシュタイン統計であり、これは ボソン と呼ばれる整数スピン(0、1、2など)を持つ同一かつ区別不可能な粒子に適用されます 。古典物理学では、 マクスウェル=ボルツマン統計は、 同一でありながら区別可能とみなされる粒子を記述するために使用されます。ボーズ=アインシュタイン統計とマクスウェル=ボルツマン統計の両方において、フェルミ=ディラック統計とは異なり、複数の粒子が同じ状態を占めることができます。
熱分布 整数スピンを持つ粒子(ボソン、赤)、半整数スピンを持つ粒子(フェルミオン、青)、そして古典的(スピンレス)粒子(緑)の平衡熱分布。平均占有率は、 系の化学ポテンシャルに対する エネルギーに対して示されている 。ここで 、は系の温度、は ボルツマン定数である。 ⟨ n ⟩ {\displaystyle \langle n\rangle } ϵ {\displaystyle \epsilon } μ {\displaystyle \mu } T {\displaystyle T} k B {\displaystyle k_{B}}
歴史 1926年にフェルミ=ディラック統計が導入される以前は、一見矛盾する現象のために、電子の挙動のいくつかの側面を理解することは困難でした。例えば、 室温 での金属の電子 熱容量は、 電流 中の 電子数 の100分の1の電子から生じているように見えました 。 [3]また、室温で金属に高電界を印加することで発生する 放出電流が 温度にほとんど依存しない理由も理解が困難でした 。
当時の金属電子論である ドルーデ模型 が直面した困難は、電子が(古典的な統計理論によれば)すべて等価であるとみなしていたことによるものでした。言い換えれば、各電子は ボルツマン定数 k B のオーダーの比熱に寄与すると考えられていました。この問題は、フェルミ=ディラック統計が開発されるまで未解決のままでした。
フェルミ・ディラック統計は、1926年に エンリコ・フェルミ [1] と ポール・ディラック [2] によって初めて発表されました。 マックス・ボルン によると 、 パスクアル・ジョルダンは 1925年に同じ統計を開発し、 パウリ 統計 と名付けましたが、タイムリーに発表されませんでした。 [4] [5] [6] ディラックによると、この統計はフェルミによって初めて研究され、ディラックはそれを「フェルミ統計」と呼び、対応する粒子を「フェルミオン」と呼びました。 [7]
フェルミ・ディラック統計は1926年に ラルフ・ファウラーによって 白色矮星 への 恒星 の崩壊を記述するために応用されました 。 [8] 1927年に アーノルド・ゾンマーフェルトは それを金属中の電子に適用して 自由電子モデル を開発しました。 [9] そして1928年にファウラーと ローター・ノルドハイムはそれを金属からの 電界電子放出 に適用しました 。 [10] フェルミ・ディラック統計は物理学の重要な部分であり続けています。
フェルミ・ディラック分布 熱力学的平衡状態にある同一のフェルミオン系において、単一粒子状態 iにおけるフェルミオンの平均数は フェルミ・ディラック(F-D)分布 で与えられる : [11] [注1]
n ¯ i = 1 e ( ε i − μ ) / k B T + 1 , {\displaystyle {\bar {n}}_{i}={\frac {1}{e^{(\varepsilon _{i}-\mu )/k_{\text{B}}T}+1}},}
ここで、 k B は ボルツマン 定数 、 T は絶対温度 、 ε i は単一粒子状態 i のエネルギー、 μは 全化学ポテンシャル である 。分布は条件
∑ i n ¯ i = N {\displaystyle \sum _{i}{\bar {n}}_{i}=N} は正または負の値をとる可能性がある こと を表すために使用できます。 [12] μ = μ ( T , N ) {\displaystyle \mu =\mu (T,N)} μ {\displaystyle \mu }
絶対温度ゼロにおいて、 μは フェルミエネルギー とフェルミオン1個あたりのポテンシャルエネルギーの 和に等しい。ただし、正のスペクトル密度 近傍 にあることが条件となる。半導体中の電子のようにスペクトルギャップがある場合、対称点 μ は通常 フェルミ準位 、あるいは電子の場合は 電気化学ポテンシャル と呼ばれ、ギャップの中央に位置する。 [13] [14]
フェルミ・ディラック分布は、システム内のフェルミオンの数が十分に大きく、システムにフェルミオンを1つ追加しても μ への影響が無視できる場合にのみ有効です。 [15]フェルミ・ディラック分布は、最大で1つのフェルミオンが各可能な状態を占めることを許容する パウリの排他原理 を使用して導出された ため、結果は次のようになります 。 [注 2] 0 < n ¯ i < 1 {\displaystyle 0<{\bar {n}}_{i}<1}
フェルミ・ディラック分布 エネルギー依存性。T が 高いほど緩やかになる 。示されていないが、 T が高いほど減少する 。 [16] μ {\displaystyle \mu } の温度依存 性 。 ε > μ {\displaystyle \varepsilon >\mu } 状態i における粒子数の 分散 は 、上記の式から計算できる : [17] [18] n ¯ i {\displaystyle {\bar {n}}_{i}}
V ( n i ) = k B T ∂ ∂ μ n ¯ i = n ¯ i ( 1 − n ¯ i ) . {\displaystyle V(n_{i})=k_{\text{B}}T{\frac {\partial }{\partial \mu }}{\bar {n}}_{i}={\bar {n}}_{i}(1-{\bar {n}}_{i}).}
エネルギーに対する粒子の分布 フェルミ関数は 、範囲の様々な温度 に対して F ( ϵ ) {\displaystyle F(\epsilon )} μ = 0.55 eV {\displaystyle \mu =0.55~{\text{eV}}} 2 K ≤ T ≤ 375 K {\displaystyle 2~{\text{K}}\leq T\leq 375~{\text{K}}} フェルミ・ディラック分布から、エネルギーに対する粒子の分布を求めることができる。 [注 3] エネルギーを持つフェルミオンの平均数は、フェルミ・ディラック分布 に 縮退 (つまり、エネルギーを持つ状態の数 ) を掛けることで求められる。 [19] ε i {\displaystyle \varepsilon _{i}} n ¯ i {\displaystyle {\bar {n}}_{i}} g i {\displaystyle g_{i}} ε i {\displaystyle \varepsilon _{i}}
n ¯ ( ε i ) = g i n ¯ i = g i e ( ε i − μ ) / k B T + 1 . {\displaystyle {\begin{aligned}{\bar {n}}(\varepsilon _{i})&=g_{i}{\bar {n}}_{i}\\&={\frac {g_{i}}{e^{(\varepsilon _{i}-\mu )/k_{\text{B}}T}+1}}.\end{aligned}}} のとき、 同じエネルギーを持つフェルミオンが占めることができる状態が複数存在するため、 となる可能性があります 。 g i ≥ 2 {\displaystyle g_{i}\geq 2} n ¯ ( ε i ) > 1 {\displaystyle {\bar {n}}(\varepsilon _{i})>1} ε i {\displaystyle \varepsilon _{i}}
エネルギーの準連続体が 状態 密度 (つまり単位エネルギー範囲当たり単位体積当たりの状態数 [20] )を持つ場合、単位エネルギー範囲当たり単位体積当たりのフェルミオンの平均数は ε {\displaystyle \varepsilon } g ( ε ) {\displaystyle g(\varepsilon )}
N ¯ ( ε ) = g ( ε ) F ( ε ) , {\displaystyle {\bar {\mathcal {N}}}(\varepsilon )=g(\varepsilon )F(\varepsilon ),} ここで フェルミ関数 と呼ばれ、 フェルミ・ディラック分布で使用される 関数 と同じである : [21] F ( ε ) {\displaystyle F(\varepsilon )} n ¯ i {\displaystyle {\bar {n}}_{i}}
F ( ε ) = 1 e ( ε − μ ) / k B T + 1 , {\displaystyle F(\varepsilon )={\frac {1}{e^{(\varepsilon -\mu )/k_{\text{B}}T}+1}},} となることによって
N ¯ ( ε ) = g ( ε ) e ( ε − μ ) / k B T + 1 . {\displaystyle {\bar {\mathcal {N}}}(\varepsilon )={\frac {g(\varepsilon )}{e^{(\varepsilon -\mu )/k_{\text{B}}T}+1}}.}
量子と古典の領域 フェルミ・ディラック分布は、 高温・低粒子密度の極限において、特別な仮定を必要とせずに マクスウェル・ボルツマン分布に近づきます。
粒子密度が低い極限では、 となるため、 あるいはそれと同値となる 。その場合、 となり、これはマクスウェル・ボルツマン統計から得られる結果である。 n ¯ i = 1 e ( ε i − μ ) / k B T + 1 ≪ 1 {\displaystyle {\bar {n}}_{i}={\frac {1}{e^{(\varepsilon _{i}-\mu )/k_{\rm {B}}T}+1}}\ll 1} e ( ε i − μ ) / k B T + 1 ≫ 1 {\displaystyle e^{(\varepsilon _{i}-\mu )/k_{\rm {B}}T}+1\gg 1} e ( ε i − μ ) / k B T ≫ 1 {\displaystyle e^{(\varepsilon _{i}-\mu )/k_{\rm {B}}T}\gg 1} n ¯ i ≈ 1 e ( ε i − μ ) / k B T = N Z e − ε i / k B T {\displaystyle {\bar {n}}_{i}\approx {\frac {1}{e^{(\varepsilon _{i}-\mu )/k_{\rm {B}}T}}}={\frac {N}{Z}}e^{-\varepsilon _{i}/k_{\rm {B}}T}} 高温限界においては、粒子は広いエネルギー範囲にわたって分布するため、各状態(特に を持つ高エネルギー状態 )における占有率は再び非常に小さくなります 。これもまた、マクスウェル・ボルツマン統計に帰着します。 ε i − μ ≫ k B T {\displaystyle \varepsilon _{i}-\mu \gg k_{\rm {B}}T} n ¯ i = 1 e ( ε i − μ ) / k B T + 1 ≪ 1 {\displaystyle {\bar {n}}_{i}={\frac {1}{e^{(\varepsilon _{i}-\mu )/k_{\rm {B}}T}+1}}\ll 1} マクスウェル・ボルツマン統計を フェルミ・ディラック統計の近似として使用できる 古典的領域は、 粒子の位置と 運動量に対する ハイゼンベルクの不確定性原理 によって課せられた限界からはるかに離れた状況を考慮することによって見つけられます。たとえば、半導体物理学では、伝導帯の状態密度がドーピング濃度よりもはるかに高い場合、伝導帯とフェルミ準位間のエネルギーギャップはマクスウェル・ボルツマン統計を使用して計算できます。そうでない場合、ドーピング濃度が伝導帯の状態密度と比較して無視できない場合は、正確な計算のために代わりにフェルミ・ディラック分布を使用する必要があります。そうすると、粒子の 濃度が 、粒子の平均 ド・ブロイ波長 よりもはるかに大きい 平均粒子間隔 に対応するときに、古典的な状況が優勢であることが示されます。 [22] R ¯ {\displaystyle {\bar {R}}} λ ¯ {\displaystyle {\bar {\lambda }}}
R ¯ ≫ λ ¯ ≈ h 3 m k B T , {\displaystyle {\bar {R}}\gg {\bar {\lambda }}\approx {\frac {h}{\sqrt {3mk_{\rm {B}}T}}},} ここで、 h は プランク定数 、 mは 粒子の質量 です 。
T = 300 K (すなわち室温付近) における典型的な金属中の伝導電子の場合、系は古典的領域から大きく離れている 。これは、電子の質量が小さく、金属中の伝導電子の濃度が高い(すなわち、小さい)ためである 。したがって、典型的な金属中の伝導電子にはフェルミ・ディラック統計が必要である。 [22] R ¯ ≈ λ ¯ / 25 {\displaystyle {\bar {R}}\approx {\bar {\lambda }}/25} R ¯ {\displaystyle {\bar {R}}}
古典的領域に属さない系のもう一つの例は、白色矮星へと崩壊した恒星の電子からなる系である。白色矮星の温度は高い(典型的には T = 表面温度は 10000 Kに も 達するが[23] 、電子濃度が高く各電子の質量が小さいため古典近似は使えず、やはりフェルミ・ディラック統計が必要となる。 [8]
派生
グランドカノニカルアンサンブル フェルミ・ディラック分布は相互作用しないフェルミオンの量子系にのみ適用され、 グランドカノニカル集団 から簡単に導出できる。 [24] この集団では、系はエネルギーを交換し、リザーバーと粒子を交換することができる(温度 T と化学ポテンシャル μは リザーバーによって固定される)。
相互作用しない性質のため、利用可能な各単粒子準位(エネルギー準位 ϵ )は、リザーバーと接触する独立した熱力学系を形成します。言い換えれば、各単粒子準位は独立した小さなグランドカノニカル集団です。パウリの排他原理により、 単粒子準位には2つの ミクロ状態、すなわち粒子なし(エネルギー E = 0)と粒子1個(エネルギー E = ε )しかありません。したがって、この単粒子準位の 分配関数は、 以下の2つの項のみで構成されます。
Z = exp ( 0 ( μ − ε ) / k B T ) + exp ( 1 ( μ − ε ) / k B T ) = 1 + exp ( ( μ − ε ) / k B T ) , {\displaystyle {\begin{aligned}{\mathcal {Z}}&=\exp {\big (}0(\mu -\varepsilon )/k_{\rm {B}}T{\big )}+\exp {\big (}1(\mu -\varepsilon )/k_{\rm {B}}T{\big )}\\&=1+\exp {\big (}(\mu -\varepsilon )/k_{\rm {B}}T{\big )},\end{aligned}}} そして、その単一粒子レベルのサブステートの平均粒子数は次のように与えられる。
⟨ N ⟩ = k B T 1 Z ( ∂ Z ∂ μ ) V , T = 1 exp ( ( ε − μ ) / k B T ) + 1 . {\displaystyle \langle N\rangle =k_{\rm {B}}T{\frac {1}{\mathcal {Z}}}\left({\frac {\partial {\mathcal {Z}}}{\partial \mu }}\right)_{V,T}={\frac {1}{\exp {\big (}(\varepsilon -\mu )/k_{\rm {B}}T{\big )}+1}}.} この結果は各単一粒子レベルに適用され、システム全体の状態に対するフェルミ・ディラック分布を与える。 [24]
粒子数の分散( 熱変動 による)も導出できる(粒子数は単純な ベルヌーイ分布 に従う)。
⟨ ( Δ N ) 2 ⟩ = k B T ( d ⟨ N ⟩ d μ ) V , T = ⟨ N ⟩ ( 1 − ⟨ N ⟩ ) . {\displaystyle {\big \langle }(\Delta N)^{2}{\big \rangle }=k_{\rm {B}}T\left({\frac {d\langle N\rangle }{d\mu }}\right)_{V,T}=\langle N\rangle {\big (}1-\langle N\rangle {\big )}.} この量は、電子ガス の 電気伝導率や 熱電係数の モットの関係式 などの輸送現象において重要であり 、 [25] エネルギー準位が輸送現象に寄与する能力はに比例します 。 ⟨ ( Δ N ) 2 ⟩ {\displaystyle {\big \langle }(\Delta N)^{2}{\big \rangle }}
標準的なアンサンブル 正準集団 においてフェルミ・ディラック統計を導くことも可能である。 相互相互作用が無視でき、熱平衡状態にある N個 の同一のフェルミオンからなる多粒子系を考える。 [15] フェルミオン間の相互作用は無視できるため、多粒子系の状態 エネルギーは、 単一粒子エネルギーの和として表すことができる。 E R {\displaystyle E_{R}} R {\displaystyle R}
E R = ∑ r n r ε r , {\displaystyle E_{R}=\sum _{r}n_{r}\varepsilon _{r},} ここでは占有数と呼ばれ、 は エネルギー を持つ 一粒子状態にある粒子の数です 。この和は、すべての可能な一粒子状態にわたって求められます 。 n r {\displaystyle n_{r}} r {\displaystyle r} ε r {\displaystyle \varepsilon _{r}} r {\displaystyle r}
多粒子系が状態にある確率は正規化された 正準分布 によって与えられる : [26] R {\displaystyle R}
P R = e − β E R ∑ R ′ e − β E R ′ , {\displaystyle P_{R}={\frac {e^{-\beta E_{R}}}{\displaystyle \sum _{R'}e^{-\beta E_{R'}}}},} ここで 、は ボルツマン因子 と呼ばれ 、その和は 多粒子系のすべての可能な状態にわたっている。占有数の平均値は [26] である。 β = 1 / k B T {\displaystyle \beta =1/k_{\text{B}}T} e − β E R {\displaystyle e^{-\beta E_{R}}} R ′ {\displaystyle R'} n i {\displaystyle n_{i}}
n ¯ i = ∑ R n i P R . {\displaystyle {\bar {n}}_{i}=\sum _{R}n_{i}P_{R}.} 多粒子系の状態は 、単一粒子状態の粒子占有率によって指定できることに注意する。 つまり、 R {\displaystyle R} n 1 , n 2 , … , {\displaystyle n_{1},n_{2},\ldots ,}
P R = P n 1 , n 2 , … = e − β ( n 1 ε 1 + n 2 ε 2 + ⋯ ) ∑ n 1 ′ , n 2 ′ , … e − β ( n 1 ′ ε 1 + n 2 ′ ε 2 + ⋯ ) , {\displaystyle P_{R}=P_{n_{1},n_{2},\ldots }={\frac {e^{-\beta (n_{1}\varepsilon _{1}+n_{2}\varepsilon _{2}+\cdots )}}{\displaystyle \sum _{{n_{1}}',{n_{2}}',\ldots }e^{-\beta (n_{1}'\varepsilon _{1}+n_{2}'\varepsilon _{2}+\cdots )}}},} そして、の式は次の ようになる。 n ¯ i {\displaystyle {\bar {n}}_{i}}
n ¯ i = ∑ n 1 , n 2 , … n i P n 1 , n 2 , … = ∑ n 1 , n 2 , … n i e − β ( n 1 ε 1 + n 2 ε 2 + ⋯ + n i ε i + ⋯ ) ∑ n 1 , n 2 , … e − β ( n 1 ε 1 + n 2 ε 2 + ⋯ + n i ε i + ⋯ ) , {\displaystyle {\begin{aligned}{\bar {n}}_{i}&=\sum _{n_{1},n_{2},\dots }n_{i}P_{n_{1},n_{2},\dots }\\&={\frac {\displaystyle \sum _{n_{1},n_{2},\dots }n_{i}e^{-\beta (n_{1}\varepsilon _{1}+n_{2}\varepsilon _{2}+\cdots +n_{i}\varepsilon _{i}+\cdots )}}{\displaystyle \sum _{n_{1},n_{2},\dots }e^{-\beta (n_{1}\varepsilon _{1}+n_{2}\varepsilon _{2}+\cdots +n_{i}\varepsilon _{i}+\cdots )}}},\end{aligned}}} ここで、和は パウリの排他原理に従うすべての値の組み合わせについてであり、 各 に対して = 0 または となる 。さらに、 の値の各組み合わせは、 粒子の総数が であるという制約を満たす 。 n 1 , n 2 , … {\displaystyle n_{1},n_{2},\ldots } n r = 0 {\displaystyle n_{r}=0} 1 {\displaystyle 1} r {\displaystyle r} n 1 , n 2 , … {\displaystyle n_{1},n_{2},\ldots } N {\displaystyle N}
∑ r n r = N . {\displaystyle \sum _{r}n_{r}=N.} 合計を整理すると、
n ¯ i = ∑ n i = 0 1 n i e − β ( n i ε i ) ∑ ∑ ( i ) n 1 , n 2 , … e − β ( n 1 ε 1 + n 2 ε 2 + ⋯ ) ∑ n i = 0 1 e − β ( n i ε i ) ∑ ∑ ( i ) n 1 , n 2 , … e − β ( n 1 ε 1 + n 2 ε 2 + ⋯ ) , {\displaystyle {\bar {n}}_{i}={\frac {\displaystyle \sum _{n_{i}=0}^{1}n_{i}e^{-\beta (n_{i}\varepsilon _{i})}\sideset {}{^{(i)}}\sum _{n_{1},n_{2},\dots }e^{-\beta (n_{1}\varepsilon _{1}+n_{2}\varepsilon _{2}+\cdots )}}{\displaystyle \sum _{n_{i}=0}^{1}e^{-\beta (n_{i}\varepsilon _{i})}\sideset {}{^{(i)}}\sum _{n_{1},n_{2},\dots }e^{-\beta (n_{1}\varepsilon _{1}+n_{2}\varepsilon _{2}+\cdots )}}},} ここで、和の符号の上の添え字は 、和が を超えていないことを示し 、和に関連する粒子の総数が であるという制約に従う 。制約 を通してが 依然として に依存することに注意する。 なぜなら、一方で は と が で評価されるのに対し、 他方で はと が で評価されるからである。 表記を簡略化し、 が依然として を通して に依存することを明確に示すために、 定義 する 。 ( i ) {\displaystyle (i)} n i {\displaystyle n_{i}} N i = N − n i {\displaystyle N_{i}=N-n_{i}} ∑ ( i ) {\displaystyle \textstyle \sum ^{(i)}} n i {\displaystyle n_{i}} N i {\displaystyle N_{i}} n i = 0 {\displaystyle n_{i}=0} ∑ ( i ) {\displaystyle \textstyle \sum ^{(i)}} N i = N , {\displaystyle N_{i}=N,} n i = 1 , {\displaystyle n_{i}=1,} ∑ ( i ) {\displaystyle \textstyle \sum ^{(i)}} N i = N − 1. {\displaystyle N_{i}=N-1.} ∑ ( i ) {\displaystyle \textstyle \sum ^{(i)}} n i {\displaystyle n_{i}} N − n i , {\displaystyle N-n_{i},}
Z i ( N − n i ) ≡ ∑ ∑ ( i ) n 1 , n 2 , … e − β ( n 1 ε 1 + n 2 ε 2 + ⋯ ) , {\displaystyle Z_{i}(N-n_{i})\equiv \sideset {}{^{(i)}}\sum _{n_{1},n_{2},\ldots }e^{-\beta (n_{1}\varepsilon _{1}+n_{2}\varepsilon _{2}+\cdots )},} となるので、前の式は について書き直して評価することができます 。 n ¯ i {\displaystyle {\bar {n}}_{i}} Z i {\displaystyle Z_{i}}
n ¯ i = ∑ n i = 0 1 n i e − β ( n i ε i ) Z i ( N − n i ) ∑ n i = 0 1 e − β ( n i ε i ) Z i ( N − n i ) = 0 + e − β ε i Z i ( N − 1 ) Z i ( N ) + e − β ε i Z i ( N − 1 ) = 1 [ Z i ( N ) / Z i ( N − 1 ) ] e β ε i + 1 . {\displaystyle {\begin{aligned}{\bar {n}}_{i}&={\frac {\displaystyle \sum _{n_{i}=0}^{1}n_{i}e^{-\beta (n_{i}\varepsilon _{i})}\,Z_{i}(N-n_{i})}{\displaystyle \sum _{n_{i}=0}^{1}e^{-\beta (n_{i}\varepsilon _{i})}\,Z_{i}(N-n_{i})}}\\&={\frac {0+e^{-\beta \varepsilon _{i}}\,Z_{i}(N-1)}{Z_{i}(N)+e^{-\beta \varepsilon _{i}}\,Z_{i}(N-1)}}\\&={\frac {1}{[Z_{i}(N)/Z_{i}(N-1)]\,e^{\beta \varepsilon _{i}}+1}}.\end{aligned}}} を代入する式を求めるために、 次の近似式 [27] が用いられる。 Z i ( N ) / Z i ( N − 1 ) {\displaystyle Z_{i}(N)/Z_{i}(N-1)}
ln Z i ( N − 1 ) ≃ ln Z i ( N ) − ∂ ln Z i ( N ) ∂ N = ln Z i ( N ) − α i , {\displaystyle {\begin{aligned}\ln Z_{i}(N-1)&\simeq \ln Z_{i}(N)-{\frac {\partial \ln Z_{i}(N)}{\partial N}}\\&=\ln Z_{i}(N)-\alpha _{i},\end{aligned}}} どこ α i ≡ ∂ ln Z i ( N ) ∂ N . {\displaystyle \alpha _{i}\equiv {\frac {\partial \ln Z_{i}(N)}{\partial N}}.}
粒子の数が十分大きく、 粒子がシステムに加えられたときに 化学ポテンシャルの変化が非常に小さい場合、 [28] 両辺に指数関数を適用し 、 N {\displaystyle N} μ {\displaystyle \mu } α i ≃ − μ / k B T . {\displaystyle \alpha _{i}\simeq -\mu /k_{\text{B}}T.} α i {\displaystyle \alpha _{i}}
Z i ( N ) / Z i ( N − 1 ) = e − μ / k B T . {\displaystyle Z_{i}(N)/Z_{i}(N-1)=e^{-\mu /k_{\text{B}}T}.} 上記を の式に代入し 、 の以前の定義を使用して を に 代入すると 、フェルミ・ディラック分布が得られます。 n ¯ i {\displaystyle {\bar {n}}_{i}} β {\displaystyle \beta } 1 / k B T {\displaystyle 1/k_{\text{B}}T} β {\displaystyle \beta }
n ¯ i = 1 e ( ε i − μ ) / k B T + 1 . {\displaystyle {\bar {n}}_{i}={\frac {1}{e^{(\varepsilon _{i}-\mu )/k_{\text{B}}T}+1}}.} マクスウェル・ボルツマン分布 や ボーズ・アインシュタイン分布 と同様に、フェルミ・ディラック分布も 平均値の ダーウィン・ファウラー法 によって導くことができる。 [29]
ミクロカノニカルアンサンブル システムの多重度を直接解析し、 ラグランジュ乗数 を使用することで結果が得られます。 [30]
指数i でラベル付けされた複数のエネルギー準位があり、 各準位のエネルギーは ε i で、合計 n i 個の粒子が含まれているとします。各準位には g i 個の 異なるサブ準位が含まれており、それらはすべて同じエネルギーを持ち、区別可能です。例えば、2つの粒子が異なる運動量(つまり、運動量が異なる方向を向いている)を持つ場合、それらは互いに区別可能ですが、それでも同じエネルギーを持つことがあります。準位 i に関連付けられた g i の値は、そのエネルギー準位の「縮退」と呼ばれます。 パウリの排他原理に よれば、そのようなサブ準位には1つのフェルミオンしか存在できません。
エネルギーレベルのg i 個のサブレベルにn i 個の区別できない粒子を 分配する方法の数は 、サブレベルごとに最大1個の粒子で構成され、 二項係数の 組み合わせ解釈 を使用して次のように与えられます 。
w ( n i , g i ) = g i ! n i ! ( g i − n i ) ! . {\displaystyle w(n_{i},g_{i})={\frac {g_{i}!}{n_{i}!(g_{i}-n_{i})!}}.} たとえば、2 つの粒子を 3 つのサブレベルに分配すると、人口数は 110、101、または 011 となり、合計 3 通りとなり、3!/(2!1!) になります。
占有数n i の集合を実現できる方法の数は、 個々のエネルギーレベルを占有できる方法の積です。
W = ∏ i w ( n i , g i ) = ∏ i g i ! n i ! ( g i − n i ) ! . {\displaystyle W=\prod _{i}w(n_{i},g_{i})=\prod _{i}{\frac {g_{i}!}{n_{i}!(g_{i}-n_{i})!}}.} マクスウェル・ボルツマン統計を 導出する際に用いたのと同じ手順に従い、粒子数とエネルギーが一定であるという制約のもとで、 W が最大となる n i の集合を求めたい。 ラグランジュ乗数 を用いて解を制約し、以下の 関数を形成する。
f ( n i ) = ln W + α ( N − ∑ n i ) + β ( E − ∑ n i ε i ) . {\displaystyle f(n_{i})=\ln W+\alpha \left(N-\sum n_{i}\right)+\beta \left(E-\sum n_{i}\varepsilon _{i}\right).} 階乗のスターリング近似 を用いて、 n i に関して微分し 、結果をゼロに設定し、 n i について解くと、フェルミ・ディラックの人口数が得られます。
n i = g i e α + β ε i + 1 . {\displaystyle n_{i}={\frac {g_{i}}{e^{\alpha +\beta \varepsilon _{i}}+1}}.} マクスウェル・ボルツマン統計の 記事で概説したのと同様の手順で 、熱力学的に と が示されるため 、 最終的に状態が占有される確率は β = 1 k B T {\displaystyle \beta ={\tfrac {1}{k_{\text{B}}T}}} α = − μ k B T {\displaystyle \alpha =-{\tfrac {\mu }{k_{\text{B}}T}}}
n ¯ i = n i g i = 1 e ( ε i − μ ) / k B T + 1 . {\displaystyle {\bar {n}}_{i}={\frac {n_{i}}{g_{i}}}={\frac {1}{e^{(\varepsilon _{i}-\mu )/k_{\text{B}}T}+1}}.}
参照 ウィキメディア コモンズには、フェルミ-ディラック分布 に関連するメディアがあります 。
注記 ^ F-D分布は ロジスティック関数 または シグモイド関数 と呼ばれる数学関数の一種です 。 ^ 同時に 1 つのフェルミオンしか同じ状態を占有できないため、 は 状態が占有されている確率でもある ことに注意してください 。 n ¯ i {\displaystyle {\bar {n}}_{i}} i {\displaystyle i} 0 < n ¯ i < 1 {\displaystyle 0<{\bar {n}}_{i}<1} ^ 状態ではなくエネルギーに対するこれらの分布はフェルミ・ディラック分布と呼ばれることもありますが、この記事ではその用語は使用しません。
参考文献 ^ ab フェルミ、エンリコ (1926)。 「スッラ・クォンティッツァツィオーネ・デル・ガス・ペルフェット・モノアトミコ」。 レンディコンティ・リンセイ (イタリア語)。 3 : 145 – 9. 、 Zannoni, Alberto (1999-12-14). 「単原子理想気体の量子化について」 arXiv : cond-mat/9912229 として翻訳 。 ^ ab ディラック, ポール AM (1926). 「量子力学の理論について」. Proceedings of the Royal Society A. 112 ( 762): 661–77 . Bibcode :1926RSPSA.112..661D. doi : 10.1098/rspa.1926.0133 . JSTOR 94692. ^ (キッテル 1971、249~250ページ) ^ 「科学史:ボーア・ハイゼンベルク・コペンハーゲン会合の謎」 Science-Week 4 (20) . 2000年5月19日. OCLC 43626035. 2009年4月11日時点のオリジナルよりアーカイブ。 2009年1月20日 閲覧 。 ^ Schücking (1999). 「ジョーダン、パウリ、ポリティクス、ブレヒト、そして可変重力定数」. Physics Today . 52 (10): 26. Bibcode :1999PhT....52j..26S. doi : 10.1063/1.882858 . ^ ^ ディラック, ポール・A・M (1967). 『量子力学の原理』(改訂第4版). ロンドン: オックスフォード大学出版局. pp. 210–1 . ISBN 978-0-19-852011-5 。 ^ ab ファウラー, ラルフ・H. (1926年12月). 「高密度物質について」. 王立天文学会月報 . 87 (2): 114–22 . Bibcode :1926MNRAS..87..114F. doi : 10.1093/mnras/87.2.114 . ^ アーノルド・ゾンマーフェルト (1927-10-14)。 「Zur Elektronentheorie der Metalle」[金属の電子理論について]。 Naturwissenschaften (ドイツ語)。 15 (41): 824–32 。 ビブコード :1927NW....15....825S。 土井 :10.1007/BF01505083。 S2CID 39403393。 ^ Fowler, Ralph H. ; Nordheim, Lothar W. (1928-05-01). 「強電界下における電子放出」 Proceedings of the Royal Society A . 119 (781): 173– 81. Bibcode :1928RSPSA.119..173F. doi : 10.1098/rspa.1928.0091 . JSTOR 95023. ^ ライフ1965年、341ページ。 ^ Landau, LD, Lifshitz, EM (2013). 統計物理学: 第5巻 (第5巻). エルゼビア. ^ ブレイクモア 2002、11ページ。 ^ キッテル, チャールズ ; クローマー, ハーバート (1980). 熱物理学(第2版). サンフランシスコ: WHフリーマン. p. 357. ISBN 978-0-7167-1088-2 。 ^ ab Reif 1965、340–342 ページ。 ^ Kittel 1971、p.245、図4および5。 ^ ピアソール、トーマス (2020). 『量子フォトニクス』第2版. 物理学大学院テキスト. シュプリンガー. doi :10.1007/978-3-030-47325-9. ISBN 978-3-030-47324-2 。 ^ (Reif 1965, p. 351) 式9.7.7、ここで 。 β = 1 / k B T , α = − μ / k B T , ∂ n ¯ i ∂ ϵ i = − ∂ n ¯ i ∂ μ {\displaystyle \beta =1/k_{\text{B}}T,\quad \alpha =-\mu /k_{\text{B}}T,\quad {\frac {\partial {\bar {n}}_{i}}{\partial \epsilon _{i}}}=-{\frac {\partial {\bar {n}}_{i}}{\partial \mu }}} ^ レイトン、ロバート・B. (1959). 『現代物理学の原理』 マグロウヒル. p. 340. ISBN 978-0-07-037130-9 。 式(1)において、 およびは それぞれ本論文の およびに対応することに注意してください 。339ページの式(32)も参照してください。 n ( ε ) {\displaystyle n(\varepsilon )} n s {\displaystyle n_{s}} n ¯ i {\displaystyle {\bar {n}}_{i}} n ¯ ( ε i ) {\displaystyle {\bar {n}}(\varepsilon _{i})} ^ ブレイクモア 2002、8ページ。 ^ ライフ1965年、389ページ。 ^ ab (ライフ 1965、246~248ページ) ^ 向井浩二、ジム・ロクナー (1997). 「天体物理学者に聞く」 NASAの宇宙を想像せよ . NASAゴダード宇宙飛行センター. 2009年1月18日時点のオリジナルよりアーカイブ。 ^ ab Srivastava, RK; Ashok, J. (2005). 「第6章」. 統計力学 . ニューデリー : PHI Learning Pvt. Ltd. ISBN 9788120327825 。 ^ Cutler, M.; Mott, N. (1969). 「電子ガスにおけるアンダーソン局在の観測」. Physical Review . 181 (3): 1336. Bibcode :1969PhRv..181.1336C. doi :10.1103/PhysRev.181.1336. ^ ab Reif 1965、203–206 ページ。 ^ 例えば、 微分§差分商による定義を 参照。これは近似値を与える。 f ( a + h ) ≈ f ( a ) + f ′ ( a ) h . {\displaystyle f(a+h)\approx f(a)+f'(a)h.} ^ Reif 1965, pp. 341–342. 式9.3.17と 近似値の妥当性に関する注釈を 参照。 ^ Müller-Kirsten, HJW (2013). 統計物理学の基礎 (第2版). World Scientific. ISBN 978-981-4449-53-3 。 ^ ブレイクモア 2002年、343-534頁。
さらに読む