求助：VESTA是如何将CHGCAR中的charge density换算成cube文件格式的

lijiaming111

CHGCAR文件部分内容：
O atom in a box                        
   1.0000000000000000     
     7.500000    0.000000    0.000000
     0.000000    8.000000    0.000000
     0.000000    0.000000    8.900000
   O
     1
Direct
  0.000000  0.000000  0.000000

   80   80   96
0.41894056779E+04 0.40241460934E+04 0.35604175570E+04 0.28913321847E+04 0.21519258957E+04
0.14755073267E+04 0.94625540123E+03 0.58147085239E+03 0.34995433362E+03 0.20891981396E+03
0.12631211025E+03 0.77852246539E+02 0.48190146336E+02 0.29379850041E+02 0.17875964654E+02

VESTA转换后的cube文件部分内容：

O atom in a box                        
Created by VESTA v3.5.8
    1    0.000000    0.000000    0.000000
   80    0.177162    0.000000    0.000000
   80    0.000000    0.188973    0.000000
   96    0.000000    0.000000    0.175193
    8    8.000000    0.000000    0.000000    0.000000
1.16256E+00  1.16796E+00  1.15349E+00  1.06132E+00  8.71425E-01  6.28309E-01
4.05333E-01  2.48268E-01  1.52835E-01  9.41869E-02  5.92514E-02  3.83458E-02
2.50756E-02  1.62092E-02  1.04392E-02  6.88702E-03  4.61872E-03  3.07339E-03

按照VESTA Manual中写的
”CHG stores lattice vectors, atomic coordinates, and total charge densities multiplied by the unit-cell volume, V. ... When these files are read in to visualize isosurfaces and sections, data values are divided by V in the unit of bohr^3. The unit of charge densities input by VESTA is, therefore, bohr^-3."CHGCAR文件中的electron density单位是eV*angstrom^3, 我的试着换算了一下CHGCAR中的electron density的第一个值0.41894056779E+04

file:///C:/Users/LIJIAM~1/AppData/Local/Temp/msohtmlclip1/01/clip_image002.png

这个值和VESTA的输出cube文件的 第一个值明显对不上，有老师知道我的换算中少算了什么值吗

hebrewsnabla

电子密度的单位是angstrom^-3或者bohr^-3，肯定不含hartree。

lijiaming111

但是第一个值0.2883是按照henkelman组vtstscripts里的chg2cube.pl方法算的 用chg2cube.pl输出的cube文件能对上第一个值 所以我不太清楚VESTA是怎么换算出1.16256的

978142355

lijiaming111 发表于 2023-9-25 22:56
但是第一个值0.2883是按照henkelman组vtstscripts里的chg2cube.pl方法算的 用chg2cube.pl输出的cube文件能 ...

2楼的提示你可以仔细看下。0.41894056779e4/7.5/8/8.9*(Bohr->Ang)**3≈1.163

王秀

想请教一下，转化为cube文件中的网格数后面的坐标分量是怎么转化的？