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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » CHGCAR内容是什么？
221/3返回列表123下一页
查看: 5001  |  回复: 21
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

anuil

新虫 (初入文坛)

[求助] CHGCAR内容是什么？已有2人参与

如题，计算得出了CHGCAR后，发现其前面一部分就是CONTCAR部分，而且后面的这一部分应该是NGXC、NGYC、NGZC网格划分的电荷密度了，如下所示部分
240  224   80
0.18511728512E+04 0.18229386693E+04 0.17438297464E+04 0.16376639284E+04 0.15329692660E+04
0.14510883238E+04 0.13981166764E+04 0.13642473664E+04 0.13304094796E+04 0.12785076342E+04
0.12000131235E+04 0.10989454800E+04 0.98856182977E+03 0.88433220587E+03 0.79719326094E+03
0.73014407380E+03 0.67940384258E+03 0.64210109804E+03 0.61600722859E+03 0.59775670218E+03
0.58417421527E+03 0.57297401667E+03 0.56292156894E+03 0.55364714257E+03 0.54529644367E+03
0.53818252828E+03 0.53254965674E+03 0.52848444053E+03 0.52594736613E+03 0.52486889660E+03
0.52525262751E+03 0.52723900229E+03 0.53110065741E+03 0.53716455389E+03 0.54568458795E+03
。。。。。。
但是，最后面的部分是什么呢？
augmentation occupancies   1 138
  0.1273849E+01  0.4212516E-04  0.5766929E-04  0.3175189E-02 -0.3700677E-04
-0.4378747E-02  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
-0.1825212E+01 -0.5402954E-04 -0.5657891E-04 -0.2565814E-02  0.2314995E-04
  0.1623432E-01  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
  0.9393123E-04  0.7156361E-04  0.3750706E-02  0.4667899E-04 -0.1894632E-01
-0.6099121E-04 -0.4466665E-04 -0.1996873E-02 -0.9221010E-05  0.7958825E-03
  0.2372896E-02 -0.9057500E-04  0.1635742E-03  0.0000000E+00  0.0000000E+00
  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
-0.2363277E-03  0.1020300E-04 -0.3149818E-04  0.0000000E+00  0.0000000E+00
  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
  0.7994550E+00  0.2954665E-04  0.3424317E-04  0.7731032E-03 -0.9565602E-05
-0.7246242E-02  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
。。。。。
发现，我计算的体系是128个原子就会有augmentation occupancies   1 138，augmentation occupancies   2 138.。。。。。一直到augmentation occupancies   128 138，且后面带着的数值是什么意思啊？
回复此楼

» 收录本帖的淘贴专辑推荐

第一性原理计算-tg

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
一定不要浪费了好资源
1楼2013-12-16 20:08:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

heechao

至尊木虫 (著名写手)

【答案】应助回帖

★
anuil: 金币+1, ★★★很有帮助 2013-12-21 21:30:41
CHGCAR 里面前面是POSCAR的内容，后面是Ng=NGXF*NGYF*NGZF个数值为rowi*Vcell的格点。格点i的电荷数为mi=rowi*Vcell/Ng。rowi为i点的电荷密度，Vcell是晶胞总的体积。
赞一下(1人) 回复此楼
4楼2013-12-20 17:17:19
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

anuil

新虫 (初入文坛)
其实我主要是想做差分电荷密度的，结果我用VESTA相减做出来的图形，与我利用那个小插件chgdiff.pl相减得出的图形不一样，所以想着仔细分析一下CHGCAR到底是怎么一回事情啊
赞一下 回复此楼
一定不要浪费了好资源
2楼2013-12-17 10:26:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

JJ20

金虫 (小有名气)
我前段时间用过CHGCAR不过只是前面那一段，后面的我也不清楚，做差分电荷密度似乎只用到前面的一段，可能是你的小插件在读取有点问题，没有将自旋给区分好。同求后面解释````
赞一下 回复此楼
选择放弃只为保留尊严
3楼2013-12-20 15:20:39
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
anuil: 金币+5, ★★★很有帮助, 其实依然不太明白，augmentation electrons 的charge这部分？charge不应该给出一个数值就可以吗？ 2013-12-21 21:27:18
不確定什麼時候會用到這部份
但這指的應該是 augmentation electrons 的 charge
你的系統應該是有 128 顆原子
每一顆的 augmentation occupancy 是 138

詳情請看 Richard M. Martin 的書
Electronic Structure. Basic Theory and Practical Methods
或閱讀以下文獻 (http://cms.mpi.univie.ac.at/vasp ... ewtopic.php?4.1042)
Kresse , Hafner J.Phys Cond Mat. 6, 8245 (1994)
Kresse,Hafner , PRB49, 14251 (1994)
Furthmueller, Hafner,Kresse, PRB50,15606 (1994)
or Vanderbilt, PRB41, 7892 (1990) (US-PPs, in general)
赞一下 回复此楼
5楼2013-12-20 20:11:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anuil

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by JJ20 at 2013-12-20 15:20:39
我前段时间用过CHGCAR不过只是前面那一段，后面的我也不清楚，做差分电荷密度似乎只用到前面的一段，可能是你的小插件在读取有点问题，没有将自旋给区分好。同求后面解释````

嗯，对的，我查看了计算差分电荷密度结束后的文件，确实是只与前面的一部分有关，而且仔细检查之后，发现插件进行减法运算时还出现了错误，所以就是插件的问题了；谢谢提醒……
赞一下 回复此楼
一定不要浪费了好资源
6楼2013-12-21 21:30:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anuil

新虫 (初入文坛)
引用回帖:
4楼: Originally posted by heechao at 2013-12-20 17:17:19
CHGCAR 里面前面是POSCAR的内容，后面是Ng=NGXF*NGYF*NGZF个数值为rowi*Vcell的格点。格点i的电荷数为mi=rowi*Vcell/Ng。rowi为i点的电荷密度，Vcell是晶胞总的体积。

谢谢你对电荷密度这一部分的详细讲解，涨姿势了……
赞一下 回复此楼
一定不要浪费了好资源
7楼2013-12-21 21:32:02
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anuil

新虫 (初入文坛)
引用回帖:
5楼: Originally posted by jpchou at 2013-12-20 20:11:21
不確定什麼時候會用到這部份
但這指的應該是 augmentation electrons 的 charge
你的系統應該是有 128 顆原子
每一顆的 augmentation occupancy 是 138

詳情請看 Richard M. Martin 的書
Electronic Struct ...

对的，我的系统是128个原子的，虽然还是不太明白你所说的“每一顆的 augmentation occupancy 是 138”的含义，但是谢谢你给出的资料，看来还得多读些基础文献才可以啊，多谢哈
赞一下 回复此楼
一定不要浪费了好资源
8楼2013-12-21 21:34:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

JJ20

金虫 (小有名气)
引用回帖:
6楼: Originally posted by anuil at 2013-12-21 21:30:31
嗯，对的，我查看了计算差分电荷密度结束后的文件，确实是只与前面的一部分有关，而且仔细检查之后，发现插件进行减法运算时还出现了错误，所以就是插件的问题了；谢谢提醒……...

因为我前段时间也是做做电荷密度的事，所以知道那么点，新手一枚，同求教
赞一下 回复此楼
选择放弃只为保留尊严
9楼2013-12-22 09:00:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

问天麦饼

捐助贵宾 (小有名气)
引用回帖:
9楼: Originally posted by JJ20 at 2013-12-22 09:00:21
因为我前段时间也是做做电荷密度的事，所以知道那么点，新手一枚，同求教...

请问，CHGCAR那5列数据每一列的含义分别是什么？
赞一下 回复此楼
10楼2014-04-21 19:02:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 anuil 的主题更新
221/3返回列表123下一页
普通表情 龙 兔 虎 猫 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 硕博巨婴，也许才刚刚开始 +30 SNaiL1995 2024-06-12 74/3700 2024-06-16 17:08 by SNaiL1995
[找工作] 初始合伙人来啦！（生物试剂耗材标准品） +8 欢快的小科研人 2024-06-15 14/700 2024-06-16 16:00 by 欢快的小科研人
[基金申请] 希望今年自己国自然面上项目和老婆青年项目能中！ +7 恐龙爸爸 2024-06-14 7/350 2024-06-16 14:48 by redfish105
[教师之家] 请问事业编制和年薪制冲突吗？ +7 ZHONGWU_U 2024-06-14 7/350 2024-06-16 14:27 by 听风Ming
[基金申请] 博后基金，博管会会提前知道消息吗？ +4 yuyiang 2024-06-13 4/200 2024-06-16 11:40 by yangyuzhong4
[博后之家] 山东大学(青岛)“天然药物生物智造”课题组 招聘“博士后”(年薪20.4-55.6万元) +5 第二种态度 2024-06-11 8/400 2024-06-16 10:14 by 午睡未进行
[考博] 2025考博 +5 自强不息a?a 2024-06-15 7/350 2024-06-16 10:07 by wang333666
[考博] 34岁读博士晚吗 +26 emitdne 2024-06-13 26/1300 2024-06-16 07:16 by liyeqik
[论文投稿] 投稿被一个审稿人恶意评审了怎么样？ +5 1chen 2024-06-14 7/350 2024-06-15 23:15 by xy66xy
[基金申请] Nature 11日发文，中国著名学者们称造假迫不得已 +5 babu2015 2024-06-14 5/250 2024-06-15 20:40 by lc231001
[基金申请] BO4的YQ答辩通知发布了吗？ +6 博学笃行 2024-06-11 6/300 2024-06-15 16:04 by 悲催科研狗
[教师之家] 饶议：什么制度能保障大学普通教师不用为领导拎包，不用看领导脸色 +8 zju2000 2024-06-12 14/700 2024-06-15 13:59 by chemhua
[基金申请] 有没有机械的前辈分享一下评上海优都是什么成果啊 +7 wulala800 2024-06-10 7/350 2024-06-15 09:33 by 晓目崇
[食品] 食品博士导师 +6 小李醒yy 2024-06-11 9/450 2024-06-14 23:37 by 小李醒yy
[基金申请] 面上基金有一个（两个）C是不是就没戏了.... 5+3 zzzm116 2024-06-13 17/850 2024-06-14 22:26 by lzt8076
[考研] 物理化学一对一辅导 +3 林大diao 2024-06-12 5/250 2024-06-14 20:57 by 林大diao
[有机交流] 旋蒸请教 +6 办公室摸鱼 2024-06-09 9/450 2024-06-13 21:31 by wuxiaominggo
[基金申请] 连续两年医学口青年项目初审体会 +11 进击的荣耀 2024-06-09 18/900 2024-06-13 17:27 by 进击的荣耀
[硕博家园] 科研求助 +5 杲www 2024-06-12 6/300 2024-06-13 16:16 by 姓李名明
[论文投稿] with editor日期变更 +3 慎独的小花卷 2024-06-12 8/400 2024-06-13 11:00 by 慎独的小花卷
信息提示
关闭
请填处理意见
关闭