计算化学公社

标题: 关于Ni(111)晶面加了被吸附分子之后难收敛的问题 [打印本页]

作者
Author: ddddnight 时间: 2024-6-22 15:29
标题: 关于Ni(111)晶面加了被吸附分子之后难收敛的问题
做了Ni(111)晶面结构优化，很快就收敛了，但是加了溴化环氧树脂单元在晶面上面，按照老师课上讲的那些改进SCF收敛的方法去做，包括加给Ni加1.0的自选磁矩（本来是0.6，老实说要加的稍微大一点），都不能做到收敛，老师有什么好的建议吗？因为做了很多体系都这样，包括Fe（110）晶面也是，很容易收敛，加了被吸附分子之后就不行了。

作者
Author: 乐平 时间: 2024-6-22 22:21
本帖最后由乐平于 2024-6-22 22:27 编辑

不收敛大部分情况是结构的问题。
大概看了一下你的模型，由于你的吸附分子的长度和表面的尺寸比较接近，所以你需要更宽的表面来避免 x, y 方向的镜像相互作用。

你的真空层厚度设置得足够大了，说明你知道为了避免周期性边界条件带来的上下表面引起的相互作用。但是你却忽视了水平方向上的相互作用……

另外，我不清楚 MULTIPLICITY 145 是怎么得到的？

作者
Author: sobereva 时间: 2024-6-22 23:34
你对Ni(111)的计算本来就不对，磁性体系错当成了闭壳层计算。先把这个算对再说。不恰当设置自旋多重度或者初猜自旋磁矩/组态，光写个UKS毫无意义，alpha和beta电子的分布还是精确对称的

作者
Author: ddddnight 时间: 2024-6-22 23:48

sobereva 发表于 2024-6-22 23:34
你对Ni(111)的计算本来就不对，磁性体系错当成了闭壳层计算。先把这个算对再说。不恰当设置自旋多重度或者 ...

老师，您在课上说的Fe（bbc）原胞的时候是只给Fe原子加了数值为3的磁矩，这个体系是不是给Ni的磁矩设置为1可以吗，我加了溴化环氧树脂分子之后就是这么算的，但是就是不收敛。还需要别的设置吗

作者
Author: sobereva 时间: 2024-6-23 00:24

ddddnight 发表于 2024-6-22 23:48
老师，您在课上说的Fe（bbc）原胞的时候是只给Fe原子加了数值为3的磁矩，这个体系是不是给Ni的磁矩设置 ...

我强调的是，不要以为不收敛是因为带了被吸附物，考虑磁性的情况下单独算Ni(111)本来就可能已经SCF很难收敛了

作者
Author: ddddnight 时间: 2024-6-23 01:49

sobereva 发表于 2024-6-23 00:24
我强调的是，不要以为不收敛是因为带了被吸附物，考虑磁性的情况下单独算Ni(111)本来就可能已经SCF很难收 ...

好的老师，那如果是要普通结构优化，算吸附能这些的时候，是应该考虑磁性的吧？还是说只是单纯为了SCF收敛呢

作者
Author: sobereva 时间: 2024-6-23 07:00

ddddnight 发表于 2024-6-23 01:49
好的老师，那如果是要普通结构优化，算吸附能这些的时候，是应该考虑磁性的吧？还是说只是单纯为了SCF收 ...

只要是磁性体系，不管算什么都要考虑磁性问题，否则结果毫无意义

作者
Author: ddddnight 时间: 2024-6-23 10:10

sobereva 发表于 2024-6-23 07:00
只要是磁性体系，不管算什么都要考虑磁性问题，否则结果毫无意义

好的老师，谢谢您的回答，我再去尝试一下

作者
Author: ddddnight 时间: 2024-6-24 23:01

乐平发表于 2024-6-22 22:21
不收敛大部分情况是结构的问题。
大概看了一下你的模型，由于你的吸附分子的长度和表面的尺寸比较接近，所 ...

谢谢老师的回答，我就算是扩胞到了足够大小的Ni（111）表面还是不收敛，MULTIPLICITY 145 是我用Multiwfn生成inp文件的时候，当我设置MAGNETIZATION为1之后，Multiwfn自己生成的MULTIPLICITY 145。之前用gaussian做计算的时候自旋多重度还会判断一下，这种晶体表面的判断就有点懵，所以Multiwfn生成的是多少就直接用了，还烦请老师指导一下，感谢！

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