24小时热门版块排行榜    

CyRhmU.jpeg
查看: 8619  |  回复: 22
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主

我用C原子做了一个测试,如果是中性体系,没有报错,当进行载流子掺杂时,则出现报错,没有找到合适的办法消除这个报错,所以我怀疑对于载流子掺杂体系可能都是这样的报错。详细参数如下所示:
CODE:
刷屏信息如下:

       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1     0.109380701932E+02    0.10938E+02   -0.22772E+03    16   0.489E+02
DAV:   2    -0.166606784476E+02   -0.27599E+02   -0.27599E+02    32   0.745E+01
DAV:   3    -0.173066078911E+02   -0.64593E+00   -0.64593E+00    16   0.166E+01
DAV:   4    -0.173119060752E+02   -0.52982E-02   -0.52982E-02    16   0.176E+00
DAV:   5    -0.173119752199E+02   -0.69145E-04   -0.69145E-04    16   0.193E-01BRMIX: very serious problems
the old and the new charge density differ
old charge density:     4.00000 new    8.00000
    0.156E+01
DAV:   6    -0.144552175833E+02    0.28568E+01   -0.23607E+02    16   0.466E+01    0.103E+01
DAV:   7    -0.109009210256E+02    0.35543E+01   -0.16980E+01    16   0.214E+01    0.474E+00
DAV:   8    -0.104594812883E+02    0.44144E+00   -0.46597E+00    16   0.949E+00    0.145E+00
DAV:   9    -0.104517588973E+02    0.77224E-02   -0.84618E-01    16   0.435E+00    0.613E-01
DAV:  10    -0.104461949914E+02    0.55639E-02   -0.32455E-02    16   0.109E+00    0.448E-01
DAV:  11    -0.104401271472E+02    0.60678E-02   -0.26793E-02    16   0.672E-01    0.168E-01
DAV:  12    -0.104399995611E+02    0.12759E-03   -0.89320E-03    16   0.419E-01    0.558E-02
DAV:  13    -0.104400102221E+02   -0.10661E-04   -0.83164E-04    24   0.190E-01
   1 F= -.10440010E+02 E0= -.10362879E+02  d E =-.154262E+00

CODE:
NELECT = 8
IALGO=38
NELM=60
NELMIN=4
EDIFF=1E-04
ISMEAR=0
LREAL=.FALSE.
ISYM=0
SIGMA=0.1
NSW=0
EDIFFG=-1E-02
LWAVE  =F
LCHARG = F

[/quote]
[quote]
K-Points
0
G
1 1 1
0  0  0

CODE:
OSZICAR
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1     0.109380701932E+02    0.10938E+02   -0.22772E+03    16   0.489E+02
DAV:   2    -0.166606784476E+02   -0.27599E+02   -0.27599E+02    32   0.745E+01
DAV:   3    -0.173066078911E+02   -0.64593E+00   -0.64593E+00    16   0.166E+01
DAV:   4    -0.173119060752E+02   -0.52982E-02   -0.52982E-02    16   0.176E+00
DAV:   5    -0.173119752199E+02   -0.69145E-04   -0.69145E-04    16   0.193E-01    0.156E+01
DAV:   6    -0.144552175833E+02    0.28568E+01   -0.23607E+02    16   0.466E+01    0.103E+01
DAV:   7    -0.109009210256E+02    0.35543E+01   -0.16980E+01    16   0.214E+01    0.474E+00
DAV:   8    -0.104594812883E+02    0.44144E+00   -0.46597E+00    16   0.949E+00    0.145E+00
DAV:   9    -0.104517588973E+02    0.77224E-02   -0.84618E-01    16   0.435E+00    0.613E-01
DAV:  10    -0.104461949914E+02    0.55639E-02   -0.32455E-02    16   0.109E+00    0.448E-01
DAV:  11    -0.104401271472E+02    0.60678E-02   -0.26793E-02    16   0.672E-01    0.168E-01
DAV:  12    -0.104399995611E+02    0.12759E-03   -0.89320E-03    16   0.419E-01    0.558E-02
DAV:  13    -0.104400102221E+02   -0.10661E-04   -0.83164E-04    24   0.190E-01
   1 F= -.10440010E+02 E0= -.10362879E+02  d E =-.154262E+00

[ Last edited by uuv2010 on 2012-2-13 at 18:54 ]
11楼2012-02-13 16:03:55
已阅   回复此楼   关注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的回帖
信息提示
请填处理意见