orca BSSE 校正TRAH-SCF procedure 不收敛

hzfish · 发表于 Post on 2024-11-11 19:36:43

在进行BSSE校正时，TRAH-SCF procedure 不收敛，计算过程中有如下WARNING：
WARNING: 54 diagonal Hessian elements are negative!
WARNING : negative HOMO - LUMO gap : -1.109443
WARNING : negative HOMO - LUMO gap : -1.109443

WARNING: 24 diagonal Hessian elements are negative!
WARNING : negative HOMO - LUMO gap : -0.876783
WARNING : negative HOMO - LUMO gap : -0.876783

有很多。

使用的orca6.0.1
inp文件文件使用Multiwfn生成。
! wB97M-V def2-TZVP def2/J RIJCOSX strongSCF noautostart miniprint nopop

$new_job
! wB97M-V def2-TZVP def2/J RIJCOSX strongSCF noautostart miniprint nopop Pmodel

wzkchem5 · 发表于 Post on 2024-11-11 19:52:34

加NoTRAH试试，或者加defgrid3

hzfish · 发表于 Post on 2024-11-11 20:02:34

wzkchem5 发表于 2024-11-11 19:52
加NoTRAH试试，或者加defgrid3

谢谢。

很奇怪，我单独把这部分取出计算单点，同样的设置，计算没有问题。在BSSE校正计算时，就不收敛。

ultramanm87 · 发表于 Post on 2024-11-12 00:18:04

我在我的体系中也发现wB97M-V在BSSE的时候难收敛，最后换了个泛函比如wB97X-V

hzfish · 发表于 Post on 2024-11-12 08:58:42

ultramanm87 发表于 2024-11-12 00:18
我在我的体系中也发现wB97M-V在BSSE的时候难收敛，最后换了个泛函比如wB97X-V

估计还是有bug，6.0.0版本BSSE结果更离谱。

wzkchem5 · 发表于 Post on 2024-11-12 09:55:55

ultramanm87 发表于 2024-11-11 17:18
我在我的体系中也发现wB97M-V在BSSE的时候难收敛，最后换了个泛函比如wB97X-V

Counterpoise计算因为需要在几乎没有电子密度的地方放基函数，会比较考验泛函在电子密度极低时的数值稳定性和连续性。其实如果用def2-TZVP的话，可以考虑改用gCP校正，gCP校正是支持def2-TZVP的

ultramanm87 · 发表于 Post on 2024-11-12 11:01:15

wzkchem5 发表于 2024-11-12 09:55
Counterpoise计算因为需要在几乎没有电子密度的地方放基函数，会比较考验泛函在电子密度极低时的数值稳定 ...

感谢建议，但是当时我的体系中，存在阴离子而且原子数比较大，需要在部分原子上使用ma-def2-TZVP（全部加内存不够，用更大的基组直到能够忽视BSSE内存也不够），此时gCP是不是就不合适了

sobereva · 发表于 Post on 2024-11-12 12:50:08

ultramanm87 发表于 2024-11-12 11:01
感谢建议，但是当时我的体系中，存在阴离子而且原子数比较大，需要在部分原子上使用ma-def2-TZVP（全部加 ...

是的

如果负电荷分布不集中，或者和相互作用关键区域离得较远，不带弥散也可以接受

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

[ORCA] orca BSSE 校正TRAH-SCF procedure 不收敛