MnO团簇模型scf不收敛

cxwudi · 发表于 Post on 2023-2-22 09:52:57

gog 发表于 2023-2-22 09:20
电子结构解析，用vasp、castep、qe、siesta，任选一个就可以。
分子动力学，半导体和绝缘体，选择cp2k。导 ...

感谢大佬提示，我也正在忧虑这件事。如果用VASP是不是要把文件给转成VASP用的形式。还是拿最后收敛的结构在用VASP算一次OPT然后分析？

卡开发发 · 发表于 Post on 2023-2-22 09:54:24

cxwudi 发表于 2023-2-22 08:53
谢谢卡卡老师，已经算完超胞吸附小分子的OPT了。用对角化SCF不收敛，然后改为OT的DIIS收敛了。这个有影响 ...

1、还是得看体系是否有gap，没gap的情况OT慎用。
2、OT算法具体实施我了解有限。原则上说，就算同样都是有gap的体系，本身按道理不同的对角化不应该得到不同的结果，但实际有些程序会这样，具体原因我说不清。
3、我不知道你指的奇怪是哪方面的问题。
4、原则上振动分析是需要的，但是体系一旦比较大其实愿意做的人不那么多。

卡开发发 · 发表于 Post on 2023-2-22 10:07:48

gog 发表于 2023-2-22 09:20
电子结构解析，用vasp、castep、qe、siesta，任选一个就可以。
分子动力学，半导体和绝缘体，选择cp2k。导 ...

个人观点，siesta的速度其实也很快，如果做电子输运可以配合transiesta做，算法相对还算稳定，可惜赝势基组也不太行，功能也偏少，我自己用的还算比较多但很少给人推荐。castep除了能在ms快速入手，迄今为止算法和功能上可圈可点的地方不太多了。

cxwudi · 发表于 Post on 2023-2-22 10:33:25

卡开发发发表于 2023-2-22 09:54
1、还是得看体系是否有gap，没gap的情况OT慎用。
2、OT算法具体实施我了解有限。原则上说，就算同样都是 ...

谢谢卡卡老师不吝指教。又得问问了这
1.OT方法还收敛了，有可能也依然不合理吗？
2.GAP具体指的是什么（能带吗）？我有些小白
3.吸附结构优化完，吸附区域晶体GV判断不成键了，然后有O原子下沉超出盒子了。
楼下附上图

cxwudi · 发表于 Post on 2023-2-22 10:33:50

卡开发发 · 发表于 Post on 2023-2-22 10:40:50

cxwudi 发表于 2023-2-22 10:33
谢谢卡卡老师不吝指教。又得问问了这
1.OT方法还收敛了，有可能也依然不合理吗？
2.GAP具体指的是什么 ...

1、具体情况我不好说，暂时判断不了，因为其他参数我也不知道，比如磁矩精度是否合理。
2、指的你价带顶和导带底的能隙。
3、O下陷的结构我暂时也判断不了原因，可能是模型不合理（比方这种单层就是稳定不了），也可能是计算参数的问题，等等。

cxwudi · 发表于 Post on 2023-2-22 11:37:38

卡开发发发表于 2023-2-22 10:40
1、具体情况我不好说，暂时判断不了，因为其他参数我也不知道，比如磁矩精度是否合理。
2、指的你价带顶 ...

1.之前扫帖的时候，看您说过先猜一个磁矩，算的时候软件会自动调整到合适的磁矩。您说的意思是要算一个精确的磁矩出来么？我之前算的时候都是盲猜一个（）
2.OT算法开始迭代的时候会给一个energy_gap : 0.08000000。是这个吗。还是要单独算
3.我附上我的输入文件

cxwudi · 发表于 Post on 2023-2-22 11:40:19

opt.inp (3.81 KB, 下载次数 Times of downloads: 4)

卡开发发 · 发表于 Post on 2023-2-22 11:53:53

cxwudi 发表于 2023-2-22 11:40

1、这个猜测不是盲猜，mp-1238773在MP网站上给的参考是铁磁性，且晶胞总磁矩是5，你可以按照这个猜，不一定要精确，但不能差太多。
2、OT可以简单看这个。
3、我随手DMol3测试了6x6x1（图方便没做晶胞优化），没出现O原子下陷的问题，结构几乎不变。

卡开发发 · 发表于 Post on 2023-2-22 14:00:36

cxwudi 发表于 2023-2-22 10:33
谢谢卡卡老师不吝指教。又得问问了这
1.OT方法还收敛了，有可能也依然不合理吗？
2.GAP具体指的是什么 ...

我用DMol3随手算了一下6x6的，结构没啥显著变化。

cxwudi · 发表于 Post on 2023-2-22 14:26:10

卡开发发发表于 2023-2-22 14:00
我用DMol3随手算了一下6x6的，结构没啥显著变化。

啊这，是我设置的参数有问题吗？或者这个用OT可能不合适

卡开发发 · 发表于 Post on 2023-2-22 14:55:49

cxwudi 发表于 2023-2-22 14:26
啊这，是我设置的参数有问题吗？或者这个用OT可能不合适

我没算吸附的，你的截图是吸附的还是未吸附的结构？

cxwudi · 发表于 Post on 2023-2-22 15:08:35

卡开发发发表于 2023-2-22 14:55
我没算吸附的，你的截图是吸附的还是未吸附的结构？

吸附的。我算未吸附的也基本没变化啊。吸附后就有个O下沉了。未吸附的对角化就可以收敛。

卡开发发 · 发表于 Post on 2023-2-22 17:24:58

cxwudi 发表于 2023-2-22 15:08
吸附的。我算未吸附的也基本没变化啊。吸附后就有个O下沉了。未吸附的对角化就可以收敛。

也许初始结构放的太糟糕了，我稍微测试了一下=O吸附的可能性也许大一些。

cxwudi · 发表于 Post on 2023-2-22 21:09:21

本帖最后由 cxwudi 于 2023-2-22 21:16 编辑

卡开发发发表于 2023-2-22 17:24
也许初始结构放的太糟糕了，我稍微测试了一下=O吸附的可能性也许大一些。

谢谢卡卡老师。

就是我那个输入文件有没有诡异的地方，FM的体系我直接把自选多重都设置成了所有单电子的加和。也没有smearing。磁矩根据MP给的4.605给Mn设置了5，别的元素都没有设置。就设置成这样原胞晶胞都能收敛就一直这样设置了，但是还是有些慌，自己到处看学了一点CP2K之前一直都是在用传统量化程序，现在还没看明白哪块设置才合理（）

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

[新手求助] MnO团簇模型scf不收敛

浏览过的版块