24小时热门版块排行榜    

查看: 5008  |  回复: 21
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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,且后面带着的数值是什么意思啊?
回复此楼
一定不要浪费了好资源
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

问天麦饼

捐助贵宾 (小有名气)

引用回帖:
14楼: Originally posted by JJ20 at 2014-04-24 08:42:27
WRITE(IU,FORM) (((C(NX,NY,NZ),NX=1,NGXC),NY=1,NGYZ),NZ=1,NGZC)
这个是说明书里面的,就是按照这个顺序写的

我就是不懂这句话到底是先z自增,还是x自增……
15楼2014-04-24 10:06:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 22 个回答

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的回帖

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是晶胞总的体积。
4楼2013-12-20 17:17:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然青年基金,1A4B能上会吗?青年和面上的上会标准是一样的吗? +6 今晚推荐22 2024-06-20 9/450 2024-06-20 23:03 by 今晚推荐22
[基金申请] 工材口青年基金上会可能性 +9 今晚推荐22 2024-06-19 11/550 2024-06-20 22:30 by andywei1028
[找工作] 药学硕士找不到工作,打算去做科研助理了 +11 pom戴墨镜 2024-06-14 24/1200 2024-06-20 21:59 by 海阔天空4924
[硕博家园] 关于硕博连读的一些疑问? +9 Lwenter 2024-06-14 11/550 2024-06-20 21:20 by rhhua11111
[论文投稿] ACS AMI 返回审稿意见,一个大修,两个据稿,编辑给的修改重投 +7 智商已更新 2024-06-19 7/350 2024-06-20 19:54 by kinlin13
[论文投稿] 求机械类四区sci推荐 5+4 迷茫小旷 2024-06-14 7/350 2024-06-20 17:42 by 不一样烟火12345
[论文投稿] 投稿求助 6+3 小短腿安琪拉 2024-06-19 4/200 2024-06-20 17:40 by 不一样烟火12345
[论文投稿] 投稿求助 +4 平凡的日子 2024-06-19 5/250 2024-06-20 16:24 by yueyueyue@
[教师之家] 每次骚扰女学生的都是院系领导,而不是普通教师,小编们要注意措辞正确 +11 zju2000 2024-06-15 13/650 2024-06-20 15:19 by ase123456
[硕博家园] 豫北虫友互识 +11 xuhongli903 2024-06-18 12/600 2024-06-20 13:32 by 1429524838
[考博] 2025年博士申请,可先做科研助理 +4 limit888 2024-06-18 7/350 2024-06-20 10:23 by 安塔瓦拉多
[考博] 有机化学迷茫学生 +4 佛系摸鱼5 2024-06-18 6/300 2024-06-20 10:05 by 295143924
[公派出国] 去英国的小伙伴儿都在哪儿租的房子呀? +7 65syn 2024-06-14 11/550 2024-06-19 10:19 by 65syn
[基金申请] 博后基金刷到的BUG,图片来的更直观 +15 carolloo 2024-06-17 16/800 2024-06-19 09:42 by msjy
[考博] 2025考博 +8 自强不息a?a 2024-06-15 13/650 2024-06-18 18:12 by 投必得科研顾问
[基金申请] 有人中过人文社科类的博后特助吗? +3 outsider1986 2024-06-16 5/250 2024-06-18 11:10 by 袁天未然
[有机交流] 求教最后一步问题出在哪 20+3 苏州浩恒 2024-06-15 4/200 2024-06-18 08:33 by like_2005
[硕博家园] 博士毕业高校和就业的相关问题 +7 SCITOPPP 2024-06-14 11/550 2024-06-18 07:51 by yinxing1995
[论文投稿] 投稿被一个审稿人恶意评审了怎么样? +5 1chen 2024-06-14 7/350 2024-06-15 23:15 by xy66xy
[基金申请] 博后基金,以往的结果点不开,怎么回事呢?最后一次机会了,两次都没中前面。 +7 kyukitu 2024-06-14 13/650 2024-06-15 06:46 by 我是王小帅
信息提示
请填处理意见