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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 求助大神， VASP 使用杂化泛函结构优化一直报错
161/2返回列表12下一页
查看: 3767  |  回复: 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下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭