24小时热门版块排行榜    

查看: 5004  |  回复: 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,且后面带着的数值是什么意思啊?
回复此楼

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

第一性原理计算-tg

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

一定不要浪费了好资源
已阅   回复此楼   关注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的回帖
查看全部 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 饶议:什么制度能保障大学普通教师不用为领导拎包,不用看领导脸色 +9 zju2000 2024-06-12 15/750 2024-06-16 20:45 by 北溟鱼1318
[基金申请] 最新消息:2024国自然下载文件名变了 +6 dream200606 2024-06-16 9/450 2024-06-16 20:28 by fhx76
[教师之家] 关于2023的收入 +35 小龙虾2008 2024-06-10 36/1800 2024-06-16 19:30 by guoguo0119
[考博] 34岁读博士晚吗 +27 emitdne 2024-06-13 27/1350 2024-06-16 18:57 by yugege2009
[硕博家园] 硕博巨婴,也许才刚刚开始 +30 SNaiL1995 2024-06-12 74/3700 2024-06-16 17:08 by SNaiL1995
[博后之家] 山东大学(青岛)“天然药物生物智造”课题组 招聘“博士后”(年薪20.4-55.6万元) +5 第二种态度 2024-06-11 8/400 2024-06-16 10:14 by 午睡未进行
[基金申请] 为什么我的博后基金还在流动站审核中?不会是学院给我卡了吧? +14 王凯12 2024-06-13 26/1300 2024-06-15 15:22 by 好人与坏人
[基金申请] 博士后基金需要结题吗? +8 zhouchuck 2024-06-13 8/400 2024-06-14 17:27 by liuyupu132
[基金申请] E12面上申请 +4 汉风之遗 2024-06-13 4/200 2024-06-14 15:28 by 天外飞去来
[考博] 申博找导师 +4 疏影横斜水清浅3 2024-06-13 6/300 2024-06-14 14:31 by zxl_1105
[有机交流] ππ堆积会发生在有机溶剂中吗 5+3 zibuyu0420 2024-06-13 4/200 2024-06-14 14:17 by 小肉干
[基金申请] 国自然基金公布的时候基金号有吗 +8 潇洒怡惜 2024-06-13 11/550 2024-06-14 11:24 by JRfei
[基金申请] 化学会年会改成一年开一次了? +5 babu2015 2024-06-12 13/650 2024-06-14 10:46 by babu2015
[基金申请] 75批博后基金 +10 kyukitu 2024-06-13 13/650 2024-06-14 10:31 by kyukitu
[论文投稿] 文章proof要求使用机构的邮箱 5+3 不可不信缘 2024-06-11 11/550 2024-06-14 07:00 by 3001160025
[论文投稿] 投稿后发现其他作者的邮箱填错了该怎么办呀 10+4 在飞的猪 2024-06-13 6/300 2024-06-14 04:45 by 小虫子咔咔
[硕博家园] 机械研究生如何拿到年薪40+w +13 阿巴阿巴哦哦 2024-06-11 15/750 2024-06-13 15:40 by 113745685
[有机交流] 原料反应完了,怎么知道是产物还是中间体 +6 小胡在努力 2024-06-11 8/400 2024-06-13 13:33 by 091602
[论文投稿] 摩擦磨损论文投稿 +3 jmysan 2024-06-12 3/150 2024-06-13 08:36 by 莱茵润色
[论文投稿] water research状态咨询 5+3 Flyyawa 2024-06-10 6/300 2024-06-11 09:45 by bobvan
信息提示
请填处理意见