计算化学公社

标题: 请教几何优化MOF晶胞中截取的一个重复片段时SCF不收敛的问题 [打印本页]

作者
Author: fdy哈哈哈 时间: 2022-3-30 12:13
标题: 请教几何优化MOF晶胞中截取的一个重复片段时SCF不收敛的问题
请教各位老师，我在（G09.Linux）用#p opt b3lyp/6-31G** em=gd3bj test 进行几何结构优化后，出现报错：
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(RB3LYP) =  -12698.4901134    A.U. after  129 cycles
         NFock=128  Conv=0.73D-04    -V/T= 2.0031
KE= 1.265961672921D+04 PE=-5.174699784148D+04 EE= 1.564442944061D+04
Convergence failure -- run terminated.
Error termination via Lnk1e in /gdata/lun7/software/G09E01/g09/l502.exe at Wed Mar 30 05:59:21 2022.
Job cpu time:    0 days 20 hours 49 minutes 28.5 seconds.
File lengths (MBytes):  RWF= 684 Int=    0 D2E=    0 Chk=    19 Scr=    2

原因是SCF不收敛，我看了sob老师的这篇贴子解决SCF不收敛问题的方法 - 思想家公社的门口：量子化学·分子模拟·二次元 (sobereva.com)
我的体系是从MOF晶胞中截取的一个重复片段，含N、Zn、C、H。目前想到的解决方案是：#p opt b3lyp/6-31G** em=gd3bj  SCF=vshift=500。
首先是这么处理是否合适，其次是在解决这种情况时，是在已报错的文件上加上命令继续算下去还是在原始文件上加上命令重新算呢？
感谢各位老师的解答。

作者
Author: niobium 时间: 2022-3-30 12:29
配体上少氢原子吧

作者
Author: 红米饭1234 时间: 2022-3-30 15:25
先用xtb软件里面预优化一下，读取最后一步放在Gaussian里面优化

作者
Author: fdy哈哈哈 时间: 2022-3-30 15:39

niobium 发表于 2022-3-30 12:29
配体上少氢原子吧

感谢答复，咪唑环上少了两个氢，已经改过来了。

作者
Author: fdy哈哈哈 时间: 2022-3-30 15:41

红米饭1234 发表于 2022-3-30 15:25
先用xtb软件里面预优化一下，读取最后一步放在Gaussian里面优化

感谢答复，请问xtb软件是什么嘞？之前没有接触过，可否分享一下

作者
Author: 红米饭1234 时间: 2022-3-30 15:53

fdy哈哈哈发表于 2022-3-30 15:41
感谢答复，请问xtb软件是什么嘞？之前没有接触过，可否分享一下

此论坛有，自己在主页搜

作者
Author: sobereva 时间: 2022-3-31 15:32

fdy哈哈哈发表于 2022-3-30 15:41
感谢答复，请问xtb软件是什么嘞？之前没有接触过，可否分享一下

看里面的简介
将Gaussian与Grimme的xtb程序联用搜索过渡态、产生IRC、做振动分析
http://sobereva.com/421（http://bbs.keinsci.com/thread-10106-1-1.html）
盘点Grimme迄今对理论化学的贡献
http://sobereva.com/388（http://bbs.keinsci.com/thread-6662-1-1.html）

作者
Author: jy1165254033 时间: 2024-8-23 12:11
请问电荷和自旋设置为多少呢

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