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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】新旧电荷密度有差异？
91/1返回列表
查看: 2157  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

keitapsp

新虫 (初入文坛)

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

我用并行版的VASP算，4个node。跑的时候开始没事，后来每一步都有这个提示：
very serious problems the old and the new charge density differ
这个是怎么回事啊？有人遇到过类似的情形吗？
谢谢
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-05-27 17:12:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
gavinliu7390(金币+1):谢谢交流！ 2010-05-28 22:32:48
计算的过程中Fermi level确定可能有问题，试试调整ISMEAR参数以及SIGMA值，增加k点。
赞一下(2人) 回复此楼
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))
赞一下(2人) 回复此楼
一定要精通MaterialsStudios,VASP.
3楼2010-05-27 20:25:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

163.com

金虫 (著名写手)

实习版主
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
还是参数设置问题，INCAR
赞一下(1人) 回复此楼
阿什顿联
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
赞一下(1人) 回复此楼
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 ，可能会避免这个问题。
赞一下(2人) 回复此楼
一定要精通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了，不知道这种想法是否有误区？
赞一下(2人) 回复此楼
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 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭