# Electron_Helmholtz_energy.py Documentation

## 1. プログラムの動作

このプログラムは、VASP（Vienna Ab initio Simulation Package）の`DOSCAR`ファイルから得られる電子状態密度（DOS）を用いて、有限温度における電子のエネルギーとエントロピーの寄与を計算します。具体的には、以下のステップを実行します：

- `DOSCAR`ファイルからエネルギーとDOSデータを読み込みます。
- フェルミ・ディラック分布を用いて、温度に依存する電子の数やエネルギーを計算します。
- ニュートン法を用いて、温度におけるフェルミエネルギーを計算します。
- 各温度における内部エネルギー、エントロピー、自由エネルギーを計算し、結果を出力します。
- 計算結果をプロットし、グラフとして表示します。

## 2. 必要な非標準ライブラリとインストールコマンドとインストール方法

このプログラムは以下の非標準ライブラリを使用します：

- `numpy`
- `scipy`
- `matplotlib`

これらのライブラリは、以下のコマンドでインストールできます：

```bash
pip install numpy scipy matplotlib
```

## 3. 必要な入力ファイル

プログラムの実行には、VASPの`DOSCAR`ファイルが必要です。このファイルには、エネルギーと対応する状態密度が含まれています。

## 4. 実行後に生成される出力ファイル

このプログラムは、実行後にファイルを生成しませんが、計算結果を標準出力に表示し、以下のグラフをプロットします：

- 温度に対する内部エネルギー、エントロピー、自由エネルギーの変化
- 温度に対するフェルミエネルギーの変化
- エネルギーに対する状態密度（DOS）

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

プログラムはコマンドラインから以下のように実行します：

```bash
python Electron_Helmholtz_energy.py
```

このコマンドを実行する前に、`DOSCAR`ファイルが同じディレクトリに存在することを確認してください。プログラムは、計算結果をコンソールに出力し、グラフを表示します。グラフを閉じるには、Enterキーを押してください。