Gibbs2 两相比较与相界

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

1. 相的能量参考必须一致#

同一化学组成的两相使用同一电子结构设置、同一每单位归一以及可比较的能量参考。各相有自己的E(V)与声子输入，不能只比较零压能量，也不能将任意积分常数当作真实相对能量。

本例使用官方MgO B1/B2输入及其中声明的同一经验修正；它用于演示准谐相比较，不是新的集群DFT或实验相界。

运行与查看#

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

module load gibbs2/1.0.68
bash run.sh official-MgO-B1-B2
bash run.sh MgO-phase-window-p5
bash run.sh MgO-phase-window-p1
ls work/official-MgO-B1-B2
tail -n 12 work/official-MgO-B1-B2/stdout.log

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

inputs/MgO-phase-window-p1/in.ing的运行输入：

root result
mm 40.3044
vfree 2
pressure 450 1 550
temperature 0 50 1000

phase B1 file ../../data/upstream/mgo_pbe/mgo.res tmodel qha \
                 units energy ry freq cm-1 \
                 prefix ../../data/upstream/mgo_pbe/ \
                 eec bpscal 126.025 161.3
phase B2 file ../../data/upstream/mgob2_pbe/mgob2.res tmodel qha \
                 units energy ry freq cm-1 \
                 prefix ../../data/upstream/mgob2_pbe/ \
                 eec use 1
end

2. 看哪些文件#

head -n 12 work/MgO-phase-window-p1/result.dgtp
head -n 12 work/MgO-phase-window-p1/result.tpstab
head -n 12 work/MgO-phase-window-p1/result.ptrans

.dgtp给相对G，.tpstab给网格上稳定相及其G/V/B，.ptrans给根据采样G差线性插值得到的相界。它不是无网格误差的非线性根。

3. 实际压力步长对照#

窄窗口只实际算450–550 GPa；本版本ptrans文本的左端点仍可能打印0，不能把这行当作计算过0至450 GPa的证据。画图按真实压力网格限定范围，再逐步细化所关注区间。

本集群的Si五点固定形状分支也随包提供，但体积区间有限，不用它直接声称收敛的高压相变压力。

真实G差与三套压力网格相界，保留实际小差异。

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

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

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

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

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