请问CP2K中分别输出价电子密度和原子实，原子实电荷为0

赵小壮

我的是 纯水体系，泛函revPBE，基组TZV2P-GTH，赝势GTH-PBE-qX。这里的"价电子"的定义是，每个水分子的两个成键电子对、O上的孤对电子共八个。
CP2K中可以在DFT模块中TOT_DENSITY_CUBE 输出全电子+核的电荷密度，纯粹的价电子密度应该是E_DENSITY_CUBE，
然而我的inp这样写之后：

1. &PRINT
   
2.   &E_DENSITY_CUBE
   
3.   FILENAME = density.dat
   
4.   STRIDE 1
   
5.   &END E_DENSITY_CUBE
   
6.      &END PRINT
   

复制代码

输出的的文件中，读取所有格点后进行运算后，确实是平均每个水分子分到-8|e|，
但是，按照sob老师在http://sobereva.com/125 中所言 第二列应该是是原子实的电荷， 但是我的文件中原子实的部分全都是0 请问这是怎么回事?

1. -Quickstep-
   
2. ELECTRON DENSITY
   
3.   192    0.000000    0.000000    0.000000
   
4.   150    0.156434    0.000000    0.000000
   
5.   150    0.000000    0.156434    0.000000
   
6.   150    0.000000    0.000000    0.156434
   
7.     8    0.000000   -6.270898   10.019415   10.675371
   
8.     1    0.000000   -6.949351   10.874765  -11.271695
   
9.     1    0.000000   -7.384103   10.539591    9.225827
   
10.     8    0.000000    5.299113    3.775028    0.705205
    
11.     1    0.000000    3.699462    3.701815    1.575235
    

复制代码

sobereva

不用管那个，不影响你分析电子密度。非要说的话算是CP2K产生的cube文件不规矩