求助：文献里水分子解离的势能曲线是怎么算的

木贞

如图是反应性力场水模型和量化计算的水分子解离曲线对比，他使用的是X3LYP/ 6-311G**，软件是Jaguar 7.0 program package，解离能大概是115kcal/mol, 我用高斯同样的基组方法扫出来相差很大是为什么，不管柔性还是刚性扫描都差很多，大概有170kcal/mol，这是我的设置文件：
%chk=h2o_scan.chk%nosave
#p X3LYP/6-311G** Opt=ModRedundant nosymm

H2O dissociation potential energy scan

0 1
O                  0.00000000    0.00000000    0.00000000
H                  0.00000000    0.00000000    0.76000000
H                  0.90493583    0.00000000   -0.32045459

1 2 S 20 0.1

请问水分子解离曲线应该怎么扫呢，看了几篇文献解离能都是115kcal/mol左右

sobereva

你这么算根本没法得到对称破缺态，解离的结构下依然还是闭壳层，明显势能曲线是错的。仔细看下文
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82
以及北京科音中级量子化学培训班（http://www.keinsci.com/KBQC）这页幻灯片：


PS：X3LYP/6-311G**已经完全过时了，数据如今没有参考意义。现在起码应当用UCCSD(T)/cc-pVQZ（刚性扫描，或柔性扫描后算最终轨迹的各个点的单点能时）

木贞

sobereva 发表于 2025-2-25 11:45
你这么算根本没法得到对称破缺态，解离的结构下依然还是闭壳层，明显势能曲线是错的。仔细看下文
谈谈片段 ...

感谢老师