| 查看: 2773 | 回复: 7 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
[求助]
opt=TS 循环步骤问题
|
|||||
|
命令这么写的，但是走完第一圈就9999出错了，请高人指教！ %mem=15000MB %nprocshared=8 #p opt(maxcyc=300,TS,calcfc) freq b3lyp/6-311++g(2d,2p) scf(maxcyc=300,tight) 该怎么弄呢？ |
回复此楼» 收录本帖的淘帖专辑推荐
 科研资料 | computational chemistry （gaussian） | Gaussian |
» 猜你喜欢
求助电催化剂制备，怎么在碳布表面生长PdCo合金
已经有0人回复
-
求助电催化剂制备，怎么在碳布表面生长PdCo合金
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有137人回复
-
304求调剂
已经有5人回复
-
上海华东理工大学机械与动力工程学院先进材料与制造研究组招聘公告
已经有0人回复
-
304求调剂
已经有6人回复
-
304求调剂
已经有1人回复
-
一志愿985中药学306，求调剂
已经有0人回复
-
河南一本高校化学，材料，化工专业招收研究生
已经有0人回复
-
大连海事大学船舶洁净能源研究中心2026年博士研究生招生启事
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于磁力循环泵和流量计的问题
已经有4人回复
- 请教关于循环伏安图中的峰电流的问题
已经有5人回复
- 循环伏安法测有机聚合物 得到聚机物能级 的相关问题
已经有8人回复
- 【整理】今天整理硬盘上的资料，偶然发现的好玩的小东西就直接发上来吧
已经有105人回复
- 【讨论】这个循环伏安图说明什么问题
已经有8人回复
- 【求助】关于循环伏安峰电流值问题
已经有10人回复
- 【求助】极化曲线与循环伏安的扫速问题
已经有7人回复
- 【求助】循环伏安的几个小问题
已经有20人回复
【答案】应助回帖
★ ★
亍寞(金币+2): 2011-06-23 14:56:44
gmy1990(金币+2): 非常感谢你对论坛的贡献，奖励一下！ 2011-06-23 15:34:43
亍寞(金币+2): 2011-06-23 14:56:44
gmy1990(金币+2): 非常感谢你对论坛的贡献，奖励一下！ 2011-06-23 15:34:43
 L9999的错误的引发一般是因为Optimization stopped.这个可以从输出文件中看到，这个产生的原因有很多，其中一个是-- Number of steps exceeded, NStep= N （N=当前优化的步数） 另一个是-- Wrong number of Negative eigenvalues: Desired= 1 Actual= N （N=当前Hessian矩阵的本征值数，找过渡态需要该值为1，通过添加noeigentest参数可以取消该检查，虽然有时也能找到过渡态，但一般情况下这说明过渡态的初始猜测并不是很合理）还有一个是 Optimization aborted. -- No acceptable step. 这个很明显，构象优化出错了，解决方法调构象或者更改计算方法或条件…… 还有很多其他的原因，我这就不一一列举了。所以解决l9999的错误最好的方法是在你输出文件中仔细找一下到底是什么原因导致计算终止，然后对症下药 盲目的加maxcyc=N opt=loose opt=gdiis等等是没什么意义的 |
6楼2011-06-23 14:50:55
gmy1990
荣誉版主 (著名写手)
- QC强帖: 1
- 应助: 163 (高中生)
- 贵宾: 0.782
- 金币: 13422.7
- 散金: 526
- 红花: 42
- 沙发: 4
- 帖子: 1627
- 在线: 1019.7小时
- 虫号: 745920
- 注册: 2009-04-11
- 专业: 理论和计算化学
- 管辖: 量子化学
【答案】应助回帖
亍寞(金币+2): 2011-06-23 14:51:46
|
看到好几个帖子说L999问题了，我把量化版的FAQ上面的贴上来 L9999错误是属于构型优化没有找到最后结果，即在规定的优化步骤内没有找到最终的合适的结构。（L9999错误）。对于优化不收敛,即L9999错误，实际上是在规定的步数内没有完成优化，即还没有找到极小值点。（或者对于过渡态优化，还没有找到过渡态），可以增加优化步骤，或是拿最后的结构 接着进行优化。此外这有几种可能性：1. 看一下能量的收敛的情况，可能正在单调减小，眼看有收敛的趋势，这样的情况下，只要加大循环的步数opt(maxcycle=200)，可能就可以解决问题了。2. 加大循环步数还不能解决的（循环步数有人说超过200再不收敛，再加也不会有用了，这虽然不一定绝对正确，但200步应该也差不多了），有两种可能。一是查看能量，发现能量在振荡了，且变化已经很小了，这时可能重新算一下，或者构型稍微变一下，继续优化，就可以得到收敛的结果（当然也有麻烦的，看运气和经验了）；二是构型变化太大，和你预计的差别过大，这很可能是你的初始构型太差了，优化不知道到哪里去了，这时最好检查一下初始构型，再从头优化。3. 对于L9999快达到收敛时,考虑减小优化步长有时对于能量振荡的情况也是有用的，opt(maxstep=1).一个建议是，对于大体系，难收敛体系，先用小基组，低精度算法优化一下，以得到较好的初始构型，再用高精度的计算接着算。如果前面的方法保留了chk文件，重新计算时需要使用 geom=allcheck 读入构型（就不必麻烦地写构型了）， guess=read（读入初始波函数，可以加快第一步SCF收敛）。 |
2楼2011-06-23 12:03:30
fatpig8832
铁杆木虫 (著名写手)
- 应助: 16 (小学生)
- 金币: 7420.2
- 散金: 31
- 红花: 28
- 帖子: 1438
- 在线: 610.9小时
- 虫号: 770403
- 注册: 2009-05-14
3楼2011-06-23 14:17:39
4楼2011-06-23 14:48:07
