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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】令人头疼的relax，大家帮忙分析一下问题出在那里？
101/1返回列表
查看: 1558  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学调剂 ，六级已过，有科研经历 +11 曦熙兮 2026-03-15 11/550 2026-03-19 17:13 by yywzz
[考研] 材料专硕274一志愿陕西师范大学求调剂 +8 薛云鹏 2026-03-13 8/400 2026-03-19 15:36 by haoshis
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[教师之家] 焦虑 +9 水冰月月野兔 2026-03-13 13/650 2026-03-19 09:50 by otani
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 【同济软件】软件（085405）考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 一志愿武理材料305分求调剂 +5 想上岸的鲤鱼 2026-03-18 6/300 2026-03-18 17:53 by 无际的草原
[考研] 0854可跨调剂，一作一项核心论文五项专利，省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 297求调剂 +8 戏精丹丹丹 2026-03-17 8/400 2026-03-18 14:30 by laoshidan
[考研] 331求调剂（0703有机化学 +7 ZY-05 2026-03-13 8/400 2026-03-18 14:13 by 007_lilei
[考研] 304求调剂 +12 小熊joy 2026-03-14 13/650 2026-03-18 12:34 by Linda Hu
[考研] 334求调剂 +3 志存高远意在机� 2026-03-16 3/150 2026-03-18 08:34 by lm4875102
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 26考研求调剂 +6 丶宏Sir 2026-03-13 6/300 2026-03-17 16:13 by 醉在风里
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
信息提示
关闭
请填处理意见
关闭