24小时热门版块排行榜    

查看: 2094  |  回复: 8

keitapsp

新虫 (初入文坛)

[交流] 【求助】新旧电荷密度有差异? 已有6人参与

我用并行版的VASP算,4个node。跑的时候开始没事,后来每一步都有这个提示:
very serious problems the old and the new charge density differ
这个是怎么回事啊?有人遇到过类似的情形吗?
谢谢
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
gavinliu7390(金币+1):谢谢交流! 2010-05-28 22:32:48
计算的过程中Fermi level确定可能有问题,试试调整ISMEAR参数以及SIGMA值,增加k点。
2楼2010-05-27 19:01:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vasp001

金虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+1):谢谢 2010-05-27 21:52:29
引用回帖:
Originally posted by keitapsp at 2010-05-27 17:12:35:
我用并行版的VASP算,4个node。跑的时候开始没事,后来每一步都有这个提示:
very serious problems the old and the new charge density differ
这个是怎么回事啊?有人遇到过类似的情形吗?
谢谢

! Check whether charge density changed
     io_begin
      IF (ABS(TOTNEL-TOTOLD)>1E-5_q*ABS(TOTNEL) .AND. LWARN) WRITE(*,11) TOTOLD,TOTNEL
   11 FORMAT('BRMIX: very serious problems',/ &
             ' the old and the new charge density differ',/ &
             ' old charge density: ',F11.5,' new',F11.5)
TOTNEL=RHO0(GRIDC,CHTOT (1,1))
TOTOLD=RHO0(GRIDC,CHTOTL(1,1))
一定要精通MaterialsStudios,VASP.
3楼2010-05-27 20:25:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

163.com

金虫 (著名写手)

实习版主


小木虫(金币+0.5):给个红包,谢谢回帖交流
还是参数设置问题,INCAR
阿什顿联
4楼2010-05-27 21:59:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

keitapsp

新虫 (初入文坛)

谢谢楼上几位的回答,现在还在算,明天去看结果,再回来讨论。
5楼2010-05-27 22:53:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫于

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by vasp001 at 2010-05-27 20:25:29:

! Check whether charge density changed
     io_begin
      IF (ABS(TOTNEL-TOTOLD)>1E-5_q*ABS(TOTNEL) .AND. LWARN) WRITE(*,11) TOTOLD,TOTNEL
   11 FORMAT('BRMIX: very serious problems',/ &am ...

请问这是一个程序用来检测新旧电荷的么?我记得好像算过一个体系,在小k点的时候没有这个错误,增加k点反而出现了这个错误,这是怎么回事呢?我的 ISMEAR=-5 SIGMA=0.2
6楼2010-06-01 09:42:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vasp001

金虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+1):谢谢交流 2010-06-01 13:59:30
引用回帖:
Originally posted by 小木虫于 at 2010-06-01 09:42:39:

请问这是一个程序用来检测新旧电荷的么?我记得好像算过一个体系,在小k点的时候没有这个错误,增加k点反而出现了这个错误,这是怎么回事呢?我的 ISMEAR=-5 SIGMA=0.2

是用来检测新旧电荷差异的,是不是与算法有关系,你用四面体矫正的方法,换用ISMEAR =0 ,可能会避免这个问题。
一定要精通MaterialsStudios,VASP.
7楼2010-06-01 10:28:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫于

金虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
xiaohunhun(金币+1):xiexie 2010-06-01 22:58:23
引用回帖:
Originally posted by vasp001 at 2010-06-01 10:28:55:


是用来检测新旧电荷差异的,是不是与算法有关系,你用四面体矫正的方法,换用ISMEAR =0 ,可能会避免这个问题。

我想得到精确的能量,但ismear=0得到的能量应该不是很精确吧。我心里有一种阴影了,凡是想要精确能量的,我都有点不敢用ismear=0了,不知道这种想法是否有误区?
8楼2010-06-01 13:33:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

error9221

银虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
可以试试这个帖子:https://tpmk.github.io/2022/01/2 ... rge-density-differ/

在INCAR里增加
LSCALAPACK = .FALSE.

或者提交作业脚本里设置参数

或者调整核数为整数
9楼2024-12-24 16:29:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 keitapsp 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见