计算化学公社

标题: 构建Fe110晶面的slab，并放置有机分子，计算吸附能。Fe超胞开启smearing也难收敛 [打印本页]

作者
Author: 故1996 时间: 2026-3-18 22:33
标题: 构建Fe110晶面的slab，并放置有机分子，计算吸附能。Fe超胞开启smearing也难收敛
老师，你好，我是做金属表面缓蚀剂吸附的，金属主要就是研究的Fe，一般做吸附能计算
一般研究的都是Fe的110晶面，我构建了一个正交的110晶面，并在上面放置了有机分子，如图，开启了smearing，开了色散校正，没有开偶极校正。
发现难以收敛，然后我把有机分子去掉，直接对Fe进行优化，也难以收敛。（smearing）

1.Fe超胞无法收敛跟盒子有关吗，我并没有在Multiwfn设置周期性，是通过高斯view加的盒子信息，参考啦讲义中的Fe_bluk.inp这个输入文件，这是我设置的输入文件，自选多重度，不在该如何设置，我记得说smearing开启后，并不需要设置电子数问题，是吗？
2.是一定要设置MAGNETIZATION吗？3 可以吗？
3.超胞和单一的晶胞优化是设置不同吗，在自旋多重度上？
4.如果后续在添加有机分子，是否要考虑偶极校正问题，盒子尺寸高35，应该是不够大？

作者
Author: sobereva 时间: 2026-3-18 23:37
该说的这里都说了：
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

开smearing后自旋多重度设置无效。但是原子初猜磁矩影响收敛难易。北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）讲磁性体系计算的一页幻灯片里展示了表面和内部的Fe的自旋磁矩的显著差异，初猜自旋磁矩尽可能考虑这一点。CP2K算Fe slab就是不容易收敛，死活解决不了的话可以改用Quantum ESPRESSO等基于平面波的免费程序

是否考虑偶极校正首先取决于用的周期性。XY二维周期性没偶极校正的事。盒子设置问题仔细复习培训里 “能量的计算及相关问题” 那一节中的 “低维体系的计算”部分

作者
Author: 故1996 时间: 2026-3-19 16:10
本帖最后由故1996 于 2026-3-19 16:13 编辑

sobereva 发表于 2026-3-18 23:37
该说的这里都说了：
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com ...

老师你好，我懂啦，那要是用QE软件，那我后续的计算吸附能是不也同样要用这个软件就行计算呢？，这个是还需要完全从新上手是吧。
老师，如果我采用Fe2O3作为slab，这样做有机分子的吸附能计算，Fe2O3的slab能够达到收敛吗，我看讲义的理解是能顾收敛，因为还是主要想在cp2k上实现吸附能计算功能。

作者
Author: sobereva 时间: 2026-3-19 16:17

故1996 发表于 2026-3-19 16:10
老师你好，我懂啦，那要是用QE软件，那我后续的计算吸附能是不也同样要用这个软件就行计算呢？，这个是还 ...

是

实测为准

作者
Author: 故1996 时间: 2026-3-19 16:20

sobereva 发表于 2026-3-18 23:37
该说的这里都说了：
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com ...

老师你好，我尝试啦建立更少原子的晶面，设置了自旋磁矩，但是依然很难达到收敛，如果要是用QE，是后续的吸附能计算都要用这个软件吗，只跟讲义课程学习过cp2k，不知重新学习，上手难度如何。
老师，我还是希望一直使用cp2k进行计算，毕竟速度优异，那我要是把slab换成Fe2O3，是否可以在cp2k上达到收敛呢？

作者
Author: northern_mar 时间: 2026-3-20 09:33
我之前发现OT算这种体系会更好收敛，但是不知道合不合理，所以也抛砖引玉跟大家讨论一下。
你可以尝试OT+pob-DZVP-rev2+PBEsol/revPBE试试，应该会比对角化算法好收敛很多。

作者
Author: 故1996 时间: 2026-3-23 16:22

northern_mar 发表于 2026-3-20 09:33
我之前发现OT算这种体系会更好收敛，但是不知道合不合理，所以也抛砖引玉跟大家讨论一下。
你可以尝试OT+p ...

感谢，我这就试试，采用全电子可能是个好办法

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