应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【调查】GaussianFAQ第三次征集答案（30金币奖励）
61/1返回列表
查看: 528  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

[交流] 【调查】GaussianFAQ第三次征集答案（30金币奖励）

已完成

[ Last edited by yjcmwgk on 2009-10-7 at 17:12 ]
回复此楼

» 猜你喜欢
1楼 2009-10-07 16:05:35
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

a349857460

木虫 (著名写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+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 ]
赞一下(21人) 回复此楼
2楼2009-10-07 16:12:50
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

flash8407

木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+0,VIP+0):请参考http://emuch.net/bbs/viewthread.php?tid=1566146的回答方式。您这样回答，不可能是最佳答案 10-7 16:27
L9999应该是步数不够，把结构保存再继续算下去就可以了
赞一下(11人) 回复此楼
3楼2009-10-07 16:23:18
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

liuchun521

铁杆木虫 (著名写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+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收敛）。
赞一下(11人) 回复此楼
4楼2009-10-07 16:28:40
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

liuchun521

铁杆木虫 (著名写手)
（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
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

liuchun521

铁杆木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
3  答
  是内存问题  可以增加内存  或是调整机器的内存。
  或是硬盘问题，两种情况：
一种是你使用的是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/
赞一下(1人) 回复此楼
6楼2009-10-07 16:45:31
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yjcmwgk 的主题更新
61/1返回列表
信息提示
关闭
请填处理意见
关闭