版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

ApolloYang

木虫 (正式写手)

应助: 15 (小学生)
金币: 2209.6
散金: 114
红花: 4
帖子: 390
在线: 236.9小时
虫号: 582657
注册: 2008-07-25
性别: GG
专业: 凝聚态物性 II ：电子结构

[求助] 求助：vasp结构优化不能收敛

优化材料是只含C H O的晶体
POSCAR
$i
8.710999488830570    0.000000000000000    0.000000000000000
0.000000000000000    5.169000148773190    0.000000000000000
0.000000000000000    0.000000000000000    10.962000846862800
H  C  N
8  20 20
Direct
0.2440000176429750  -0.1260000020265580 0.1790000051259990
0.0150000005960464 0.4250000119209290 0.1700000017881390
  -0.2440000176429750 0.1260000020265580 0.6790000051259990
  -0.0150000005960464  -0.4250000119209290 0.6700000017881390
0.7440000176429750 0.1260000020265580 0.1790000051259990
0.5150000005960465  -0.4250000119209290 0.1700000017881390
0.2559999823570250  -0.1260000020265580 0.6790000051259990
……
  ……
……

驰豫计算的INCAR
ENCUT = 400
ISTART = 0
ICHARGE = 2
NBANDS=300
ISMEAR = 0
SIGMA = 0.2
NSW = 100
IBRION = 3
ISIF = 5
IALGO=48
#ALGO= Fast
#ISYM = 0
POTIM = 0.2
EDIFF = 1E-4
EDIFFG = -1E-3
PREC = Accurate

问题出现在：
驰豫步数为100步，但是算完几圈之后，每次驰豫迭代圈数就很小，而且不收敛。
请问这样都优化是否有效？有什么解决方法：
（IBRION＝1，3等都试过，均出现这种情况）
……
……
RMM: 8 -0.371093682477E+03 -0.10907E-06 -0.36825E-05  1557 0.774E-03
4 F= -.37109368E+03 E0= -.37100276E+03  d E =-.236048E+00
   N    E                   dE          d eps    ncg    rms       rms(c)
RMM: 1 -0.371095287506E+03 -0.16051E-02 -0.22090E-02  2432 0.209E-01 0.128E-01
RMM: 2 -0.371095076068E+03 0.21144E-03 -0.17902E-04  2225 0.178E-02 0.141E-01
RMM: 3 -0.371094841359E+03 0.23471E-03 -0.61025E-05  1821 0.918E-03 0.557E-02
RMM: 4 -0.371094754099E+03 0.87260E-04 -0.70320E-05  1476 0.105E-02
5 F= -.37109475E+03 E0= -.37100357E+03  d E =-.107162E-02
   N    E                   dE          d eps    ncg    rms       rms(c)
RMM: 1 -0.371096827810E+03 -0.19865E-02 -0.91344E-02  2432 0.423E-01 0.254E-01
RMM: 2 -0.371095937496E+03 0.89031E-03 -0.74723E-04  2432 0.364E-02 0.279E-01
RMM: 3 -0.371095012793E+03 0.92470E-03 -0.23466E-04  2240 0.172E-02 0.112E-01
RMM: 4 -0.371094658377E+03 0.35442E-03 -0.30030E-04  2426 0.213E-02 0.300E-02
RMM: 5 -0.371094641050E+03 0.17327E-04 -0.15466E-05  1451 0.539E-03
6 F= -.37109464E+03 E0= -.37100313E+03  d E =0.113049E-03
   N    E                   dE          d eps    ncg    rms       rms(c)
RMM: 1 -0.371095354007E+03 -0.69563E-03 -0.97841E-03  2432 0.141E-01 0.862E-02
RMM: 2 -0.371095293345E+03 0.60662E-04 -0.91573E-05  2055 0.126E-02
7 F= -.37109529E+03 E0= -.37100411E+03  d E =-.652295E-03
   N    E                   dE          d eps    ncg    rms       rms(c)
RMM: 1 -0.370481562311E+03 -0.69726E-06 -0.74509E-05  1330 0.179E-02 0.875E-03
RMM: 2 -0.370481560887E+03 0.14238E-05 -0.66862E-06 896 0.735E-03
8F= -.37048156E+03 E0= -.67048156E+03  d E =-.719179E+00
   N    E                   dE          d eps    ncg    rms       rms(c)
RMM: 1 -0.370481561444E+03 0.86728E-06 -0.31566E-06 813 0.553E-03 0.318E-03
RMM: 2 -0.370481561069E+03 0.37484E-06 -0.50675E-07 480 0.372E-03
  9 F= -0.37048156E+03 E0= -.67048156E+03  d E =-.719179E+00
…………
如这次OSZCAR到第6次驰豫时候，迭代圈数就只有两步了，算到100步也未收敛

[ Last edited by ApolloYang on 2011-10-10 at 16:32 ]

回复此楼

» 猜你喜欢

售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有3人回复
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有3人回复
售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有3人回复
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急已经有3人回复
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急已经有3人回复
售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有3人回复
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有3人回复
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急已经有3人回复
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有3人回复
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急已经有4人回复

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

» 本主题相关价值贴推荐，对您同样有帮助:

关于vasp收敛和优化已经有24人回复
关于VASP一些琐事的求助（K点，吸附优化等）已经有12人回复
vasp计算：优化能收敛，但同样的设置静态计算不收敛已经有3人回复
vasp5.2 优化vasp4.6 优化过的结构不收敛已经有5人回复
【求助完结】vasp优化结构运行一段时间后没有出现错误，但是停止了！已经有20人回复
【求助】K点取得越大越不收敛，什么原因！已经有20人回复
【求助完毕】vasp-slab-边界附近原子顶位吸附小分子的结构优化很难收敛吗？已经有5人回复
【求助】vasp静态计算能量值错误已经有21人回复
【求助】vasp不同磁结构最后的收敛问题已经有6人回复
【求助】结构优化后，静态计算竟然失败？已经有8人回复
【求助】vasp 离子驰豫与电子自旋的关系已经有11人回复
【求助】VASP结构优化不收敛怎么办？已经有5人回复
【求助】帮忙编个用于VASP小程序已经有22人回复

1楼 2011-10-10 16:08:32

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

贺仪

铁杆木虫 (著名写手)

1ST强帖: 5
应助: 93 (初中生)
贵宾: 0.636
金币: 5937
散金: 1063
红花: 38
沙发: 2
帖子: 2499
在线: 588.7小时
虫号: 995125
注册: 2010-04-12
性别: GG
专业: 凝聚态物性 II ：电子结构

★
zzy870720z(金币+1): 鼓励讨论 2011-10-12 10:04:16

引用回帖:

9楼: Originally posted by uuv2010 at 2011-10-10 22:37:29:
为什么？这个算法有重大缺陷？您给科普一下吧

如下是vasp手册的说明：
A faily robust mixture of both algorithm is selected for ALGO = Fast. In this case, Davidson (IALGO = 38) is used for the in ...

不是这个意思，各个算法有自己的优缺点，在一种方法不能得出结果就换一种。我不要用这个是在前面让他去掉IALGO=48计算，结果他又加上了ALGO=48，所以让他都去点掉使用DAV算

赞一下(1人)

回复此楼

嘟啦啦嘟啦啦嘟

10楼2011-10-11 09:14:19

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

查看全部 19 个回答

贺仪

铁杆木虫 (著名写手)

1ST强帖: 5
应助: 93 (初中生)
贵宾: 0.636
金币: 5937
散金: 1063
红花: 38
沙发: 2
帖子: 2499
在线: 588.7小时
虫号: 995125
注册: 2010-04-12
性别: GG
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★
uuv2010(金币+1): 多谢提示 2011-10-10 18:18:58

去掉IALGO=48，
用IBRION = 2 跑吧

赞一下(1人)

回复此楼

嘟啦啦嘟啦啦嘟

2楼2011-10-10 16:34:03

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

fzx2008

荣誉版主 (著名写手)

1ST强帖: 9
应助: 272 (大学生)
贵宾: 1.045
金币: 21814.1
散金: 710
红花: 33
沙发: 1
帖子: 1808
在线: 1601小时
虫号: 621341
注册: 2008-10-09
性别: GG
专业: 凝聚态物性 II ：电子结构
管辖: 第一性原理

【答案】应助回帖

★
uuv2010(金币+1): 多谢提示 2011-10-10 18:18:34
ApolloYang(金币+100): 10 2011-10-13 09:39:09

我的建议是：EDIFF = 1E-4改成1E-5或者EDIFF = 1E-6
试试！

赞一下(1人)

回复此楼

3楼2011-10-10 16:39:47

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

ApolloYang

木虫 (正式写手)

应助: 15 (小学生)
金币: 2209.6
散金: 114
红花: 4
帖子: 390
在线: 236.9小时
虫号: 582657
注册: 2008-07-25
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

2楼: Originally posted by 贺仪 at 2011-10-10 16:34:03:
去掉IALGO=48，
用IBRION = 2 跑吧

用IBRION＝2
ALGO＝FAST 跑过，
问题依然存在。

赞一下

回复此楼

4楼2011-10-10 16:57:06

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

查看全部 19 个回答

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3	E4V9bALS54sj 2026-09-19	3/150	2026-09-20 01:39 by 9i7PCSZyMzot
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3	E4V9bALS54sj 2026-09-19	3/150	2026-09-20 01:27 by 9i7PCSZyMzot
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3	E4V9bALS54sj 2026-09-19	3/150	2026-09-20 01:27 by 9i7PCSZyMzot
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3	E4V9bALS54sj 2026-09-19	3/150	2026-09-20 01:15 by 9i7PCSZyMzot
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3	E4V9bALS54sj 2026-09-19	3/150	2026-09-20 01:15 by 9i7PCSZyMzot
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3	E4V9bALS54sj 2026-09-19	3/150	2026-09-20 01:03 by 9i7PCSZyMzot
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3	E4V9bALS54sj 2026-09-19	3/150	2026-09-20 00:51 by 9i7PCSZyMzot
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3	4pRyVGCVerAA 2026-09-19	3/150	2026-09-20 00:03 by 9i7PCSZyMzot
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4	4pRyVGCVerAA 2026-09-19	4/200	2026-09-20 00:03 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3	4pRyVGCVerAA 2026-09-19	4/200	2026-09-19 23:51 by 9i7PCSZyMzot
[基金申请] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3	4pRyVGCVerAA 2026-09-19	4/200	2026-09-19 23:39 by 9i7PCSZyMzot
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3	4pRyVGCVerAA 2026-09-19	4/200	2026-09-19 23:27 by 9i7PCSZyMzot
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3	4pRyVGCVerAA 2026-09-19	4/200	2026-09-19 23:27 by 9i7PCSZyMzot
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3	4pRyVGCVerAA 2026-09-19	4/200	2026-09-19 23:03 by 9i7PCSZyMzot
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3	4pRyVGCVerAA 2026-09-18	7/350	2026-09-19 21:03 by 9i7PCSZyMzot
[教师之家] 现代”学阀”该如何界定 +3	iaeyuan 2026-09-19	3/150	2026-09-19 20:52 by 没昵称呀
[基金申请] 国社科系统bug了，是不是要放榜了？ +7	kynobel 2026-09-16	9/450	2026-09-19 14:52 by 接好运婷
[考博] 申博发邮件 +8	Lmengk 2026-09-15	9/450	2026-09-19 09:27 by 中间体超市
[有机交流] 各位大神，目前国内有哪些比较好用的逆合成软件？ 20+3	zapen 2026-09-15	9/450	2026-09-18 18:17 by zapen
[考博] 上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告 +7	水士口 2026-09-16	8/400	2026-09-16 13:47 by kkklens