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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】令人头疼的relax，大家帮忙分析一下问题出在那里？
101/1返回列表
查看: 1875  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

vaspstudy

银虫 (小有名气)

[交流] 【求助】令人头疼的relax，大家帮忙分析一下问题出在那里？ 已有3人参与

做的ZnO中N间隙，relax离子步不收敛，大家帮忙分析一下问题出在那里：
INCAR
ENCUT    =   520.0
ISTART   =   0   
ICHARG   =   2
ISMEAR   =   -5
NSW      =   200
IBRION   =   2
ISIF     =   2
POTIM    =   0.2
EDIFF    =   1E-4
EDIFFG   =   -0.05
PREC     =   Accurate
LREAL    =   .TRUE.
LWAVE    =   .FALSE.


运行信息（目前还在运行，一点收敛的意思都没有；由于帖子长度限制，省去了电子迭代信息）：
   1 F= -.17031400E+03 E0= -.17031400E+03  d E =-.170314E+03                           
   
   2 F= -.17062041E+03 E0= -.17062041E+03  d E =-.306407E+00                                 
   
   3 F= -.17083441E+03 E0= -.17083441E+03  d E =-.520408E+00                                 

   4 F= -.17090323E+03 E0= -.17090323E+03  d E =-.688234E-01                                 

   5 F= -.17096062E+03 E0= -.17096062E+03  d E =-.126209E+00                                 

   6 F= -.17100118E+03 E0= -.17100118E+03  d E =-.405607E-01                                 

   7 F= -.17103422E+03 E0= -.17103422E+03  d E =-.736048E-01                                 

   8 F= -.17105799E+03 E0= -.17105799E+03  d E =-.237658E-01                                 

   9 F= -.17105953E+03 E0= -.17105953E+03  d E =-.253081E-01                                 

  10 F= -.17105694E+03 E0= -.17105694E+03  d E =-.227136E-01                                 

  11 F= -.17105972E+03 E0= -.17105972E+03  d E =-.255003E-01                                 

  12 F= -.17107533E+03 E0= -.17107533E+03  d E =-.156108E-01                                 

  13 F= -.17108304E+03 E0= -.17108304E+03  d E =-.233180E-01                                 

  14 F= -.17108466E+03 E0= -.17108466E+03  d E =-.249336E-01                                 

  15 F= -.17109913E+03 E0= -.17109913E+03  d E =-.144775E-01                                 

  16 F= -.17110209E+03 E0= -.17110209E+03  d E =-.174314E-01                                 

  17 F= -.17110262E+03 E0= -.17110262E+03  d E =-.179646E-01                                 

  18 F= -.17111000E+03 E0= -.17111000E+03  d E =-.737632E-02                                 

  19 F= -.17110972E+03 E0= -.17110972E+03  d E =-.709768E-02                                 

  20 F= -.17111338E+03 E0= -.17111338E+03  d E =-.365639E-02                                 

  21 F= -.17111439E+03 E0= -.17111439E+03  d E =-.466673E-02                                 

  22 F= -.17112106E+03 E0= -.17112106E+03  d E =-.667709E-02                                 

  23 F= -.17112291E+03 E0= -.17112291E+03  d E =-.852443E-02                                 

  24 F= -.17112826E+03 E0= -.17112826E+03  d E =-.534966E-02                                 

  25 F= -.17113908E+03 E0= -.17113908E+03  d E =-.161667E-01                                 

  26 F= -.17113840E+03 E0= -.17113840E+03  d E =-.154921E-01                                 

  27 F= -.17114328E+03 E0= -.17114328E+03  d E =-.203736E-01                                 

  28 F= -.17114359E+03 E0= -.17114359E+03  d E =-.206772E-01                                 

  29 F= -.17114438E+03 E0= -.17114438E+03  d E =-.214651E-01                                 

  30 F= -.17114639E+03 E0= -.17114639E+03  d E =-.201855E-02                                 

  31 F= -.17115177E+03 E0= -.17115177E+03  d E =-.739082E-02                                 

  32 F= -.17115726E+03 E0= -.17115726E+03  d E =-.128875E-01                                 

  33 F= -.17116538E+03 E0= -.17116538E+03  d E =-.811287E-02                                 

  34 F= -.17118634E+03 E0= -.17118634E+03  d E =-.290753E-01                                 

  35 F= -.17120964E+03 E0= -.17120964E+03  d E =-.523775E-01                                 

  36 F= -.17124173E+03 E0= -.17124173E+03  d E =-.320904E-01                                 

  37 F= -.17133467E+03 E0= -.17133467E+03  d E =-.125029E+00                                 

  38 F= -.17147823E+03 E0= -.17147823E+03  d E =-.268590E+00                                 

  41 F= -.17149772E+03 E0= -.17149772E+03  d E =-.288078E+00                                 

  42 F= -.17150862E+03 E0= -.17150862E+03  d E =-.109035E-01                                 

  43 F= -.17153884E+03 E0= -.17153884E+03  d E =-.411173E-01                                 

  44 F= -.17158220E+03 E0= -.17158220E+03  d E =-.844868E-01                                 

  45 F= -.17154379E+03 E0= -.17154379E+03  d E =-.460673E-01                                 

  46 F= -.17159143E+03 E0= -.17159143E+03  d E =-.937155E-01                                 

  47 F= -.17164903E+03 E0= -.17164903E+03  d E =-.575946E-01                                 

  48 F= -.17179773E+03 E0= -.17179773E+03  d E =-.206294E+00
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-09-27 16:19:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
vaspstudy(金币+2): 2010-09-27 17:38:06
建议：
PREC=NORMAL
ISMEAR =0,SIGMA=0.5
ENCUT=400
POTIM=0.4
EDIFFG=-0.1
NSW=30
收敛后再来详细优化。
赞一下(1人) 回复此楼
2楼2010-09-27 17:05:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vaspstudy

银虫 (小有名气)
引用回帖:
Originally posted by cenwanglai at 2010-09-27 17:05:21:
建议：
PREC=NORMAL
ISMEAR =0,SIGMA=0.5
ENCUT=400
POTIM=0.4
EDIFFG=-0.1
NSW=30
收敛后再来详细优化。

看您的参数设置应该是进行了一个很粗的relax。但您提到的“收敛后再来详细优化”不是很懂，是指用以上粗优的离子位置，然后按照我上面提供的INCAR再细优一遍吗？这样优化有什么好处呢？还望指点！
赞一下 回复此楼
3楼2010-09-27 17:13:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
mazuju028(金币+2):谢谢交流 2010-09-27 17:39:52
离子步收敛困难大部分是因为初始结构不合理。所以先降低精度优化，在更宽的范围内较快的找到一个比较合理的结构。然后再继续提高精度优化。

这就像是要去浙大图书馆找书，先找到校区，再找图书馆，再找楼层和架位。
赞一下(2人) 回复此楼
4楼2010-09-27 17:23:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
vaspstudy(金币+2): 2010-09-27 17:38:20
cenwanglai(金币+1):呵呵，静态计算我也常用ISMEAR=0.半导体超胞。 2010-09-27 17:43:47
离子驰豫为什么用ISMEAR=-5呢？
这很费时，修正的四面体方法一般用于静态计算。
建议设置ISMEAR=0.
赞一下(2人) 回复此楼
5楼2010-09-27 17:33:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hughluhan

金虫 (小有名气)
★
vaspstudy(金币+1):谢谢参与
引用回帖:
Originally posted by vaspstudy at 2010-09-27 17:13:03:

看您的参数设置应该是进行了一个很粗的relax。但您提到的“收敛后再来详细优化”不是很懂，是指用以上粗优的离子位置，然后按照我上面提供的INCAR再细优一遍吗？这样优化有什么好处呢？还望指点！

一般你开始时候不收敛是因为你的初始结构不太合理，所以才会出现不收敛的情况，所以你要先找出一个较为粗的的范围，然后再在这个范围周围寻找精确的位置。所以就是这个意思。
赞一下(1人) 回复此楼
6楼2010-09-27 22:44:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hughluhan

金虫 (小有名气)
vaspstudy(金币+1): 2010-09-28 15:29:06
引用回帖:
Originally posted by stractor at 2010-09-27 17:33:46:
离子驰豫为什么用ISMEAR=-5呢？
这很费时，修正的四面体方法一般用于静态计算。
建议设置ISMEAR=0.

这是你总结出来的，还是有什么依据吗？请教一下！
赞一下 回复此楼
7楼2010-09-27 22:46:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)
★
vaspstudy(金币+1):谢谢参与
引用回帖:
Originally posted by hughluhan at 2010-09-27 22:46:59:

这是你总结出来的，还是有什么依据吗？请教一下！

rtfm  
回复此楼
8楼2010-09-28 16:24:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
★
vaspstudy(金币+1): 2010-09-28 19:34:04
ym23(金币+1):欢迎常来交流 2010-09-28 19:44:07
引用回帖:
Originally posted by hughluhan at 2010-09-27 22:46:59:

这是你总结出来的，还是有什么依据吗？请教一下！

不是我总结出来的。VASP手册上和文献上都这么做的。ISMEAR=-5是修正的四面体方法，这个方法用于静态计算，主要是为了得到较准的能量。由于计算费时，所以不推荐用于离子驰豫。此外，ZnO是半导体，半导体用ISMEAR=0(VASP 手册）.
赞一下(1人) 回复此楼
9楼2010-09-28 18:41:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hughluhan

金虫 (小有名气)
引用回帖:
Originally posted by stractor at 2010-09-28 18:41:18:

不是我总结出来的。VASP手册上和文献上都这么做的。ISMEAR=-5是修正的四面体方法，这个方法用于静态计算，主要是为了得到较准的能量。由于计算费时，所以不推荐用于离子驰豫。此外，ZnO是半导体，半导体用ISMEA ...

哦，是我忽略了，谢了。
赞一下 回复此楼
10楼2010-09-28 18:47:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 vaspstudy 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +17 爱看书的可乐 2026-08-26 19/950 2026-08-30 22:24 by xiangy672
[基金申请] 面上意见出来了 +9 黄鸟于飞Chao 2026-08-29 18/900 2026-08-30 16:47 by 黄鸟于飞Chao
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
关闭
请填处理意见
关闭