| 查看: 1356 | 回复: 6 | |||
ww1987铁杆木虫 (小有名气)
|
[交流]
【求助】计算错误 帮忙看看 已有5人参与
|
|
我计算C80（D5d）构型的分子 计算方法用B3LYP/6-31G(d)水平 算到一半出现这种现象 请帮忙！ >>>>>>>>>> Convergence criterion not met. SCF Done: E(RB+HF-LYP) = -3048.38575494 A.U. after 129 cycles Convg = 0.7361D-06 -V/T = 2.0094 S**2 = 0.0000 KE= 3.020073824779D+03 PE=-3.325628686430D+04 EE= 1.408728660716D+04 Convergence failure -- run terminated. Error termination via Lnk1e in C:\G03W\l502.exe at Sat Jul 03 19:20:23 2010. Job cpu time: 0 days 10 hours 15 minutes 11.0 seconds. File lengths (MBytes): RWF= 897 Int= 0 D2E= 0 Chk= 98 Scr= 1 |
回复此楼» 猜你喜欢
中南大学化学化工学院招收电池储能、理论计算方向博士生
已经有16人回复
-
请问哪里可以代烧带压力的炉子？
已经有6人回复
-
物理化学论文润色/翻译怎么收费?
已经有85人回复
-
关于合成MIL-125的求助
已经有0人回复
-
双非本硕光催化方向想读博，请各位老师看看我！
已经有2人回复
-
瑞典林雪平大学范立洲课题组现招聘博士后、岗位制博士生
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
寻求物理化学学报最新的投稿模板！
已经有4人回复
-
投稿咨询
已经有7人回复
-
自科面上 散金币 祈福
已经有61人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 高斯irc计算出错，请高手帮忙看一下怎么解决
已经有8人回复
- 得不到纯产物，哪出错了，亲们帮忙看看啊
已经有3人回复
- 求两个B级错误的处理方法,非常感谢！
已经有6人回复
- vasp自洽计算K点不匹配问题，急！！
已经有15人回复
- G09集群计算结果的chk文件转换问题
已经有7人回复
- 一般英文稿写好之后 怎么改语法和一些语句的错误 大家一般都是怎么走这个过程的
已经有18人回复
- 1000金币悬赏，做过回收率和线性范围的同学进
已经有32人回复
- 关于Dictra的计算问题（主要是关于计算过程出错的问题）
已经有8人回复
- 【求助完成】提交vasp时出现以下错误，请大家帮忙看下
已经有5人回复
- 【求助】含时密度泛函计算金团簇的激发态输出出错
已经有8人回复
sunhaitao
木虫 (小有名气)
- QC强帖: 1
- 应助: 12 (小学生)
- 金币: 3138.4
- 散金: 11
- 红花: 2
- 帖子: 281
- 在线: 416.9小时
- 虫号: 718371
- 注册: 2009-03-09
- 性别: GG
- 专业: 理论和计算化学
2楼2010-07-05 15:15:55
3楼2010-07-05 15:43:43
ter20
铁杆木虫 (正式写手)
- QC强帖: 1
- 应助: 28 (小学生)
- 金币: 6026.5
- 散金: 99
- 红花: 6
- 帖子: 426
- 在线: 119.8小时
- 虫号: 541263
- 注册: 2008-04-07
- 性别: GG
- 专业: 理论和计算化学
4楼2010-07-05 17:21:54
ter20
铁杆木虫 (正式写手)
- QC强帖: 1
- 应助: 28 (小学生)
- 金币: 6026.5
- 散金: 99
- 红花: 6
- 帖子: 426
- 在线: 119.8小时
- 虫号: 541263
- 注册: 2008-04-07
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
heyo_123(金币+2):鼓励交流！！ 2010-07-05 17:34:41
小木虫(金币+0.5):给个红包，谢谢回帖交流
heyo_123(金币+2):鼓励交流！！ 2010-07-05 17:34:41
|
网上的建议： 1． 检查是否有初始文件错误 在命令行中加入 %kJob L301 or %kJob L302 如果通过则一般初始文件ok。 常见初级错误： a. 自旋多重度错误 b. 变量赋值为整数 c. 变量没有赋值或多重赋值 d. 键角小于等于0度，大于等于180度 e. 分子描述后面没有空行 f. 二面角判断错误，造成两个原子距离过近 g. 分子描述一行内两次参考同一原子,或参考原子共线 2． SCF(自洽场)不收敛 则一般是L502错误 省却情况做64个cycle迭代(G03缺省128 cycles) a. 修改坐标，使之合理 b. 改变初始猜 Guess=Huckel 或其他的，看Guess关键词。 c. 增加叠代次数 SCFCYC=N （对小分子作计算时最好不要增加，很可能结构不合理） d. iop（5/13=1）这样忽略不收敛，继续往下做。 3． 分子对称性改变 a. 修改坐标，强制高对称性或放松对称性 b. 给出精确的、对称性确定的角度和二面角。 如CH4的角度给到109.47122 c. 放松对称性判据 Symm=loose d. 不做对称性检查 iop(2/16=1) （最好加这个选项） iop(2/16=2) 则保持新的对称性来计算 4． Opt时收敛的问题 a. 修改坐标，使之合理 b. 增加叠代次数optcyc=N 5． 优化过渡态，若势能面太平缓，则不好找到。 iop(1/8=10) 默认30（下一个结构和该结构的差别0.3Å），可改成10。如果每一步都要用到小的步长，应该加opt(notrustupdate) 6． 在CI（组态）方法中 如QCISD(T),CCSD(T),CID方法中，省却最大循环50，若出错（L913错误） 解决方法：#P QCISD(maxcyc=N) 注：N≤512 7． 优化过渡态 opt=TS (给出过渡态) opt=qst2 (给出反应物和产物) opt=qst3 (给出反应物和产物和过渡态) a. 用G03时的出错 opt=ts 必须加FC (force constant) 写法：opt=(TS, calcFc) or opt=(TS,calchffc)计算HF力常数，对QCISD，CCSD等方法用； or opt=(TS,modRedundant) (最好写这个) b. 如果计算采用QCISD计算（不好计算FC） 则写为QCISD opt=(TS, calcHFFC) (用HF计算FC) 8. 无法写大的Scratch文件RWF a. 劈裂RWF文件 %rwf=loc1,size1,loc2,size2,……..,locN,-1 b. 改变计算方法 MP2=Direct可以少占硬盘空间 c. 限制最大硬盘 maxdisk=N GB，****MB,有些系统写2GB会出错，可以写2000MB 9. FOPT出错 原因是变量数与分子自由度数不相等。 可用POPT 或直接用OPT 10. 优化过渡态只能做一个STEP 原因是负本征数目不对 添加 iop（1/11）=1或者noeigentest |
5楼2010-07-05 17:24:04
6楼2010-07-05 18:35:05
beyondboyy
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 82
- 散金: 24
- 红花: 2
- 帖子: 361
- 在线: 67.2小时
- 虫号: 713701
- 注册: 2009-03-03
- 性别: GG
- 专业: 功能陶瓷
7楼2010-07-13 19:08:56
20