EF-T-metal.py 技術ドキュメント

プログラムの動作

プログラム EF-T-metal.py は、金属のフェルミ準位（EF）の温度依存性を数値的に計算し、近似式による結果と比較するPythonスクリプトです。このプログラムは、自由電子モデルに基づき、フェルミ・ディラック統計に従う電子系の振る舞いを解析します。

主な機能:

フェルミ準位の計算: 数値積分とニュートン法を組み合わせて、指定された電子数密度を満たすフェルミ準位を正確に計算します。
温度依存性の評価: $0 \text{ K}$ から指定された最高温度までの範囲で、フェルミ準位の温度変化を追跡します。
近似式との比較: 温度依存性フェルミ準位の一般的な近似式による結果も計算し、数値解と比較します。
結果の可視化: Matplotlibライブラリを用いて、計算されたフェルミ準位（数値解と近似解）の温度依存性をグラフとして表示します。
コマンドライン引数による設定: 電子数密度や有効質量をコマンドライン引数で変更できるため、異なる材料条件での解析が可能です。

このプログラムは、高温におけるフェルミ準位の正確な振る舞いを評価し、近似式の適用限界を理解するのに役立ちます。

原理

このプログラムは、金属中の電子が自由電子モデルとフェルミ・ディラック統計に従うという仮定に基づいています。

状態密度 (Density of States, DOS) 3次元の自由電子モデルにおいて、電子の状態密度 $D(E)$ はエネルギー $E$ の平方根に比例します。 $$D(E) = D_0 \sqrt{E}$$ ここで $D_0$ は以下の式で与えられる定数です。 $$D_0 = \frac{(2S+1)2\pi(2m^*)^{1.5}}{h^3}$$ ここで $S$ はスピン（電子の場合は $1/2$）、$m^*$ は有効質量、$h$ はプランク定数です。プログラム中では、単位を $cm^{-3}eV^{-1.5}$ に変換しています。
フェルミ・ディラック分布 (Fermi-Dirac Distribution) 電子が特定のエネルギー準位 $E$ を占有する確率は、フェルミ・ディラック分布関数 $f(E, T, E_F)$ で与えられます。 $$f(E, T, E_F) = \frac{1}{\exp\left(\frac{E - E_F}{k_B T}\right) + 1}$$ ここで $T$ は絶対温度、$E_F$ はフェルミ準位、$k_B$ はボルツマン定数です。
電子数密度とフェルミ準位の決定 全電子数密度 $N$ は、状態密度とフェルミ・ディラック分布関数の積を全エネルギー範囲で積分することによって得られます。 $$N = \int_0^\infty D(E) f(E, T, E_F) dE$$ このプログラムでは、与えられた電子数密度 $N$ を満たす $E_F$ を見つけるためにニュートン法を使用します。積分はフェルミ準位近傍で数値的に (`scipy.integrate.quad`)、それ以外の領域では解析的に行われます。具体的には、関数 $g(E_F) = N_e(T, E_F, dE) - N$ を定義し、$g(E_F) = 0$ となる $E_F$ をニュートン法で探索します。ここで $N_e(T, E_F, dE)$ は現在の $E_F$ で計算される電子数密度です。
$0 \text{ K}$ におけるフェルミ準位 ($E_{F0}$) $T=0 \text{ K}$ では、フェルミ・ディラック分布関数は $E < E_F$ で $1$、$E > E_F$ で $0$ となるステップ関数になります。この場合のフェルミ準位 $E_{F0}$ は解析的に計算できます。 $$N = \int_0^{E_{F0}} D_0 \sqrt{E} dE = \frac{2}{3} D_0 E_{F0}^{3/2}$$ したがって、 $$E_{F0} = \left(\frac{1.5 N}{D_0}\right)^{2/3}$$
フェルミ準位の近似式 高温におけるフェルミ準位の温度依存性に対する近似式も計算されます。これは、低い温度範囲で有効な近似です。 $$E_F(T) \approx E_{F0} - \frac{\pi^2}{6} (k_B T)^2 \left. \frac{1}{D(E)} \frac{dD(E)}{dE} \right|_{E=E_{F0}}$$ プログラム中では、$k_B T$ はeV単位に変換され、$dD(E)/dE$ の部分は $E_{F0}$ で評価されます。

必要な非標準ライブラリとインストール方法

本プログラムの実行には、以下のPython非標準ライブラリが必要です。

numpy: 数値計算を効率的に行うためのライブラリ。
scipy: 科学技術計算ライブラリ。特に数値積分 (scipy.integrate.quad) およびニュートン法 (scipy.optimize.newton) を使用します。
matplotlib: グラフ描画ライブラリ。結果を可視化するために使用します。

これらのライブラリは pip コマンドを使用してインストールできます。

pip install numpy scipy matplotlib

必要な入力ファイル

このプログラムは、実行に必要な入力ファイルを特に必要としません。電子数密度と有効質量はプログラム内の定数として設定されており、コマンドライン引数を通じて変更することが可能です。

生成される出力ファイル

このプログラムは、ファイル出力を生成しません。計算結果は標準出力に表形式で表示され、フェルミ準位の温度依存性を示すグラフがインタラクティブに表示されます。

標準出力の内容:

使用された有効質量 (m) の値。
$0 \text{ K}$ におけるフェルミ準位 (EF at 0K) とその際の電子数密度 (Ne at 0K)。
以下の列を持つ表形式のデータ:
- T(K): 絶対温度 (ケルビン)
- EF(eV): 数値積分とニュートン法で計算されたフェルミ準位 (電子ボルト)
- EFapprox(eV): 近似式で計算されたフェルミ準位 (電子ボルト)
- Ncheck(cm-3): 計算された EF における電子数密度。入力 N と一致しているはずです。

グラフ出力:

温度 (T (K)) をX軸に、フェルミ準位 (EF (eV)) をY軸にとったグラフが表示されます。
「EF(exact)」として数値解が、「EF(approx)」として近似解がプロットされます。
グラフウィンドウを閉じるか、ターミナルでEnterキーを押すまで表示が維持されます。

コマンドラインでの使用例 (Usage)

基本的な実行コマンドは以下の形式です。

python EF-T-metal.py [N_electron] [effective_mass]

N_electron: 電子数密度を指定します (単位: $cm^{-3}$)。デフォルト値は 5.0e22 です。
effective_mass: 有効質量を電子質量 $m_e$ の倍数で指定します。デフォルト値は 1.0 です。

引数を指定しない場合、デフォルト値が使用されます。

コマンドラインでの具体的な使用例

1. デフォルト設定での実行

電子数密度 $5.0 \times 10^{22} \text{ cm}^{-3}$、有効質量 $1.0 \, m_e$ のデフォルト設定でプログラムを実行します。

python EF-T-metal.py

実行結果の説明:

標準出力に $0 \text{ K}$ でのフェルミ準位と、指定された温度範囲（$0 \text{ K}$ から $4000 \text{ K}$ まで $200 \text{ K}$ ステップ）におけるフェルミ準位の数値解および近似解が表形式で出力されます。
最終的に、これらのデータポイントをプロットしたグラフウィンドウが表示されます。グラフには、数値的に計算された厳密なフェルミ準位と、近似式によるフェルミ準位が異なる色でプロットされ、温度に対するフェルミ準位の変化が視覚的に確認できます。グラフはユーザーがEnterキーを押すまで開いたままになります。

2. 電子数密度を変更して実行

電子数密度を $8.0 \times 10^{22} \text{ cm}^{-3}$ に変更して実行します。有効質量はデフォルトの $1.0 \, m_e$ です。

python EF-T-metal.py 8.0e22

実行結果の説明:

デフォルト設定の実行例と同様に、表形式のデータとグラフが表示されますが、電子数密度が高くなったため、フェルミ準位の全体的な値が上昇します。

3. 電子数密度と有効質量の両方を変更して実行

電子数密度を $5.0 \times 10^{22} \text{ cm}^{-3}$ に、有効質量を $0.5 \, m_e$ に変更して実行します。

python EF-T-metal.py 5.0e22 0.5

実行結果の説明:

この例では、有効質量が電子質量の半分になるため、$D_0$ の値が変化し、その結果としてフェルミ準位の具体的な値も変更されます。表形式データとグラフに反映された、新たな条件でのフェルミ準位の温度依存性が示されます。