| 查看: 864 | 回复: 3 | |||
| 当前主题已经存档。 | |||
liumiao05至尊木虫 (知名作家)
秦武汉风
|
[交流]
【求助】请问表面能计算时slab如何建立的问题，急
|
||
|
先说一下我想解决的问题 我想计算一个二元合金的表面结构，同一层上肯定是同种原子，它们相间排列，即: ---A原子层---B原子层---A原子层---B原子层---A原子层---B原子层--- 现在有两种可能，一是表面层全是A原子，另外就是表面层全是B原子（实验证明是其中一种） 我的想法是取相同层数的A，B两种原子，按两种方式排列: (体相)A--B--A--B.........--A--B(表面) (体相)B--A--B--A.........--B--A(表面) 固定下面若干层作体相然后对上面数层做弛豫，分别计算两种情况的表面能，表面能小的更稳定,跟这个帖子差不多 http://muchong.com/bbs/viewthread.php?tid=1199160&fpage=4 举个例子，模拟的时候我选择一个晶胞比如有六层A原子和六层B原子，全部固定不动作静态计算得到一个Eb, 然后固定四层A原子和四层B原子作为体相，弛豫上面两层A,B原子（上面还有真空层）得到一个Ea，这个时候只有一个表面，就是最上层原子与真空接触的面积。那么该怎么计算表面能呢？是不是公式中就不用除以2，还是说这种想法不可行？ 我看了下很多帖子及一些文献，表面能的公式是 （Eslab-N*Ebulk）/2A ，这个2是因为slab有两个表面，而且可以层数从1开始递增，我不太明白这时候slab是怎样建立的，是不是没有固定层，而是直接取1，2，3，……层原子，然后全部弛豫，看表面能变化趋势。如果是，这个弛豫方式在vasp里怎么设置呢，还是表面弛豫方式吗，把所有的坐标分量都设置为T？ 另外问题在于，我如果用这种模型，A原子和B原子的层数就不一样了，如下 (1)(表面)A--B--A--B--A(表面) (2)(表面)B--A--B--A--B(表面) 如果将两种结果得到的表面能进行比较能说明我最开始提出的问题吗？我个人觉得系统组分变了，比较不同排列的表面能没有说服力吧 谢谢，请大家指点 [ Last edited by liumiao05 on 2009-5-13 at 03:53 ] |
回复此楼» 猜你喜欢
面上可以超过30页吧？
已经有13人回复
-
网上报道青年教师午睡中猝死、熬夜猝死的越来越多，主要哪些原因引起的？
已经有10人回复
-
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有13人回复
-
什么是人一生最重要的？
已经有8人回复
-
版面费该交吗
已经有17人回复
-
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有19人回复
-
【博士招生】太原理工大学2026化工博士
已经有8人回复
-
280求调剂
已经有4人回复
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
qasd(金币+2,VIP+0):谢谢~ 5-13 19:56
liumiao05(金币+4,VIP+0):谢谢；明天俺和老师讨论下，今天看的头昏脑胀的才有了点头绪 5-14 01:42
liumiao05(金币+6,VIP+0):问题还是没那么简单啊，两种元素的表面能耦合在一起，还得引入化学势，光靠选取结构没法消去未知数 5-16 02:57
qasd(金币+2,VIP+0):谢谢~ 5-13 19:56
liumiao05(金币+4,VIP+0):谢谢；明天俺和老师讨论下，今天看的头昏脑胀的才有了点头绪 5-14 01:42
liumiao05(金币+6,VIP+0):问题还是没那么简单啊，两种元素的表面能耦合在一起，还得引入化学势，光靠选取结构没法消去未知数 5-16 02:57
|
我猜猜，（Eslab-N*Ebulk）/2A，中的2 是因为表面有上下之分 只要你的slab中的原子和bulk中总的原子数是一致的 然后A取成层的面积即可 驰豫时还是固定底层，放开上面两层 第二个问题 个人看法是你不妨取 ABABAB 和BABABA，这样原子数就一致了 |
2楼2009-05-13 09:19:18
3楼2009-05-16 09:33:48
liumiao05
至尊木虫 (知名作家)
秦武汉风
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 金币: 20224.3
- 红花: 1
- 帖子: 6935
- 在线: 61.1小时
- 虫号: 697035
- 注册: 2009-02-05
- 性别: GG
- 专业: 无机合成和制备化学
4楼2009-05-20 00:31:18