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

不会扣篮的后卫

在使用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的数据？
还是我的输入文件有问题？
（附上修改后的输入文件和新的输出报错文件）

sobereva

建议先考虑直接支持的GFN1-xTB，省得折腾参数

卡开发发

CP2K就没内置Cu-Cu的Slater-Koster参数，事实上DFTB官网也没Cu-Fe等这些参数，要么自己拟合要么换方法。
需要注意的是xTB对磁性体系可能不友好，建议慎重测试。

不会扣篮的后卫

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

谢谢卢老师！

不会扣篮的后卫

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

好的，谢谢老师！

不会扣篮的后卫

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

请问老师，为什么需要Cu-Cu，Cu-Fe的Slater-Koster参数（如下图所示，在CuFeO2的结构中，并没有发现Cu-Cu bond，Cu-Fe bond)？
图片来自于CuFeO2文章，The effect of oxygen vacancy on magnetism of geometrically frustrated
triangular lattice CuFeO2: Ab initio study

卡开发发

不会扣篮的后卫 发表于 2022-2-18 02:12
请问老师，为什么需要Cu-Cu，Cu-Fe的Slater-Koster参数（如下图所示，在CuFeO2的结构中，并没有发现Cu-Cu ...

你给的图只是显示的时候根据距离和元素判断Cu-Cu和Cu-Fe没有形成化学键，例如从上往下第七个原子层的Cu之间光看距离也是有作用，实际上的计算在处理Cu和Fe对应的矩阵元还是需要相关的参数，判断能否形成化学键是通过这样的计算后做电子结构的分析得到的。

不会扣篮的后卫

卡开发发 发表于 2022-2-18 09:06
你给的图只是显示的时候根据距离和元素判断Cu-Cu和Cu-Fe没有形成化学键，例如从上往下第七个原子层的Cu之 ...

谢谢老师的解答。我有三个问题麻烦老师再解答一下，
1.根据老师之前的解答，是不是可以说第7个原子层的Cu之间，存在非成键的弱相互作用，需要用Cu-Cu的Slater-Koster参数文件来针对Cu之间的非成键的弱相互作用进行计算？

-------------------------------------------------------------------------------------------------------
还有我在DFTB官网提供的文件夹，来自于Compressed archive with unpacked SK-sets:
slako-unpacked\slako\matsci\matsci-0-3里找到了
Cu-Cu，Cu-O的Slater-Koster参数文件
在slako-unpacked\slako\trans3d\trans3d-0-1里找到了
Fe-Fe，Fe-O，O-Fe的Slater-Koster参数文件

2. 请问老师，虽然还是没有Cu-Fe的Slater-Koster参数文件，
但有没有可能在有Cu-Cu，Cu-O, Fe-Fe，Fe-O，O-FeSlater-Koster参数文件情况下，外加Cu和Fe距离较远的情况下，
可以在输入文件中使用DFTB方法来获得合理的输出文件结果？

--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
3.请问老师，为什么Fe-O，O-Fe的Slater-Koster参数文件会有差异？（感觉两个Slater-Koster参数文件里的数据不一样，但pair interaction between the atom types是一样的，对吗？）

（附上Fe-O，O-Fe的Slater-Koster参数文件）

卡开发发

不会扣篮的后卫 发表于 2022-2-18 10:28
谢谢老师的解答。我有三个问题麻烦老师再解答一下，
1.根据老师之前的解答，是不是可以说第7个原子层的C ...

1、具体是什么样的作用应该要看电子结构分析完的结果，不一定是弱的作用。
2、即便Cu和Fe距离比较远，这个部分的参数应该要拟合或者去构造，应该并不能直接按照0来处理，更何况你的距离也不够长。
3、我对这个文件格式定义不大清楚，这方面我帮不上什么忙，没有构造这个方法参数的经验，你可以翻找一下原文看看，或者是等有DFTB资深使用经验的人来给予一些指导。

不会扣篮的后卫

卡开发发 发表于 2022-2-18 12:37
1、具体是什么样的作用应该要看电子结构分析完的结果，不一定是弱的作用。
2、即便Cu和Fe距离比较远，这 ...

好的，谢谢老师！

不会扣篮的后卫

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

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

卡开发发

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

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

不会扣篮的后卫

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

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

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

Extended tight-binding quantum chemistry methods

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

卡开发发

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

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

对。

不会扣篮的后卫

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

谢谢老师！