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,且后面带着的数值是什么意思啊? 返回小木虫查看更多
因为我前段时间也是做做电荷密度的事,所以知道那么点,新手一枚,同求教
,
请问,CHGCAR那5列数据每一列的含义分别是什么?
那五列数据就是代表个点上的电荷密度的具体数值啊,总共有就是你取的的格点总数个值啊
那么第一个数字是(0,0,0)点的密度,第二个是(0,0,1)还是(1,0,0)点的密度呢?
WRITE(IU,FORM) (((C(NX,NY,NZ),NX=1,NGXC),NY=1,NGYZ),NZ=1,NGZC)
这个是说明书里面的,就是按照这个顺序写的
我就是不懂这句话到底是先z自增,还是x自增……
这个是先x递增再y再到z,你可以自己写个小的fortran程序验证一下