| 查看: 1896 | 回复: 2 | |||||||||||
honey宝贝金虫 (初入文坛)
|
[交流]
【转载】Gaussian入门技巧三（关于收敛问题） 已有2人参与
|
|
关于收敛问题 (L502, L508, L9999) 对于一个优化计算，它的过程是先做一个SCF计算，得到这个构型下的能量，然后优化构型，再做SCF，然后再优化构型。。。因此，会有两种不收敛的情况：一是在某一步的SCF不收敛(L502错误)，或者构型优化没有找到最后结果（L9999错误）。 预备知识:计算时保存chk文件,可以在后续计算中使用guess=read读初始猜测. 对于SCF不收敛，通常有以下的解决方法： 1. 使用小基组，或低级算法计算，得到scf收敛的波函数，用guess=read读初始波函数。 2. 使用scf=qc，这个计算会慢，而且需要用stable关键字来测试结果是否波函数稳定。如果这个还不收敛，会提示L508错误。 3. 改变键长，一般是缩小一点，有时会有用。 4. 计算相同体系的其他电子态，比如相应的阴离子、阳离子体系或单重态体系，得到的收敛波函数作为初始猜测进行计算。 5. 待补充. 对于优化不收敛,即L9999错误，实际上是在规定的步数内没有完成优化，即还没有找到极小值点。（或者对于过渡态优化，还没有找到过渡态） 这有几种可能性： 1. 看一下能量的收敛的情况，可能正在单调减小，眼看有收敛的趋势，这样的情况下，只要加大循环的步数（opt(maxcycle=200)），可能就可以解决问题了。 2. 加大循环步数还不能解决的（循环步数有人说超过200再不收敛，再加也不会有用了，这虽然不一定绝对正确，但200步应该也差不多了），有两种可能。一是查看能量，发现能量在振荡了，且变化已经很小了，这时可能重新算一下，或者构型稍微变一下，继续优化，就可以得到收敛的结果（当然也有麻烦的，看运气和经验了）；二是构型变化太大，和你预计的差别过大，这很可能是你的初始构型太差了，优化不知道到哪里去了，这时最好检查一下初始构型，再从头优化。 3. 对于L9999快达到收敛时,考虑减小优化步长有时对于能量振荡的情况也是有用的，opt(maxstep=1).(flyingheart ) 一个建议是，对于大体系，难收敛体系，先用小基组，低精度算法优化一下，以得到较好的初始构型，再用高精度的计算接着算。如果前面的方法保留了chk文件，重新计算时需要使用 geom=allcheck 读入构型（就不必麻烦地写构型了）， guess=read（读入初始波函数，可以加快第一步SCF收敛）。 |
回复此楼» 收录本帖的淘帖专辑推荐
 Materials-studio | 量子化学资料 | 高斯/MOPAC | 计算化学自用 |
| 量子化学 | Gaussian | 高斯计算 | Photonic Crystal |
» 猜你喜欢
美国杰克逊州立大学周剑歌教授的化学计算与AI实验室诚聘博士后研究员
已经有4人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有261人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
Tc数据拟合
已经有1人回复
-
寻求物理化学学报最新的投稿模板！
已经有4人回复
-
韩国釜山大学储能电池与功能材料实验室诚招博后（电池/电化学/器件制备/MOF合成）
已经有187人回复
-
祈福！
已经有79人回复
-
投稿咨询
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Gaussian为何对cpu和内存使用率很低？
已经有13人回复
- 高斯03软件基础运用（gaussian03）
已经有7人回复
- 【活动征集帖子】征集"从入门到精通”的DFT软件的入门帖子
已经有77人回复
- 求一个可以在i686，也就是ia32的CUP上使用的gaussian09
已经有10人回复
- 【转载】Gaussian入门技巧六（关于虚频）
已经有14人回复
- 【转载】Gaussian入门技巧四（关于对称性）
已经有7人回复
- 【转载】Gaussian入门技巧二（开壳层和闭壳层）
已经有12人回复
- 【转载】Gaussian入门技巧一（自旋多重度）
已经有8人回复
- 【其他】扫描了一份Gaussian 09使用手册（但进来以后，你会很失望）
已经有34人回复
- 【分享】Gaussian 09 使用手冊（中文繁体）
已经有32人回复
hx7111
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 1488.9
- 帖子: 310
- 在线: 226.6小时
- 虫号: 1038865
- 注册: 2010-06-09
- 性别: GG
- 专业: 化学工程及工业化学
2楼2011-04-17 15:53:42
xinglaidor
金虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 1301.7
- 散金: 366
- 红花: 4
- 帖子: 423
- 在线: 177.1小时
- 虫号: 938587
- 注册: 2010-01-07
- 性别: MM
- 专业: 理论和计算化学
3楼2011-04-17 16:21:55
