24小时热门版块排行榜    

查看: 4461  |  回复: 15

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 论文编辑与投稿

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

golddoushi

木虫 (正式写手)

ALGO=Normal
就行了
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个核。。。
瞅瞅是不是你这个参数设置过大,导致内存不够引起的错误
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
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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-28 05:13 by winsaint
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +11 晴天加油 2026-08-26 12/600 2026-08-28 01:22 by li09465290@
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 XLGfUIdOmEAO 2026-08-27 3/150 2026-08-28 00:49 by 3RAFyYXYSwGS
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 XLGfUIdOmEAO 2026-08-27 4/200 2026-08-27 22:26 by 3RAFyYXYSwGS
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 出来了 +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
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +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 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见