计算化学公社

标题: 高斯优化scf不收敛，尝试过多次，谢谢指点 [打印本页]

作者
Author: 小范范1989 时间: 2018-1-30 17:33
标题: 高斯优化scf不收敛，尝试过多次，谢谢指点
优化计算了如下两个分子，一直scf不收敛，我按照sob老师的帖子http://sobereva.com/61都尝试了。而且我计算是高斯16.一直是不收敛，各位老师还有什么其他好办法吗？谢谢指点。

RMSDP=1.89D-03 MaxDP=2.01D-01 DE=-2.00D-01 OVMax= 4.00D-02

>>>>>>>>>> Convergence criterion not met.
SCF Done: E(RPBE-PBE) = -1431.71077994 A.U. after 129 cycles
NFock=128 Conv=0.19D-02 -V/T= 2.0091
KE= 1.418836190477D+03 PE=-9.325551999246D+03 EE= 3.479196941588D+03
Convergence failure -- run terminated.

作者
Author: 小范范1989 时间: 2018-1-30 19:52
难道是结构问题？需要打破这个平面结构吗？还会自旋多重读问题？

作者
Author: xyyan1207 时间: 2018-1-30 20:13
再难的情况依照sob老师博文中第10种方法基本都能收敛，不过一定要记得多做一步波函数稳定性检测，详情可以查看http://bbs.keinsci.com/thread-9012-1-1.html，这个我之前咨询的帖子。
如果这个方法还解决不了，可以试一下调整结构了

作者
Author: sobereva 时间: 2018-1-30 21:23
你确定边缘的N不需要加氢么？

作者
Author: 小范范1989 时间: 2018-1-30 21:37
本帖最后由小范范1989 于 2018-1-30 21:52 编辑

sobereva 发表于 2018-1-30 21:23
你确定边缘的N不需要加氢么？

对，谢谢sob老师，我弄错了。谢谢老师指点。这个边缘的B是不是也要加H饱和？

那中间破缺的部分，也用H饱和，还是说这样露着？多重读是不是应该改变？

作者
Author: 小范范1989 时间: 2018-1-30 21:38
本帖最后由小范范1989 于 2018-1-30 21:44 编辑

xyyan1207 发表于 2018-1-30 20:13
再难的情况依照sob老师博文中第10种方法基本都能收敛，不过一定要记得多做一步波函数稳定性检测，详情可以 ...

谢谢指点，我看一下这个帖子，按照sob老师说的，我改一下结构，我的结构有问题，之前太粗心了，人为的认为了。
同时按照方法在调节一下，
谢谢指点

作者
Author: xyyan1207 时间: 2018-1-30 21:51

小范范1989 发表于 2018-1-30 21:38
谢谢指点，我看一下这个帖子，按照sob老师说的，我改一下结构，我的结构有问题，之前太粗心了，人为的认 ...

嗯嗯，多试几次就行了，确实感觉边缘的B和N应该都加H饱和才合理

作者
Author: 小范范1989 时间: 2018-1-30 21:59

xyyan1207 发表于 2018-1-30 21:51
嗯嗯，多试几次就行了，确实感觉边缘的B和N应该都加H饱和才合理

恩，谢谢指点。
不知道中间那部分用不用饱和，这个是不是类似石墨烯破缺？中间扣掉的部分，用H饱和吗？还是说露着，然后多重读会改变呢？
第一次计算这样的，谢谢指点。

作者
Author: xyyan1207 时间: 2018-1-30 22:05

小范范1989 发表于 2018-1-30 21:59
恩，谢谢指点。
不知道中间那部分用不用饱和，这个是不是类似石墨烯破缺？中间扣掉的部分，用H饱和吗？ ...

我也不是这个方向的，这个问题我不太确定，感觉中间如果加了H的话，位阻会很大吧。

我之前听过清华大学李隽教授的报告，他做的类似体系有很多，你可以去网上找几篇相关文献读一读，看看能否解决这个问题。

作者
Author: 小范范1989 时间: 2018-1-30 22:09

xyyan1207 发表于 2018-1-30 22:05
我也不是这个方向的，这个问题我不太确定，感觉中间如果加了H的话，位阻会很大吧。

我之前听过清华大 ...

好的，非常感谢老师指点，谢谢老师。

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