| 查看: 529 | 回复: 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
» 猜你喜欢
0831生医工第一轮调剂失败求助
已经有17人回复
-
291求调剂
已经有11人回复
-
297，工科调剂?河南农业大学本科
已经有14人回复
-
0854求调剂
已经有20人回复
-
26药学专硕105500求调剂
已经有7人回复
-
一志愿华中农业071010，320求调剂
已经有15人回复
-
2026 WR青拔
已经有7人回复
-
一志愿中科大材料与化工，353分还有调剂学校吗
已经有9人回复
-
271求调剂
已经有37人回复
-
申博自荐
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
