求助，在CP2K中使用DFTB方法时，出现了报错

不会扣篮的后卫 · 发表于 Post on 2022-2-11 14:09:36

本帖最后由不会扣篮的后卫于 2022-2-12 01:02 编辑

在使用CP2K时，我在3个chemical unit CuFeO2的结构优化输入文件中加入了DFTB方法，想得到优化结构后的最低能量时，然后与同学采用的DFT+U方法进行数据比较。结果运行CP2K后出现了报错。
Failure in assigning KINDS <CU> and <CU> to a DFTB interaction pair!

请教老师，同学该怎么解决这个报错呢？
（附上输入文件和输出的报错文件）

更新于（02/11/2022）尝试加入原子间的interaction，于是在输入文件中加了以下内容：
      SK_FILE  Cu    Cu       cucu
      SK_FILE  Fe    Fe       fefe
      SK_FILE  O       O       oo
      SK_FILE  Cu    Fe       cufe
      SK_FILE  Fe    Cu       fecu
      SK_FILE  Cu    O       cuo
      SK_FILE  O       Cu       ocu
      SK_FILE  Fe    O       feo
      SK_FILE  O       Fe       ofe

但是又出现新的报错，The specified OLD file <DFTB/nonscc/cucu> cannot be opened. It does not exist.
Data directory path:/util/academic/cp2k/8.1-sse/cp2k/data

请教老师，同学是DFTB方法的data文件夹没更新以至于没有cucu的数据？
还是我的输入文件有问题？
（附上修改后的输入文件和新的输出报错文件）

yjb · 发表于 Post on 2022-5-21 19:33:21

你好，我的体系里面有FE
出现报错： Failure in assigning KINDS <FE> and <FE> to a DFTB interaction pair!

想请教一下，您是怎么解决的这个问题。

不会扣篮的后卫 · 发表于 Post on 2022-3-5 03:43:41

本帖最后由不会扣篮的后卫于 2022-3-5 03:44 编辑

sobereva 发表于 2022-3-4 18:50
你这个体系也不大，完全没必要用GFN-xTB
本来GFN-xTB是对于分子体系训练的参数，用于这种过渡金属氧化物 ...

谢谢卢老师！我尝试一下其他方法去。

不会扣篮的后卫 · 发表于 Post on 2022-3-5 03:41:08

Eudaimonia 发表于 2022-3-4 14:23
你使用OT方法时，普通SCF迭代到步数上限未收敛会进入outer SCF迭代，这个是outer SCF最大步数和收敛限设 ...

谢谢老师！

sobereva · 发表于 Post on 2022-3-4 18:50:05

不会扣篮的后卫发表于 2022-3-4 14:06
请教卢老师，我在用CP2K跑6个chemical unit CuFeO2组成的优化结构输入文件里，选择用GFN1-xTB方法[/backc ...

你这个体系也不大，完全没必要用GFN-xTB
本来GFN-xTB是对于分子体系训练的参数，用于这种过渡金属氧化物，又牵扯不同自旋的复杂情况，可靠度很有限。

Eudaimonia · 发表于 Post on 2022-3-4 14:23:39

不会扣篮的后卫发表于 2022-3-4 14:06
请教卢老师，我在用CP2K跑6个chemical unit CuFeO2组成的优化结构输入文件里，选择用GFN1-xTB方法[/backc ...

你使用OT方法时，普通SCF迭代到步数上限未收敛会进入outer SCF迭代，这个是outer SCF最大步数和收敛限设置的标准

不会扣篮的后卫 · 发表于 Post on 2022-3-4 14:06:06

sobereva 发表于 2022-2-24 11:27
不用管
GFN-xTB方法自身定义了基函数

请教卢老师，我在用CP2K跑6个chemical unit CuFeO2组成的优化结构输入文件里，选择用GFN1-xTB方法跑run_type:cell_opt，自旋多重度分别设为1（3个Fe1 alpha 电子数为5，3个Fe2 beta 电子数为5）
和31（6个Fe1 alpha 电子数为 5）.
在跑自旋多重度分别设为1的输入文件时，
得到理想的输出结果，优化结构成功收敛，

但在跑自旋多重度分别设为31的输入文件时，
优化结构没有成功收敛。

尤其是压力变的很不稳定。
第一个问题，请问卢老师为什么在高自旋的情况下，压力变化这么大？有什么输入文件的参数设置，可以使压力稳定和优化结构完成收敛吗？
（我尝试了选择run_type:Optimizing structure(cell is fixed)跑自旋多重度分别设为1（3个Fe1 alpha 电子数为5，3个Fe2 beta 电子数为5）
和31（6个Fe1 alpha 电子数为 5）的输入文件，这样可以避免产生关于压力
Pressure Deviation [bar] =    38288.6656706373
  Pressure Tolerance [bar] =    100.0000000000
  Conv. for  PRESSURE       =                NO    的信息报告，但是在run_type:Optimizing structure(cell is fixed)这种情况下，自旋多重度分别设为1和31都没有完成收敛。）

第二个问题，

请问卢老师，这个&OUTER_SCF参数设置有什么作用吗？
（附上相关的文件）

不会扣篮的后卫 · 发表于 Post on 2022-2-24 13:17:04

sobereva 发表于 2022-2-24 11:27
不用管
GFN-xTB方法自身定义了基函数

好的，谢谢卢老师！

sobereva · 发表于 Post on 2022-2-24 11:27:48

不会扣篮的后卫发表于 2022-2-24 07:08
卢老师，您好。我在使用GFN1-xTB方法计算结构的能量时（该结构由 6个chemical unit的CuFeO2组成），
出 ...

不用管
GFN-xTB方法自身定义了基函数

不会扣篮的后卫 · 发表于 Post on 2022-2-24 07:08:00

sobereva 发表于 2022-2-12 03:07
建议先考虑直接支持的GFN1-xTB，省得折腾参数

卢老师，您好。我在使用GFN1-xTB方法计算结构的能量时（该结构由 6个chemical unit的CuFeO2组成），
出现了以下的warning，
Cu和O原子的信息都被忽略了。
请问卢老师，这个warning该怎么解决呢？
（附上相关的输入和输出文件）

*** WARNING in qs_kind_types.F:2043 :: Information provided in the input ***
*** file regarding POTENTIAL for KIND <Cu> will be ignored!             ***

*** WARNING in qs_kind_types.F:2049 :: Information provided in the input ***
*** file regarding BASIS for KIND <Cu> will be ignored!                ***

*** WARNING in qs_kind_types.F:2043 :: Information provided in the input ***
*** file regarding POTENTIAL for KIND <O> will be ignored!             ***

*** WARNING in qs_kind_types.F:2049 :: Information provided in the input ***
*** file regarding BASIS for KIND <O> will be ignored!                ***

不会扣篮的后卫 · 发表于 Post on 2022-2-20 12:51:22

卡开发发发表于 2022-2-20 12:19
对。

谢谢老师！

卡开发发 · 发表于 Post on 2022-2-20 12:19:07

不会扣篮的后卫发表于 2022-2-20 09:52
请问老师，还记得文章题目的关键词吗？

我试着谷歌搜了一下关键词“grimme xtb high multiplicity to ...

对。

不会扣篮的后卫 · 发表于 Post on 2022-2-20 09:52:13

卡开发发发表于 2022-2-20 08:12
这个在Grimme本人的文章里面有提到，就是会出现高自旋算成低自旋的，低自旋的可能算出来就没磁性了。

请问老师，还记得文章题目的关键词吗？

我试着谷歌搜了一下关键词“grimme xtb high multiplicity to low multiplicity”出现了这篇文章

Extended tight-binding quantum chemistry methods

请问老师，是这篇文章吗？
（附上文章和感觉相关的文章内容）

卡开发发 · 发表于 Post on 2022-2-20 08:12:23

不会扣篮的后卫发表于 2022-2-20 07:44
请问老师，“xTB对磁性体系可能不友好”这个观点，可不可以简单讲一讲（我在网上没有找到类似的结论）？

这个在Grimme本人的文章里面有提到，就是会出现高自旋算成低自旋的，低自旋的可能算出来就没磁性了。

不会扣篮的后卫 · 发表于 Post on 2022-2-20 07:44:56

卡开发发发表于 2022-2-12 09:26
CP2K就没内置Cu-Cu的Slater-Koster参数，事实上DFTB官网也没Cu-Fe等这些参数，要么自己拟合要么换方法。
...

请问老师，“xTB对磁性体系可能不友好”这个观点，可不可以简单讲一讲（我在网上没有找到类似的结论）？

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 求助，在CP2K中使用DFTB方法时，出现了报错

浏览过的版块