| 查看: 526 | 回复: 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
» 猜你喜欢
291求调剂
已经有7人回复
-
300求调剂
已经有8人回复
-
求调剂推荐
已经有7人回复
-
289 分105500药学专硕求调剂(找B区学校)
已经有4人回复
-
0854求调剂
已经有13人回复
-
初试324 中药学 一志愿天中医 求调剂
已经有3人回复
-
药学求调剂
已经有14人回复
-
327求调剂
已经有27人回复
-
急需调剂
已经有5人回复
-
271求调剂
已经有33人回复
» 本主题相关商家推荐: (我也要在这里推广)
