| 查看: 502 | 回复: 0 | ||
[求助]
几何优化计算中 nstep 设定失效，是个Bug？
|
|
昨天刚刚在virtualbox 虚拟机（ubuntu 10.04）上并行编译好QE ，今天开始跑下简单的几何优化任务（就只有一个CO2和一个H2O 分子，算下最小能量构型），问题出现了： 当我提高能量和力收敛阈值的时候，例如： etot_conv_thr=1.0D-10, forc_conv_thr=1.0D-9, 计算结果里显示： The maximum number of steps has been reached. End of BFGS Geometry Optimization 这个不难理解，就是更难收敛以后，需要的步数增多，超过了 nstep的值，但是非常奇怪，我仔细查看了out文件后发现， 一共仅仅算了10步就结束了。用xcrysden看结果文件也是， 只有10帧。后来我把nstep的值设定为500了，结果没有任何 变化，依然是仅仅算10步。 in文件和out 文件都在附件里，希望大家能帮忙看看是什么 原因导致了仅仅计算10步，还提示我超过最大步数限制。 这会是PWscf的bug吗？In 文件内容如下： &control calculation = "relax", prefix = 'Si_paw2Long4', etot_conv_thr=1.0D-10, forc_conv_thr=1.0D-9, nstep=500 / &system ibrav = 1, celldm(1) = 47.2431497, nat = 6, ntyp = 3, ecutwfc = 20 / &electrons conv_thr = 1.D-9, mixing_beta = 0.7, electron_maxstep = 1000 / &IONS pot_extrapolation = "second_order", wfc_extrapolation = "second_order", / ATOMIC_SPECIES H 1.00784 H.pbe-kjpaw_psl.0.1.UPF C 12.0096 C.pbe-n-kjpaw_psl.0.1.UPF O 15.99903 O.pbe-n-kjpaw_psl.0.1.UPF ATOMIC_POSITIONS {crystal} H 0.4687179890644891 0.4999999999999999 0.3488926097145480 H -0.4687179890644891 -0.4999999999999999 0.3488926097145479 C 0.5000000000000000 -0.5000000000000000 0.5442073777314737 O 0.4503999948501588 -0.5000000000000001 0.5442073777314737 O -0.4503999948501588 0.5000000000000001 0.5442073777314737 O 0.5000000000000000 0.5000000000000000 0.3732574409310301 K_POINTS (gamma) 我知道我的etot_conv_thr和 forc_conv_thr设置的特别高，但是这个不能成为仅仅算了几步就 说超过最大步数限制的理由。希望大家看看问题在哪里？ [ Last edited by xiaohang07 on 2013-6-21 at 14:00 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : test.in
- 附件 2 : test.out
2013-06-19 21:29:09, 1.02 K
2013-06-19 21:29:22, 63.35 K
» 猜你喜欢
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有5人回复
-
博士读完未来一定会好吗
已经有29人回复
-
博士申请都是内定的吗？
已经有5人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有12人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
