计算化学公社

标题: 求助：做H3体系单点能扫描时，出现l502，把不收敛那一帧单独计算单点能又没问题 [打印本页]

作者
Author: 木贞 时间: 2023-2-27 16:48
标题: 求助：做H3体系单点能扫描时，出现l502，把不收敛那一帧单独计算单点能又没问题
本帖最后由木贞于 2023-2-27 16:50 编辑

想算一下H3的势能面看一下，固定了H-H-H键角为60°，然后改变两个H-H键的距离进行扫描。当两个键分别等于1.1和1.2Å时出现l502，把这一帧单独拿出来计算单点能没有问题。想问下怎么解决扫描时出现的问题。com文件上传了，log上传不了截图了最后一段，目前试过添加SCF=NoVarAcc还是不行，期待大佬回复。 (, 下载次数 Times of downloads: 3) (, 下载次数 Times of downloads: 4)

作者
Author: wzkchem5 时间: 2023-2-27 17:55
H-H-H正三角形的结构的电子态是二重简并的，单参考态方法第一收敛困难，第二即使能收敛结果也是定性错误的。应该用多参考态方法，考虑到体系很小，直接用full CI就行

作者
Author: 喵星大佬 时间: 2023-2-27 20:36
显然这玩意用单参考的方法结果是没有意义的，体系只有3个电子，直接CAS(3,3)等价于极小基下的FCI

这边随手扫了一下，没几分钟的事情，固定正三角形的H3，从0.55埃扫描到1.23埃
(, 下载次数 Times of downloads: 5)

作者
Author: 木贞 时间: 2023-2-28 10:26

喵星大佬发表于 2023-2-27 20:36
显然这玩意用单参考的方法结果是没有意义的，体系只有3个电子，直接CAS(3,3)等价于极小基下的FCI

这边随 ...

%chk=H3_scan
#p casscf(3,3)/aug-cc-pvqz scan geom=NoCrowd

H3

0 2
H
H                1          B1
H                2          B2 1          A1

B1 0.6
B2 0.6
A1 60.0

我输入文件这样写的，跑出来结果不对，显式是没有收敛

作者
Author: 木贞 时间: 2023-2-28 10:29

wzkchem5 发表于 2023-2-27 17:55
H-H-H正三角形的结构的电子态是二重简并的，单参考态方法第一收敛困难，第二即使能收敛结果也是定性错误的 ...

%chk=H3_scan
#p FullCI/aug-cc-pvqz scan geom=NoCrowd

H3

0 2
H
H                1          B1
H                2          B2 1          A1

B1 0.6 16 0.05
B2 0.6 16 0.05
A1 60.0

请问输入文件是这样写吗，但我跑的时候报错，显示No executable for file l921.exe，好像是缺少个文件，我去那个文件夹下看了确实没有

作者
Author: 喵星大佬 时间: 2023-2-28 10:43
本帖最后由喵星大佬于 2023-2-28 10:50 编辑

瓜丝没有FullCI
而且你的扫法不对，保持正三角形要用广义内坐标

作者
Author: 喵星大佬 时间: 2023-2-28 10:45

木贞发表于 2023-2-28 10:26
%chk=H3_scan
#p casscf(3,3)/aug-cc-pvqz scan geom=NoCrowd

瓜丝的cas也就图一乐，收敛性并不咋样

你先摆个初始结构，然后单点算对了再读波函数开始扫描

作者
Author: yflchx 时间: 2023-2-28 11:54
单独算没问题可能是初始猜测的缘故。

这种体系高精度的势能面应该早就有了，就像前面已回复的说法：低level的计算没有什么实在的意义。

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