计算化学公社

标题: 反磁性物质表面的吸附 [打印本页]

作者
Author: 来巡老和山 时间: 2016-12-27 16:38
标题: 反磁性物质表面的吸附
各位大神，我想请问下：我在对反磁性物质进行计算时，优化、静态和切表面优化时得到的mag=0，但是Hg吸附在上面时，总磁矩不在为0。想咨询下这是为什么，是要对Hg原子进行磁矩的设置吗（我计算时对Hg设为0）？恳请各位大神予以指导，或者给点意见也好。

作者
Author: 卡开发发 时间: 2016-12-27 21:58
给个LORIBIT=11看看各自的磁矩是多少，同时检查一下是否每一步的SCF都是收敛的。实在搞不定把输入输出发来看看。

作者
Author: 来巡老和山 时间: 2016-12-28 10:01
本帖最后由来巡老和山于 2016-12-28 18:16 编辑

卡开发发发表于 2016-12-27 21:58
给个LORIBIT=11看看各自的磁矩是多少，同时检查一下是否每一步的SCF都是收敛的。实在搞不定把输入输出发来 ...

那我先加LORBIT进行计算！如果解决不掉，再向您请教吧！

作者
Author: 来巡老和山 时间: 2016-12-28 20:19

卡开发发发表于 2016-12-27 21:58
给个LORIBIT=11看看各自的磁矩是多少，同时检查一下是否每一步的SCF都是收敛的。实在搞不定把输入输出发来 ...

大神您好，这是我今天计算的Co3O4（110）表面的优化。A和B为两种不同的切法。设置Co2+磁矩为3，一正一负。得到的结果中A的mag=2，而B的mag=0。

作者
Author: 卡开发发 时间: 2016-12-28 21:03

来巡老和山发表于 2016-12-28 20:19
大神您好，这是我今天计算的Co3O4（110）表面的优化。A和B为两种不同的切法。设置Co2+磁矩为3，一正一负 ...

1、两个表面原子数不同；
2、两个表面均不符合化学计量比。

作者
Author: 来巡老和山 时间: 2016-12-28 21:22

卡开发发发表于 2016-12-28 21:03
1、两个表面原子数不同；
2、两个表面均不符合化学计量比。

切得平面是按照文献中来的。

作者
Author: 卡开发发 时间: 2016-12-29 00:42

来巡老和山发表于 2016-12-28 21:22
切得平面是按照文献中来的。

注意你建立的模型顶上还有原子，我建议你在z方向建立个超胞之后再和文献对比。

作者
Author: 来巡老和山 时间: 2016-12-29 12:30

卡开发发发表于 2016-12-29 00:42
注意你建立的模型顶上还有原子，我建议你在z方向建立个超胞之后再和文献对比。

我上传的两个图就是建立晶胞之后的，只是我画图的时候把边界去掉了。但是我建立的超晶胞上方并不存在原子啊。而且昨天我切了四层，两个结构都满足化学计量数之比吗，然后进行计算，两个表面的总磁矩都不为0了。

作者
Author: 卡开发发 时间: 2016-12-29 13:43

来巡老和山发表于 2016-12-29 12:30
我上传的两个图就是建立晶胞之后的，只是我画图的时候把边界去掉了。但是我建立的超晶胞上方并不存在原子 ...

哦，可能vesta里面看起来比较奇怪。
1、一般切面thickness 是整数才可以保证化学计量比，至于是否符合化学计量比，你可以看一下poscar 文件的各个原子数。
2、bulk反铁磁系统也不见得能够保证表面总磁矩是0，上下表面化学环境不对称的话可能有些许偏差，但不会是2这么大。

作者
Author: 来巡老和山 时间: 2016-12-29 13:48

卡开发发发表于 2016-12-29 13:43
哦，可能vesta里面看起来比较奇怪。
1、一般切面thickness 是整数才可以保证化学计量比，至于是否符合化 ...

今天刚上传的两个图是我用ms切得的，图片也是在ms中截屏得到的。
切表面时确实thickness不为整数。
所以现在很想知道到底哪里出错了，感觉设置的参数并没有出现问题啊

作者
Author: 卡开发发 时间: 2016-12-29 16:56

来巡老和山发表于 2016-12-29 13:48
今天刚上传的两个图是我用ms切得的，图片也是在ms中截屏得到的。
切表面时确实thickness不为整数。
所 ...

不符合化学计量比的话，多出的原子当然会引入额外磁矩。

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