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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助成功】BRMIX:very serious problems the old and the new charge density
51/1返回列表
查看: 9657  |  回复: 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的回帖

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

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的回帖
查看全部 23 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +16 爱看书的可乐 2026-08-26 18/900 2026-08-28 16:06 by winsaint
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +4 gdd2018 2026-08-28 9/450 2026-08-28 14:22 by jklily
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +12 晴天加油 2026-08-26 13/650 2026-08-28 11:32 by skyjudy
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+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
信息提示
关闭
请填处理意见
关闭