| 查看: 2325 | 回复: 2 | ||
| 【悬赏金币】回答本帖问题，作者呆瓜呀呀将赠送您 100 个金币 | ||
[求助]
[VASP]请教一下结构优化能量一直下降但力不收敛INCAR是否有问题 已有1人参与
|
||
|
老师们好，我最近在计算两层石墨烯转角体系，目前计算的是28个原子，但是已经计算了一周，能量一直在下降，但是力一直不收敛。初始结构是根据文献构建的。用的16个核。因为最大离子步设100，所以这是接着算的第三次了。一个电子步大概是两百多秒。是INCAR文件有问题吗？因为原子体系不是很大。INCAR每一个参数都是对照着别的帖子设的。 通过网盘分享的文件：INCAR等4个文件 链接: https://pan.baidu.com/s/1rwqodS8ED1eBJq33-MShFA?pwd=ex2q 提取码: ex2q 发自小木虫手机客户端 |
回复此楼» 猜你喜欢
跳槽后在研项目怎么办？
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有9人回复
» 本主题相关商家推荐: (我也要在这里推广)
rlafite
至尊木虫 (正式写手)
- 应助: 199 (高中生)
- 金币: 11352.5
- 红花: 29
- 沙发: 1
- 帖子: 956
- 在线: 361.9小时
- 虫号: 4563202
- 注册: 2016-04-02
- 性别: GG
- 专业: 同步辐射技术及其应用
【答案】应助回帖
感谢参与，应助指数 +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
|
谢谢老师！但是老师我刚刚试了一下，使用PBE泛函的DFT-D3一个电子步花费500s左右，为什么时间更长呢。用vdw-DF2-B86r泛函花费300s左右，这个是我最初用的泛函，算了一个周还没算完。 发自小木虫手机客户端 |
3楼2025-07-09 11:22:07
40