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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于vasp CHGCAR获取体系电荷的计算问题
51/1返回列表
查看: 1901  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

refors

木虫 (小有名气)

不坚定者

[求助] 关于vasp CHGCAR获取体系电荷的计算问题 已有1人参与

从vasp的手册上可以知道，CHGCAR里面的数据是指 p(r)*Vcell     前面的P(r)是指电荷密度, 后面的Vcell是指体系的原胞体积——我看到一种方法是将CHGCAR内的电荷加和到一起后除以三个方向上网格的分割数目即可获得体系的价电子数，我做了一个12个C和16个N的体系，三个方向分割分别为96 96 280，体系体积为880.3埃的三次方，我的电荷文件是关于一条能带的partial charge（charge文件格式与CHGCAR相同），也是说体系的电荷应该是2个电子，其中charge部分的积分为512000e，但是按照这个计算方法我得到的体系电荷是0.2个电子，刚好差了一个数量级。请教大神错漏之处
回复此楼

» 猜你喜欢
磨心砺性
1楼 2019-02-27 17:07:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jmlv

金虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
refors: 金币+5, ★有帮助 2019-03-01 08:53:33
VASP-VTST工具包可自动算出体系中各个原子的电子数目，无需自己算。
赞一下 回复此楼
2楼2019-02-28 16:36:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

refors

木虫 (小有名气)

不坚定者
引用回帖:
2楼: Originally posted by jmlv at 2019-02-28 16:36:04
VASP-VTST工具包可自动算出体系中各个原子的电子数目，无需自己算。

可以说的稍微详细点么，我vtst官网上没有看到相关的功能，谢谢
赞一下 回复此楼
磨心砺性
3楼2019-02-28 18:48:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jmlv

金虫 (著名写手)
引用回帖:
3楼: Originally posted by refors at 2019-02-28 18:48:06
可以说的稍微详细点么，我vtst官网上没有看到相关的功能，谢谢...

VASP-VTST有个Bader小程序，下载网页：http://theory.cm.utexas.edu/henkelman/code/bader/
赞一下 回复此楼
4楼2019-03-01 09:54:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

refors

木虫 (小有名气)

不坚定者
引用回帖:
4楼: Originally posted by jmlv at 2019-03-01 09:54:35
VASP-VTST有个Bader小程序，下载网页：http://theory.cm.utexas.edu/henkelman/code/bader/...

bader我有用过，但是这次我需要积分的是partial charge的电荷，包含单个k点和整个能带两种情况，最后一步调用bader脚本的时候好像是有问题的，所以不太能用bader直接去算
赞一下 回复此楼
磨心砺性
5楼2019-03-01 17:57:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 refors 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 11:47 by xSrQNSSjcgoX
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 10:52 by xSrQNSSjcgoX
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 09:46 by xSrQNSSjcgoX
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 SX5mJmPWDfOO 2026-11-13 3/150 2026-11-15 08:29 by xSrQNSSjcgoX
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 SX5mJmPWDfOO 2026-11-13 3/150 2026-11-15 07:56 by xSrQNSSjcgoX
[基金申请] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 XJXMuRjD7cK7 2026-11-13 4/200 2026-11-15 07:12 by xSrQNSSjcgoX
[基金申请] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 XJXMuRjD7cK7 2026-11-13 3/150 2026-11-15 05:44 by xSrQNSSjcgoX
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 YHiw1gou1Eda 2026-11-10 7/350 2026-11-14 22:44 by Equinoxhua
[基金申请] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +6 C6CRn66RYYXS 2026-11-10 7/350 2026-11-13 10:55 by Equinoxhua
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 TgvkbZZl9tqB 2026-11-11 6/300 2026-11-12 12:42 by ECLTXnJnRSXf
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:31 by HajghKOMynmz
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 6/300 2026-11-12 11:22 by HajghKOMynmz
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 11:11 by HajghKOMynmz
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 C6CRn66RYYXS 2026-11-10 7/350 2026-11-12 10:16 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 C6CRn66RYYXS 2026-11-10 8/400 2026-11-12 10:16 by HajghKOMynmz
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 YHiw1gou1Eda 2026-11-10 6/300 2026-11-12 10:14 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +8 zSbvVX9JFl9J 2026-11-09 12/600 2026-11-12 09:21 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 08:15 by HajghKOMynmz
[基金申请] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 5/250 2026-11-12 08:15 by HajghKOMynmz
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 05:19 by 9fOTqj074tIg
信息提示
关闭
请填处理意见
关闭