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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 过渡态
51/1返回列表
查看: 1434  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lm214019

银虫 (小有名气)

[求助] 过渡态

计算过渡态时，一直出现9999错误（超过5个），应该加个什么关键词可以避免？
回复此楼

» 猜你喜欢
1楼 2012-05-16 11:13:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fishdewow

铜虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
lm214019: 金币+5, ★★★很有帮助 2012-05-18 10:41:37
本班有帖子说过，optcyc设置为200，优化不出来的话，盲目加大循环次数意义不大，几何构型优化不收敛的解决方案。引起这个问题的原因比较多，可能也说不全，尽量做到全面。
　　　　　　(1) 如果是很小的分子（10原子以内），初始结构可能离平衡结构比较远，又在输出文件优化的最后一步判断是否收敛的位置（可以通过查找Threshold字段找到，它下面有四个判断项，都是YES才代表优化收敛）看到了“-- Number of steps exceeded,  NStep=xxx”，可以通过加大优化的循环次数来解决问题。使用关键词opt=(maxcycle=n).
　　　　　　(2) 优化到如果四个收敛标准中前两项早就收敛了，而后两项尤其是第三项不收敛，这需要判断原因。如果每一步优化过程中的能量始终在下降，那么可以继续让他算，超过了最大步数停掉了的话把结构拿出来，重新提交就行了。如果能量忽大忽小出现跳跃，说明遇到了比较平坦的势能面，或者优化的算法不好。此时建议按如下顺序处理：
　　　　　　　　　a. 使用opt=maxstep=n来缩小最大步长，即原子移动的距离，n的默认值是30，只能设置整数。
　　　　　　　　　b. 使用关键词opt=gdiis，很多时候可以很快收敛。
　　　　　　　　　c. 使用opt=calcfc关键词，在几何优化的第一步计算力常数，为优化定一个大方向，类似在第一个三岔路口选一条路。
　　　　　　　　　d. 使用opt=calcall关键词，在几何优化计算的每一步都计算力常数，类似于在每一个三岔路口都来选一条路。这个关键词不到万不得已不要使用，非常耗时。
　　　　　　(3) 四个判断标准的后两项早就收敛了，但前两项似乎固定在某个数值上了。遇到这种情况，首先去检查你优化第一步的时候给出的每个轨道得本征值部分，也就是轨道的能量，如果HOMO和LUMO的能量相同或者非常接近，那么你的设置有问题，即自旋多重度不对。如果HOMO-LUMO能量没问题，那么将目前的计算停掉，提出结构，使用opt=gdiis关键词，很多时候管用。
　　　　　　(4) 如果是在几何构型优化的过程中出现因为scf不收敛导致无法优化，可以先用上面scf不收敛的解决方案来处理，如果都不行，可以通过在#后面添加p输出更多信息，如果看到每次scf能量都能收敛到10的-6次方以下然后开始出现能量的跳跃，就是收敛不到10的-8次方（这种情况还比较常见），可以使用scf=(conver=6)来让几何优化继续下去，等优化好了，再将此关键词去掉，我的经验，优化的结构一般不会有太大的差异。如果不到10的-6能量就开始跳越，不适合此方法。
　　　　　　(5) 如果几何优化过程中发现键变得异常短，而且此时几何优化总是不收敛。遇到这种情况，问题出在你的基组使用，肯定是使用了赝势基组却没有读入赝势。由于赝势基组将内层电子和原子核用一个有效核心势来表示，你没有读入赝势就说明没有读入这个有效核心势，那么相当于把原子的内层电子都给除掉了，巨大的核引力把原子之间的距离变得异常短也不足为奇。


上述解决方案转自本版 http://muchong.com/bbs/viewthread.php?tid=1566146
赞一下 回复此楼
8楼2012-05-18 10:10:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

gujing136

禁虫 (文坛精英)

【答案】应助回帖

感谢参与，应助指数 +1
呃，虽然我没有做过过渡态，但是发帖的时候注意一下要求，把输入文件，和错误文件贴出来一些吧……
赞一下 回复此楼
叫老猫的兔子
2楼2012-05-16 14:35:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lifeycc

铁虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
将Maxcyc=20变成Maxcyc=50
赞一下 回复此楼
superman爱喝摩卡
3楼2012-05-16 14:47:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lm214019

银虫 (小有名气)
引用回帖:
2楼: Originally posted by gujing136 at 2012-05-16 14:35:01:
呃，虽然我没有做过过渡态，但是发帖的时候注意一下要求，把输入文件，和错误文件贴出来一些吧……

我的输入文件是这样的：
%chk=5Tcuzsm-5.chk
%mem=6MW
%nproc=1
#N opt=(ts,z-matrix,noeigentest,nrscale) freq b3lyp/3-21g
optcyc=300

Title Card Required

2 1
Si                 0.00000000    0.00000000    0.00000000
Si                 0.00000000    0.00000000    6.41770817
Al                 1.29801685    0.00000000    3.21107229
O                  0.33237084    0.00010657    4.73012718
O                  0.34620934   -0.04882568    1.68408882
O                  2.76575769   -1.21674610    3.16045999
Si                 2.92082130   -2.93769984    3.07576355
Si                 2.33863405    3.04541412    3.53816768
O                  2.62399645    1.35919585    3.27713846
H                  0.75407383   -1.08709958   -0.66331527
H                  0.42704873    1.30726459   -0.54689002
H                 -1.45070396   -0.19388262   -0.21969879
H                  1.00301304    0.85066069    7.09648737
H                  0.09712464   -1.38610032    6.92730870
H                 -1.35773233    0.53666734    6.66050660
H                  4.34792661   -3.31939068    2.98586063
H                  2.41839893   -3.48644048    1.79636703
H                  2.54807112   -3.59559693    4.34801696
H                  3.30496835    3.73939673    2.65783471
H                  2.57380310    3.28988594    4.97876793
H                  0.92877746    3.22178395    3.12392369
Al                 4.09876604    0.15116763    3.15956939
Cl                 6.01828139    0.19570053    4.68140303
Cl                 4.46132208    0.12056725    1.01014697
Mg                -1.43180142   -0.02717497    3.22700087
Si                 8.26631441   -0.28422264    4.12834776
C                  6.81483360   -0.25970718    2.59567709
H                  6.13745546   -1.00011244    2.96697067
H                  6.35381359    0.70439162    2.64928705
H                  7.06218171   -0.48070527    1.57838717
Cl                 9.47291120    1.07354540    3.07450236
Cl                 8.49773348   -0.03973570    6.20118913
Cl                 9.01568130   -2.18821704    3.65584951
赞一下 回复此楼
4楼2012-05-16 15:00:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭