24小时热门版块排行榜    

查看: 2744  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 相丽温 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +4 Ldrop2023 2026-08-13 4/200 2026-08-14 14:50 by Tide man
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[基金申请] 好奇怪的filecode +6 布布和一二 2026-08-08 7/350 2026-08-14 13:38 by Equinoxhua
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +19 医学老男孩 2026-08-13 43/2150 2026-08-14 13:34 by 医学老男孩
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +5 工厂打螺丝 2026-08-08 5/250 2026-08-14 11:03 by 玖戈弋
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
信息提示
请填处理意见