计算化学公社

标题: orca BSSE 校正TRAH-SCF procedure 不收敛 [打印本页]

作者
Author: hzfish 时间: 2024-11-11 19:36
标题: orca BSSE 校正TRAH-SCF procedure 不收敛
在进行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

作者
Author: wzkchem5 时间: 2024-11-11 19:52
加NoTRAH试试，或者加defgrid3

作者
Author: hzfish 时间: 2024-11-11 20:02

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

谢谢。

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

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

作者
Author: hzfish 时间: 2024-11-12 08:58

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

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

作者
Author: wzkchem5 时间: 2024-11-12 09:55

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

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

作者
Author: ultramanm87 时间: 2024-11-12 11:01

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

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

作者
Author: sobereva 时间: 2024-11-12 12:50

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

是的

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

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