计算化学公社

标题: 求助：VESTA是如何将CHGCAR中的charge density换算成cube文件格式的 [打印本页]

作者
Author: lijiaming111 时间: 2023-9-24 18:16
标题: 求助：VESTA是如何将CHGCAR中的charge density换算成cube文件格式的
本帖最后由 lijiaming111 于 2023-9-24 18:50 编辑

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

(, 下载次数 Times of downloads: 3)

(, 下载次数 Times of downloads: 2)

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

作者
Author: hebrewsnabla 时间: 2023-9-24 19:03
电子密度的单位是angstrom^-3或者bohr^-3，肯定不含hartree。

作者
Author: lijiaming111 时间: 2023-9-25 22:56
但是第一个值0.2883是按照henkelman组vtstscripts里的chg2cube.pl方法算的用chg2cube.pl输出的cube文件能对上第一个值所以我不太清楚VESTA是怎么换算出1.16256的

作者
Author: 978142355 时间: 2023-9-26 11:17

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

作者
Author: 王秀 时间: 2025-10-23 09:51
本帖最后由王秀于 2025-11-6 10:02 编辑

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)