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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 怎么把VASP计算得到的两个CHGCAR文件相加作为一个新的电荷密度文件
152/2返回列表上一页12
查看: 2068  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

恒静无言

金虫 (正式写手)
引用回帖:
10楼: Originally posted by neweroica at 2016-01-29 22:01:13
ISTART应该设置为0，如果是1则表示初始時读取WAVECAR...

改为ISTART=0,仍然是不对的，错误还是出现在CHGCAR上
WARNING: chargedensity file is incomplete
ERROR: charge density could not be read from file CHGCAR for ICHARG>10
赞一下 回复此楼
11楼2016-01-30 08:52:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

neweroica

木虫 (著名写手)

【答案】应助回帖

★ ★
franch: 金币+2, 谢谢回帖交流。。。。 2016-02-16 06:54:22
引用回帖:
11楼: Originally posted by 恒静无言 at 2016-01-30 08:52:47
改为ISTART=0,仍然是不对的，错误还是出现在CHGCAR上
WARNING: chargedensity file is incomplete
ERROR: charge density could not be read from file CHGCAR for ICHARG>10...

确保你的PREC和ENCUT在加和前和后来的非自洽计算都是一致的。另外非自洽计算的IBRION取-1
赞一下 回复此楼
Simpler, stronger, more tolerant
12楼2016-02-01 01:46:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

恒静无言

金虫 (正式写手)
引用回帖:
8楼: Originally posted by neweroica at 2016-01-29 20:59:53
用chgsum.pl脚本就行。见http://theory.cm.utexas.edu/vtsttools/scripts.html

虽然我的问题没有解决，但是非常感谢你
赞一下 回复此楼
13楼2016-02-26 10:59:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gao.f

新虫 (初入文坛)
问题解决了吗，我也遇到了同样的问题，刚开始也想用vesta，然后输出数据保存成.3ed的格式，但仔细看数据，发现和原始CHGCAR相差很大，我觉得可能vesta读取数据的时候已经除以FFT-grid的体积了（参见vasp guide 42 页），所以和原来的CHGCAR没有可比性

后来我用了一个软件叫vasptools，里面有个功能可以把CHGCAR里的数据加和，但是输出格式要自己规范一下，然后我把ICHARG设成11，进行DFT运算，发现两个铝原子能量居然是-108eV. 我觉得把电子密度直接相加肯定有不合理的地方
赞一下 回复此楼
14楼2016-03-17 06:02:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gao.f

新虫 (初入文坛)
引用回帖:
11楼: Originally posted by 恒静无言 at 2016-01-30 08:52:47
改为ISTART=0,仍然是不对的，错误还是出现在CHGCAR上
WARNING: chargedensity file is incomplete
ERROR: charge density could not be read from file CHGCAR for ICHARG>10...

我也遇到过这个提示
WARNING: chargedensity file is incomplete
ERROR: charge density could not be read from file CHGCAR for ICHARG>10
我的问题是忽略了CHGCAR最后还有一部分，叫aumentation occupacies， 逐个拷贝进去就好了，可以看看常规计算的CHGCAR的格式是怎样的
赞一下 回复此楼
15楼2016-03-17 06:06:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 恒静无言 的主题更新
152/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 2026国自然放榜时间 +3 布布和一二 2026-08-08 3/150 2026-08-08 18:15 by Lightingo
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:22 by oEVWOejN9taj
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:19 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:19 by oEVWOejN9taj
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:02 by oEVWOejN9taj
[论文投稿] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 15:07 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:42 by oEVWOejN9taj
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[基金申请] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 14:27 by oEVWOejN9taj
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 有没有H口的？有收到消息的吗？ +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 面上提前没消息，有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
信息提示
关闭
请填处理意见
关闭