| 查看: 1875 | 回复: 5 | ||
| 【悬赏金币】回答本帖问题，作者maxuexiang将赠送您 5 个金币 | ||
maxuexiang木虫 (正式写手)
|
[求助]
结构优化求助已有1人参与
|
|
|
在优化结构时出现能量不收敛的情况，也试了几种方法仍未解决，恳请大家帮忙看看是什么情况，不胜感激，谢谢啦！！！ 输入文件收尾贴出如下： %nprocshared=12 %mem=10GB # M06/gen opt(maxstep=1) freq pseudo=read opt 0 1 Ni -0.00706200 -0.00440700 -0.00376200 P 0.92149800 0.65951800 1.84502200 …… …… H -2.84472700 0.88597100 -2.13012300 H -2.29412300 2.37536800 -2.95207900 H -2.20024300 3.77058900 -0.72978100 H -2.68876800 2.42748600 0.35141400 H -1.23080000 3.36549500 0.71345100 C H P 0 6-31G** **** Ni 0 LANL2DZ **** Ni 0 f 1 1.0 3.130 1.0 **** Ni 0 LANL2DZ  QQ截图20190325210559.png |
回复此楼» 猜你喜欢
各位佬，合成GO（氧化石墨），哪家的石墨粉好一些啊，求告知！
已经有5人回复
-
OER的i-t稳定性测试
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有110人回复
-
PMS活化
已经有11人回复
-
k型热电偶的测温线出现两个地方短接，而只在一端有高温对测温有影响吗？
已经有1人回复
-
制备比表面积2平米每克的氧化锡载体
已经有1人回复
-
华算科技诚聘理论计算工程师（兼职，有限元方向）
已经有0人回复
lkk1213
新虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 2552.1
- 散金: 5
- 红花: 10
- 帖子: 441
- 在线: 65.3小时
- 虫号: 4861545
- 注册: 2016-07-22
- 性别: GG
- 专业: 金属有机化学
2楼2019-03-26 10:48:47
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 564 (博士)
- 贵宾: 0.01
- 金币: 24596.4
- 红花: 338
- 帖子: 1969
- 在线: 1334.1小时
- 虫号: 2466443
- 注册: 2013-05-16
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
maxuexiang: 金币+15, ★★★很有帮助, 非常感谢，我都试试，弄好了来跟帖 2019-03-26 15:18:49
小红豆: 金币+5, 回帖奖励 2019-04-05 09:02:01
感谢参与，应助指数 +1
maxuexiang: 金币+15, ★★★很有帮助, 非常感谢，我都试试，弄好了来跟帖 2019-03-26 15:18:49
小红豆: 金币+5, 回帖奖励 2019-04-05 09:02:01
|
也不知楼主试过哪些方法了。就贴出来的这个input文件来看，命令行里已经用了maxstep=1 把步长已经收缩到很小了。已经是想解决这个几何优化震荡问题。不过呢 楼主这样的写法只是让opt的第一步步长是1 之后又会恢复到默认值30，想要每一步都步长减小，要多加一个 notrust 关键词。除此之外还可以试试其它的方法帮助收敛。可以试试以下的关键词 1、 M06 泛函对格点比较敏感。如果楼主用的G09 可以加入int=ultrafine 命令。这个基本是m06系列泛函的计算 在g09中都要加入的关键词。 如果是用的G16 就不必了，因为已经是默认的了。 2. 换一种opt模式。大体系可以尝试用opt=GDIIS。综合1和2， 关键词改写成这样 opt(GDIIS,maxstep=1,notrust) int=ultrafine 3 Opt=calcall 这个可以每一步都计算Hessian ，会帮助收敛。这样之后的freq命令就需要了 calcall已经包含了。 这个关键词这样会大幅减慢计算速度 要是体系很大 会时间比较久。 考虑把基组改成 6-31G*， 这个多加第二个*，也就是 对H加极化意义不大的 毕竟就是个DZ基组。如果用这个 就不需要再用maxstep了 因为这样已经回每一步都比较精准了 减小步长反而不好。 4. 如果上面都做了 还不行 就尝试降低收敛标准了。不过这个有风险 有可能优化完有虚频什么的。关键词是 opt=loose |
3楼2019-03-26 11:24:12
wzhsun
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 281
- 散金: 50
- 帖子: 140
- 在线: 20小时
- 虫号: 13171742
- 注册: 2018-12-19
- 性别: GG
- 专业: 有机合成
4楼2019-04-16 14:16:56
maxuexiang
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2960.9
- 散金: 313
- 红花: 4
- 帖子: 522
- 在线: 155.3小时
- 虫号: 3011865
- 注册: 2014-03-03
- 专业: 理论和计算化学
5楼2020-05-06 22:55:59
paramecium86
版主 (著名写手)
  |
6楼2020-05-07 01:26:39
10