计算化学公社

标题: 用MP2、MP3计算氢分子解离时有突变点，如何消除？ [打印本页]

作者
Author: Mulliken 时间: 2022-10-7 13:38
标题: 用MP2、MP3计算氢分子解离时有突变点，如何消除？
本帖最后由 Mulliken 于 2022-10-7 14:49 编辑

用 Gaussian 中的 MP2、MP3 方法计算氢分子解离过程时，发现能量曲线上有一个奇怪的突变点

(, 下载次数 Times of downloads: 15)

所用输入文件为

#p ump3/sto-3g scan unit=bohr geom=nocrowd guess(mix,always)
H2 scan
0 1
H
H,1,r
r .6 72 .05

复制代码

请问是什么原因？用的基组很小，请轻喷... 主要是想复现一下 Szabo 那本 Modern Quantum Chemistry 书里的图 6.2

更新：
突变点位于 Coulson-Fischer point，谢谢 @wzkchem5 的提醒！

(, 下载次数 Times of downloads: 15)

作者
Author: granvia 时间: 2022-10-7 13:51
先检查下对应的UHF的曲线如何？

作者
Author: Mulliken 时间: 2022-10-7 13:55

granvia 发表于 2022-10-7 13:51
先检查下对应的UHF的曲线如何？

对应的 UHF 曲线也画在上面了，看起来似乎没问题

作者
Author: scf 时间: 2022-10-7 14:02
印象里高斯的scan每一个点的默认初始波函数是上一个点收敛的波函数, 不确定guess(mix,always)是不是更优先. 每一个单点用一个独立的文件试试?

作者
Author: hebrewsnabla 时间: 2022-10-7 14:17
不妨先展示一下你画出来的曲线

作者
Author: Mulliken 时间: 2022-10-7 14:19

hebrewsnabla 发表于 2022-10-7 14:17
不妨先展示一下你画出来的曲线

整半天图没贴成功...
重贴了一下，应该能看见了

作者
Author: wzkchem5 时间: 2022-10-7 14:25
突变点是不是位于Coulson-Fischer点？

作者
Author: Mulliken 时间: 2022-10-7 14:28

scf 发表于 2022-10-7 14:02
印象里高斯的scan每一个点的默认初始波函数是上一个点收敛的波函数, 不确定guess(mix,always)是不是更优先. ...

guess(always) 应该能保证每一次都重新做初猜，我等一下也试试单个单个地算

作者
Author: hebrewsnabla 时间: 2022-10-7 14:45
本帖最后由 hebrewsnabla 于 2022-10-7 14:50 编辑

如果你对某个点单独算一个uhf guess=mix stable=opt --link1-- ump2的结果（能量和S^2)和你扫描的结果一致，就说明没有问题。

我看了一下以前算的结果，确实UMP2在FC点处会有一个不易察觉的一阶导不连续（无论什么基组）（如果扫描步长大一点就更不易察觉了

作者
Author: Mulliken 时间: 2022-10-7 14:47
本帖最后由 Mulliken 于 2022-10-7 14:52 编辑

wzkchem5 发表于 2022-10-7 14:25
突变点是不是位于Coulson-Fischer点？

是的！谢谢提醒！如一楼更新的图， CF 点在 2.173 Bohr 附近

作者
Author: Mulliken 时间: 2022-10-7 14:54

hebrewsnabla 发表于 2022-10-7 14:45
如果你对某个点单独算一个uhf guess=mix stable=opt --link1-- ump2的结果（能量和S^2)和你扫描的结果一致 ...

我对这个FC点不是很熟悉，今天第一次真的碰到，感谢确认

作者
Author: Mulliken 时间: 2022-10-7 15:06

hebrewsnabla 发表于 2022-10-7 14:45
如果你对某个点单独算一个uhf guess=mix stable=opt --link1-- ump2的结果（能量和S^2)和你扫描的结果一致 ...

那 Szabo 的图 6.2 画得不对呀...

作者
Author: hebrewsnabla 时间: 2022-10-7 15:37

Mulliken 发表于 2022-10-7 15:06
那 Szabo 的图 6.2 画得不对呀...

我觉得是对的，只不过画得比较粗糙

作者
Author: zjxitcc 时间: 2022-10-8 10:43
你可以试试OO-UMP2，说不定没有这个问题

作者
Author: zjxitcc 时间: 2022-10-8 10:44
你可以试试OO-UMP2，说不定没有这个问题

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