计算化学公社

标题: 请问CP2K中分别输出价电子密度和原子实，原子实电荷为0 [打印本页]
作者
Author: 赵小壮    时间: 2023-11-1 13:20
标题: 请问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
复制代码


作者
Author: sobereva    时间: 2023-11-1 16:08
不用管那个，不影响你分析电子密度。非要说的话算是CP2K产生的cube文件不规矩



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3