求助关于pms分子吸附在TiO2(101)表面的bader电荷问题

yueyueawyl

各位老师好，我最近在利用cp2k计算pms（HSO5-）吸附在在TiO2(101)表面遇到了关于bader电荷的问题。首先模型建立几何优化是slab+pms（charge是-1 自旋是1），然后拿到几何优化模型后，对slab+pms模型单点能计算（charge是-1 自旋是1），对slab模型（charge是0，自旋是1）和pms模型（charge是-1，自旋是1）此时没有做几何优化，直接从总的模型中删除不要的原子。然后对相应的输出文件的electron density 做bader分析，得到结果是pms向slab模型转移了0.24个电子，可是几乎所有的文章都是slab向pms转移电子，我就觉得很奇怪。下面是部分的输出参数文件。感谢大家提供的宝贵意见


      ! ------------------------------------------------------
      ! 电子密度 cube
      ! 做差分电荷时，system / surface / adsorbate 必须使用
      ! 相同晶胞、相同网格、相同 STRIDE、相同坐标参考。
      ! ------------------------------------------------------
      &E_DENSITY_CUBE
        FILENAME electron_density.cube
        STRIDE 1 1 1
        ADD_LAST NUMERIC
                 APPEND F
        &EACH
          QS_SCF 0
        &END EACH
      &END E_DENSITY_CUBE

reid

很正常，因为PMS的O与Ti配位，会提供2个电子给Ti共用，所以，你计算电子转移量，必然会显示PMS向基底转移电子

yueyueawyl

reid 发表于 2026-6-17 19:41
很正常，因为PMS的O与Ti配位，会提供2个电子给Ti共用，所以，你计算电子转移量，必然会显示PMS向基底转移电 ...

是的，这就是我所困惑的地方，然而大多数的文章都是显示pms是得电子的一方，除非他们一开始认为pms是不带电

vasp冲冲

yueyueawyl 发表于 2026-6-17 19:44
是的，这就是我所困惑的地方，然而大多数的文章都是显示pms是得电子的一方，除非他们一开始认为pms是不带 ...

你的判断可能是正确的的

yueyueawyl

vasp冲冲 发表于 2026-6-17 20:13
你的判断可能是正确的的

非常感谢您的回复，可是如果不带电，pms就应该是自由基（·HSO5而不是HSO5-),可是如果按照自由基来判断pms就不太会是一个稳定的结构可能会发生断键。

reid

用vasp可以不用设电荷。但是，算完大概率还是显示PMS向基底转移电子

yueyueawyl

reid 发表于 2026-6-17 22:09
用vasp可以不用设电荷。但是，算完大概率还是显示PMS向基底转移电子

非常感谢您的回复，我也看到有人这么设置过（bilibili上一个up主）。他的确没有设置电荷，但是谁知道呢