| 查看: 1086 | 回复: 8 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】读CHGCAR里的数据
|
|||
|
虽然用vesta能很好的看电荷密度，但是不能得到几个相减的电荷密度差。网上有小软件做这个，但我的老是出问题，所以特别想自己拿最原始的数据，自己动手在origin做。 问题是chgcar里面的数据根本读不同，有没有哪位高人能指点下？ 不胜感激。 |
回复此楼» 猜你喜欢
博士延得我，科研能力直往上蹿
已经有7人回复
-
退学或坚持读
已经有27人回复
-
面上基金申报没有其他的参与者成吗
已经有5人回复
-
有70后还继续奋斗在职场上的吗？
已经有5人回复
-
遇见不省心的家人很难过
已经有22人回复
-
多组分精馏求助
已经有6人回复
2楼2009-12-04 13:24:41
3楼2009-12-04 17:50:41
4楼2009-12-04 17:51:21
|
能不能解释一下"WRITE(IU,FORM) (((C(NX,NY,NZ),NX=1,NGXC),NY=1,NGYZ),NZ=1,NGZC) "? 因为我看chgcar里面电荷密度部分每行有五个数据，按照（x,y,z,density）理解好像不对啊。 还有，上面的语句是Fortran写的吗？ 多谢了！ CHGCAR 如下： Fe-CO 1.00000000000000 0.000000 5.732800 0.000000 0.000000 0.000000 5.732800 17.166000 0.000000 0.000000 1 1 24 Direct 0.500000 0.464705 0.451881 0.500000 0.641670 0.501367 0.000000 0.000000 0.000000 0.000000 0.500000 0.000000 0.500000 0.000000 0.000000 ... 48 48 144 0.11710375410E+05 0.11082654420E+05 0.93465338201E+04 0.69185564327E+04 0.43845496481E+04 0.23036430807E+04 0.98892296509E+03 0.40174174331E+03 0.23369168341E+03 0.16858018791E+03 0.13671233690E+03 0.12337598318E+03 0.11990845078E+03 0.12335148597E+03 0.13665683923E+03 0.16846542258E+03 0.23345974810E+03 0.40125016805E+03 0.98777940755E+03 ... augmentation occupancies 25 138 -0.8339124E-02 0.1781598E-03 -0.2282241E-03 -0.1983454E-03 -0.3674109E-04 0.4917278E-03 -0.1912815E-01 0.8767631E-03 -0.7149171E-03 -0.8271628E-03 -0.2565947E-03 0.2019992E-02 0.9772558E-02 -0.5522550E-02 0.3923657E-01 ... augmentation occupancies 26 138 -0.9033073E-02 0.1112419E-03 0.3308886E-03 -0.2789434E-03 0.9789807E-04 0.5204789E-03 -0.1939311E-01 0.4184915E-03 0.1347132E-02 -0.1239935E-02 0.2499834E-03 0.2262646E-02 0.5939990E-02 0.3275988E-02 0.3513974E-01 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 ... |
5楼2009-12-05 09:24:11
6楼2009-12-05 09:30:26
7楼2009-12-05 10:50:00
8楼2009-12-05 12:34:51
9楼2009-12-05 13:30:22