计算化学公社

 找回密码 Forget password
 注册 Register

NiO表面能计算问题

查看数: 14534 | 评论数: 12 | 收藏 Add to favorites 3
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
sky
发布时间: 2016-12-27 14:50

正文摘要:

我在做NiO表面能计算的时候遇到了些问题,NiO的表面有重构现象,要拿掉最表层的四分之三的原子和第二层的四分之一的原子,体系才能稳定存在。如图: 求大神指点一下,谢谢!

回复 Reply

卡开发发 发表于 Post on 2018-1-17 13:22:59
幸福每个人 发表于 2018-1-17 11:25
OK,多谢,还有个问题想请教下,就是在计算表面吸附时,要确定表面层数、表面尺寸大小、以及收敛精度,kp ...

一般优先测试SCF的参数,比如动能截断、kpoints,然后再是层高和真空,最后才考虑覆盖度。多个变量需要测试的时候可以把其中几个设置到充分大去测试其中一个,比如处理层高的时候可以先取个足够大的真空层。

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
sobereva + 1

查看全部评分 View all ratings

幸福每个人 发表于 Post on 2018-1-17 11:25:19
卡开发发 发表于 2018-1-16 17:26
有文章是按照先从bulk->slab->缺陷slab这两步做热化学循环来算的,这样第一步的表面能就能算出来。不过真 ...

OK,多谢,还有个问题想请教下,就是在计算表面吸附时,要确定表面层数、表面尺寸大小、以及收敛精度,kpoint测试,这些参数确定有先后顺序吗,如果有,应该遵循什么样的顺序较为合理呢
卡开发发 发表于 Post on 2018-1-16 17:26:58
幸福每个人 发表于 2018-1-16 17:03
所以你在导入NiO晶胞时 是导入的它的原胞,从而得到的能量是计算的单胞的能量,也就是Ebulk,不知道之前 ...

有文章是按照先从bulk->slab->缺陷slab这两步做热化学循环来算的,这样第一步的表面能就能算出来。不过真按照表面能定义较真一下,bulk->slab->缺陷slab还是bulk->缺陷bulk->缺陷slab结果其实并不一样,因此缺陷的表面这样处理表面能我觉得不是well-defined。

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
sobereva + 1

查看全部评分 View all ratings

幸福每个人 发表于 Post on 2018-1-16 17:03:53
zorow 发表于 2018-1-15 22:09
这个N是什么?是表面能计算公式中的N吗,如果是,这里的N等于slab中含有的原胞数量,具体的解释请看这帖 ...

所以你在导入NiO晶胞时 是导入的它的原胞,从而得到的能量是计算的单胞的能量,也就是Ebulk,不知道之前你有没有切出这样的表面, 就是ABC原子数与原胞原子数不成比例,我有切过化学式是Ca2CuO3的物质,切出来的表面没法算n,不知道对于这种情况,要怎么解决n的问题
zorow 发表于 Post on 2018-1-15 22:09:16
幸福每个人 发表于 2018-1-15 21:16
那这样的话在切完平面后,n怎么确定呢

这个N是什么?是表面能计算公式中的N吗,如果是,这里的N等于slab中含有的原胞数量,具体的解释请看这帖子中卡开发发老师的回答:http://bbs.keinsci.com/forum.php ... &extra=page%3D1
幸福每个人 发表于 Post on 2018-1-15 21:16:43
zorow 发表于 2018-1-15 11:28
先从MS中导入NiO的晶胞,经几何优化后得到的结构其能量就是Ebulk吧

那这样的话在切完平面后,n怎么确定呢
zorow 发表于 Post on 2018-1-15 11:28:56
幸福每个人 发表于 2018-1-15 11:14
想问下对于AB化合物那个Ebulk,在MS里您是如何确定的呢

先从MS中导入NiO的晶胞,经几何优化后得到的结构其能量就是Ebulk吧
幸福每个人 发表于 Post on 2018-1-15 11:14:39
zorow 发表于 2018-1-12 20:22
话说楼主是用什么软件做的啊,也是MS吗

想问下对于AB化合物那个Ebulk,在MS里您是如何确定的呢
itpfeng 发表于 Post on 2018-1-12 20:49:41
想知道实验上NiO的表面是什么样的,建的模是不是符合实际。
zorow 发表于 Post on 2018-1-12 20:22:16
话说楼主是用什么软件做的啊,也是MS吗
zorow 发表于 Post on 2018-1-12 20:21:34
追枫少年 发表于 2018-1-12 11:12
楼主的问题解决了吗  我这几天也在琢磨表面能
我认为计算的话γ=[Eslab-nEbulk] /2A这个公式合理一点吧,n ...

是的,听了卡开发发老师的指导,我也觉得应该是用这个公式
追枫少年 发表于 Post on 2018-1-12 11:12:37
楼主的问题解决了吗  我这几天也在琢磨表面能
我认为计算的话γ=[Eslab-nEbulk] /2A这个公式合理一点吧,n是slab模型的原子格式,Ebulk是体材料的平均原子能量

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-18 00:49 , Processed in 0.187082 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list