24小时热门版块排行榜    

查看: 2488  |  回复: 14
当前主题已经存档。

sunhg1983

银虫 (小有名气)

[交流] 【求助】表面与表面能的计算

最近想研究一下晶体的表面问题,看了些资料,还是不太明白,想请教各位大侠指点一二。困惑在于收敛测试和表面能的计算。请问:
1.收敛测试有原子层数、真空层厚度以及relax层数,那么先进行哪种测试呢?
2.表面能的计算公式:E=E(slab)-nE(bulk)/2A,我师妹认为只有当上下两个面一致的时候才除以2,是这样的么?可我的晶体没法做到上下两个面一样啊,这怎么办?
3.收敛测试是对优化好的构型进行的么?例如要测试3.6.9.12等等,E(slab)是优化后的模型能量还是仅仅切出来的原始构型能量?
4.表面能的数量级是多大呢?文献上有一点几和二点几的,我算了个表面能都3到4了,是不是很大?
5.我是学化学的,看了几篇PRB关于表面的文章,很晕。能不能帮忙给篇偏化学的文献呢,多谢了。
我的分就这么多,还请大侠们不吝赐教。

[ Last edited by aylayl08 on 2010-1-27 at 19:27 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

liumiao05

至尊木虫 (知名作家)

秦武汉风

★ ★
qasd(金币+2):谢谢 2010-01-30 17:31
引用回帖:
Originally posted by wei.ji at 2010-01-29 13:31:23:
再加一层Zn或者一层O很难么?

加了的话不能保证原子比和体相一样了
比如原子数相同的CsCl结构的AB合金,(100)方向上是交替的A原子层和B原子层,你要让两个面都是A的话比如ABABABABA, 就没法保证原子比例,你要保证原子比例ABABABAB的话用上面那个方法得到的就只是两种原子表面能之和
天下英雄出我辈,一入江湖岁月催。鸿图霸业谈笑间,不胜人生一场醉。提剑跨骑挥鬼雨,白骨如山鸟惊飞。尘世如潮人如水,只叹江湖几人回。
10楼2010-01-29 21:11:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

wei.ji

木虫 (正式写手)

★ ★
ice_rain(金币+2):谢谢指导 1-27 18:23
sunhg1983(金币+1):沙发给分,谢谢啦 2010-01-30 13:58
1.测试层数收敛。真空层不知道怎么收敛。
2.没道理不可能
3.有些测试不用,有些测试需要。
4.1、2、3、4本来就没太大区别
5.还是看物理的吧。
http://sim.phys.ruc.edu.cn
2楼2010-01-27 18:10:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hgye

金虫 (小有名气)

aylayl08(金币+2):谢谢提示,那对于上下表面不同的话该怎样计算呢? 2010-01-28 19:10
sunhg1983(金币+1):多谢啦 2010-01-30 14:00
只回答第二个问题,你师妹的理解是正确的。 当slab上下两个面不一样时,只能算出两个面表面能之和,无法算出每个表面能的绝对值。
3楼2010-01-28 18:04:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aylayl08

荣誉版主 (文坛精英)

至尊天蝎

优秀版主


zdhlover(金币+1):谢谢 2010-01-29 04:30
sunhg1983(金币+1):是不一样啊 2010-01-30 14:02
想想slab上下怎么可能不是一样的呢?毕竟你那个是算的周期性的结构,只是你看到角度不同罢了
珍惜拥有,快乐人生!欢迎光临计算模拟区(*^__^*)。。。。。。。。。。。。。。
4楼2010-01-28 19:24:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hgye

金虫 (小有名气)


ice_rain(金币+1):xiexie~ 2010-01-29 20:13
引用回帖:
Originally posted by aylayl08 at 2010-01-28 19:24:08:
想想slab上下怎么可能不是一样的呢?毕竟你那个是算的周期性的结构,只是你看到角度不同罢了

slab上下两个面不一样的情况很多了。比如钎锌矿结构ZnO垂直于c轴的面,slab上下一个是Zn在最外,一个是O在最外。
5楼2010-01-29 12:08:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wei.ji

木虫 (正式写手)

再加一层Zn或者一层O很难么?
引用回帖:
Originally posted by hgye at 2010-01-29 12:08:30:

slab上下两个面不一样的情况很多了。比如钎锌矿结构ZnO垂直于c轴的面,slab上下一个是Zn在最外,一个是O在最外。

http://sim.phys.ruc.edu.cn
6楼2010-01-29 13:31:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aylayl08

荣誉版主 (文坛精英)

至尊天蝎

优秀版主


zdhlover(金币+1):谢谢 2010-01-29 04:30
引用回帖:
Originally posted by hgye at 2010-01-29 12:08:30:

slab上下两个面不一样的情况很多了。比如钎锌矿结构ZnO垂直于c轴的面,slab上下一个是Zn在最外,一个是O在最外。

但我看有些文献并没有提这个上下表面不同的情况,即使是看到的模型是这样的
珍惜拥有,快乐人生!欢迎光临计算模拟区(*^__^*)。。。。。。。。。。。。。。
7楼2010-01-29 13:54:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liumiao05

至尊木虫 (知名作家)

秦武汉风

★ ★ ★ ★ ★ ★
aylayl08(金币+6):谢谢详细指点 2010-01-29 21:06
sunhg1983(金币+3):谢谢啦 2010-01-30 13:57
1.先真空层吧,这个比较简单,一般来说大于10A应该够了,slab稍微大点看看真空到多少精度足够,还可以检查真空层中间是不是没有电荷分布了。另外两个比较麻烦,因为计算量要大些。
2.的确有这种情况。那个公式对单质成立,对化合物的话不是每种表面都可以,因为选slab的时候不能同时兼顾原子比例和表面组成(就是说两个上下两个面都是你计算的表面)
举个例子来说,对于CsCl结构的B2晶体来说,如果你计算(100)面的表面能的话,你选的slab就不可能保证原子比1:1,外层原子肯定要多
所以需要计算化学势,具体的话你可以看看关于point defect(点缺陷)的文献
比如这篇
Concentrations of Atomic Defects in B2FexAl1-x. An Ab-Initio Study
J. MAYER, ELSASSER, M. FAHNLE
phys. stat. sol. (b) 191, 283 (1995)
3.用原始的就可以了,比如晶格常数,结构什么的可以用实验值,精度够了才谈的上优化,当然你的初始结构不能太不靠谱了。
4.这个看你的体系了吧,数量级肯定没问题,不同的体系,计算方法和软件肯定会有差别。当然如果你是纯粹重复别人的工作作为检验,那差别是有点大。
5.上面给了篇,你可以查查参考文献和引用文献

[ Last edited by liumiao05 on 2010-1-29 at 21:05 ]
天下英雄出我辈,一入江湖岁月催。鸿图霸业谈笑间,不胜人生一场醉。提剑跨骑挥鬼雨,白骨如山鸟惊飞。尘世如潮人如水,只叹江湖几人回。
8楼2010-01-29 21:03:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aylayl08

荣誉版主 (文坛精英)

至尊天蝎

优秀版主


zdhlover(金币+1):谢谢 2010-01-29 04:31
引用回帖:
Originally posted by liumiao05 at 2010-01-29 21:03:16:
1.先真空层吧,这个比较简单,一般来说大于10A应该够了,slab稍微大点看看真空到多少精度足够,还可以检查真空层中间是不是没有电荷分布了。另外两个比较麻烦,因为计算量要大些。
2.的确有这种情况。那个公式对 ...

还想问一句,人为地使slab上下一致可以吗?比如说添加原子层或是剔除,是否对模型的合理性有影响呢?回头再好好研究你给的文献和其他文献,谢谢
珍惜拥有,快乐人生!欢迎光临计算模拟区(*^__^*)。。。。。。。。。。。。。。
9楼2010-01-29 21:10:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sunhg1983 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 朋友圈看到的 +4 wangzilk 2026-08-18 5/250 2026-08-18 19:25 by wangzilk
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +51 医学老男孩 2026-08-13 110/5500 2026-08-18 19:07 by icerianxia
[基金申请] 时间戳变了,能看出什么问题? +8 基诺咪客 2026-08-17 11/550 2026-08-18 17:13 by 悠然熊猫
[基金申请] 重要消息,中午系统在维护 +9 yuleib84 2026-08-18 10/500 2026-08-18 16:52 by archvillain
[教师之家] 两姐妹,一个辛辛苦苦30岁工作当大学教师,另一个23岁就工作当开车教练 +3 瞬息宇宙 2026-08-11 7/350 2026-08-18 16:42 by LNP@mRNA
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 8/400 2026-08-18 11:47 by gltch
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
请填处理意见