24小时热门版块排行榜    

查看: 1782  |  回复: 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的回帖

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的回帖
查看全部 6 个回答

oupengfei717

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
楼主这个是静态计算吗?怎么只有一步离子迭代计算呢?我猜想可能是该结构在初期全弛豫的时候没有能找到平衡结构?不知道楼主是否进行了弛豫得到了平衡结构?楼主的输入文件是否都检查过不存在问题呢?
2楼2013-08-17 10:31:52
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-13 7/350 2026-09-16 03:26 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-16 03:23 by DgNGHc3h5tPl
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-16 03:14 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-16 03:02 by DgNGHc3h5tPl
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-16 02:38 by DgNGHc3h5tPl
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 6/300 2026-09-16 02:02 by DgNGHc3h5tPl
[考博] 申博发邮件 +4 Lmengk 2026-09-15 4/200 2026-09-16 00:05 by 木心想继续深造
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:47 by zNcEhTcH7ihq
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:26 by zNcEhTcH7ihq
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:39 by CgyNCDVNhGVg
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 09:38 by VMDqgnlKimtZ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:34 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 4/200 2026-09-14 23:24 by vZfe6xYu34yj
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
信息提示
请填处理意见