| 查看: 2051 | 回复: 19 | ||
| 【悬赏金币】回答本帖问题，作者sky423将赠送您 50 个金币 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
sky423新虫 (著名写手)
|
[求助]
高斯不收敛，如何优化已有1人参与
|
|
|
#p b3lyp/lanl2dz opt pop=(full,npa) >>>>>>>>>> Convergence criterion not met. SCF Done: E(UB+HF-LYP) = -1714.45115901 A.U. after 129 cycles Convg = 0.8068D-05 -V/T = 2.0529 S**2 = 0.9172 KE= 1.628346096785D+03 PE=-1.275059909267D+04 EE= 5.005782636124D+03 Annihilation of the first spin contaminant: S**2 before annihilation 0.9172, after 0.7568 Convergence failure -- run terminated. Error termination via Lnk1e in d:program filesg03l502.exe at Fri Jun 18 18:59:12 2021. Job cpu time: 1 days 5 hours 55 minutes 14.0 seconds. File lengths (MBytes): RWF= 206 Int= 0 D2E= 0 Chk= 18 Scr= 1 请问高手，我这个从报错看应该是没有收敛，但是应该如何优化或者修改？救救孩子吧。。。。@jiewei |
回复此楼» 猜你喜欢
中国科学技术大学苏州高研院外籍院士工作站招聘科研助理1名
已经有0人回复
-
美国夏威夷大学生物胶体与界面实验室诚招全奖博士生和博士后
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有158人回复
-
求助共轭分子的荧光强度？
已经有0人回复
-
1,4-萘醌的苯环上引入一个溴
已经有1人回复
-
CO-SCR
已经有0人回复
-
XPS峰位移求助
已经有6人回复
-
拉曼光谱。。。。。。
已经有1人回复
-
一篇一作AM，找西安985博导推荐
已经有1人回复
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 566 (博士)
- 贵宾: 0.01
- 金币: 24942.9
- 红花: 338
- 帖子: 1972
- 在线: 1338.6小时
- 虫号: 2466443
- 注册: 2013-05-16
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
4楼2021-06-22 05:30:15
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 566 (博士)
- 贵宾: 0.01
- 金币: 24942.9
- 红花: 338
- 帖子: 1972
- 在线: 1338.6小时
- 虫号: 2466443
- 注册: 2013-05-16
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
【答案】应助回帖
感谢参与，应助指数 +1
|
导致scf不收敛的原因很多。通常默认的128步是够用的。如果不收敛。可能有其他一些原因。直接增加步数用处可能会很小。有一些比较通用的做法可以看我这个教程里第2.2节scf不收敛的应对方法。 http://muchong.com/t-14242735-1 因为没法发截图了 那就把内容复制到这里，不过排版就没有那么好了。 如下： a. 首先检查输入结构是不是合理。绝大部分情况下都是因为结构跟 “真实” 结构差的很远，导致的scf难以收敛。尽量给一个能给出的最好结构 （比如先用uff粗略优化一下，再用pm6 在优化一下。然后再用dft计算之类的） b. 检查体系的自旋多重度是否设置的合理。尤其是过渡金属体系。如果多重度设置的偏离 “真实” 情况很远，则非常难以收敛 c. 检查基组是否合理。尤其是包含重原子，是否合理的设置了混合基组之类的都会影响。具体参考本文1.1节设置混合基组和赝势的方法。这里你可以多看看 我看你用的lanl2dz如果你这个体系不是一个重金属的团簇 对整个体系都用这个赝势基组未必是好选择 d. 检查当前给的初始结构是否真实存在于 “气相” 中。因为有时候气相和溶液相结构很不一样。这就会导致初始结构以为很准其实是很差的情况。如果是溶液中的结构则需要加上溶剂模型scrf，见本文1.2 .有时候加了溶剂模型就可以帮助收敛。 e. 以上确定是没有问题的，这时候就考虑加入一些关键词帮助收敛 scf=(xqc,maxconventionalcycles=64) 这个关键词可以应用在绝大部分情形下都有效。但是绝不要一上来就用，一定要检查之前的步骤一都没问题了再去使用。g09的话就用scf=xqc就好。 g16可以写这一串。 scf=novaracc 高斯在计算scf过程中，最开始的几步为了加速会降低积分精度。这样有时候就会让波函数不小心进入了一个错误的区间。很难再优化下去。所以可以用这个关键词，取消这种加速降低积分精度的方法。 scf(vshift=500) 这个一般可以用在过渡金属体系。这种体系，可能homo-lumo间距过小导致波函数不容易收敛。可以计算中人为太高这个值。这个不会影响结果只会影响收敛过程。一般设定在500上下比较合理。 int(acc2e=12, ultrafine) 如果用g09进行计算，也可以考虑这个关键词提高了dft的积分精度和格点精度。这个仅仅适用在某些dft计算中。比如m062x 这类对于格点要求比较高的泛函中。g16中这个已经是默认值不需要写。 实际计算可以用以上提到的一个或者几个关键词一起帮助收敛。 如果体系实在是太过巨大，这些都做了还是很难收敛。可以使用 scf(conver=6) 这个关键词降低scf收敛限到10^-6， 而默认是10^-8。但是慎用。这样的结构可能得到的根本不准。最好是优化完成后，读取波函数和结构恢复到默认的收敛限再进行一次优化。 |
» 本帖已获得的红花（最新10朵）
sky4235楼2021-06-22 05:39:13
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 566 (博士)
- 贵宾: 0.01
- 金币: 24942.9
- 红花: 338
- 帖子: 1972
- 在线: 1338.6小时
- 虫号: 2466443
- 注册: 2013-05-16
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
10楼2021-06-23 01:09:12
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 566 (博士)
- 贵宾: 0.01
- 金币: 24942.9
- 红花: 338
- 帖子: 1972
- 在线: 1338.6小时
- 虫号: 2466443
- 注册: 2013-05-16
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
11楼2021-06-23 01:12:18
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 566 (博士)
- 贵宾: 0.01
- 金币: 24942.9
- 红花: 338
- 帖子: 1972
- 在线: 1338.6小时
- 虫号: 2466443
- 注册: 2013-05-16
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
14楼2021-06-25 15:36:55
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 566 (博士)
- 贵宾: 0.01
- 金币: 24942.9
- 红花: 338
- 帖子: 1972
- 在线: 1338.6小时
- 虫号: 2466443
- 注册: 2013-05-16
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
|
如果基组设置没有问题。电荷和多重度也设置的没有什么问题。结构本身原子间的距离并没有太大异常而仅仅是gaussview不显示某些键那一般不会有什么问题。 见我这个教程3.3节 http://muchong.com/t-14242735-1 3.3 化学键的显示问题。 诸如 为什么优化完成后“缺了什么什么键“，“哪个键变成了单/双/三键”，“那个键不见了”等问题。 这类问题的回答是：没有关系。gaussview怎么显示的，都没什么所谓。键显示的符合预期也好或者不符合预期也罢都不能说明计算结果本身是对还是错。 用高斯做的量化计算需要注意的是量化计算中，或者说薛定谔方程中根本没有“键”这种东西。一切都是波函数。也就是一个一个的分子轨道。优化出来的就是原子间的相对位置，因为薛定谔方程需要的只是原子的位置信息。 而gaussview在显示分子的时候会根据一个内在的键长标准（gaussview内置了一套原子的半径信息，会拿原子间的距离相应的两个原子的半径之和作比较，比如键长超过了原子半径之和的百分之多少就属于单键，之类之类的。），完全是根据原子间的距离来显示键的种类的。跟电子真正在分子中的排布一毛钱关系都没有。所以要分析原子间的成键需要进行波函数分析。算算键级之类的去判断。gaussview显示的完全就是为了好看，但是没有任何化学意义的。并且也是可以自己手动更改键的种类的。 |
16楼2021-06-28 13:59:22
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 566 (博士)
- 贵宾: 0.01
- 金币: 24942.9
- 红花: 338
- 帖子: 1972
- 在线: 1338.6小时
- 虫号: 2466443
- 注册: 2013-05-16
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
|
如我所说 原子的距离如果没有什么大问题 仅仅不显示键那 就不影响。你判断有没有问题不能看有没有键，而是看比如你觉得原子间距离就是离得太远了 不和常规。 那可能就是计算有问题。而且为了美观，高斯view的键也是可以自己手动自行添加上去的 ，它是没有任何化学含义的 发自小木虫IOS客户端 |
18楼2021-06-28 16:09:45
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 566 (博士)
- 贵宾: 0.01
- 金币: 24942.9
- 红花: 338
- 帖子: 1972
- 在线: 1338.6小时
- 虫号: 2466443
- 注册: 2013-05-16
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
20楼2021-06-28 16:50:42
5