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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】令人头疼的relax，大家帮忙分析一下问题出在那里？
北京石油化工学院2026年研究生招生接收调剂公告
101/1返回列表
查看: 1587  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕调剂 +5 椰椰。 2026-03-29 5/250 2026-03-29 20:50 by 唐沐儿
[考研] 总分293求调剂 +8 加一一九 2026-03-25 11/550 2026-03-29 19:53 by 无际的草原
[考研] 309求调剂 +7 谁不是少年 2026-03-29 7/350 2026-03-29 13:17 by mumin1990
[考研] 0856材料化工调剂 总分330 +14 zhubinhao 2026-03-27 14/700 2026-03-29 10:01 by Sjndkwm
[考研] 330分求调剂 +5 qzenlc 2026-03-29 5/250 2026-03-29 07:37 by 无际的草原
[考研] 289求调剂 +13 新时代材料 2026-03-27 13/650 2026-03-29 01:16 by 544594351
[考研] 11408软件工程求调剂 +3 Qiu学ing 2026-03-28 3/150 2026-03-28 21:50 by zhq0425
[考研] 食品工程专硕一志愿中海洋309求调剂 +4 小张zxy张 2026-03-26 8/400 2026-03-28 19:25 by lbsjt
[考研] 一志愿华理，数一英一285求A区调剂 +8 AZMK 2026-03-25 12/600 2026-03-28 18:15 by AZMK
[考研] 求化学调剂 +4 wulanna 2026-03-28 4/200 2026-03-28 13:37 by 唐沐儿
[考研] 085701环境工程，267求调剂 +16 minht 2026-03-26 16/800 2026-03-28 12:16 by zllcz
[考研] 调剂 +3 好好读书。 2026-03-28 3/150 2026-03-28 12:04 by 王保杰33
[考研] 材料与化工考研调剂 +17 孅華 2026-03-22 17/850 2026-03-28 08:35 by WYUMater
[考研] 一志愿南师大0703化学 275求调剂 +4 Ripcord上岸 2026-03-27 4/200 2026-03-27 17:00 by zhyzzh
[考研] 308求调剂 +7 墨墨漠 2026-03-25 7/350 2026-03-27 14:47 by 狂炫麦当当
[考研] 324求调剂 +8 hanamiko 2026-03-26 10/500 2026-03-27 08:06 by hypershenger
[考研] 293求调剂 +7 加一一九 2026-03-24 7/350 2026-03-25 12:02 by userper
[考研] 0854电子信息求调剂 324 +4 Promise-jyl 2026-03-23 4/200 2026-03-25 11:36 by Sugarlight
[考研] 340求调剂 +5 话梅糖111 2026-03-24 5/250 2026-03-25 06:53 by ilovexiaobin
[论文投稿] 急发核心期刊论文 +3 贤达问津 2026-03-23 5/250 2026-03-23 17:13 by 妹子不好惹
信息提示
关闭
请填处理意见
关闭