| 查看: 1503 | 回复: 7 | |||
lm214019实习版主
|
[求助]
过渡态
|
| 计算过渡态时，一直出现9999错误（超过5个），应该加个什么关键词可以避免？ |
回复此楼» 猜你喜欢
深圳信息职业技术学院2026年教师招聘公告（最新）
已经有0人回复
-
诚邀有志之士加入江苏大学环境学院合成生物学课题组！(长期接收学硕专硕、硕博连读)
已经有16人回复
-
物理化学论文润色/翻译怎么收费?
已经有71人回复
-
研究生招生
已经有0人回复
-
一志愿南大理论与计算化学方向（070300）342分，求调剂
已经有1人回复
-
欢迎化学、化工、物理、材料、或计算机背景的考研咨询
已经有0人回复
-
山西大学资源与环境工程研究所化工分离/环境功能材料小组招收2026年博士生
已经有3人回复
-
一志愿南师大0703化学 275求调剂
已经有5人回复
-
A 区过线就可申！重庆双非公办高校 0856调剂，接受跨专业！纳米材料 / 润滑热门方向
已经有0人回复
-
欢迎调剂长江师范学院汤老师课题组的硕士研究生
已经有0人回复
-
调剂招生
已经有3人回复
2楼2012-05-16 14:35:01
lifeycc
实习版主
- 应助: 6 (幼儿园)
- 金币: 519.9
- 散金: 1430
- 红花: 2
- 帖子: 111
- 在线: 57小时
- 虫号: 1145832
- 注册: 2010-11-13
- 性别: GG
- 专业: 理论和计算化学
3楼2012-05-16 14:47:44
lm214019
专家顾问
- 应助: 0 (幼儿园)
- 金币: 389
- 散金: 20
- 帖子: 143
- 在线: 74.5小时
- 虫号: 1155115
- 注册: 2010-11-24
- 性别: MM
- 专业: 理论和计算化学
|
我的输入文件是这样的： %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
lm214019
超级版主
- 应助: 0 (幼儿园)
- 金币: 389
- 散金: 20
- 帖子: 143
- 在线: 74.5小时
- 虫号: 1155115
- 注册: 2010-11-24
- 性别: MM
- 专业: 理论和计算化学
|
出错部分： Sum of Mulliken charges= 2.00000 Electronic spatial extent (au): 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
6楼2012-05-17 08:42:16
meteoric30
兑换贵宾
- 应助: 43 (小学生)
- 金币: 1719.5
- 散金: 12
- 红花: 2
- 帖子: 436
- 在线: 320.7小时
- 虫号: 1046846
- 注册: 2010-06-24
- 性别: GG
- 专业: 催化化学
7楼2012-05-18 09:33:45
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
lm214019: 金币+5, ★★★很有帮助 2012-05-18 10:41:37
感谢参与，应助指数 +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
130