求助：做H3体系单点能扫描时，出现l502，把不收敛那一帧单独计算单点能又没问题

木贞 · 发表于 Post on 2023-2-27 16:48:32

本帖最后由木贞于 2023-2-27 16:50 编辑

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

H3_scan.com (223 Bytes, 下载次数 Times of downloads: 3)

wzkchem5 · 发表于 Post on 2023-2-27 17:55:44

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

喵星大佬 · 发表于 Post on 2023-2-27 20:36:02

显然这玩意用单参考的方法结果是没有意义的，体系只有3个电子，直接CAS(3,3)等价于极小基下的FCI

这边随手扫了一下，没几分钟的事情，固定正三角形的H3，从0.55埃扫描到1.23埃

h3cas_scan.rar (27.28 KB, 下载次数 Times of downloads: 5)

木贞 · 发表于 Post on 2023-2-28 10:26:12

喵星大佬发表于 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

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

木贞 · 发表于 Post on 2023-2-28 10:29:07

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，好像是缺少个文件，我去那个文件夹下看了确实没有

喵星大佬 · 发表于 Post on 2023-2-28 10:43:50

本帖最后由喵星大佬于 2023-2-28 10:50 编辑

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

喵星大佬 · 发表于 Post on 2023-2-28 10:45:59

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

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

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

yflchx · 发表于 Post on 2023-2-28 11:54:48

单独算没问题可能是初始猜测的缘故。

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

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

浏览过的版块