请教几何优化MOF晶胞中截取的一个重复片段时SCF不收敛的问题

fdy哈哈哈 · 发表于 Post on 2022-3-30 12:13:34

请教各位老师，我在（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。
首先是这么处理是否合适，其次是在解决这种情况时，是在已报错的文件上加上命令继续算下去还是在原始文件上加上命令重新算呢？
感谢各位老师的解答。

niobium · 发表于 Post on 2022-3-30 12:29:32

配体上少氢原子吧

红米饭1234 · 发表于 Post on 2022-3-30 15:25:22

先用xtb软件里面预优化一下，读取最后一步放在Gaussian里面优化

fdy哈哈哈 · 发表于 Post on 2022-3-30 15:39:52

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

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

fdy哈哈哈 · 发表于 Post on 2022-3-30 15:41:10

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

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

红米饭1234 · 发表于 Post on 2022-3-30 15:53:57

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

此论坛有，自己在主页搜

sobereva · 发表于 Post on 2022-3-31 15:32:08

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）

jy1165254033 · 发表于 Post on 2024-8-23 12:11:28

请问电荷和自旋设置为多少呢

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

[新手求助] 请教几何优化MOF晶胞中截取的一个重复片段时SCF不收敛的问题