计算化学公社

标题: 求助：CP2K使用k点出现Method not implemented for k-points问题 [打印本页]

作者
Author: XWH 时间: 2021-10-23 14:48
标题: 求助：CP2K使用k点出现Method not implemented for k-points问题
本帖最后由 XWH 于 2021-10-23 14:48 编辑

请问一下，CP2K运用K点结构优化一个晶胞（12个原子），出现这个Method not implemented for k-points问题要怎么解决？感谢！输入文件如下：

作者
Author: 卡开发发 时间: 2021-10-23 18:34

Method not implemented for k-points

这话你细品，这方法对多k点不支持。从你的输入上看，可能的原因时他们的DFT+U还没支持多k点。

作者
Author: XWH 时间: 2021-10-23 18:47

卡开发发发表于 2021-10-23 18:34
这话你细品，这方法对多k点不支持。从你的输入上看，可能的原因时他们的DFT+U还没支持多k点。

感谢，可以请问下对于这种过渡组金属元素，结构优化可以先不加U吗？后期计算加U可以吗？不用k点的话，后面切面，体系可能会有点大

作者
Author: 卡开发发 时间: 2021-10-23 18:52

XWH 发表于 2021-10-23 18:47
感谢，可以请问下对于这种过渡组金属元素，结构优化可以先不加U吗？后期计算加U可以吗？不用k点的话，后 ...

U会影响势能面，但是否是有利于计算结果，我不好说。

作者
Author: XWH 时间: 2021-10-23 18:55

卡开发发发表于 2021-10-23 18:52
U会影响势能面，但是否是有利于计算结果，我不好说。

感谢！

作者
Author: haisenburg 时间: 2021-10-27 21:50
我也遇到同样问题

作者
Author: XWH 时间: 2021-10-28 01:26

haisenburg 发表于 2021-10-27 21:50
我也遇到同样问题

CP2K中k点设置目前不适用DFT+U方法

作者
Author: sigma-jlu 时间: 2023-6-18 14:36

XWH 发表于 2021-10-28 01:26
CP2K中k点设置目前不适用DFT+U方法

http://bbs.keinsci.com/thread-27607-3-1.html
这个帖子中用了PLUS_U_METHOD MULLIKEN关键段，未见不适用或not implemented的报错。

作者
Author: sobereva 时间: 2023-6-18 15:20

sigma-jlu 发表于 2023-6-18 14:36
http://bbs.keinsci.com/thread-27607-3-1.html
这个帖子中用了PLUS_U_METHOD MULLIKEN关键段，未见不适 ...

目前最新的CP2K 2023.1仍不支持DFT+U考虑k点，直接就会报错
输入输出例子：
(, 下载次数 Times of downloads: 34)
(, 下载次数 Times of downloads: 15)

作者
Author: sigma-jlu 时间: 2023-6-18 16:25

sobereva 发表于 2023-6-18 15:20
目前最新的CP2K 2023.1仍不支持DFT+U考虑k点，直接就会报错
输入输出例子：

多谢卢老师提供补充输入输出文件，明白前面所说的DFT+U的含义，指的是Kind部分的&DFT_PLUS_U字段，我刚才将PLUS_U_METHOD MULLIKEN误以为跟&DFT_PLUS_U是一回事。
&DFT_PLUS_U字段跟KPOINTS互斥的问题，目前在cp2k的确未见解决。

作者
Author: sigma-jlu 时间: 2023-6-18 16:59

sobereva 发表于 2023-6-18 15:20
目前最新的CP2K 2023.1仍不支持DFT+U考虑k点，直接就会报错
输入输出例子：

在卢老师的例子文件中，删除&DFT_PLUS_U字段，并在points字段中增加symmtry控制，发现一个费解的问题：
k点设成2 2 2，symmetry无论T还是F，都能正常运行；但
k点设成3 3 5，symmetry设成F能够正常输出，但设成T就出现"CPASSERT failed"错误闪退。
不知何故？恳请明白的老师解答。

相关4轮测试的输入输出已附。

作者
Author: sobereva 时间: 2023-6-19 11:17

sigma-jlu 发表于 2023-6-18 16:59
在卢老师的例子文件中，删除&DFT_PLUS_U字段，并在points字段中增加symmtry控制，发现一个费解的问题：
...

在我来看CP2K这方面代码没完全搞好，至少没个正经的报错提示。
可以问开发者，或者如果有钻研精神，可以去按照提示信息右下角说的，去看源代码文件cryssym.F的361行前后的信息试图理解

作者
Author: sigma-jlu 时间: 2023-7-20 16:00
本帖最后由 sigma-jlu 于 2023-7-20 16:08 编辑

sobereva 发表于 2023-6-19 11:17
在我来看CP2K这方面代码没完全搞好，至少没个正经的报错提示。
可以问开发者，或者如果有钻研精神，可以 ...

开发者juerghutter坦言了: “Unfortunately, there are no new developments available. Two projects that have been started, but have not yet finished are related with an OT type optimization for k-points and HFX for k-points.”
具体来说，kpoints的若干问题本质上难搞：

1.1 OT for non-metals: using the same algorithms but a lot of programming needed
1.2 OT with smearing: must include subspace rotations, even more to program
2. Would need totally new approaches
3. Should be rather straightforward, only acts on density
4. Should be rather straightforward, updates Hamiltonian with a density dependent term
5. NEW k-point symmetrization
6. NEW Wannier functions and interface to Wannier90
7. NEW (ADMM)-HFX

ref： https://github.com/cp2k/cp2k/issues/225

不过，kpoints的自洽场迭代速度问题，开发者aussy正在通过qs_diis.F源程序的commit在改进：

https://github.com/cp2k/cp2k/pull/2721

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