计算化学公社

标题: 使用CP2K优化亚铁磁尖晶石晶胞结构时报错怎么解决？ [打印本页]

作者
Author: zhangtp 时间: 2024-9-30 10:34
标题: 使用CP2K优化亚铁磁尖晶石晶胞结构时报错怎么解决？
我使用使用ms建模Fe2CoO4尖晶石晶胞，导入到multiwfn编辑cp2k输入文件，进行结构优化，输入指令是mpirun -np 4 cp2k.popt Fe2CoO4.inp |tee test.out，出现如下图报错，但在计算sob老师提供的例子时COF Cu SiC的计算时未出现报错，请问这个问题要怎么排查和解决呢？下面会把输入inp文件和testout发到下面。

作者
Author: HAC 时间: 2024-9-30 12:22
看out文件的提示，应该让你改自旋多重度，这个sob老师讲过，你可以看看http://bbs.keinsci.com/thread-4836-1-1.html

作者
Author: zhangtp 时间: 2024-9-30 13:49

HAC 发表于 2024-9-30 12:22
看out文件的提示，应该让你改自旋多重度，这个sob老师讲过，你可以看看http://bbs.keinsci.com/thread-4836 ...

好的，非常感谢，我先看一下帖子信息。

作者
Author: 1372133091 时间: 2024-9-30 17:50
本帖最后由 1372133091 于 2024-9-30 18:03 编辑

看了下你的输入文件，有一些问题你可能会遇到，
比如“&DFT”部分的设置，你可能需要添加 BASIS_SET_FILE_NAME BASIS_MOLOPT_UCL 来解决报错问题（我测试了一下并不需要，但是以防万一说一下），相关解答在此贴中有简要概述，可以看一下（【新提醒】请教关于CP2K找不到势文件的解决方(法 - 第一性原理 (First Principle) - 计算化学公社 (keinsci.com)），此外，这个网站也可以帮助你选择合适的基组和势文件。（cp2k-basis • Get basis set and potential (pierrebeaujean.net)）
关于图中LSD的报错，建议去Materials Project（MP）网站看一下Fe和Co对应的磁矩（如图1），然后在Multiwfn中得到输出文件的过程中进行如下操作修改你的自旋多重度：
在主菜单中输入cp2k后，输入-9，你可以观察到 2 Set spin multiplicity, current: 2 此时说明自旋多重度默认为2，但是不建议去手动修改自旋多重度，原因下面会说。
然后输入9，定义Fe和Co的初始磁矩，根据我的查找：Fe为4.144，Co为2.799。如果你也是在MP上下载的这个晶胞，其磁矩你直接套用即可。
输入完毕后你会观察到2 Set spin multiplicity, current: 89 ，是因为Multiwfn自动帮你根据你所输入的磁矩计算了对应的自旋多重度。
连贯起来
cp2k
-9
9
*Fe
4.144
*Co
2.799
q

此外，如果你输入了这些磁矩，得到的自旋多重度仍然报错，那么最好是调整磁矩，使得自选多重度不报错，例如从奇数(89)到偶数(88)，因为很多时候可能单纯只是奇偶数的问题（比如其实你设置MULTIPLICITY 1 也能运行不报错，见图2），至于其深层原因，就需要做一些功夫去理解了。
以上仅为我自己解决自旋多重度过程中所总结的一些经验，如有出错还望见谅。
最后祝你科研顺利

作者
Author: zhangtp 时间: 2024-9-30 22:49

1372133091 发表于 2024-9-30 17:50
看了下你的输入文件，有一些问题你可能会遇到，
比如“&DFT”部分的设置，你可能需要添加 BASIS_SET_FILE_ ...

真的非常非常感谢你，我跨专业初学第一性原理计算软件，学起来很吃力。还想请教你一下，比如我的结构是亚铁磁尖晶石结构，应该怎么在Multiwfn设置Fe的反磁呢

作者
Author: zhangtp 时间: 2024-9-30 23:34

1372133091 发表于 2024-9-30 17:50
看了下你的输入文件，有一些问题你可能会遇到，
比如“&DFT”部分的设置，你可能需要添加 BASIS_SET_FILE_ ...

抱歉我上个问题问的不清楚，我看文献里Fe的反磁是磁矩大小相等方向相反，但具体到结构里就不知道哪一个Fe的磁矩设为正，哪一个设为负，请问设定规则是什么呢？下面是我用ms建模设置的结果，你可以看一下设置的有什么问题吗？

作者
Author: sobereva 时间: 2024-10-1 03:12

zhangtp 发表于 2024-9-30 22:49
真的非常非常感谢你，我跨专业初学第一性原理计算软件，学起来很吃力。还想请教你一下，比如我的结构是亚 ...

十分建议通过北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）完整、系统学一遍，专门有磁性体系计算的一节，详细说了怎么设置初始磁矩并给了具体完整的计算例子。本来CP2K就是用起来相对复杂的程序，没第一性原理计算经验的人从零摸索不仅效率极其低下还容易各种踩坑、算出错误结果，犯了致命问题还浑然不知。像是当前居然设REL_CUTOFF 200，完全是胡来，计算效率奇低，令CP2K的multigrid策略起不到效果。

作者
Author: 1372133091 时间: 2024-10-1 12:18
赞成社长的话，参加北京科音的第一性原理计算培训班是及其有必要的，并且是通过及其正确的方法来开始你的第一性原理计算之路，此外，如果学校里面有开设类似于“材料物理计算”这类的课程的话，建议可以去旁听或者选课学一学，亦是受益匪浅。

作者
Author: zhangtp 时间: 2024-10-1 20:45

sobereva 发表于 2024-10-1 03:12
十分建议通过北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）完 ...

sob老师，不知道明年三月份的cp2k培训会有新内容更新吗，我现在在犹豫是买今年三月份的cp2k培训资料，还是等着明年的三月份去现场参加培训

作者
Author: sobereva 时间: 2024-10-1 23:35

zhangtp 发表于 2024-10-1 20:45
sob老师，不知道明年三月份的cp2k培训会有新内容更新吗，我现在在犹豫是买今年三月份的cp2k培训资料，还 ...

所有培训每届都会有更新，正如同手机每一代都有新改进
早学早受益。自行决定

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