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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 过渡态
81/1返回列表
查看: 1433  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

lm214019

银虫 (小有名气)

[求助] 过渡态

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

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

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的回帖

lm214019

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

出错部分：
Sum of Mulliken charges=   2.00000
Electronic spatial extent (au):  = 12022.6191
Charge=     2.0000 electrons
Dipole moment (field-independent basis, Debye):
    X=     8.7970    Y=     5.8364    Z=    -1.2550  Tot=    10.6314
Quadrupole moment (field-independent basis, Debye-Ang):
   XX=  -101.3779   YY=  -152.1076   ZZ=  -129.0925
   XY=    -2.8804   XZ=   -40.2458   YZ=     1.6128
Traceless Quadrupole moment (field-independent basis, Debye-Ang):
   XX=    26.1481   YY=   -24.5816   ZZ=    -1.5665
   XY=    -2.8804   XZ=   -40.2458   YZ=     1.6128
Octapole moment (field-independent basis, Debye-Ang**2):
  XXX=    11.3080  YYY=    60.5999  ZZZ=   -58.2902  XYY=    72.0253
  XXY=    39.9780  XXZ=   110.9472  XZZ=    74.3849  YZZ=    21.3764
  YYZ=   -39.7445  XYZ=   -15.4819
Hexadecapole moment (field-independent basis, Debye-Ang**3):
XXXX= -9084.0365 YYYY= -2720.6517 ZZZZ= -2253.0993 XXXY=  -192.6681
XXXZ=  -447.6416 YYYX=    80.3685 YYYZ=     0.6904 ZZZX=  -245.0311
ZZZY=   -26.2051 XXYY= -1943.6951 XXZZ= -1936.7095 YYZZ=  -889.3074
XXYZ=   -51.2733 YYXZ=  -120.6486 ZZXY=    14.7472
N-N= 3.193312440948D+03 E-N=-1.564641506277D+04  KE= 3.902425821285D+03


THERE IS NO SUBJECT, HOWEVER COMPLEX,
WHICH, IF STUDIED WITH PATIENCE AND INTELLIGIENCE
WILL NOT BECOME
MORE COMPLEX
                        QUOTED BY D. GORDON ROHMAN
Error termination request processed by link 9999.
Error termination via Lnk1e in /home/g03/l9999.exe at Tue May 15 07:16:07 2012.
Job cpu time:  0 days 10 hours 40 minutes 22.0 seconds.
File lengths (MBytes):  RWF=     55 Int=      0 D2E=      0 Chk=     21 Scr=      1
赞一下 回复此楼
5楼2012-05-16 15:05:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baoa9992

木虫 (初入文坛)

【答案】应助回帖

感谢参与，应助指数 +1
计算的步数不够，用最后的结构接着计算
赞一下 回复此楼
开启教师新路程
6楼2012-05-17 08:42:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

木虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
引用回帖:
5楼: Originally posted by lm214019 at 2012-05-16 15:05:01:
出错部分：
Sum of Mulliken charges=   2.00000
Electronic spatial extent (au):  <R**2>= 12022.6191
Charge=     2.0000 electrons
Dipole moment (field-independent basis, Debye):
    X= ...

opt=（maxcycle=1000）
赞一下 回复此楼
7楼2012-05-18 09:33:45
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 lm214019 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭