请教一个castep计算计算SO在a-Fe2O3表面的吸附的问题

永远都要学会

新人小白，最近刚开始学castep，按照文献的方法在计算SO在a-Fe2O3表面的吸附。参数和模型都和文献设置的一样的，为什么软件跑不下去，有时只跑一步就失败了。文献和计算参数、模型都截图放在下面了，恳请各位老师指点，谢谢谢谢！
氧化铁脱除单质汞的机理和实验研究_刘亭 氧化铁脱除单质汞的机理和实验研究_刘亭

卡开发发

0、问题是自洽场（SCF）没收敛引起的
1、目测初始磁矩没设置，a-Fe2O3应该是反铁磁的。
2、动能截断也太低了，建议好好做个测试。
3、SCF如果还是难收敛试试调整density mixing的charge。

永远都要学会

卡开发发 发表于 2021-1-25 10:02
0、问题是自洽场（SCF）没收敛引起的
1、目测初始磁矩没设置，a-Fe2O3应该是反铁磁的。
2、动能截断也太 ...

谢谢卡神的回复！！
1、初始磁矩我是设置了，看文献设置的初始磁矩为5，不知道是不是这样设置：
2、截断能和k点我没有测试，但看的几篇文献都是测试350ev时能量最低：
3、请问截断能测试时候测试原胞和测试切面扩胞后的晶体得到的结论是一样的吗？
谢谢卡神老师！！

卡开发发

永远都要学会 发表于 2021-1-25 15:43
谢谢卡神的回复！！
1、初始磁矩我是设置了，看文献设置的初始磁矩为5，不知道是不是这样设置：
2、截 ...

1、如果按反铁磁结构设置，那应该没关系。
2、收敛性测试是看能量随动能截断增加是否收敛（比如10meV/atom），并非取能量最低的那个动能截断。
3、原则上表面会高一些，所以块体测试的情况最好预留大一些。
4、如果还不能收敛，调整density mixing是有效的办法。

永远都要学会

卡开发发 发表于 2021-1-25 16:49
1、如果按反铁磁结构设置，那应该没关系。
2、收敛性测试是看能量随动能截断增加是否收敛（比如10meV/at ...

谢谢卡神！祝工作顺利！！

永远都要学会

卡开发发 发表于 2021-1-25 16:49
1、如果按反铁磁结构设置，那应该没关系。
2、收敛性测试是看能量随动能截断增加是否收敛（比如10meV/at ...

卡开发发老师，density mixing的charge值已经调到0.1了还是无法收敛怎么办？一步都计算不下去

卡开发发

永远都要学会 发表于 2021-1-27 09:31
卡开发发老师，density mixing的charge值已经调到0.1了还是无法收敛怎么办？一步都计算不下去

块体和空表面可以收敛嘛？

永远都要学会

卡开发发 发表于 2021-1-27 10:43
块体和空表面可以收敛嘛？

可以，就是吸附不收敛

卡开发发

永远都要学会 发表于 2021-1-27 11:31
可以，就是吸附不收敛

吸附的结构和磁矩也许需要稍作调整，但是具体也得逐步做尝试。

永远都要学会

卡开发发 发表于 2021-1-27 13:20
吸附的结构和磁矩也许需要稍作调整，但是具体也得逐步做尝试。

卡开发发老师，您看我初始磁矩这样设置对吗？初始磁矩为5

卡开发发

永远都要学会 发表于 2021-1-27 14:35
卡开发发老师，您看我初始磁矩这样设置对吗？初始磁矩为5

按照块体优化的进行切面，应该问题不大。

Zhangchenyuan

你好楼主可以加你个联系方式吗？我也遇到了这种问题