24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

爱惜花草

木虫 (著名写手)

应助: 7 (幼儿园)
金币: 3380.2
散金: 762
红花: 9
沙发: 1
帖子: 1173
在线: 196.5小时
虫号: 805026
注册: 2009-07-07
专业: 半导体材料

[交流] 【讨论】我的这个任务结构优化时间过长是不是体系设置有问题

我INCAR里面设置的NSW=600（当然一般运动不了这么多步数）
做的体系有点大，大概一个cell里有近60个离子
开始做结构优化，现在已经跑了200多圈了还没有达到平衡位置
不过输出文件里仔细检查了能量值，发现总体趋势是在下降的，只是变化不大，大概在0.0001这个数量级左右，也没有出现那种振荡的现象
这种情况是否应该将任务kill？初始结构给的有问题？

PS：精度设置也不是很高
EDIFF = 1E-3
EDIFFG = -1E-2

IBRION=2
POTIM=0.2

[ Last edited by 爱惜花草 on 2009-10-16 at 10:12 ]

回复此楼

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

VASP

» 猜你喜欢

博士读完未来一定会好吗已经有21人回复
导师想让我从独立一作变成了共一第一已经有5人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有11人回复
读博已经有4人回复
JMPT 期刊投稿流程已经有4人回复
心脉受损已经有5人回复
Springer期刊投稿求助已经有4人回复
小论文投稿已经有3人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有9人回复
申请2026年博士已经有6人回复

一步一步去实现自己的目标，不要让自己成为自己的绊脚石，克制自己！

1楼2009-10-16 08:14:54

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

nkleof

银虫 (正式写手)

应助: 0 (幼儿园)
金币: 2093.2
散金: 649
红花: 5
帖子: 813
在线: 165.9小时
虫号: 152581
注册: 2006-01-02
专业: 凝聚态物性 II ：电子结构

★ ★ ★ ★
aylayl08(金币+3,VIP+0):感谢你的解答，欢迎常来 10-16 17:35
爱惜花草(金币+1,VIP+0): 10-17 10:14

引用回帖:

Originally posted by 爱惜花草 at 2009-10-16 10:11:

设置的IBRION=2，
看了2楼和3楼的回复，POTIM的概念把我弄糊涂了

设置IBRION=2的时候，可以不用设置POTIM，用vasp的默认值，它会在优化结束后的任务报告中给出一个合适的trialstep（在任务报告的最末），下一次用IBRION=1优化的时候就可以采用这个trialstep*0.5作为POTIM的值。在vasp的说明中有介绍。
http://cms.mpi.univie.ac.at/vasp/vasp/node115.html

Within all relaxation algorithms ( IBRION=1,2 and 3) the parameter POTIM should be supplied in the INCAR file. For IBRION$>$0, 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.

In this case, the OUTCAR file and stdout will contain a line indicating a reliable POTIM. For IBRION=2, the following lines will be written to stdout after each corrector step (usually each odd step):

trial: gam= .00000 g(F)= .152E+01 g(S)= .000E+00 ort = .000E+00
(trialstep = .82)

The quantity gam is the conjugation parameter to the previous step, g(F) and g(S) are the norm of the force respectively the norm of the stress tensor. The quantity ort is an indicator whether this search direction is orthogonal to the last search direction (for an optimal step this quantity should be much smaller than ( g(F) + g(S)). The quantity trialstep is the size of the current trialstep. This value is the average step size leading to a line minimization in the previous ionic step. An optimal POTIM can be determined, by multiplying the current POTIM with the quantity trialstep.

[ Last edited by nkleof on 2009-10-16 at 10:33 ]