24小时热门版块排行榜    

查看: 1513  |  回复: 5

whw19850730

铁杆木虫 (小有名气)

[求助] VASP 电子迭代发散 诡异问题

我做一个六角晶体LuMnO3,电子迭代发散,遇到不收敛问题可以增加 K点,改变
电荷密度混合参数等很多办法,可是我的不是迭代 震荡,而是迭代能量越来越大最后崩溃掉,我调了很久都解决不掉,希望高手大牛给些意见。  
     N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.242426654415E+03   -0.24243E+03   -0.49552E+02 12636   0.152E+01
DAV:   2    -0.244734802042E+03   -0.23081E+01   -0.13019E+01  8820   0.106E+00
DAV:   3    -0.244911974268E+03   -0.17717E+00   -0.17362E+00 12996   0.458E-01
DAV:   4    -0.244917066602E+03   -0.50923E-02   -0.45943E-02  8136   0.601E-02
DAV:   5    -0.244917150324E+03   -0.83722E-04   -0.83640E-04  5904   0.106E-02
DAV:   6    -0.247540272456E+03   -0.26231E+01   -0.31613E+00  5976   0.408E-03
DAV:   7    -0.246111205273E+03    0.14291E+01   -0.36981E+01  5760   0.395E-01
DAV:   8    -0.260796035821E+03   -0.14685E+02   -0.25748E+01  5760   0.277E-01    0.148E+00
DAV:   9    -0.244897893351E+03    0.15898E+02   -0.23689E+01 12600   0.638E+00    0.122E+00
DAV:  10    -0.244879497047E+03    0.18396E-01   -0.12455E-01 12564   0.182E+01    0.761E-01
DAV:  11    -0.244883124169E+03   -0.36271E-02   -0.56369E-02 12564   0.110E+01    0.688E-01
DAV:  12    -0.244879821921E+03    0.33022E-02   -0.94784E-03 13140   0.449E+00    0.367E-01
DAV:  13    -0.245058699872E+03   -0.17888E+00   -0.16888E+00 13032   0.277E+00    0.294E-01
DAV:  14    -0.252709857733E+03   -0.76512E+01   -0.84904E+01 12960   0.465E+00    0.152E-01
DAV:  15    -0.245505026950E+03    0.72048E+01   -0.31400E+01 13068   0.124E+00    0.641E-02
DAV:  16    -0.245391127504E+03    0.11390E+00   -0.13041E+01 12888   0.127E+00    0.531E-02
DAV:  17    -0.244879335901E+03    0.51179E+00   -0.69883E+00 12348   0.371E-01    0.707E-02
DAV:  18    -0.267404373109E+03   -0.22525E+02   -0.22226E+02 12996   0.465E-01    0.757E-01
DAV:  19    -0.244846617464E+03    0.22558E+02   -0.92693E+01 10080   0.445E+00    0.452E-02
DAV:  20    -0.317293990069E+03   -0.72447E+02   -0.71922E+02 13284   0.191E-01    0.566E-01
DAV:  21    -0.520146965043E+05   -0.51697E+05   -0.51779E+05  9288   0.219E+00    0.187E+00
DAV:  22    -0.248042514172E+03    0.51767E+05   -0.54223E+02  8028   0.761E+00    0.460E-01
DAV:  23    -0.543114596657E+04   -0.51831E+04   -0.63755E+02 10764   0.247E+00    0.164E+00
DAV:  24    -0.589055666359E+03    0.48421E+04   -0.36204E+03 10296   0.657E+00    0.713E-01
DAV:  25    -0.308703180751E+06   -0.30811E+06   -0.24597E+04 10044   0.731E+00
   1 F= -.30870318E+06 E0= -.30870318E+06  d E =-.308703E+06  mag=     0.0000
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

oupengfei717

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
楼主这个是静态计算吗?怎么只有一步离子迭代计算呢?我猜想可能是该结构在初期全弛豫的时候没有能找到平衡结构?不知道楼主是否进行了弛豫得到了平衡结构?楼主的输入文件是否都检查过不存在问题呢?
2楼2013-08-17 10:31:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whw19850730

铁杆木虫 (小有名气)

引用回帖:
2楼: Originally posted by oupengfei717 at 2013-08-17 10:31:52
楼主这个是静态计算吗?怎么只有一步离子迭代计算呢?我猜想可能是该结构在初期全弛豫的时候没有能找到平衡结构?不知道楼主是否进行了弛豫得到了平衡结构?楼主的输入文件是否都检查过不存在问题呢?

我打算做几何优化的,第一步电子迭代就发散了,之所以迭代25步,是因为我做测试设定的NELM =25。

我是输入文件

SYSTEM = LuMnO3 calculation
ISTART = 0
ICHARG = 2
ISPIN = 2
LREAL = F
#PREC= Accurate
#ADDGRID=.TRUE.

Electronic Relaxation
ENCUT = 500
EDIFF = 1E-07
GGA=PS
NELM =40
NELMDL=8
TIME = 0.1

Ionic relaxation
NSW= 2
IBRION = 1
ISIF = 3
EDIFFG = -0.001
MAGMOM= 6*0 4 4 4 -4 -4 -4 18*0

LDAU= .TRUE.
LDAUTYPE= 2
LDAUL= -1  2  -1
LDAUU=  0  4  0
LDAUJ=  0  1  0
LMAXMIX= 4

AMIX=0.1
BMIX=0.0001
AMIX_MAG=0.4
BMIX_MAG=0.0001

DOS related parameters
ISMEAR = 0
SIGMA = 0.05

OUTPUT OPTION
LWAVE =  True
LCHARG = True
3楼2013-08-17 11:18:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoluoj

金虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
NELM确定的是一个离子步中电子步迭代的最大步数,25设定的太小了
4楼2013-08-17 12:47:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whw19850730

铁杆木虫 (小有名气)

引用回帖:
4楼: Originally posted by xiaoluoj at 2013-08-17 12:47:07
NELM确定的是一个离子步中电子步迭代的最大步数,25设定的太小了

25步已经发现-0.30811E+06了  在增加迭代步数没有作用,我已经试过了,最后会越来越大。
5楼2013-08-18 02:28:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

oupengfei717

金虫 (小有名气)

【答案】应助回帖

楼主将IBRION改为2试试,感觉LuMnO3的晶体结构应该比较复杂吧?是多少个原子的晶胞?如果开始的晶体结构不是很确定的话,建议IBRION改为2。另外,你的EDIFFG,EDIFF和SIGMA值都设置的太小了,建议将这些参数的数值增大,比如EDIFFG=-0.02,EDIFF=1E-6,SIGMA=0.2应该已经足以,有些参数我在计算的时候也不常用,所以没有太多的经验。
如果这样还是不能收敛的话,第一步试着将ISIF改为7试试,只弛豫晶胞体积,然后第二步改为5,弛豫晶胞形状,最后再全部弛豫找到平衡结构。
6楼2013-08-18 11:10:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 whw19850730 的主题更新
信息提示
请填处理意见