计算化学公社

标题: 使用CP2K进行CDFT计算时怎么判断收敛趋势 [打印本页]

作者
Author: ljh123 时间: 2025-2-8 16:14
标题: 使用CP2K进行CDFT计算时怎么判断收敛趋势
(, 下载次数 Times of downloads: 9) 输入文件和输出文件已经上传。进行了很多尝试一直无法收敛到想要的电子态。CDFT的收敛似乎和常规的收敛过程不太一样，有几步好像已经到收敛限了但是还在计算，请教下各位该怎么看CDFT过程的收敛趋势呢

作者
Author: sobereva 时间: 2025-2-8 18:16
参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）讲CDFT部分的这页幻灯片：
(, 下载次数 Times of downloads: 1)

作者
Author: ljh123 时间: 2025-2-8 19:02

sobereva 发表于 2025-2-8 18:16
参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）讲CDFT部分的这页幻灯片：

sob老师，请问CDFT除了SCF收敛判断，是不是还有约束值收敛的判断呢（也就是约束的当前值和目标值的偏离）？如果有的话，约束值收敛限该如何调整呢？

作者
Author: sobereva 时间: 2025-2-8 20:02

ljh123 发表于 2025-2-8 19:02
sob老师，请问CDFT除了SCF收敛判断，是不是还有约束值收敛的判断呢（也就是约束的当前值和目标值的偏离） ...

(, 下载次数 Times of downloads: 1)

作者
Author: ljh123 时间: 2025-2-8 20:20
本帖最后由 ljh123 于 2025-2-8 20:27 编辑

sobereva 发表于 2025-2-8 20:02

sob老师，看您上一张图里显示的，lamda指的是strength of constrant吧？这一张图中EPS_SCF设置lamda收敛限为1E-3，但是上一张图收敛时，lamda的值为0.22736。我看了下官方手册，请问是因为lamda这个值还会进行迭代吗？但是我在输出文件里，似乎没看到lamda值得迭代过程，而是直接给出了lamda值。

作者
Author: sobereva 时间: 2025-2-9 00:07

ljh123 发表于 2025-2-8 20:20
sob老师，看您上一张图里显示的，lamda指的是strength of constrant吧？这一张图中EPS_SCF设置lamda收敛 ...

lamda是拉格朗日乘子，收敛限是指其变化量。原理在培训里有非常详细的讲授

作者
Author: ljh123 时间: 2025-2-10 18:20

sobereva 发表于 2025-2-9 00:07
lamda是拉格朗日乘子，收敛限是指其变化量。原理在培训里有非常详细的讲授

sob老师，我用CDFT的方法计算了电子转移积分，我对两个电子态的约束定义为一个约束片段自旋多重度，另一个约束片段电子数。但是里面没有rotation的那一项电子转移积分，请问这样是存在问题的吗。
  Atomic group:                                                             1
  Strength of constraint I:                                     0.000000000000
  Strength of constraint J:                                     0.127112735436
  Final value of constraint I:                                  41.972530249564
  Final value of constraint J:                                  -1.089728454315

  Overlap between states I and J:                               0.334300203940
  Charge transfer energy (J-I) (Hartree):                      0.065222477238

  Diabatic electronic coupling (Lowdin, mHartree):             26.102118274764
  ------------------------------------------------------------------------------
NO FORCE_EVAL section calculated the dipole

ENERGY| Total FORCE_EVAL ( MIXED ) energy [a.u.]:       -2725.829405303261410

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