计算化学公社

标题: cp2k进行CI-NEB计算时，中间结构无法进行scf计算 [打印本页]

作者
Author: Shana 时间: 2025-5-5 20:24
标题: cp2k进行CI-NEB计算时，中间结构无法进行scf计算
各位老师好，
我在使用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还没开始电子步。我等了一个晚上，也是一点信息都不输出，一直卡在这个界面。
(, 下载次数 Times of downloads: 2)

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

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

作者
Author: Stardust0831 时间: 2025-5-7 09:39
NEB是研究化学反应而非研究电子转移的。
研究电荷转移的多少可以用原子电荷；研究电荷转移速率可以了解一下Marcus theory、算算电子耦合之类的。

作者
Author: Shana 时间: 2025-5-7 15:41

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

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

作者
Author: Stardust0831 时间: 2025-5-7 16:00

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

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

作者
Author: Shana 时间: 2025-5-7 17:02

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

好的。谢谢

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