计算化学公社

 找回密码 Forget password
 注册 Register
Views: 756|回复 Reply: 3
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] LiNiO2的slab模型掺杂W原子后做结构优化,W周围的氧原子受力不收敛怎么办?

[复制链接 Copy URL]

1

帖子

0

威望

73

eV
积分
74

Level 2 能力者

本帖最后由 1275716025 于 2024-3-20 15:54 编辑

参照了这个帖子修改了各种参数后力都不能收敛http://bbs.keinsci.com/forum.php?mod=viewthread&tid=13065
粗优化时使用EDIFFG=-0.1是可以收敛的,但是改到EDIFFG=-0.05就不能收敛了,用脚本查了一下发现是掺杂的W原子右侧的氧原子受力始终降不到0.05以下。求教各位老师,这种情况应该修改什么来让力收敛?

不收敛的INCAR如下:
ISTART=0
ICHARG=2
PREC=N
ISPIN=2
MAGMOM=18*1  36*2  17*2  2
LDAU=.T.
LDAUL=-1  -1  2  -1
LDAUU=0.0  0.0   6.4   0.0
LDAUJ=0  0  0  0
                                    
EDIFF=1E-7
ENCUT=450
IBRION=1
AMIX=0.1
MAXMIX=50
POTIM=0.3


NELMDL=-20
NELM=200
NELMIN=6
NSW=300


ISIF=2
EDIFFG=-0.05
ISMEAR=0
SIGMA=0.05


LWAVE=.F.
LCHARG=.F.
LREAL=Auto
LDIPOL=.TRUE.
IDIPOL=3

QQ截图20240320150422.jpg (47.25 KB, 下载次数 Times of downloads: 25)

灰黑色为掺杂W原子

灰黑色为掺杂W原子

3

帖子

0

威望

137

eV
积分
140

Level 2 能力者

2#
发表于 Post on 2024-4-24 09:01:51 | 只看该作者 Only view this author
出现类似问题,请问楼主解决了吗?也是按照这个帖子办法逐一尝试,尝试将IBRION=2修改为1,POTIM=0.5修改为0.1,但力却始终在0.07附近波动,无法收敛

1

帖子

0

威望

181

eV
积分
182

Level 3 能力者

3#
发表于 Post on 2024-8-3 02:13:16 | 只看该作者 Only view this author
请问是否可以分享下查找达不到力收敛标准原子的脚本吗?

7

帖子

0

威望

41

eV
积分
48

Level 2 能力者

4#
发表于 Post on 2024-11-6 10:29:41 | 只看该作者 Only view this author
后来解决了吗,可否分享一下办法

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-26 23:53 , Processed in 0.566894 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list