计算化学公社

标题: Fe(110)结构优化怎么一直不收敛? [打印本页]

作者
Author:
gch123    时间: 2019-9-25 21:40
标题: Fe(110)结构优化怎么一直不收敛?
如题,切了个表面一直优化不成功,应该怎么才能使表面收敛呢?

作者
Author:
卡开发发    时间: 2019-9-25 22:47
1、设置一下磁矩。
2、减小density mixing-charge。
作者
Author:
gch123    时间: 2019-9-28 09:30
卡开发发 发表于 2019-9-25 22:47
1、设置一下磁矩。
2、减小density mixing-charge。

感谢老师的回复,我charge从0.5减小到0.3后任然不收敛,但我把泛函换成RPBE之后收敛了,另外请教一下这个磁矩一般是根据什么原理去设置呢?
作者
Author:
卡开发发    时间: 2019-9-28 09:36
本帖最后由 卡开发发 于 2019-9-28 09:44 编辑
gch123 发表于 2019-9-28 09:30
感谢老师的回复,我charge从0.5减小到0.3后任然不收敛,但我把泛函换成RPBE之后收敛了,另外请教一下这个 ...

对Fe这个肯定能断定是铁磁性的,具体磁矩的具体大小在这个体系不是很重要,磁序调整正确就行。如果不设置磁矩可能结构会被算到没磁性或者反铁磁的状态,所以才导致构型优化难收敛。按道理,使用RPBE和PW91都能算出来。

作者
Author:
gch123    时间: 2019-9-28 10:20
卡开发发 发表于 2019-9-28 09:36
对Fe这个肯定能断定是铁磁性的,具体磁矩的具体大小在这个体系不是很重要,磁序调整正确就行。如果不设置 ...

老师,我想计算的是Fe(110)-石墨烯界面界面结构,搭建的结构进行结构优化一直不收敛,搭建的结构如图,K点选择为13*4*1,减小了density mix charge也不收敛,这应该怎么修改参数呢?
作者
Author:
卡开发发    时间: 2019-9-28 10:30
gch123 发表于 2019-9-28 10:20
老师,我想计算的是Fe(110)-石墨烯界面界面结构,搭建的结构进行结构优化一直不收敛,搭建的结构如图,K ...

density mixing其实是解决SCF难收敛的。如果磁序设置正确,且计算参数合理(比如收敛性测试问题),就得考虑模型是否合理,比如里面两个界面是否匹配比较好。
作者
Author:
gch123    时间: 2019-9-28 15:39
卡开发发 发表于 2019-9-28 10:30
density mixing其实是解决SCF难收敛的。如果磁序设置正确,且计算参数合理(比如收敛性测试问题),就得 ...

晶面匹配度是根据MS里面的限制来设置的,失配率应该是小于5%,另外设置磁矩是在 initial spin这个位置设置吗?
作者
Author:
卡开发发    时间: 2019-9-28 18:03
gch123 发表于 2019-9-28 15:39
晶面匹配度是根据MS里面的限制来设置的,失配率应该是小于5%,另外设置磁矩是在 initial spin这个位置设 ...

modify->electronic configuration




欢迎光临 计算化学公社 (http://ccc.keinsci.com/) Powered by Discuz! X3.3