Gaussian水二聚体与BSSE修正

下载小体系输入、真实结果与绘图数据，在自己的目录解压：

unzip small-examples-*.zip
cd small-examples-20261027

1 准备个人环境和算例#

Gaussian 16 C.01使用厂商提供的AVX二进制，按用户独立初始化，不使用公共Gaussian module。先运行 yskit gaussian-init；安装及配额要求见Gaussian环境说明。

yskit gaussian-init
mkdir -p work/water-dimer-bsse
cp gaussian/water-dimer-bsse/input.gjf work/water-dimer-bsse/
cd work/water-dimer-bsse

2 查看完整输入#

下列参考输入对应四线程CPU计算。资源行与计算方法分开：%Mem是主机内存，坐标单位为Å，0 1表示中性单重态。

%chk=result.chk
%NProcShared=4
%Mem=2GB
#p B3LYP/6-31G(d) Counterpoise=2 Int=UltraFine SCF=Tight

water-dimer-bsse

0 1 0 1 0 1
O(Fragment=1) 0.00000000 0.00000000 0.00000000
H(Fragment=1) 0.00000000 0.00000000 0.96000000
H(Fragment=1) 0.93000000 0.00000000 -0.24000000
O(Fragment=2) 0.00000000 0.00000000 2.85000000
H(Fragment=2) 0.75700000 0.00000000 3.43700000
H(Fragment=2) -0.75700000 0.00000000 3.43700000

Counterpoise=2启用两个片段；Fragment=1、Fragment=2标明归属，电荷行依次为总体系及两个片段的电荷和多重度。修正相互作用能=Counterpoise corrected energy−sum of fragments。这里片段保持二聚体中的几何，没有加入单体形变能，也不是含热力学校正的结合自由能。

3 提交并检查结果#

公共8173m CPU模板按56核申请，提交前将输入线程数改为56。实际参考结果使用四线程CPU后端；这里只比较科学结果，不提供CPU节点性能结论。请在CPU分区可用时提交，不在登录终端直接启动正式计算。

sed -i 's/%NProcShared=4/%NProcShared=56/' input.gjf
yskit gaussian-init
sbatch sub_gaussian_8173m.sh input.gjf
squeue -u "$USER"
grep 'Normal termination' input.log
grep 'SCF Done' input.log

4 读取真实参考结果#

grep -E 'Counterpoise corrected|BSSE energy|sum of fragments' input.log

5 换成自己的体系#

修改坐标、电荷和自旋多重度，并检查方法与基组是否适合元素及研究问题；增加原子数前先用小模型确认收敛。优化、频率和能量对照保持可比的方法设置，不混用不同几何或基组下的总能量。

重新绘制下载包中的数据图#

回到解压包根目录执行，图片输出到新的 replot 目录：

module load ase/3.29.0
python plot.py

脚本使用NumPy和Matplotlib，参考原始数值在 results.json，轨迹和工具输出在相应软件目录。