Gibbs2 电子热贡献的使用范围

下载Gibbs2实操包。解压后进入包目录；完整输入、原始数据、提取脚本与高清图随包提供。本专题使用1.0.68固定源码版本、GNU Fortran13.3构建的双精度串行CPU入口，配套Python/Matplotlib与Gnuplot处理结果。每个流程从新的目录开始。

1. 这些关键字表示什么#

ELEC有Sommerfeld自由电子、给定费米能DOS和POL4形式的输入。本版本能读入这些关键字，但当前热力学输出中的Fel/Uel/Sel/Cv,el仍为零；不能用这些列声称真实电子贡献为零，也不能据此完成高温金属电子自由能计算。

运行与查看#

在已分配的CPU资源内，从包根目录执行：

module load gibbs2/1.0.68
bash run.sh official-electronic-Al
ls work/official-electronic-Al
tail -n 12 work/official-electronic-Al/stdout.log

程序每次1个串行进程；不需要GPU或mpirun。原始输入、.eos_static、.eos与.efit按条件隔离保存。

inputs/official-electronic-Al/in.ing的运行输入：

# Aluminium
root 001_elec
set notrans
mm 26.981538
vfree 1
pressure 0 
temperature -1
volume input

# free electron model
phase Al:free file ../../data/upstream/al_lda/al.res tmodel qha phfield 12 \
         units energy ry freq cm-1 \
         prefix ../../data/upstream/al_lda/ eec bpscal 112.04 72.7 \
         elec sommerfeld free nelec 3

# free electron model with calculated N(Ef)
phase Al:somm_nef file ../../data/upstream/al_lda/al.res tmodel qha phfield 12 \
         units energy ry freq cm-1 edos ev \
         prefix ../../data/upstream/al_lda/ eec bpscal 112.04 72.7 \
         elec sommerfeld nelec 3

# Fel 4th-order polynomial
phase Al:pol4 file ../../data/upstream/al_lda/al.res tmodel qha phfield 12 \
         units energy ry freq cm-1 \
         prefix ../../data/upstream/al_lda/ eec bpscal 112.04 72.7 \
         elec pol4 4
end

2. 实际Al例子#

三个输入在298.15 K都打印电子项0。以其中自由电子近似、3电子和实际输出体积建立独立低温Sommerfeld参考，F电子项约为-0.040 kJ/mol，说明程序零列不等于物理零修正。

head -n 20 work/official-electronic-Al/001_elec.eos
ls figures-reference/output-limits.png

3. 当前怎样使用#

本专题用此入口处理静态EOS和明确的振动项。若目标包含电子热贡献，应在能输出该贡献的电子结构/热力学工具里完成，并保持能量参考、单位和温度定义一致；不要将这里的零列直接加入发表级总自由能结果。外部振动F/S/Cv接口也不等同于自动补全电子分项。

原生电子零字段与独立模型参考分开标示；BM推荐入口也在同图说明。

从原始结果重新提取、画图#

若想先看配套数据，在另一份新解压目录执行：

bash analyze.sh
ls analysis-reference
ls figures-reference

它读取 reference/ 原始输出，重新提取表格并生成16张高清图。完整重算时，在新目录运行 bash run-all.sh，最终数据在 analysis/、图片在 figures/。同名结果目录已存在时会受到保护；只算一个条件时，不要调用需要全部条件的整包分析。

相关入口：Gibbs2专题、固定版本源码与手册。以本专题对应的集群版本与输入为准；有限体积、温压和频率网格是操作示例，真实材料计算需要检查相应收敛范围。