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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » CHGCAR内容是什么？
71/1返回列表
查看: 5018  |  回复: 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的回帖

JJ20

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

JJ20

金虫 (小有名气)
引用回帖:
10楼: Originally posted by 问天麦饼 at 2014-04-21 19:02:37
请问，CHGCAR那5列数据每一列的含义分别是什么？...

那五列数据就是代表个点上的电荷密度的具体数值啊，总共有就是你取的的格点总数个值啊
赞一下 回复此楼
选择放弃只为保留尊严
11楼2014-04-22 10:17:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

JJ20

金虫 (小有名气)
引用回帖:
12楼: Originally posted by 问天麦饼 at 2014-04-22 13:50:07
那么第一个数字是(0,0,0)点的密度，第二个是(0,0,1)还是(1,0,0)点的密度呢？...

这个说明书上有的，我也记不清了，应该是先z递增再y再x，是个三循环嵌套的

[ 发自小木虫客户端 ]
赞一下 回复此楼
选择放弃只为保留尊严
13楼2014-04-23 23:30:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

JJ20

金虫 (小有名气)
WRITE(IU,FORM) (((C(NX,NY,NZ),NX=1,NGXC),NY=1,NGYZ),NZ=1,NGZC)
这个是说明书里面的，就是按照这个顺序写的
赞一下 回复此楼
选择放弃只为保留尊严
14楼2014-04-24 08:42:27
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

JJ20

金虫 (小有名气)
引用回帖:
15楼: Originally posted by 问天麦饼 at 2014-04-24 10:06:16
我就是不懂这句话到底是先z自增，还是x自增……...

这个是先x递增再y再到z，你可以自己写个小的fortran程序验证一下
赞一下 回复此楼
选择放弃只为保留尊严
16楼2014-04-24 19:19:52
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 anuil 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 我在写关于多维空间的论文，希望能与大家交流，准备在science上发表 +3 lgf519 2024-06-24 4/200 2024-06-24 21:14 by JKO
[基金申请] 博后面上今天下午会公布吗？有无消息？ +18 hajkdfdf 2024-06-24 23/1150 2024-06-24 19:32 by wangdachui0
[硕博家园] 申博雅思求助 +4 刘嘉研究究 2024-06-24 4/200 2024-06-24 16:20 by 仙女科研
[教师之家] 试用期辞职 +12 ZHONGWU_U 2024-06-18 23/1150 2024-06-24 15:40 by fenda
[基金申请] 说博后基金7月出的真打电话了吗？ +12 antonysole 2024-06-24 14/700 2024-06-24 13:39 by sizhouyi
[基金申请] 国自然资助比率是不是要下降了？？ +8 今晚推荐22 2024-06-21 12/600 2024-06-24 11:15 by Pickfoot
[催化] 求助原位红外 +3 jfdhj 2024-06-21 7/350 2024-06-24 10:52 by qvhm2609
[基金申请] 博后面上和特助今天出吗？ +41 逗您玩 2024-06-21 77/3850 2024-06-24 09:47 by weichuny
[论文投稿] 论文提交二审还有三天就三个月了，连续问了编辑部几次 10+3 大王叫我来寻山� 2024-06-22 9/450 2024-06-24 08:50 by 大王叫我来寻山�
[教师之家] 有没有今年的影响因子？ +6 jurkat.1640 2024-06-22 8/400 2024-06-24 07:25 by HiRocLiN丶
[有机交流] 生成亚胺的反应怎么能进行完全 +3 1369836 2024-06-23 3/150 2024-06-23 18:44 by hwqMSE
[基金申请] 请问三个字的中文名字，在文中到底怎么表示？ 5+3 zhongyuan764 2024-06-22 6/300 2024-06-23 17:42 by 鱼翔浅底1
[基金申请] 工材口青年基金上会可能性 +10 今晚推荐22 2024-06-19 13/650 2024-06-23 15:11 by Pickfoot
[基金申请] 工材口青年基金大概什么样能上会？ +15 今晚推荐22 2024-06-20 21/1050 2024-06-22 23:04 by qbn0326
[基金申请] 江南大学到瑞士招聘，称取消非升即走，改预聘+长聘 +24 babu2015 2024-06-18 29/1450 2024-06-22 16:29 by 风今25
[有机交流] 间氨基苯甲酸甲酯合成 5+3 pengdatao 2024-06-19 6/300 2024-06-22 10:58 by zyp0009928
[基金申请] 教育部基金 +5 m1393 2024-06-21 5/250 2024-06-21 21:13 by odes
[论文投稿] 审稿 +8 香瓜木香 2024-06-19 9/450 2024-06-21 15:51 by season7865
[考博] 2025年博士申请，可先做科研助理 +5 limit888 2024-06-18 8/400 2024-06-21 12:45 by 半简体
[海外博后] 墨尔本大学博后offer要不要接 +3 kyxblmm 2024-06-18 3/150 2024-06-19 22:39 by blake1111
信息提示
关闭
请填处理意见
关闭