| 查看: 2032 | 回复: 4 | |||||
| 当前主题已经存档。 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
[交流]
[求助]关于opt(maxcycle=100)的设置问题。
|
|||||
|
问问， 优化过程中出现link9999错误后，用opt(maxcycle=100)，是不是非得优化100次？我设置了后，优化到47次（快一天的时间了，而且我发现从结构/键长/键角上看没有什么大的变化了），实在没有耐心了，我强行退出了，呵呵 大侠们给指点下吧？ |
回复此楼» 收录本帖的淘帖专辑推荐
 第一原理资料汇编 | 好的资料 |
» 猜你喜欢
求助火焰封管的时候管子炸了
已经有1人回复
-
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有11人回复
-
物理化学论文润色/翻译怎么收费?
已经有214人回复
-
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有10人回复
-
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
-
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？
已经有2人回复
-
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
-
半透性容器进行高湿或者低湿研究
已经有0人回复
-
Chemical Bonding at Surfaces and Interfaces,最经典的一本表面上化学相互作用教材
已经有0人回复
5楼2008-03-24 09:35:42
cuihang
铁杆木虫 (正式写手)
VASPer
- 应助: 0 (幼儿园)
- 金币: 7507.1
- 散金: 17
- 红花: 2
- 帖子: 697
- 在线: 78.5小时
- 虫号: 447551
- 注册: 2007-11-01
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2008-03-19 16:49:12
★ ★ ★ ★
lei0736(金币+4,VIP+0):呵呵 真详细 鼓励啊
lei0736(金币+4,VIP+0):呵呵 真详细 鼓励啊
|
你应该看看每步优化后的能量，如果能量持续下降，那你继续优化是能的得到优化结果，如果发现能量的变化忽大忽小，就可以选用Geom=(checkpoint, step=n)，从能量最小的那一步为起点，重新开始优化构型。 这里有一些关于优化困难的方法，希望有用： 对于一个优化计算，它的过程是先做一个SCF计算，得到这个构型下的能量，然后优化构型，再做SCF，然后再优化构型。。。因此，会有两种不收敛的情况：一是在某一步的SCF不收敛(L502错误)，或者构型优化没有找到最后结果（L9999错误）。 预备知识:计算时保存chk文件,可以在后续计算中使用guess=read读初始猜测. 对于SCF不收敛，通常有以下的解决方法： 1. 使用小基组，或低级算法计算，得到scf收敛的波函数，用guess=read读初始波函数。 2. 使用scf=qc，这个计算会慢，而且需要用stable关键字来测试结果是否波函数稳定。如果这个还不收敛，会提示L508错误。 3. 改变键长，一般是缩小一点，有时会有用。 4. 计算相同体系的其他电子态，比如相应的阴离子、阳离子体系或单重态体系，得到的收敛波函数作为初始猜测进行计算。 对于优化不收敛,即L9999错误，实际上是在规定的步数内没有完成优化，即还没有找到极小值点。（或者对于过渡态优化，还没有找到过渡态） 这有几种可能性： 1. 看一下能量的收敛的情况，可能正在单调减小，眼看有收敛的趋势，这样的情况下，只要加大循环的步数（opt(maxcycle=200)），可能就可以解决问题了。 2. 加大循环步数还不能解决的（循环步数有人说超过200再不收敛，再加也不会有用了，这虽然不一定绝对正确，但200步应该也差不多了），有两种可能。一是查看能量，发现能量在振荡了，且变化已经很小了，这时可能重新算一下，或者构型稍微变一下，继续优化，就可以得到收敛的结果（当然也有麻烦的，看运气和经验了）；二是构型变化太大，和你预计的差别过大，这很可能是你的初始构型太差了，优化不知道到哪里去了，这时最好检查一下初始构型，再从头优化。 3. 对于L9999快达到收敛时,考虑减小优化步长有时对于能量振荡的情况也是有用的，opt(maxstep=1). 一个建议是，对于大体系，难收敛体系，先用小基组，低精度算法优化一下，以得到较好的初始构型，再用高精度的计算接着算。如果前面的方法保留了chk文件，重新计算时需要使用 geom=allcheck 读入构型（就不必麻烦地写构型了）， guess=read（读入初始波函数，可以加快第一步SCF收敛）。 [ Last edited by liliracial on 2008-3-19 at 17:51 ] |
3楼2008-03-19 17:50:41
4楼2008-03-20 08:53:30