|
|
本帖最后由 Hamilton1758 于 2025-7-12 23:14 编辑
如题，在UM06-2X/6-311++G(d,p)水平下对反应H+O2X (单重态氧气)中的反应物和过渡态进行几何优化。
参考Zou老师之前回复过的帖子http://bbs.keinsci.com/thread-22039-1-1.html，计算得到CASPT2(1,1)/cc-pvtz下H为-0.49980981, CASPT2(12,12)/cc-pvtz下单重态氧气为-150.09597685，三重态氧气为-150.13322099。
然而当我将H和O2X放在一起（距离15A）作为reactant，自旋多重度设为2下计算CASPT2(13,13)/cc-pvtz单点能-150.63304444，自旋多重度设为4下计算CASPT2(13,13)/cc-pvtz单点能150.63303710。接着计算过渡态单点能得到：自旋多重度设为2下，CASPT2(13,13)/cc-pvtz单点能-150.63514444，加上关键词（ist=2,readuhf='TS_uhf.fch', crazywfn）后单点能-150.63514429，自旋多重度设为4下，CASPT2(13,13)/cc-pvtz单点能-150.61677315。
综上，我认为和之前帖子中出现的情况有类似之处，将H和O2X放在一起的时候，无法正确处理结构中单重态氧气。
附件是所有的MOKIT输入输出文件 (氢原子, 三重态氧气，单重态氧气，H+O2X反应物，过渡态)，恳请各位大佬指教。@zjxitcc
TS_CASPT2_cc-pvtz_13_13_multi4.zip
(1014.27 KB, 下载次数 Times of downloads: 1)
TS_CASPT2_cc-pvtz_13_13_readuhf.zip
(638.45 KB, 下载次数 Times of downloads: 0)
TS_CASPT2_cc-pvtz_13_13.zip
(1007.56 KB, 下载次数 Times of downloads: 0)
reactant_CASPT2_cc-pvtz_13_multi4.zip
(869.29 KB, 下载次数 Times of downloads: 0)
O2_CASPT2_cc-pvtz_12.zip
(3.56 MB, 下载次数 Times of downloads: 0)
reactant_CASPT2_cc-pvtz_13.zip
(919.04 KB, 下载次数 Times of downloads: 0)
O2X_CASPT2_cc-pvtz_12.zip
(2.5 MB, 下载次数 Times of downloads: 0)
H_CASPT2_cc-pvtz_1.zip
(82.46 KB, 下载次数 Times of downloads: 1)
|
|
发表于 Post on 2025-7-12 23:13:42
收藏 Add to favorites