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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助成功】BRMIX:very serious problems the old and the new charge density
231/3返回列表123下一页
查看: 8622  |  回复: 22
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

uuv2010

荣誉版主 (职业作家)

优秀版主

[求助] 【求助成功】BRMIX:very serious problems the old and the new charge density

电子掺杂体系计算出现如下报错：
BRMIX: very serious problems  the old and the new charge density differ

INCAR相关参数如下：
EDIFF=1E-06
ISMEAR=0
LREAL=Auto
ISYM=0
VOSKOWN = 1
AMIX = 0.2
BMIX = 0.00001
AMIX_MAG = 0.8
BMIX_MAG = 0.00001

其中用O1和O2的编译版本计算均出现如上报错，请问如何解决？多谢！

[ Last edited by liliangfang on 2012-2-25 at 21:10 ]
回复此楼

» 收录本帖的淘帖专辑推荐

dft_vasp

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-02-13 08:35:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

mylifecg

木虫 (正式写手)
引用回帖:
19楼: Originally posted by 浴雪飘零_ at 2017-12-15 08:27:53
老哥我也出现这个问题了，你解决了吗...

没有，我只是换了台服务器计算，就没有报错。
赞一下(2人) 回复此楼
20楼2017-12-15 10:09:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

error9221

银虫 (初入文坛)
漫天飘雪: 回帖置顶 2022-02-17 15:25:44
或许可以参考这个：
https://755452800.github.io/2022 ... rge-density-differ/
在提交作业的脚本中加入：
export I_MPI_ADJUST_REDUCE=3
或者在INCAR中添加：
LSCALAPACK = .FALSE.
亲测都可以消除该错误
赞一下(1人) 回复此楼
22楼2022-02-08 01:19:25
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

likai0106

金虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
liliangfang(金币+2): 多谢交流 2012-02-13 09:37:53
uuv2010(金币+3): 多谢提示 2012-02-13 09:54:26
这个问题，看到过很多一直也受争议。
有说是编译器的问题，有的说是优化级别的问题。
楼主可以试试把KPOINTS中改成Gamma的形式进行计算看看能不能解决上述问题。
赞一下 回复此楼
2楼2012-02-13 09:10:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主
liliangfang(金币+1): 多谢提示 2012-02-13 10:17:41
引用回帖:
2楼: Originally posted by likai0106 at 2012-02-13 09:10:15:
这个问题，看到过很多一直也受争议。
有说是编译器的问题，有的说是优化级别的问题。
楼主可以试试把KPOINTS中改成Gamma的形式进行计算看看能不能解决上述问题。

k点是包含gamma点的
编译器是ifort11
优化采用了O1，O2的编译都是一样的报错
赞一下 回复此楼
3楼2012-02-13 09:55:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主
难道带电体系都是这个样子的报错？还是有别的原因？
赞一下 回复此楼
4楼2012-02-13 09:56:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

likai0106

金虫 (小有名气)
liliangfang(金币+1): 多谢交流 2012-02-13 10:17:57
引用回帖:
3楼: Originally posted by uuv2010 at 2012-02-13 09:55:49:
k点是包含gamma点的
编译器是ifort11
优化采用了O1，O2的编译都是一样的报错

楼主有没有试试看我说的把KPOINTS中改成gamma点展开的情况呢？
赞一下 回复此楼
5楼2012-02-13 10:15:11
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主
引用回帖:
5楼: Originally posted by likai0106 at 2012-02-13 10:15:11:
楼主有没有试试看我说的把KPOINTS中改成gamma点展开的情况呢？

谢谢您！
是的，KPOINTS如下：
K-Points
0
G
4 8 1
0  0  0
赞一下 回复此楼
6楼2012-02-13 10:31:58
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
uuv2010(金币+2): 多谢提示 2012-02-13 14:12:30
ben_ladeng(专家考核): 2012-02-13 18:48:04
ifort 11应该是可以直接用O3优化的～甚至据说O4都可以～实在不成，上makefile吧（帖有用代码部分就成～那些注释就不用了～）
赞一下 回复此楼
思想重于技巧，内涵重于表象
7楼2012-02-13 12:44:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主
引用回帖:
ifort 11应该是可以直接用O3优化的～甚至据说O4都可以～实在不成，上makefile吧（帖有用代码部分就成～那些注释就不用了～）

怀疑是载流子掺杂的体系都是有这样的报错，用一个原子做了一系列测试，包括ifort10，ifort11，O1,O2,O3,ISYM=0,ADDGID=TRUE,这个报错仍然存在，所以猜测是这样的：在没有update电子密度的时候，这个电子数是中性体系的电子数；当开始update电子密度的时候，计入载流子掺杂的电子数，所以新的电子密度和开始的电子密度不一致，导致这样的报错；正式由于这个原因，计算会照常进行而不至于停止。
现在的疑问是进行载流子掺杂的体系，计算时是怎样加入相反电荷的背景的？或者这段代码在哪个子文件中？

[ Last edited by uuv2010 on 2012-2-13 at 14:23 ]
赞一下 回复此楼
8楼2012-02-13 14:19:19
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

liliangfang(金币+2, 专家考核): 多谢指教 2012-02-13 16:06:18
uuv2010(金币+2): 多谢指导！ 2012-02-13 16:10:09
引用回帖:
8楼: Originally posted by uuv2010 at 2012-02-13 14:19:19:
怀疑是载流子掺杂的体系都是有这样的报错，用一个原子做了一系列测试，包括ifort10，ifort11，O1,O2,O3,ISYM=0,ADDGID=TRUE,这个报错仍然存在，所以猜测是这样的：在没有update电子密度的时候，这个电子数是中性 ...

记得那个错误一般实在rms(c)的值小于0.01，也就是0.1E-01时开始出现的，并且那个时候显示的rms(c)是有问题的～
不知道版大尝试的单原子是怎样的，我试了一下：
CODE:
# INCAR
  SYSTEM = Ti -single
    PREC = High
   ENCUT = 500
  IBRION = -1
     NSW = 0
    ALGO = Very Fast (RMM-DIIS)
    NELM = 60
  NELMIN = 2
   EDIFF = 1.0e-05
  EDIFFG = -0.02
VOSKOWN = 1
  NBLOCK = 1
   ISPIN = 1
  INIWAV = 1
  ISTART = 0
  ICHARG = 2
   LWAVE = .FALSE.
  LCHARG = .FALSE.
ADDGRID = .FALSE.
  ISMEAR = 1
   SIGMA = 0.2
   LREAL = .FALSE.
   RWIGS = 1.45

CODE:
#KPOINTS
Automatic mesh
0
Gamma
  5  5  5
  0.  0.  0.

CODE:
#POSCAR
Ti -single
    1.0
            10.00000000 0.00000000 0.00000000
            0.00000000 10.00000000 0.00000000
            0.00000000 0.00000000 10.00000000
  1
Direct
        0.50000000 0.50000000 0.50000000

计算的过程（在60步内没有收敛，不过盒子单原子一向容易如此）：
CODE:
       N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1     0.224307323808E+03    0.22431E+03   -0.82795E+02   120   0.302E+02
RMM:   2     0.157998692312E+03   -0.66309E+02   -0.38650E+02   120   0.783E+01
RMM:   3     0.180314187439E+02   -0.13997E+03   -0.58858E+02   120   0.105E+02
RMM:   4     0.282401550681E+01   -0.15207E+02   -0.17415E+02   120   0.689E+01
RMM:   5    -0.122423688476E+00   -0.29464E+01   -0.28167E+01   120   0.331E+01
RMM:   6    -0.416626367826E+00   -0.29420E+00   -0.25995E+00   120   0.157E+01
RMM:   7    -0.721549369634E+00   -0.30492E+00   -0.72052E-01   120   0.768E+00
RMM:   8    -0.127405531025E+01   -0.55251E+00   -0.17794E+00   120   0.553E+00
RMM:   9    -0.173989953448E+01   -0.46584E+00   -0.29661E+00   300   0.540E+00
RMM:  10    -0.194337486504E+01   -0.20348E+00   -0.10627E+00   346   0.170E+00
RMM:  11    -0.195320832601E+01   -0.98335E-02   -0.15824E-01   377   0.738E-01
RMM:  12    -0.196251503403E+01   -0.93067E-02   -0.46917E-02   359   0.189E-01    0.414E+00
RMM:  13    -0.127331119427E+01    0.68920E+00   -0.26089E+00   272   0.105E+01    0.128E+00
RMM:  14    -0.129749053952E+01   -0.24179E-01   -0.10811E+00   320   0.452E+00    0.432E+00
RMM:  15    -0.153442614239E+01   -0.23694E+00   -0.18426E+00   275   0.100E+01    0.234E+00
RMM:  16    -0.136438956798E+01    0.17004E+00   -0.87094E-02   310   0.221E+00    0.166E+00
RMM:  17    -0.196578833937E+01   -0.60140E+00   -0.22319E+00   298   0.109E+01    0.110E+01
RMM:  18    -0.122712276835E+01    0.73867E+00   -0.18171E-01   310   0.374E+00    0.881E-02
RMM:  19    -0.130710771260E+01   -0.79985E-01   -0.27563E-01   293   0.351E+00    0.114E+00
RMM:  20    -0.128118685540E+01    0.25921E-01   -0.92875E-03   315   0.658E-01    0.963E-01
RMM:  21    -0.123965985128E+01    0.41527E-01   -0.66660E-02   311   0.186E+00    0.777E-01
RMM:  22    -0.123977968172E+01   -0.11983E-03   -0.46937E-03   331   0.266E-01    0.739E-01
RMM:  23    -0.123102956896E+01    0.87501E-02   -0.76543E-03   316   0.613E-01    0.383E-01
RMM:  24    -0.123664918908E+01   -0.56196E-02   -0.18193E-03   298   0.298E-01    0.644E-01
RMM:  25    -0.123857007215E+01   -0.19209E-02   -0.85204E-04   339   0.138E-01    0.678E-01
RMM:  26    -0.123632916050E+01    0.22409E-02   -0.78203E-04   324   0.139E-01    0.623E-01
RMM:  27    -0.123503597287E+01    0.12932E-02   -0.57164E-04   343   0.121E-01    0.576E-01
RMM:  28    -0.123067580224E+01    0.43602E-02   -0.30523E-03   315   0.361E-01    0.180E-01
RMM:  29    -0.123059089273E+01    0.84910E-04   -0.80113E-04   371   0.981E-02    0.112E-01
RMM:  30    -0.123063972421E+01   -0.48831E-04   -0.36072E-04   307   0.308E-02    0.120E-01
RMM:  31    -0.123079966063E+01   -0.15994E-03   -0.36917E-04   322   0.382E-02    0.167E-01
RMM:  32    -0.123074285808E+01    0.56803E-04   -0.24040E-04   292   0.265E-02    0.141E-01
RMM:  33    -0.123080161148E+01   -0.58753E-04   -0.67147E-04   283   0.186E-02    0.139E-01
RMM:  34    -0.123084324720E+01   -0.41636E-04   -0.29660E-04   263   0.308E-02    0.142E-01
RMM:  35    -0.123076521027E+01    0.78037E-04   -0.58403E-04   297   0.118E-01    0.159E-01
RMM:  36    -0.123068255251E+01    0.82658E-04   -0.35434E-04   332   0.949E-02    0.119E-01
RMM:  37    -0.123074702773E+01   -0.64475E-04   -0.26248E-04   314   0.313E-02    0.142E-01
RMM:  38    -0.123076951027E+01   -0.22483E-04   -0.21484E-04   285   0.220E-02    0.134E-01
RMM:  39    -0.123077544736E+01   -0.59371E-05   -0.30475E-04   272   0.161E-02    0.130E-01
RMM:  40    -0.123085133506E+01   -0.75888E-04   -0.89492E-04   218   0.133E-02    0.141E-01
RMM:  41    -0.123094526214E+01   -0.93927E-04   -0.28226E-04   261   0.701E-02    0.195E-01
RMM:  42    -0.123092380925E+01    0.21453E-04   -0.12559E-04   282   0.216E-02    0.187E-01
RMM:  43    -0.123080558274E+01    0.11823E-03   -0.24531E-04   278   0.353E-02    0.156E-01
RMM:  44    -0.123078735192E+01    0.18231E-04   -0.23950E-04   279   0.210E-02    0.157E-01
RMM:  45    -0.123081798046E+01   -0.30629E-04   -0.24921E-04   299   0.564E-02    0.199E-01
RMM:  46    -0.123071476809E+01    0.10321E-03   -0.19343E-04   310   0.403E-02    0.170E-01
RMM:  47    -0.123066984389E+01    0.44924E-04   -0.17490E-04   276   0.221E-02    0.164E-01
RMM:  48    -0.123068437793E+01   -0.14534E-04   -0.18366E-04   257   0.217E-02    0.174E-01
RMM:  49    -0.123069311390E+01   -0.87360E-05   -0.14480E-04   232   0.132E-02    0.175E-01
RMM:  50    -0.123070560334E+01   -0.12489E-04   -0.83441E-05   184   0.904E-03    0.177E-01
RMM:  51    -0.123075544272E+01   -0.49839E-04   -0.40682E-04   268   0.110E-01    0.221E-01
RMM:  52    -0.123072354163E+01    0.31901E-04   -0.36267E-04   271   0.110E-01    0.190E-01
RMM:  53    -0.123064857732E+01    0.74964E-04   -0.20526E-04   261   0.661E-02    0.180E-01
RMM:  54    -0.123067334688E+01   -0.24770E-04   -0.15396E-04   291   0.381E-02    0.175E-01
RMM:  55    -0.123066981523E+01    0.35316E-05   -0.12629E-04   251   0.185E-02    0.177E-01
RMM:  56    -0.123068328463E+01   -0.13469E-04   -0.16984E-04   236   0.140E-02    0.178E-01
RMM:  57    -0.123070436166E+01   -0.21077E-04   -0.14611E-04   227   0.194E-02    0.183E-01
RMM:  58    -0.123075063842E+01   -0.46277E-04   -0.20489E-04   261   0.557E-02    0.191E-01
RMM:  59    -0.123072464068E+01    0.25998E-04   -0.13786E-04   293   0.340E-02    0.185E-01
RMM:  60    -0.123074878365E+01   -0.24143E-04   -0.19776E-04   228   0.983E-03
   1 F= -.12307488E+01 E0= -.11788354E+01  d E =-.155740E+00

顺便附上使用的赝势文件和OUTCAR于附件中～
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : test-single-Ti.tar.gz
    • 2012-02-13 15:12:25, 88.36 K
思想重于技巧，内涵重于表象
9楼2012-02-13 15:12:28
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主
是的，在这些文件中没有
赞一下 回复此楼
10楼2012-02-13 16:01:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 uuv2010 的主题更新
231/3返回列表123下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭