大佬您好，我同样是想做Mn3O4对水分子的吸附能，但是现在在模型优化阶段就出现了问题，我通过castep模块测试了ecut截断能和k-point，得到了final energy差值小于0.01ev，但是在模型优化的时候出现了问题，无法跑出曲线就失败了，报错是这样的。请问这是由于我在最初设置就有问题吗还是怎么呢？这是我的一些设置

我在算PMS吸附的时候遇到过类似的问题，不知道你解决了吗？我当时是同样吸附能算的很低，但是吸附能的计算比较有技巧。吸附的过程中底物和添加剂的结构扭曲所产生的能量变化同样会被计入吸附能中，所以我猜测是否是你的结构没有优化到能量最低点导致吸附能算的不够准确？因为看你的计算参数里使用的是med精度。一般来说，使用fine精度的收敛标准更符合文献里提到的参数。

楼主，请问一下是什么原因导致这样的呢，我也遇见了，吸附之后，吸附能是负数，但是吸附分子到催化剂的距离变长，吸附分子间的距离变短

看你的初始构型图中，slab模型的底部原子不太规整，请问有没有固定底部原子再进行驰豫呢

卡开发发 发表于 2023-2-17 16:20
掺杂就是做一定替换那样。VCA应该是混合原子性质。

谢谢老师的解答

fd7 发表于 2023-2-17 15:44
谢谢老师，老师，掺杂模型是指混合原子的那种吗

掺杂就是做一定替换那样。VCA应该是混合原子性质。

卡开发发 发表于 2023-2-17 12:33
如果你实际的体系是那种有分数占位的情况，我建议处理起来慎重一些，要么按照掺杂模型做，要么按照虚晶近 ...

谢谢老师，老师，掺杂模型是指混合原子的那种吗

fd7 发表于 2023-2-17 10:33
谢谢老师。老师，我用的是mp-674482。老师，自旋极化是要对Mn和Fe分别进行设置吗？应该使用高自旋还是低 ...

如果你实际的体系是那种有分数占位的情况，我建议处理起来慎重一些，要么按照掺杂模型做，要么按照虚晶近似（VCA）做，但两者各自有各自的问题。

卡开发发 发表于 2023-2-17 10:15
乍一眼看不出你的MnFeO3是哪种相，不过Materials Project搜索得到两个结构mp-759744和mp-674482应该都是铁 ...

谢谢老师。老师，我用的是mp-674482。老师，自旋极化是要对Mn和Fe分别进行设置吗？应该使用高自旋还是低自旋呢？（老师，其实我一开始用XRD推断出是Mn1.16Fe0.84O3。根据Findit进行模型构建，我发现是混合原子，我也进行了相关的计算，计算方法跟图3一样，能够得到O-O键长变大，吸附距离变短，吸附能为-5左右。后来看别人说最好不要进行混合原子的计算，麻烦多。所以找到与XRD也能符合的MnFeO3。其实MnFeO3按照Findit构建出来也是混合原子，然后我在Material Project上找到了不是混合原子的。不知道这样做是不是合理）。辛苦老师解答。

乍一眼看不出你的MnFeO3是哪种相，不过Materials Project搜索得到两个结构mp-759744和mp-674482应该都是铁磁性的，所以计算时应该考虑自旋极化。