计算化学公社

标题: orca结构优化遇到问题 WARNING : negative HOMO - LUMO gap [打印本页]

作者
Author: lll123 时间: 2024-4-24 20:24
标题: orca结构优化遇到问题 WARNING : negative HOMO - LUMO gap
优化的是带有Zr4+的有机配合物。采用下面的方法进行结构优化
! r2SCAN-3c opt freq tightSCF noautostart miniprint nopop
%maxcore 3000
%pal nprocs 16 end
! CPCM(water)
* xyz 4 1
遇到的问题是：WARNING: 1 diagonal Hessian elements are negative! WARNING : negative HOMO - LUMO
请问各位老师如何解决。是体系的自旋多重度设置的有问题嘛。或者是Zr需要采用赝式基组吗？

作者
Author: wzkchem5 时间: 2024-4-24 21:26
分子结构极其不合理，必须做显式溶剂化，配体该脱的质子必须脱
建议仔细看文献是如何算溶液里的Zr配合物的

作者
Author: lll123 时间: 2024-4-24 21:50

wzkchem5 发表于 2024-4-24 21:26
分子结构极其不合理，必须做显式溶剂化，配体该脱的质子必须脱
建议仔细看文献是如何算溶液里的Zr配合物的

好的。谢谢老师！我去试一下

作者
Author: hzfish 时间: 2024-12-15 21:55
借楼问解决办法。
使用r2SCAN-3c进行优化结构和频率计算，然后使用wB97M-V def2-QZVP def2/J RIJCOSX tightSCF进行单点能计算。
出现多处
WARNING: 134 diagonal Hessian elements are negative!
WARNING : negative HOMO - LUMO gap : -0.502892
WARNING : negative HOMO - LUMO gap : -0.502892
SCF不收敛如何解决？

作者
Author: zjxitcc 时间: 2024-12-15 22:40
本帖最后由 zjxitcc 于 2024-12-17 21:24 编辑

hzfish 发表于 2024-12-15 21:55
借楼问解决办法。
使用r2SCAN-3c进行优化结构和频率计算，然后使用wB97M-V def2-QZVP def2/J RIJCOSX tigh ...

十分容易解决，以下技巧可联用：
（1）记得在%scf中加sthresh 1e-6关键词。Gaussian/Gamess等量化软件都默认使用这个阈值删除线性依赖的基函数，而ORCA默认这个阈值为1e-7。但平时大家进行r2SCAN-3c或wB97M-V/def2-QZVP计算时，都很容易碰到基函数线性依赖，因此这个选项建议作为默认值写进输入文件，而非算了半天检查后回来再加这个关键词。如果体系没有基函数线性依赖，这个关键词不会产生作用。
（2）如果仍在使用ORCA 5，应换成ORCA 6.0.1。
（3）可以先用ORCA做wB97M-V/def2-TZVP计算；若能收敛，轨道给def2-QZVP计算读取。
（4）可以先用Gaussian做wB97X/def2QZVP计算；若能收敛，使用MOKIT的fch2mkl小程序传轨道给ORCA的wB97M-V/def2-QZVP计算读取。若使用此技巧，上述（1）会被自动考虑。
请注意，本人不建议一个个尝试关键词（例如SlowConv KDIIS等等）；若您一定要尝试，碰到报错/不收敛请不要来问我，我只回答我上面提出的建议之疑问。

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