计算化学公社

标题: CP2k计算磁性异质结出现如下报错 [打印本页]

作者
Author: 机智的咸鱼 时间: 2024-2-3 20:32
标题: CP2k计算磁性异质结出现如下报错
本帖最后由机智的咸鱼于 2024-2-3 20:33 编辑

各位老师好，
我准备用CP2k做一个垂直的vdW异质结的结构弛豫。输入文件大部分内容由Multiwfn产生。但是每次运行计算都会出现ABORT：Invalid last k specified (A)，不知道是什么原因。

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

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

PS: Multiwfn真是良心，本来我这个结构丑陋，又要加U又要设反铁磁磁矩还要给Fe分类，本来自己挺怵的，结果很轻易的就完成了。

作者
Author: sobereva 时间: 2024-2-4 05:05
设置很不合理。MT要求非周期方向上盒子尺寸是相应方向有电子密度电子跨度的两倍以上，当前盒子明显太小，这明显不是Multiwfn用二维周期性结合MT时直接给你定义的盒子尺寸，肯定自己改过。这种slab很厚的情况应当当三维周期性算，免得用过大的真空区。

北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）的ppt里专门做了强调
(, 下载次数 Times of downloads: 30)

如果用三维周期性结合PSOLVER periodic之后还有问题，上传输出文件

作者
Author: 机智的咸鱼 时间: 2024-2-4 10:31
本帖最后由机智的咸鱼于 2024-2-4 10:32 编辑

sobereva 发表于 2024-2-4 05:05
设置很不合理。MT要求非周期方向上盒子尺寸是相应方向有电子密度电子跨度的两倍以上，当前盒子明显太小，这 ...

感谢社长的指正，Multiwfn直接产生的c方向晶格大约是111 Ang，我也用这个跑过，同样的报错。按照您说的三维周期性结合PSOLVER periodic跑了，还是一样的报错。输入输出都放在下面了。另外还附加了一个二维周期MT再将c加到120的结果，也是一样报错。
(, 下载次数 Times of downloads: 10)

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

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

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

作者
Author: sobereva 时间: 2024-2-5 04:09

机智的咸鱼发表于 2024-2-4 10:31
感谢社长的指正，Multiwfn直接产生的c方向晶格大约是111 Ang，我也用这个跑过，同样的报错。按照您说的三 ...

XYZ_PERIODIC_DIPOLE.inp第17行还是PERIODIC XY啊，改成PERIODIC XYZ。然后把smearing关了，我这里用2024.1能正常跑。不用smearing时还可以用OT，对这么大体系速度能快很多

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