| 查看: 670 | 回复: 5 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
yjcmwgk禁虫 (文坛精英)
密度泛函·小卒
|
[交流]
【调查】GaussianFAQ第三次征集答案（30金币奖励）
|
||
|
已完成 [ Last edited by yjcmwgk on 2009-10-7 at 17:12 ] |
回复此楼» 猜你喜欢
分子动力学模拟合作需求
已经有4人回复
-
中国科学院成都山地灾害与环境研究所2026年特别研究助理（博士后）招聘启事
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有52人回复
-
湖北文理学院材料与化工专硕还有调剂名额，欢迎大家推荐考生！
已经有1人回复
-
湖南大学刘巧玲课题组2026年博士研究生招生信息
已经有28人回复
-
南方科技大学表面催化反应动力学研究课题组招收2026届博士生
已经有1人回复
-
北京邮电大学物理学院老师招收博士生
已经有3人回复
-
Comsol 6.0 桶内肥料产热产湿过程模拟求助
已经有1人回复
-
微纳hub，微纳加工的宝库
已经有0人回复
-
微纳加工里的 lift-off 工艺
已经有0人回复
flash8407
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 1717.6
- 散金: 477
- 帖子: 435
- 在线: 724.1小时
- 虫号: 430876
- 注册: 2007-08-11
- 性别: GG
- 专业: 理论和计算化学
3楼2009-10-07 16:23:18
a349857460
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 2552.2
- 散金: 1066
- 红花: 13
- 帖子: 1075
- 在线: 142.3小时
- 虫号: 576761
- 注册: 2008-06-21
- 性别: GG
- 专业: 有机分子功能材料化学
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+0,VIP+0):请参考http://emuch.net/bbs/viewthread.php?tid=1566146的回答方式。您这样回答，不可能是最佳答案 10-7 16:27
yjcmwgk(金币+30,VIP+0):您二位的回答都太好了，我怎么取舍啊……我给您二位合并答案了，每人30分吧！ 10-7 16:48
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+0,VIP+0):请参考http://emuch.net/bbs/viewthread.php?tid=1566146的回答方式。您这样回答，不可能是最佳答案 10-7 16:27
yjcmwgk(金币+30,VIP+0):您二位的回答都太好了，我怎么取舍啊……我给您二位合并答案了，每人30分吧！ 10-7 16:48
|
1. 对于一个优化计算，它的过程是先做一个SCF计算，得到这个构型下的能量，然后优化构型，再做SCF，然后再优化构型。。。因此，会有两种不收敛的情况：一是在某一步的SCF不收敛(L502错误)，或者构型优化没有找到最后结果（L9999错误）。 （1）对于SCF不收敛，通常有以下的解决方法： a. 使用小基组，或低级算法计算，得到scf收敛的波函数，用guess=read读初始波函数。 b. 使用scf=qc，这个计算会慢，而且需要用stable关键字来测试结果是否波函数稳定。如果这个还不收敛，会提示L508错误。 c. 改变键长，一般是缩小一点，有时会有用。 d. 计算相同体系的其他电子态，比如相应的阴离子、阳离子体系或单重态体系，得到的收敛波函数作为初始猜测进行计算。 （2）. 对于优化不收敛,即L9999错误，实际上是在规定的步数内没有完成优化，即还没有找到极小值点。（或者对于过渡态优化，还没有找到过渡态） 这有几种可能性： a. 看一下能量的收敛的情况，可能正在单调减小，眼看有收敛的趋势，这样的情况下，只要加大循环的步数（opt(maxcycle=200)），可能就可以解决问题了。 b. 加大循环步数还不能解决的（循环步数有人说超过200再不收敛，再加也不会有用了，这虽然不一定绝对正确，但200步应该也差不多了），有两种可能。一是查看能量，发现能量在振荡了，且变化已经很小了，这时可能重新算一下，或者构型稍微变一下，继续优化，就可以得到收敛的结果当然也有麻烦的，看运气和经验了）；二是构型变化太大，和你预计的差别过大，这很可能是你的初始构型太差了，优化不知道到哪里去了，这时最好检查一下初始构型，再从头优化。 c. 对于L9999快达到收敛时,考虑减小优化步长有时对于能量振荡的情况也是有用的，opt(maxstep=1).(flyingheart )一个建议是，对于大体系，难收敛体系，先用小基组，低精度算法优化一下，以得到较好的初始构型，再用高精度的计算接着算。如果前面的方法保留了chk文件，重新计算时需要使用 geom=allcheck 读入构型（就不必麻烦地写构型了）， guess=read（读入初始波函数，可以加快第一步SCF收敛）。 引用网页：http://bbs.theochem.net/viewthread.php?tid=20 2.writwa错误突然中断 硬盘问题，两种情况： 一种是你使用的是32位的系统或Gaussian软件，那么你的计算最多只能写16GB的硬盘空间，即使你的电脑剩余空间再多，需要的空间超过了16G也要出错，如果是这样的，在命令行尝试使用maxdisk=2000mw关键词，也许会解决问题，如果不行，你只能找64位的版本来做计算了。windows下使用老版本的gaussian作计算时还会出现单个文件2G的限制，这时需要在将这16G拆分成8个左右的磁盘读写文件。在命令行的上面添加 %rwf=name1.rwf,name2.rwf,name3.rwf.......这里的文件名随便起 另一种情况是你的系统和软件是64位的，这种情况下如果你没有定义最大硬盘，gaussian就会无限的写下去，直到满足要求为止，这时经常会出现硬盘被写满的情况，也就是你的硬盘没有空间了，此时同样需要定义maxdisk关键词，容量至少要比你的硬盘剩余空间小一些。 引自：http://blog.163.com/lhwhxu/blog/static/10253754120097792051215/ [ Last edited by a349857460 on 2009-10-7 at 16:33 ] |
2楼2009-10-07 16:12:50
liuchun521
铁杆木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 5621.6
- 散金: 20
- 红花: 19
- 帖子: 1562
- 在线: 263.3小时
- 虫号: 674434
- 注册: 2008-12-14
- 专业: 理论和计算化学
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+30,VIP+0):您二位的回答都太好了，我怎么取舍啊……我给您二位合并答案了，每人30分吧！ 10-7 16:48
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+30,VIP+0):您二位的回答都太好了，我怎么取舍啊……我给您二位合并答案了，每人30分吧！ 10-7 16:48
|
1 答： 是属于构型优化没有找到最后结果，即在规定的优化步骤内没有找到最终的合适的结构。（L9999错误）。 对于优化不收敛,即L9999错误，实际上是在规定的步数内没有完成优化，即还没有找到极小值点。（或者对于过渡态优化，还没有找到过渡态），可以增加优化步骤，或是拿最后的结构 接着进行优化。 此外这有几种可能性： 1. 看一下能量的收敛的情况，可能正在单调减小，眼看有收敛的趋势，这样的情况下，只要加大循环的步数（opt(maxcycle=200)），可能就可以解决问题了。 2. 加大循环步数还不能解决的（循环步数有人说超过200再不收敛，再加也不会有用了，这虽然不一定绝对正确，但200步应该也差不多了），有两种可能。一是查看能量，发现能量在振荡了，且变化已经很小了，这时可能重新算一下，或者构型稍微变一下，继续优化，就可以得到收敛的结果（当然也有麻烦的，看运气和经验了）；二是构型变化太大，和你预计的差别过大，这很可能是你的初始构型太差了，优化不知道到哪里去了，这时最好检查一下初始构型，再从头优化。 3. 对于L9999快达到收敛时,考虑减小优化步长有时对于能量振荡的情况也是有用的，opt(maxstep=1).(flyingheart ) 一个建议是，对于大体系，难收敛体系，先用小基组，低精度算法优化一下，以得到较好的初始构型，再用高精度的计算接着算。如果前面的方法保留了chk文件，重新计算时需要使用 geom=allcheck 读入构型（就不必麻烦地写构型了）， guess=read（读入初始波函数，可以加快第一步SCF收敛）。 |
4楼2009-10-07 16:28:40
liuchun521
铁杆木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 5621.6
- 散金: 20
- 红花: 19
- 帖子: 1562
- 在线: 263.3小时
- 虫号: 674434
- 注册: 2008-12-14
- 专业: 理论和计算化学
|
（2）答 你的结构在优化到某一部的时候的SCF不收敛(L502错误) 预备知识:计算时保存chk文件,可以在后续计算中使用guess=read读初始猜测. 对于SCF不收敛，一般有一下几种的解决方法： 1. 使用小基组，或低级算法计算，得到scf收敛的波函数，用guess=read读初始波函数。 2. 使用scf=qc，这个计算会慢，而且需要用stable关键字来测试结果是否波函数稳定。如果这个还不收敛，会提示L508错误。 3. 改变结构（键角，键长等），这个很重要，需要经验的积累。 4. 计算相同体系的其他电子态，比如相应的阴离子、阳离子体系或单重态体系，得到的收敛波函数作为初始猜测进行计算。 5. 加入一些关键字。如：guess=core 等 6. 放低收敛的精度要求。 7. |
5楼2009-10-07 16:33:27
