cp2k进行CI-NEB计算时，中间结构无法进行scf计算

Shana

各位老师好，
我在使用cp2k进行CI-NEB计算时，遇到了几个问题。想请各位老师指点。
问题1：
我使用sob老师的插点工具sobNEB进行插点。我先是只插一个点，想了解一下过渡态的信息。
sobNEB.in的设置如下：
IS.xyz : 1
FS.xyz : 3
得到了中间点2.xyz。
然后设置BAND信息如下：
  &BAND
    NPROC_REP 32
    NUMBER_OF_REPLICA 3
    BAND_TYPE CI-NEB
    K_SPRING 0.05
    &CONVERGENCE_CONTROL # 跟结构优化类似
      MAX_FORCE 0.001
    &END CONVERGENCE_CONTROL
    &OPTIMIZE_BAND
      OPT_TYPE DIIS
      &DIIS
          MAX_STEPS 1000
      &END DIIS
    &END OPTIMIZE_BAND
    &REPLICA
      COORD_FILE_NAME  IS.xyz
    &END REPLICA
    &REPLICA
      COORD_FILE_NAME 2.xyz
    &END REPLICA
    &REPLICA
      COORD_FILE_NAME  FS.xyz
    &END REPLICA
    &PROGRAM_RUN_INFO
      INITIAL_CONFIGURATION_INFO
    &END
  &END BAND
96核提交计算后，生成了三个PROJECT-r.out文件，和三个PROJECT-BAND.out文件。奇怪的是，BAND1.out和BAND3.out都已经计算完成了，BAND2.out还没开始电子步。我等了一个晚上，也是一点信息都不输出，一直卡在这个界面。


问题2：
我取消中间插点后，cp2k能够正常进行计算。三个BAND.out都有信息输出。但是cp2k的插点却存在问题。我研究的是电荷转移的过程，想找到转移中的中间态。cp2k插完点后，只用到了初态结构。replica2和replica3与末态结构完全不同，过渡态的能差更是达到几个a.u.，而replica3得能量比末态高了近10个a.u.。cp2k插得点完全偏离初末态了。

现在得情况是，手动插点，cp2k不计算中间点，但是初末态是对的。cp2k自动插点，所有得replica都能正常计算，但是只有初态结构没变，最后一个replica已经偏离末态了。请各位老师指点一二。

Stardust0831

NEB是研究化学反应而非研究电子转移的。
研究电荷转移的多少可以用原子电荷；研究电荷转移速率可以了解一下Marcus theory、算算电子耦合之类的。

Shana

Stardust0831 发表于 2025-5-7 09:39
NEB是研究化学反应而非研究电子转移的。
研究电荷转移的多少可以用原子电荷；研究电荷转移速率可以了解一 ...

好的。谢谢。电子耦合有哪些理论或方法可以算呢？

Stardust0831

Shana 发表于 2025-5-7 15:41
好的。谢谢。电子耦合有哪些理论或方法可以算呢？

我觉得一般可以用CDFT或者Projection-Operator Diabatization来算。

Shana

Stardust0831 发表于 2025-5-7 16:00
我觉得一般可以用CDFT或者Projection-Operator Diabatization来算。

好的。谢谢