24小时热门版块排行榜    

查看: 1121  |  回复: 2
【悬赏金币】回答本帖问题,作者呆瓜呀呀将赠送您 100 个金币

呆瓜呀呀

新虫 (小有名气)

[求助] [VASP]请教一下结构优化能量一直下降但力不收敛INCAR是否有问题 已有1人参与

老师们好,我最近在计算两层石墨烯转角体系,目前计算的是28个原子,但是已经计算了一周,能量一直在下降,但是力一直不收敛。初始结构是根据文献构建的。用的16个核。因为最大离子步设100,所以这是接着算的第三次了。一个电子步大概是两百多秒。是INCAR文件有问题吗?因为原子体系不是很大。INCAR每一个参数都是对照着别的帖子设的。

通过网盘分享的文件:INCAR等4个文件
链接: https://pan.baidu.com/s/1rwqodS8ED1eBJq33-MShFA?pwd=ex2q 提取码: ex2q
发自小木虫手机客户端
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

rlafite

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
你可以先参考下面几个关键参数,确认你是否设置得比较合理:
1. 与力收敛相关的参数
EDIFF = 1E-6     # 能量收敛准则
EDIFFG = -0.01   # 力收敛准则(更松可设为-0.05)
建议放宽 EDIFFG,先达到一个较松弛的结构再细化。
2. 离子优化设置
IBRION = 2       # 共轭梯度算法
NSW = 100        # 离子步数
如果在弛豫后半段力下降缓慢,可以改用 IBRION = 1 或 IBRION = 3(准牛顿或BFGS)
3. 与离子弛豫质量有关:
POTIM = 0.5      # 原子的步长(可设小些,如 0.3)
若振荡严重,可减小 POTIM,更稳定。
4. vdw 相互作用(扭转双层石墨烯尤其重要)
IVDW = 12   # 使用Grimme DFT-D3校正
或者使用 optB88 或 vdW-DF 等非局域泛函(但速度会更慢)。

5. 精度控制
PREC = Accurate
ENCUT = 500      # 确保能量截止足够大(可根据 POTCAR 建议)
2楼2025-07-07 23:05:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

呆瓜呀呀

新虫 (小有名气)

引用回帖:
2楼: Originally posted by rlafite at 2025-07-07 23:05:01
你可以先参考下面几个关键参数,确认你是否设置得比较合理:
1. 与力收敛相关的参数
EDIFF = 1E-6     # 能量收敛准则
EDIFFG = -0.01   # 力收敛准则(更松可设为-0.05)
建议放宽 EDIFFG,先达到一个较松弛的 ...

谢谢老师!但是老师我刚刚试了一下,使用PBE泛函的DFT-D3一个电子步花费500s左右,为什么时间更长呢。用vdw-DF2-B86r泛函花费300s左右,这个是我最初用的泛函,算了一个周还没算完。
发自小木虫手机客户端
3楼2025-07-09 11:22:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 呆瓜呀呀 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
信息提示
请填处理意见