24小时热门版块排行榜    

查看: 2276  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085601专硕,总分342求调剂,地区不限 +4 share_joy 2026-03-16 4/200 2026-03-17 21:32 by hmn_wj
[考研] 341求调剂 +5 捣蛋猪猪 2026-03-11 7/350 2026-03-17 19:09 by 捣蛋猪猪
[考研] 344求调剂 +4 knight344 2026-03-16 4/200 2026-03-17 17:27 by ruiyingmiao
[考研] 332求调剂 +6 Zz版 2026-03-13 6/300 2026-03-17 17:03 by ruiyingmiao
[考研] 26考研求调剂 +6 丶宏Sir 2026-03-13 6/300 2026-03-17 16:13 by 醉在风里
[考研] 求调剂,总分315,考的生物医药,一志愿湖南师范大学。调剂到任何专业都可以 +4 小丁想进步 2026-03-11 5/250 2026-03-17 16:05 by 外星文明
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考博] 26申博 +4 八旬速览 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] (081700)化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
信息提示
请填处理意见