求助计算ITO吸附能时SCF不收敛

BBH · 发表于 Post on 2023-3-11 23:21:40

本帖最后由 BBH 于 2023-4-23 14:42 编辑

各位大佬好，本人刚学MS，有较多不懂的地方，还望海涵，我做的是SO2在ITO（110）面吸附能的计算，结构图如下：

建立的真空层厚度为15A，SO2在Sn原子正上方2.4A左右，计算采取P1对称性下广义梯度近似GGA/PBE泛函，kpoint2*2*1、截断能571.4 eV 、SCF cycles100

单击RUN运行后，弹出如下对话框

计算文件里的SCF结果如下：

想请问一下为什么会出现这种情况，以及如何修改。

重新进行了一次几何优化，只进行了截断能的修改，改至650eV，输出的CASTEP文件如下：

优化后失败的结构如下：

-------------问题分界线--------------------------------
In2O3：SnO2=90%：10%，选取的In2O3的空间群号为206，晶格参数a=10.117A，取（110）面进行掺杂，a=8.898786A,b=8.898786，c=21.648757，In：Sn：O =22：2：36

进行几何优化，Ecut=700eV,K点选取为2*2*2
输入文件：

1输入.txt (1022 Bytes, 下载次数 Times of downloads: 4)
输出文件：

1输出.txt (33.65 KB, 下载次数 Times of downloads: 1)
在输出文件里面SCF还是突然就断掉了，我用的是超算进行计算（不知道从哪看超算的内存），cores提交的是32，想问一下这是为什么，以及如何处理。
第二次提交只改变了cores为48，出现了以下问题：

同样也是SCF突然停止
输入文件：

2输入.txt (1022 Bytes, 下载次数 Times of downloads: 5)
输出文件：

2输出.txt (324.68 KB, 下载次数 Times of downloads: 4)
我又提交了第三次，这次是迭代了几次突然提示失败，cores这次选择的是64，去超算上查询作业，该项作业已经停止运行了

补充材料：
优化前：

优化后：

参数设置：

BBH · 发表于 Post on 2023-4-23 14:32:04

本帖最后由 BBH 于 2023-4-23 17:10 编辑

DesignChemistry 发表于 2023-3-16 20:51
在我做过类似的主族金属氧化物体系（如In2O3）中，结构优化用收敛性测试得到的值就够用了，基本只微小地 ...

老师，您好，我最近又遇到了一些问题想请教一下您，我再算ITO吸附能计算，首先构建的晶面是In2O3的（100）面，选取的是In终止面，为了不破化化学计量比，所以在切面时top值为0.033，thickness值为1，并且对In原子和O原子均加了U值，分别为3eV和5eV，切出来的面包含8层原子，为极性表面，所以我进行了偶极修正（只是在calculation里勾选了self-consistent），固定了下面的四层原子，真空层厚度为20A，但是在优化前后原子的整体位置进行了移动，想问一下我这种情况正确吗，计算时core选择了32，跑了一天了也之跑出一个点，也能否麻烦您帮我看一下我的参数设置是否正确，同时也想问一下，对于这种极性表面，若不加偶极修正，对计算吸附能、能带结构、PDOS以及mulliken有影响吗？麻烦老师了。
（材料我补充在了一楼）

BBH · 发表于 Post on 2023-3-27 19:30:38

DesignChemistry 发表于 2023-3-26 22:11
不考虑温度影响的话，算出来的电子能量（E）直接减就行了，算电子能量无法考虑温度
需要考察温度影响的 ...

好的，谢谢老师

DesignChemistry · 发表于 Post on 2023-3-26 22:11:57

BBH 发表于 2023-3-26 10:13
老师，您好，我还有一个问题想请教一下，就是ITO吸附一般是在较高温度下进行的，在用castep计算吸附能时 ...

不考虑温度影响的话，算出来的电子能量（E）直接减就行了，算电子能量无法考虑温度
需要考察温度影响的话，建议AIMD

BBH · 发表于 Post on 2023-3-26 10:13:22

DesignChemistry 发表于 2023-3-21 16:15
自旋没必要开；+U与否看你后边需不需要进行电子结构的分析，以及多找些文献看看LDA/GGA级别描述的偏差有 ...

老师，您好，我还有一个问题想请教一下，就是ITO吸附一般是在较高温度下进行的，在用castep计算吸附能时，可以从哪里修改温度吗，我在网上搜索没有搜索到

BBH · 发表于 Post on 2023-3-25 17:38:20

DesignChemistry 发表于 2023-3-21 16:15
自旋没必要开；+U与否看你后边需不需要进行电子结构的分析，以及多找些文献看看LDA/GGA级别描述的偏差有 ...

好的，谢谢老师

DesignChemistry · 发表于 Post on 2023-3-21 16:15:15

BBH 发表于 2023-3-20 16:43
好的，谢谢老师，最近在忙别的事情，就没有怎么登这个网站，不好意思这个时候才看到消息，同时我也想问一 ...

自旋没必要开；+U与否看你后边需不需要进行电子结构的分析，以及多找些文献看看LDA/GGA级别描述的偏差有多大。我当时的体系只算个吸附能和自由能台阶图，就没有加；如果后边需要DOS、能带之类的分析，应该从几何优化就开始加

BBH · 发表于 Post on 2023-3-20 16:43:44

reid 发表于 2023-3-16 23:10
初始结构中，为什么不让SO2躺下？上几天用cp2k算二水合铊，发现如果站着，容易跑散，也不容易收敛；让它们 ...

好的，谢谢老师

BBH · 发表于 Post on 2023-3-20 16:43:14

本帖最后由 BBH 于 2023-3-20 16:50 编辑

DesignChemistry 发表于 2023-3-16 20:51
在我做过类似的主族金属氧化物体系（如In2O3）中，结构优化用收敛性测试得到的值就够用了，基本只微小地 ...

好的，谢谢老师，最近在忙别的事情，就没有怎么登这个网站，不好意思这个时候才看到消息，同时我也想问一下老师一个问题，就是老师在对In2O3进行的几何优化过程中，有选择自旋或者LDA+U吗，因为我看的一些参考文献里面有用LDA+U去对SnO2进行几何优化，但是Sn和In一样都是主族元素，师兄告诉我对于主族元素不需要用这些去加长计算时间，想问一下，老师是怎么选择的

reid · 发表于 Post on 2023-3-16 23:10:57

初始结构中，为什么不让SO2躺下？上几天用cp2k算二水合铊，发现如果站着，容易跑散，也不容易收敛；让它们躺下则容易地多

DesignChemistry · 发表于 Post on 2023-3-16 20:51:59

BBH 发表于 2023-3-16 20:26
好的，谢谢老师，提高截断能是因为我去问组里面一个也用MS的师兄，师兄说提高截断能试一下，他之前就是靠 ...

在我做过类似的主族金属氧化物体系（如In2O3）中，结构优化用收敛性测试得到的值就够用了，基本只微小地影响数值，没遇到过影响SCF和几何优化收敛的，升高并不会有明显增益。

BBH · 发表于 Post on 2023-3-16 20:26:17

DesignChemistry 发表于 2023-3-16 16:58
抱歉这两天忙着干实验qwq
最后的这个新问题不像是计算本身问题，可能是超算问题，这个我没经验（一直都 ...

好的，谢谢老师，提高截断能是因为我去问组里面一个也用MS的师兄，师兄说提高截断能试一下，他之前就是靠提高截断能成功的。

DesignChemistry · 发表于 Post on 2023-3-16 16:58:34

BBH 发表于 2023-3-14 18:04
老师，您好，我这里出现了新的问题，资料已补充至一楼，能否劳烦您帮我看一眼该如何处理

抱歉这两天忙着干实验qwq
最后的这个新问题不像是计算本身问题，可能是超算问题，这个我没经验（一直都是自装机跑计算）
有个疑问是，你似乎一直执着于提高截断能呢？可以考虑k点的调整，有时用Gamma点会让受力计算粗糙、力震荡而不收敛（个人经验）

BBH · 发表于 Post on 2023-3-14 23:01:15

含光君发表于 2023-3-14 20:48
感觉是在超算使用上存在一些问题，这个就没法儿具体地判定哪个环节出错了。而且你给的输入输出文件都是同 ...

老师，我之前发现了输入输出文件的问题，后面已经修改过了，如果您有时间，还劳烦您再看一次

含光君 · 发表于 Post on 2023-3-14 20:48:46

BBH 发表于 2023-3-14 18:04
老师，您好，我这里出现了新的问题，资料已补充至一楼，能否劳烦您帮我看一眼该如何处理

感觉是在超算使用上存在一些问题，这个就没法儿具体地判定哪个环节出错了。而且你给的输入输出文件都是同一份输入文件。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CASTEP/Dmol3/MS] 求助计算ITO吸附能时SCF不收敛