24小时热门版块排行榜    

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

wangle_xq

银虫 (小有名气)

[交流] 【求助】VASP计算表面的疑问

刚刚开始使用VASP计算表面的能量。

遇到这样一个问题,文献中提到建立一个由7层原子和一个真空层组成的超晶胞,由于第一性原理的周期性的特性,这样就有两个表面。未免翻译错误,附上原文:
The calculation is performed for a fully relaxed slab with 7 layers, each layer
consisting of an (8×8) lattice sites matrix. This setup models a thin film with free surfaces at the top and the bottom.

到这儿我还是可以理解的,随后文献中又提到,为了控制计算的精确度,在处于平衡的系统下加入了扩散势垒区(diffusion barrier),扩散势垒区共有15层原子组成,其中底部的3层原子固定不动。
To control the accuracy of the calculation, the diffusion barrier for a system
in equilibrium was calculated for a slab consisting of 15 (8×8) layers, with the bottom 3 layers fixed.

我的问题是:
1.为什么要在表面计算的时候在底部加入3层固定不动的原子?

2.加入这样的原子如何控制计算的精确度?要怎么分析?文献中没有具体展开,只是说7层原子和15层原子的计算结果未见明显差别,然后就说7层原子的计算结果可靠;

3.计算表面时,对于真空层有怎样的要求?真空层需要的最小厚度如何确定?

初来乍到,还望各位多多帮忙啊!

[ Last edited by wuchenwf on 2009-6-23 at 19:51 ]
回复此楼

» 收录本帖的淘帖专辑推荐

Photochemistry

» 猜你喜欢

网络日志 http://blog.163.com/wangle_xq/ 关注计算材料科学的应用,以及相关的计算机技术。文章以转贴为主,少量原创。欢迎参观
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

yhjiangshiyu

木虫 (正式写手)

★ ★ ★ ★ ★ ★
wuli8(金币+3,VIP+0):(*^__^*)辛苦了 3-5 09:20
wangle_xq(金币+2,VIP+0):我刚来,还要写评价啊,那先说谢谢吧。 3-5 14:00
wangle_xq(金币+1,VIP+0):最后一分给沙发吧! 3-9 02:40
1 底部固定是为了节约计算时间,提高计算效率;
2 比较7层与15层的优化结果,能量,结构弛豫等,主要是所计算体系的结果差异,如无差别或差别较小(具体就要自己把握了),则用较少的原子模型;
3 真空层厚度有8,10,12,15nm等等,也是要看对你的计算体系有没有影响,如果是计算小分子吸附的话一般都要8-10nm左右,如果吸附分子较大则真空层相应增大。
流星,虽然短暂,却迸发出最强的光芒,耀眼夺目。
2楼2009-03-05 08:21:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyb4142

金虫 (正式写手)

新手初来,学习一下!
4楼2009-03-05 13:14:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

tfh00

银虫 (小有名气)

★ ★ ★ ★ ★
wuli8(金币+3,VIP+0):(*^__^*)谢谢,辛苦了。 3-5 09:20
wangle_xq(金币+2,VIP+0):金币不多了,剩下一个,看看还能给谁吧,谢谢啊 3-5 14:00
我觉得 底部固定可以这样理解:当固体被切,露出表面时,接近表面的几层原子会发生变化,但是接近内部的层就和没切相似了,在你这里可以认为是表面只对7层中的4层有明显影响,4层以下的就没什么影响了,所以4层以下的就可以固定不动。
7和15层没有明显差别,也可以做类似的理解。
真空层一般大于10埃
3楼2009-03-05 08:44:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cbFeng200881

金虫 (著名写手)

计算前期的工作都是为了找到一个最优化的模型。。。
5楼2009-03-05 13:32:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chineseli


qasd(金币+1,VIP+0):xiexie 3-11 18:07
使用这种板层模型都是在bloch定理下三维无限重复的,所以要用真空层使每个重复的层之间没有作用,这样那个“表面”就可以认为是表面了~
6楼2009-03-05 14:05:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangle_xq

银虫 (小有名气)

谢谢2,3楼,另外还有一个附加的问题,

关于原子驰豫的设置,我在论坛上看到有人讨论时提问:算表面能的时候,弛豫之后原子把真空层填满了,就变成连续的比较松散的bulk了。

而后又有人回答说:这种情况很有可能是他将ISIF设置为3来计算的,这样原子驰豫会很强,晶格也会发生很强畸变。故,做表面计算时最好固定晶格而只驰豫原子。

这个回答的意见,是不是ISIF设定为2啊。这个设定是表面计算必须的么?
网络日志 http://blog.163.com/wangle_xq/ 关注计算材料科学的应用,以及相关的计算机技术。文章以转贴为主,少量原创。欢迎参观
7楼2009-03-05 14:09:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tfh00

银虫 (小有名气)


zdhlover(金币+1,VIP+0):谢谢 3-13 20:27
可以试试  ISIF=4
8楼2009-03-05 15:43:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangdelin0000

木虫 (著名写手)

tfh00  你能具体说说ISIF=2和ISIF=4的区别吗,对计算的结果到底有什么样的影响,新手,非常感谢!
9楼2009-03-06 23:22:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)

★ ★
spur(金币+2,VIP+0): 3Q!!谢谢参与! 3-11 17:01
为什么要在表面计算的时候在底部加入3层固定不动的原子?

真实的材料的表明只有一面,也就是仅一面与真空接触。由于VASP软件处理的周期性,我们构建的超胞就有两个表面,这不符合实际情况。因此计算时就让下面几层原子固定在体材料的位置,只允许上面几层原子动。所以文献中时常有这样的话 the A atoms of the bottom two A layers were frozen during the relaxation to mimic the bulk.
10楼2009-03-11 16:42:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wangle_xq 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见