计算化学公社

标题: 求助：文献里水分子解离的势能曲线是怎么算的 [打印本页]

作者
Author: 木贞 时间: 2025-2-25 11:03
标题: 求助：文献里水分子解离的势能曲线是怎么算的
(, 下载次数 Times of downloads: 4)
如图是反应性力场水模型和量化计算的水分子解离曲线对比，他使用的是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左右

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

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

作者
Author: 木贞 时间: 2025-2-25 15:36

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

感谢老师

欢迎光临计算化学公社 (http://bbs.keinsci.com/)