CHGCAR内容是什么？
如题，计算得出了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，且后面带着的数值是什么意思啊？
返回小木虫查看更多
京公网安备 11010802022153号
其实我主要是想做差分电荷密度的，结果我用VESTA相减做出来的图形，与我利用那个小插件chgdiff.pl相减得出的图形不一样，所以想着仔细分析一下CHGCAR到底是怎么一回事情啊
我前段时间用过CHGCAR不过只是前面那一段，后面的我也不清楚，做差分电荷密度似乎只用到前面的一段，可能是你的小插件在读取有点问题，没有将自旋给区分好。同求后面解释````
CHGCAR 里面前面是POSCAR的内容，后面是Ng=NGXF*NGYF*NGZF个数值为rowi*Vcell的格点。格点i的电荷数为mi=rowi*Vcell/Ng。rowi为i点的电荷密度，Vcell是晶胞总的体积。
不確定什麼時候會用到這部份
但這指的應該是 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)，
嗯，对的，我查看了计算差分电荷密度结束后的文件，确实是只与前面的一部分有关，而且仔细检查之后，发现插件进行减法运算时还出现了错误，所以就是插件的问题了；谢谢提醒……
谢谢你对电荷密度这一部分的详细讲解，涨姿势了……
对的，我的系统是128个原子的，虽然还是不太明白你所说的“每一顆的 augmentation occupancy 是 138”的含义，但是谢谢你给出的资料，看来还得多读些基础文献才可以啊，多谢哈