| 查看: 1739 | 回复: 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 |
回复此楼» 猜你喜欢
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有72人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
-
招聘 金属平磨液，抛光液研发工程师
已经有3人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于electrostatic potential
已经有5人回复
- CFX 能量方程迭代发散
已经有5人回复
- 一个很菜的问题。。。管道流量或速度的计算？
已经有4人回复
- 在做fluent例子的时候遇到了问题，发现无法进行迭代计算，求教大家
已经有9人回复
- 请问~Fluent 中两相流如何设置各自的压力入口？另外，迭代后发现发散咋办？
已经有4人回复
oupengfei717
金虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 744.8
- 散金: 80
- 红花: 3
- 帖子: 258
- 在线: 144小时
- 虫号: 1435139
- 注册: 2011-10-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2013-08-17 10:31:52
whw19850730
铁杆木虫 (小有名气)
虫王
- 应助: 1 (幼儿园)
- 金币: 5877.3
- 红花: 1
- 帖子: 64
- 在线: 121.3小时
- 虫号: 1048864
- 注册: 2010-06-28
- 性别: GG
- 专业: 物理学I
|
我打算做几何优化的，第一步电子迭代就发散了，之所以迭代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
4楼2013-08-17 12:47:07
whw19850730
铁杆木虫 (小有名气)
虫王
- 应助: 1 (幼儿园)
- 金币: 5877.3
- 红花: 1
- 帖子: 64
- 在线: 121.3小时
- 虫号: 1048864
- 注册: 2010-06-28
- 性别: GG
- 专业: 物理学I
5楼2013-08-18 02:28:01
oupengfei717
金虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 744.8
- 散金: 80
- 红花: 3
- 帖子: 258
- 在线: 144小时
- 虫号: 1435139
- 注册: 2011-10-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2013-08-18 11:10:56
