Gaussian连续介质溶剂单点计算

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

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-pcm
cp gaussian/water-pcm/input.gjf work/water-pcm/
cd work/water-pcm

2 查看完整输入#

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

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

water-pcm

0 1
O 0.0 0.0 0.0
H 0.0 0.757 0.587
H 0.0 -0.757 0.587

SCRF=(PCM,Solvent=Water)指定连续介质与溶剂。图中气相参照取自相同初始几何的第一次SCF，不与优化后的气相能量混用。这里是电子能差，没有加入完整标准态、热力学和溶剂化自由能修正。

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 'SCF Done|Polarizable Continuum|Solvent' input.log

5 换成自己的体系#

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

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

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

module load ase/3.29.0
python plot.py

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