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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 求助大神， VASP 使用杂化泛函结构优化一直报错
161/2返回列表12下一页
查看: 3948  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

junfeiyan

新虫 (小有名气)

[求助] 求助大神， VASP 使用杂化泛函结构优化一直报错 已有5人参与

小弟最近使用杂化泛函进行结构优化，原因是我算的是含f电子稀土化合物， vasp 给的赝势用GGA+U 收敛极其困难，用杂化泛函好些，已经有发表在nature 子篇的论文，做过类似的体系，用的杂化泛函优化， 求助各位大神，看看是小弟的输入文件有问题，还是天河超算中心的VASP 编译问题，感谢大神
小弟的INCAR

Name of system
   SYSTEM    =     a

Start parameter
   ISTART    =     0
   INIWAV    =     1
   ICHARG    =     2

Write flags
   NWRITE    =     2
   LWAVE     =     T
   LCHARG    =     T

Parallelization
   NPAR      =     72
   LPLANE    =     T
   LSCALU    =     F
   LSCALAPACK =    F

Electronic Relaxation 1
   ENCUT     =     500.00  eV
   EDIFF     =     1.0E-05  eV
   PREC      =     High
   NELM      =     200
   NELMIN    =     2

Electronic Relaxation 2
   ALGO      =     Damped
   TIME      =     0.4
   LREAL     =     Auto

Ionic Relaxation
   NSW       =     150
   IBRION    =     2
   ISIF      =     3
   EDIFFG    =     -0.05
   POTIM     =     0.3
   ISYM      =     3
   IWAVPR    =     1
   LCORR     =     T

DOS related values
   ISMEAR    =     0
   SIGMA     =     0.025
   NEDOS     =     3000
#   NBANDS    =     3958

Magnetic parameter
   LORBIT    =     11
Exchange-correlation functional
   GGA       =     PS
   ISPIN = 2
  MAGMOM = 4*8 4*4 12*0
   HSE06 parameter
   LHFCALC   =     T
   HFSCREEN  =     0.2
   PRECFOCK  =     Fast
#   NKRED     =     2

还有一个问体就是我想用HSE 考虑稀土元素的f 电子， 那么我选择赝势的时候是否可以用PBE 里面不含f 电子的赝势， 不知道杂化泛函会不会自己考虑f 电子
前提是我用GGA+U 算不含f电子的赝势，很快收敛， 用含f电子的赝势极难收敛
报错一
ROT:  1   0.159495398721E+04    0.17831E+04BRMIX: very serious problems
the old and the new charge density differ
old charge density:   160.00000 new***********
    0.603E+06
ERROR FEXCF: supplied exchange-correlation table
    is too small, maximal index :      912695

这个报错有人说应该减小POTIM。 有成功的经验吗
报错2
ROT:  1  -0.155058592914E+03   -0.30390E+02    0.284E+01
   ROT:  2  -0.141575035880E+03    0.13484E+02    0.160E+01
   ROT:  3  -0.139046123851E+03    0.25289E+01    0.114E+01
   ROT:  4  -0.139495088906E+03   -0.44897E+00    0.453E+00
   ROT:  5  -0.139199449050E+03    0.29564E+00    0.133E+00
gam= 0.382 g(H,U,f)=  0.305E+03 0.218E+02 0.320E+00 ort(H,U,f) =-0.188E+03-0.426E+02-0.465E+00
DMP:   7    -0.112154452901E+03   -0.50536E+02   -0.95574E+02  2160   0.327E+03-0.231E+03
         N     E                     dE            rms(c)
   ROT:  1  -0.159501125680E+03   -0.20006E+02    0.795E+01
   ROT:  2  -0.157770040931E+03    0.17311E+01    0.495E+01
   ROT:  3  -0.155549182836E+03    0.22209E+01    0.173E+01
   ROT:  4  NaN                    NaN            NaN      
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image              PC                Routine            Line        Source            
vasp               00000000005A769A  Unknown               Unknown  Unknown
vasp               000000000070A100  Unknown               Unknown  Unknown

求助大神， VASP 使用杂化泛函结构优化一直报错
U31~GV(QLRG}V{D75RS6[]8.png
回复此楼

» 收录本帖的淘帖专辑推荐

soft study 论文编辑与投稿

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-08-28 09:52:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

golddoushi

木虫 (正式写手)
ALGO=Normal
就行了
赞一下(1人) 回复此楼
2楼2015-08-28 11:36:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lshok

铁杆木虫 (著名写手)

【答案】应助回帖

★
感谢参与，应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-09-05 18:29:29
NPAR=72?
VASP 推荐NPAR＝sqrt（cores），你设置的NPAR＝72对应5184个核。。。
瞅瞅是不是你这个参数设置过大，导致内存不够引起的错误
赞一下(1人) 回复此楼
5楼2015-08-29 16:44:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)

【答案】应助回帖

★
liliangfang: 金币+1, 谢谢交流 2015-09-05 18:29:36
引用回帖:
3楼: Originally posted by junfeiyan at 2015-08-28 15:15:50
谢谢您，麻烦在问一下，你说的这个参数，是对应哪个报错的呢，还有帮我看看其他参数有没有需要修改的呢，小弟没怎么用过杂化泛函...

上面的人肯定没算过杂化泛函。ALGO=Normal不用于杂化泛函计计算， 对于杂化泛函需要ALGO=ALL or Damped
仔细看看vasp manual http://cms.mpi.univie.ac.at/vasp ... k_calculations.html
赞一下(1人) 回复此楼
7楼2015-09-01 13:53:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

junfeiyan

新虫 (小有名气)
引用回帖:
2楼: Originally posted by golddoushi at 2015-08-28 11:36:41
ALGO=Normal
就行了

谢谢您，麻烦在问一下，你说的这个参数，是对应哪个报错的呢，还有帮我看看其他参数有没有需要修改的呢，小弟没怎么用过杂化泛函
赞一下 回复此楼
3楼2015-08-28 15:15:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zuocuiping

木虫 (职业作家)
能把你说的参考文献的题目发给我吗
我想看一下
谢谢
赞一下 回复此楼
4楼2015-08-29 12:50:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

junfeiyan

新虫 (小有名气)
引用回帖:
5楼: Originally posted by lshok at 2015-08-29 16:44:55
NPAR=72?
VASP 推荐NPAR＝sqrt（cores），你设置的NPAR＝72对应5184个核。。。
瞅瞅是不是你这个参数设置过大，导致内存不够引起的错误

谢谢，回复，我也了解前辈说的这个参数，很多人推荐设置您说的核数的平方根， 也有人说，直接设置为核数也行，我再测试测试
赞一下 回复此楼
6楼2015-08-29 17:53:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)

【答案】应助回帖

★
liliangfang: 金币+1, 谢谢交流 2015-09-05 18:29:43
引用回帖:
6楼: Originally posted by junfeiyan at 2015-08-29 17:53:56
谢谢，回复，我也了解前辈说的这个参数，很多人推荐设置您说的核数的平方根， 也有人说，直接设置为核数也行，我再测试测试...

这回复。。。。瞎说！ 杂化泛函要要求NPAR=#of cores
你的NPAR设置没问题
赞一下 回复此楼
8楼2015-09-01 13:54:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)

【答案】应助回帖

★
liliangfang: 金币+1, 谢谢交流 2015-09-05 18:29:51
你的问题我来回答一下吧

1. ERROR FEXCF: supplied exchange-correlation table
    is too small, maximal index :      912695
这个问题是由于intel编译器的问题，（我猜）天河上用intel可能是2013 or 2015。你用intel2011编译的vasp运行是不会出现这个错误的。
BTW，这个问题偶尔通过减小SIGMA 还是POTIM可以得到解决(以前试过，忘了，anyway, 这不是终极解决方法）

2. ROT:  4  NaN                    NaN            NaN      
forrtl: severe (174): SIGSEGV, segmentation fault occurred
这个错误明显是你体系energy 不收敛。你看energy 都NaN了！所以才死掉。
4f电子计算有很多tricks。你慢慢研究研究吧。
赞一下 回复此楼
9楼2015-09-01 13:59:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

junfeiyan

新虫 (小有名气)
引用回帖:
9楼: Originally posted by emilyoyang at 2015-09-01 13:59:55
你的问题我来回答一下吧

1. ERROR FEXCF: supplied exchange-correlation table
    is too small, maximal index :      912695
这个问题是由于intel编译器的问题，（我猜）天河上用intel可能是2013 or 2015。 ...

谢谢，谢谢，希望前辈介绍一下关于4f电子计算的经验，让晚生少有些弯路！

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
10楼2015-09-01 21:09:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 junfeiyan 的主题更新
161/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 314求调剂 +8 无懈可击的巨人 2026-03-12 8/400 2026-03-18 14:50 by haxia
[考研] 085601专硕，总分342求调剂，地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 化工学硕306求调剂 +10 42838695 2026-03-12 10/500 2026-03-18 14:42 by haxia
[考研] 0703化学调剂 +4 pupcoco 2026-03-17 7/350 2026-03-18 12:14 by djl2006
[考研] 材料，纺织，生物（0856、0710），化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考研] 环境工程调剂 +8 大可digkids 2026-03-16 8/400 2026-03-18 09:36 by zhukairuo
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 考研化学学硕调剂，一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 10/500 2026-03-17 15:07 by 一鸭鸭哟
[考研] 材料与化工专硕调剂 +5 heming3743 2026-03-16 5/250 2026-03-17 14:03 by 勇敢太监王公公
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 0703 物理化学调剂 +3 我可以上岸的对� 2026-03-13 5/250 2026-03-16 10:50 by 我可以上岸的对�
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 土木第一志愿276求调剂，科研和技能十分丰富，求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生，2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
信息提示
关闭
请填处理意见
关闭