24小时热门版块排行榜    

查看: 5066  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

Prettyswan

木虫 (小有名气)

[交流] 【求助】结构优化中,如何选择POTIM? 已有6人参与

说明手册上说:
引用回帖:
Within all relaxation algorithms (IBRION=1,2 and 3) the parameter POTIM should be supplied in the INCAR file. For IBRION0, the forces are scaled internally before calling the minimization routine. Therefore for relaxations, POTIM has no physical meaning and serves only as a scaling factor. For many systems, the optimal POTIM is around 0.5. Because the Quasi-Newton algorithm and the damped algorithms are sensitive to the choice of this parameter, use IBRION=2, if you are not sure how large the optimal POTIM is.

这里推荐在结构优化中使用IBRION=2和POTIM=0.5.
在侯柱峰编写的《vasp软件包使用入门指南》中《§6.2结构参数的优化》中“一次成型”一节写道:
引用回帖:
SYSTEM = Mg-hex
ENCUT = 250
ISTART = 0; ICHARG = 2
ISMEAR = 1; SIGMA = 0.2
NSW = 60; IBRION = 2
ISIF = 3
POTIM = 0.2
EDIFF = 1E-6; EDIFFG = -1E-3
PREC = Accurate

那么应该如何选择POTIM的值呢?
回复此楼

» 收录本帖的淘帖专辑推荐

VASP

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zhang668(金币+3):多谢指教! 2010-12-23 18:51:30
引用回帖:
Originally posted by swx0789 at 2010-12-23 08:24:33:
我现在也在优化晶体的结构参数。看了后博士的中文入门手册,虽然讲的很详细,当有个地方还不是很明白。就他例子的关于Mg的晶格参数的优化(P25)。
我取10个点,循环完后,拟合得到最佳的a值。他后面说“在该体积 ...

“改变$i的循环值,再运行run_cell计算一次”
此时$i只附了一个值,就是前面由EOS拟合得到的平衡体积V_0计算(V_0 = 0.5*sqrt(3.0) a^3 * c/a,  c/a的值取自run_cell所设置的,即c/a的值与前一步中run_cell中的相同。这样得到了一个平衡体积下初始的a, c/a)。这一次run_cell就是固定体积(将体积固定在所拟合出来的V_0),根据run_cell所设置的初始a, c/a进行固定体积的优化。

E(a,c) -----> E(V, a)------->E(V_0,a)------> V_0下最优的a, c的值。
6楼2010-12-23 18:47:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Prettyswan 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见