24小时热门版块排行榜    

查看: 2347  |  回复: 11

相丽温

新虫 (著名写手)

[求助] 吸附体系结构优化难收敛 已有2人参与

在优化清洁FeS表面时,可收敛,但是水分子吸附在FeS表面时难以收敛,其中INCAR,KPOINTS参数均相同,下图是我的参数设置以及输出文件OUTCAR中的POSITION和OSZICAR的部分信息。
问题:为何小分子吸附后模型难以收敛?应该如何处理模型和参数,才能使吸附模型优化尽快收敛。

望大家不吝赐教!!!谢谢!!!

吸附体系结构优化难收敛
1.png


吸附体系结构优化难收敛-1
2.png


吸附体系结构优化难收敛-2
3.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

xyj880506

银虫 (小有名气)

可以用较小的POTIM,比如0.05
还有就是电子步的精度,可以试试更高的。
也可以提高截断能,比如用500
4楼2017-04-07 11:41:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

相丽温

新虫 (著名写手)

有做表面吸附的吗?还请给点建议,困惑好久了,不知道该怎么处理?谢谢了!
2楼2017-04-06 22:36:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

前田庆次

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
相丽温: 金币+5 2017-04-08 11:06:40
你把你INCAR打成文字,不要图片,我可以给你提提意见
本来吸附分子收敛就慢,有时候经常几百步
水分子吸附,极有可能是物理吸附哦
长枪一横花飘零,松风追月伴我行。无双人间世如梦,倾奇万世永留名。
3楼2017-04-07 10:56:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

冲天飞龙

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
相丽温: 金币+5, 有帮助 2017-04-08 11:13:32
你这是电子步吧,第一个离子步的循环跑几十个电子步很正常啊,毕竟刚开始的wavefuntion是random的,要跑很多电子步来逼近求解
5楼2017-04-07 23:09:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

相丽温

新虫 (著名写手)

引用回帖:
3楼: Originally posted by 前田庆次 at 2017-04-07 10:56:53
你把你INCAR打成文字,不要图片,我可以给你提提意见
本来吸附分子收敛就慢,有时候经常几百步
水分子吸附,极有可能是物理吸附哦

首先,非常感谢您的回复!

SYSTEM = FeS001-H2O
ENCUT = 400
PREC = Normal
EDIFF = 1E-5
EDIFFG = -0.01
LREAL  = .FALSE.
ISTART = 0
ICHARG = 2
ISMEAR = 1
SIGMA = 0.2
IBRION = 2
ISIF = 2
NSW = 100
POTIM = 0.5
NELM = 60
LWAVE = .FALSE.
LCHARG = .FALSE.
ALGO = Fast

VDW_RADIUS= 50.0
VDW_S6=0.75
VDW_SR=1.00
VDW_SCALING =0.75
VDW_D = 20.0
VDW_C6=10.80 5.57  0.700  0.140
VDW_R0=1.562 1.683 1.342  1.001
LVDW=.TRUE.
这是我的INCAR,迭代100步不收敛,然后把CONTCAR拷贝成POSCAR设NSW=300,POTIM=0.2迭代300步还是不收敛。
这个是物理吸附,水分子中的H和表面的S形成弱氢键。
我还想问一下,迭代这么多步都不收敛是不是构建模型的问题,还有一般迭代多少步收敛正常?
我是新手,刚刚开始学习计算,还请多多指教,谢谢!
6楼2017-04-08 11:06:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

相丽温

新虫 (著名写手)

引用回帖:
4楼: Originally posted by xyj880506 at 2017-04-07 11:41:54
可以用较小的POTIM,比如0.05
还有就是电子步的精度,可以试试更高的。
也可以提高截断能,比如用500

首先谢谢您的回复!
我把迭代100步的CONTCAR拷贝成POSCAR,修改NSW=300,POTIM=0.2迭代300步还是不能收敛!提高电子步的精度是EDIFFG=-0.01,改为-0.02吗?
我是第一次计算吸附,不知道吸附体系迭代这么多步还不收敛正常吗?
7楼2017-04-08 11:10:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

相丽温

新虫 (著名写手)

引用回帖:
5楼: Originally posted by 冲天飞龙 at 2017-04-07 23:09:50
你这是电子步吧,第一个离子步的循环跑几十个电子步很正常啊,毕竟刚开始的wavefuntion是random的,要跑很多电子步来逼近求解

谢谢您的回复!
可是设置NSW=100迭代100次还不收敛,复制CONTCAR为POSCAR,修改NSW=300,POTIM=0.2,迭代300次还是不能收敛!
迭代这么多次不收敛正常吗?
8楼2017-04-08 11:13:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangfrank

金虫 (正式写手)

EDIFF = 1E-5
EDIFFG = -0.01这俩精度降低一些,分别为1E-4 -0.02试试吧

发自小木虫Android客户端
计算顺利点吧
9楼2017-04-08 13:16:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

相丽温

新虫 (著名写手)

引用回帖:
9楼: Originally posted by zhangfrank at 2017-04-08 13:16:16
EDIFF = 1E-5
EDIFFG = -0.01这俩精度降低一些,分别为1E-4 -0.02试试吧

谢谢您的回复,我试试看~~~

10楼2017-04-08 15:02:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 相丽温 的主题更新
信息提示
请填处理意见