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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP 电子迭代发散 诡异问题
61/1返回列表
查看: 1754  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-08-17 02:47:00
已阅   回复此楼   关注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 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 系统查不到 +9 董八千 2026-08-26 9/450 2026-08-26 12:30 by yanmofang
[基金申请] 国际合作可查了，中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 能否退出参与的面上项目解除限项 +22 koalala 2026-08-24 25/1250 2026-08-26 11:14 by dxcharlary
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 出来了 +8 trojank 2026-08-26 8/400 2026-08-26 09:00 by gxc
[基金申请] 范进中举一文的中心思想 +8 炎黄贵胄 2026-08-22 9/450 2026-08-26 07:47 by WASM
[基金申请] 2026国自然函评费到账 +21 羊腰板 2026-08-21 24/1200 2026-08-26 07:20 by finnigan
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
关闭
请填处理意见
关闭