Seebeck_scan.py 技術ドキュメント

プログラムの動作

Seebeck_scan.py は、半導体材料のゼーベック係数を、フェルミ準位 (EF)、キャリア密度 (Ne)、温度 (T)、有効質量 (\(m^*/m_0\)) といった主要な物理パラメータに対して計算し、その結果をプロット表示またはCSVファイルに出力するPythonスクリプトです。このプログラムの主な目的は、熱電材料の特性評価において、これらのパラメータがゼーベック係数に与える影響を系統的に解析することにあります。

主な機能:

  • ゼーベック係数の計算: フェルミ・ディラック積分を用いて、与えられた還元フェルミ準位 (\(\xi\)) と散乱因子 (\(r\)) からゼーベック係数を算出します。

  • パラメータスキャン:

    • EFスキャン: フェルミ準位を変化させたときのゼーベック係数を計算します。

    • Neスキャン: キャリア密度を変化させたときのゼーベック係数を計算します。

    • Tスキャン: 温度を変化させたときのゼーベック係数を計算します。この際、EFまたはキャリア密度のいずれかを固定することができます。

    • 有効質量スキャン: 有効質量を変化させたときのゼーベック係数を計算します。この際、EFまたはキャリア密度のいずれかを固定することができます。

  • 散乱因子 (\(r\)) の指定: 複数の散乱因子を一度に指定し、それぞれの条件でのゼーベック係数を比較できます。

  • プロット表示: スキャン結果をグラフとして可視化します。

  • CSV出力: 計算結果の詳細データをCSV形式で保存します。

解決する課題:

このプログラムは、熱電材料設計において、材料の組成やドーピング量、動作温度などがゼーベック係数にどのように影響するかを理論的に予測し、実験データを解釈するためのツールとして役立ちます。特に、キャリア輸送特性を支配する散乱メカニズムの違いを散乱因子 \(r\) の変化でシミュレーションできるため、様々な材料システムへの応用が可能です。

原理

本プログラムは、非縮退および縮退半導体におけるゼーベック係数の一般的な表現に基づいて計算を行います。主要な物理式とアルゴリズムは以下の通りです。

1. フェルミ・ディラック積分

ゼーベック係数やキャリア密度の計算には、フェルミ・ディラック積分が不可欠です。本プログラムでは、以下の定義に基づいて Fj(j, xi) 関数を実装しています。

\[F_j(\xi) = \int_0^\infty \frac{x^j}{1 + \exp(x - \xi)} dx\]

ここで、\(j\) は積分の次数、\(\xi = (E_F - E_c)/(k_B T)\) は還元フェルミ準位です。Fj 関数は、scipy.integrate.quad を用いた数値積分により計算されます。計算効率のため、functools.lru_cache により計算結果がメモ化されています。

2. 有効状態密度 (3次元の場合)

3次元半導体における有効状態密度 \(N_c\) は、以下の式で与えられます。

\[N_c = 2 \left( \frac{2 \pi m^* k_B T}{h^2} \right)^{3/2}\]

ここで、\(m^*\) は有効質量、\(k_B\) はボルツマン定数、\(T\) は絶対温度、\(h\) はプランク定数です。本プログラムでは、Nc_3D 関数でこの値を計算します。

3. 電子密度と還元フェルミ準位の関係

電子密度 \(n\) と還元フェルミ準位 \(\xi\) の関係は、有効状態密度 \(N_c\) とフェルミ・ディラック積分を用いて次のように表されます。

\[n = N_c \cdot \frac{F_{1/2}(\xi)}{\Gamma(3/2)}\]

ここで、\(\Gamma(3/2)\) はガンマ関数です。electron_density_from_xi 関数はこの式に基づいて電子密度を計算します。 逆に、与えられた電子密度 \(n\) から還元フェルミ準位 \(\xi\) を求めるために、xi_from_electron_density 関数では scipy.optimize.brentq アルゴリズム（ブリッジ法）を用いた数値的求根を行います。

4. ゼーベック係数

半導体のゼーベック係数 \(S\) は、還元フェルミ準位 \(\xi\) と散乱因子 \(r\) を用いて以下の式で計算されます。

\[S = -\frac{k_B}{q} \left[ \frac{(r+2) F_{r+1}(\xi)}{(r+1) F_r(\xi)} - \xi \right]\]

ここで、\(q\) はキャリアの電荷（電子の場合は \(-e\)、ホールの場合は \(+e\)）、\(k_B\) はボルツマン定数です。 散乱因子 \(r\) は、キャリアの緩和時間 \(\tau(E)\) がエネルギー \(E\) に対して \(\tau(E) \propto E^{r-1/2}\) のように依存すると仮定したときの指数です。典型的な散乱機構に対応する \(r\) の値は以下の通りです。

  • \(r = 2\): イオン化不純物散乱

  • \(r = 1\): 極性光学フォノン散乱

  • \(r = 0\): 音響フォノン散乱、非極性光学フォノン散乱

  • \(r = -0.5\): 格子欠陥散乱、中性不純物散乱

seebeck_from_xi 関数では、上記の式を用いてゼーベック係数を計算します。キャリアの種類（電子またはホール）に応じて電荷の符号を適切に処理します。

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

Seebeck_scan.py の実行には、以下のPython非標準ライブラリが必要です。

  • numpy: 数値計算

  • matplotlib: グラフ描画

  • scipy: 科学技術計算（数値積分、求根、特殊関数）

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

pip install numpy matplotlib scipy

必要な入力ファイル

Seebeck_scan.py は、コマンドライン引数のみで全てのパラメータを指定するため、特別な入力ファイルは必要ありません。

生成される出力ファイル

本プログラムは、実行時に以下の種類の出力ファイルを生成する可能性があります。

1. CSVファイル

--csv オプションが指定された場合、計算結果がCSV（カンマ区切り値）形式で保存されます。

  • ファイル名: --csv オプションで指定されたファイル名（例: result.csv）。

  • 内容: 各スキャンポイントにおける計算結果が以下のヘッダを持つ行として記録されます。

    • mode: スキャンモード (EF, Ne, T, m_eff)

    • scan_value: スキャン対象のパラメータ値

    • r: 使用された散乱因子

    • S_V_per_K: ゼーベック係数（V/K単位）

    • S_uV_per_K: ゼーベック係数（μV/K単位）

    • xi: 還元フェルミ準位

    • EF_eV: フェルミ準位（バンド端からの相対値、eV単位）

    • Ne_cm3: 電子密度（\(cm^{-3}\)単位）

    • T_K: 温度（K単位）

    • m_eff_m0: 有効質量（\(m_0\)に対する比）

2. プロット（画面表示）

--no-show オプションが指定されていない場合、計算されたゼーベック係数のプロットが matplotlib ウィンドウとして画面に表示されます。このプロットはファイルとしては保存されませんが、ユーザーはGUIを通じて保存することができます。

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

Seebeck_scan.py は、以下の基本的な形式でコマンドラインから実行します。

python Seebeck_scan.py [OPTIONS]

主なオプション:

  • --mode {EF,Ne,T,m_eff}: スキャンするパラメータを指定します。デフォルトは EF。

  • --fixed {EF,Ne,xi}: mode が T または m_eff の場合に、何らかの値を固定するかを指定します。デフォルトは EF。

  • --carrier {electron,hole}: キャリアの種類を指定します。デフォルトは electron。

  • --T FLOAT: 基準温度 [K]。デフォルトは 300.0。

  • --m_eff FLOAT: 基準有効質量 \(m^*/m_0\)。デフォルトは 1.0。

  • --EF FLOAT: 基準EF-Ec [eV]。デフォルトは 0.10。

  • --Ne FLOAT: 基準電子密度 [\(cm^{-3}\)]。デフォルトは 1e18。

  • --xi FLOAT: 基準還元フェルミ準位。デフォルトは 0.0。

  • --r-list TEXT: コンマ区切りの散乱因子 \(r\) のリスト。デフォルトは "2,1.5,1,0.5,0,-0.5"。

  • --scan-min FLOAT: スキャン範囲の最小値。モードに応じてデフォルト値が設定されます。

  • --scan-max FLOAT: スキャン範囲の最大値。モードに応じてデフォルト値が設定されます。

  • --n INT: スキャンポイントの数。デフォルトは 121。

  • --logscan: スキャン値を対数スケールで生成します。

  • --xscale {auto,linear,log}: プロットのX軸スケール。auto の場合、Ne モードや --logscan では log、それ以外では linear になります。

  • --csv FILENAME: 計算結果をCSVファイルに保存します。

  • --no-show: プロットを表示しません。

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

以下に、プログラムのコメントで提示されている具体的な使用例とその実行結果について説明します。

1. EF (フェルミ準位) スキャン

フェルミ準位 \(E_F - E_c\) を -0.2 eV から 0.5 eV まで 121 点でスキャンします。温度は300 K、有効質量は \(m^*/m_0=1.0\) に固定されます。

python Seebeck_scan.py --mode EF --T 300 --m_eff 1.0 --scan-min -0.2 --scan-max 0.5 --n 121

実行結果: X軸に \(E_F - E_c\) (eV)、Y軸に \(S\) (\(\mu\)V/K) をとったグラフが表示されます。複数の散乱因子 \(r\) に対して、ゼーベック係数がどのように変化するかが曲線としてプロットされます。通常、\(E_F - E_c\) が増加（つまりフェルミ準位がバンド端に近づき、キャリア濃度が高くなる）するにつれて、ゼーベック係数の絶対値は減少します。

2. Ne (電子密度) スキャン

電子密度を \(1 \times 10^{15} \text{ cm}^{-3}\) から \(1 \times 10^{21} \text{ cm}^{-3}\) まで121点で対数スケールでスキャンします。温度は300 K、有効質量は \(m^*/m_0=1.0\) に固定されます。

python Seebeck_scan.py --mode Ne --T 300 --m_eff 1.0 --scan-min 1e15 --scan-max 1e21 --n 121 --logscan

実行結果: X軸に \(N_e\) (\(cm^{-3}\)) を対数スケール、Y軸に \(S\) (\(\mu\)V/K) をとったグラフが表示されます。電子密度が増加するにつれて、ゼーベック係数の絶対値は減少する傾向が見られます。

3. T (温度) スキャン、EF固定

温度を 100 K から 800 K まで 121 点でスキャンします。フェルミ準位 \(E_F - E_c\) は 0.10 eV に固定され、有効質量は \(m^*/m_0=1.0\) です。

python Seebeck_scan.py --mode T --fixed EF --EF 0.10 --m_eff 1.0 --scan-min 100 --scan-max 800 --n 121

実行結果: X軸に \(T\) (K)、Y軸に \(S\) (\(\mu\)V/K) をとったグラフが表示されます。温度の上昇に伴うゼーベック係数の変化が観察できます。フェルミ準位を固定した場合、通常、温度の上昇とともにゼーベック係数の絶対値は増加します。

4. T (温度) スキャン、Ne固定

温度を 100 K から 800 K まで 121 点でスキャンします。電子密度は \(1 \times 10^{18} \text{ cm}^{-3}\) に固定され、有効質量は \(m^*/m_0=1.0\) です。

python Seebeck_scan.py --mode T --fixed Ne --Ne 1e18 --m_eff 1.0 --scan-min 100 --scan-max 800 --n 121

実行結果: X軸に \(T\) (K)、Y軸に \(S\) (\(\mu\)V/K) をとったグラフが表示されます。キャリア密度を固定した場合、通常、温度の上昇とともにゼーベック係数の絶対値は減少します。これは、温度の上昇に伴って還元フェルミ準位 \(\xi\) が減少するためです。

5. \(m_{eff}\) (有効質量) スキャン、Ne固定

有効質量を \(0.1 m_0\) から \(3.0 m_0\) まで 121 点でスキャンします。電子密度は \(1 \times 10^{18} \text{ cm}^{-3}\)、温度は 300 K に固定されます。

python Seebeck_scan.py --mode m_eff --fixed Ne --Ne 1e18 --T 300 --scan-min 0.1 --scan-max 3.0 --n 121

実行結果: X軸に \(m^*/m_0\)、Y軸に \(S\) (\(\mu\)V/K) をとったグラフが表示されます。有効質量が増加するにつれて、ゼーベック係数の絶対値が増加する傾向が見られます。

6. CSV保存

電子密度スキャンの結果を result.csv というファイルに保存します。プロットも同時に表示されます。

python Seebeck_scan.py --mode Ne --csv result.csv

実行結果: 上記の「Ne (電子密度) スキャン」と同様のプロットが表示されるとともに、result.csv というファイルが生成されます。このCSVファイルには、各スキャンポイントでの mode, scan_value, r, S_V_per_K, S_uV_per_K, xi, EF_eV, Ne_cm3, T_K, m_eff_m0 といった詳細な計算データが格納されます。