应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】VASP结构弛豫不收敛的问题
111/2返回列表12下一页
查看: 2109  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

stractor

金虫 (著名写手)

[交流] 【求助】VASP结构弛豫不收敛的问题

最近遇到一个体系，离子弛豫了200步还不收敛。我的INCAR是这样的：
SYSTEM=……
ISTART=0
ICHARG=2
ISMEAR=0
SIGMA=0.05
ENCUT = 340 eV
ISIF= 2
NSW=200
LREAL = Auto
GGA = 91
VOSKOWN = 1
MAXMIX = 45
NELMIN = 4
采用这个INCAR200步后没有收敛，但如果我降低离子弛豫收敛的标准，只增加一项EDIFFG=0.01，则到67步就收敛了。我觉得降低收敛标准不妥，试问有没有其它解决的办法？
    此外，VASP里面默认离子弛豫收敛标准时按能量判断的，但文献里大多采用力判断，我的问题是某个力判断标准是否和能量收敛标准对应，因为我在文献中看到了这样的表达：The atomic positions were relaxed until the total energy difference was converged within 0.001 eV, which results in having the residual forces below 0.03 eV/A.

[ Last edited by stractor on 2009-5-14 at 20:51 ]
回复此楼

» 猜你喜欢
1楼 2009-05-14 20:02:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leigao

金虫 (小有名气)
★ ★
stractor(金币+1,VIP+0):模型没什么问题，就是在超胞里面加了一个空位 5-14 22:10
wuchenwf(金币+1,VIP+0):谢谢 5-14 23:03
是不是模型建的有问题 超过50步收敛就有可能出问题了
赞一下(12人) 回复此楼
2楼2009-05-14 21:39:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sars518

木虫 (著名写手)

小木虫灌水大队队长
★ ★ ★
spur(金币+2,VIP+0):谢谢帮助！ 5-14 21:58
stractor(金币+1,VIP+0):模型没什么问题，就是在超胞里面加了一个空位。 5-14 22:12
呵呵，我的算了 300步 才收敛  可能是建的结构 不合理

不过，我的是默认精度

EDIFFG = -0.01  力判据  貌似就够了

[ Last edited by sars518 on 2009-5-14 at 21:50 ]
赞一下(12人) 回复此楼
人生是一场旅行，在欣赏风景的同时，别忘了你也是别人的风景。
3楼2009-05-14 21:46:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spur

金虫 (正式写手)
★ ★
stractor(金币+1,VIP+0):你一般用什么标准判断呢？如果我力的标准设置为-0.03觉得粗糙吗 5-14 22:12
wuchenwf(金币+1,VIP+0):谢谢 5-14 23:03
EDIFFG为负值则表示以受力为标准，如果为正值则表示以能量为标准。
赞一下(12人) 回复此楼
4楼2009-05-14 21:59:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangyc

木虫 (正式写手)
★ ★ ★
wuchenwf(金币+1,VIP+0):谢谢 5-14 23:02
stractor(金币+2,VIP+0):很好，谢谢你的解答 5-15 15:34
文献中的意思是文献采用的是能量收敛条件，不过结束后检查作用在各个原子上的力均小于0.03了。

个人认为0.03差不多吧，不算粗，力收敛太费时了
赞一下(11人) 回复此楼
5楼2009-05-14 22:23:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuxinxing1

铁杆木虫 (著名写手)
★
wuchenwf(金币+1,VIP+0):谢谢 5-14 23:02
接近 能量最小点附近，CG 算法有时 弛豫得比较慢.

试试 IBRION=1，准牛顿弛豫.
赞一下(10人) 回复此楼
6楼2009-05-14 22:31:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spur

金虫 (正式写手)
★ ★ ★
stractor(金币+2,VIP+0):我试了，把力的收敛标准设置为-0.05eV/A，一百步以内就完成了离子弛豫 5-15 15:36
stractor(金币+1,VIP+0): 5-16 15:23
如果体系大，原子数量多，则受力相对放宽-0.03或-0.05也可；如果原子数量少，体系小，则受力相对收紧，比如-0.01
赞一下(1人) 回复此楼
7楼2009-05-14 23:05:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
引用回帖:
Originally posted by spur at 2009-5-14 23:05:
如果体系大，原子数量多，则受力相对放宽-0.03或-0.05也可；如果原子数量少，体系小，则受力相对收紧，比如-0.01

我的体系80个原子。我试了，把力的收敛标准设置为-0.05eV/A，一百步以内就完成了离子弛豫。我也查到一篇文献使用的力的收敛标准就是-0.05eV/A.谢谢wuxinxing1  和 spur的回答。

[ Last edited by stractor on 2009-5-15 at 15:39 ]
赞一下 回复此楼
8楼2009-05-15 15:38:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hookah

木虫 (正式写手)
★ ★
zdhlover(金币+1,VIP+0):谢谢，欢迎常来 5-15 00:02
stractor(金币+1,VIP+0):谢谢 5-16 15:22
觉得采用VASP的默认标准还是比较报考的！否则还是参考文献的为准吧！
赞一下(11人) 回复此楼
9楼2009-05-15 16:31:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hlshi

金虫 (正式写手)
★ ★ ★ ★ ★ ★
fegg7502(金币+5,VIP+0):thank you very much! 5-17 17:29
stractor(金币+1,VIP+0):谢谢 5-20 22:00
用 IBRION=2 应该会快一些
引用回帖:
Originally posted by stractor at 2009-5-14 20:02:
最近遇到一个体系，离子弛豫了200步还不收敛。我的INCAR是这样的：
SYSTEM=……
ISTART=0
ICHARG=2
ISMEAR=0
SIGMA=0.05
ENCUT = 340 eV
ISIF= 2
NSW=200
LREAL = Auto
GGA = 91
VOSKOWN = 1
MAXMIX  ...

赞一下(10人) 回复此楼
10楼2009-05-17 09:59:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 stractor 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭