24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 2646  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 338求调剂 +5 小猪红色 678 2026-04-06 6/300 2026-04-07 21:18 by 乔哒哒哒
[考研] 388求调剂 +5 四川王涛 2026-04-07 7/350 2026-04-07 19:01 by dongzh2009
[考研] 081700,311,求调剂 +17 冬十三 2026-04-04 18/900 2026-04-07 12:50 by Sammy2
[考研] 0703调剂 +16 拾玖壹 2026-04-04 18/900 2026-04-07 12:49 by flydream1314
[考研] 一志愿北交大材料工程总分358求调剂 +10 cs0106 2026-04-05 12/600 2026-04-06 19:41 by 无际的草原
[考研] 求助 +3 卡卡东88 2026-04-06 4/200 2026-04-06 15:28 by going home
[考研] 考研调剂生寻找导师 +3 顾瞻考研啊 2026-04-05 3/150 2026-04-05 18:18 by 啵啵啵0119
[考研] 298求调剂 +3 manman511 2026-04-05 3/150 2026-04-05 18:09 by kk112233
[考研] 282求调剂 +7 aaa车辆 2026-04-02 11/550 2026-04-05 17:24 by yulian1987
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 22:08 by 啵啵啵0119
[考研] 294求调剂 +6 Grey_Ey 2026-04-02 9/450 2026-04-04 22:07 by hemengdong
[考研] 338求调剂 +7 晟功? 2026-04-03 7/350 2026-04-04 20:37 by 蓝云思雨
[考研] 怎么删帖子啊 +3 缝曦1000 2026-04-04 3/150 2026-04-04 14:20 by 土木硕士招生
[考研] 336求调剂 +8 kiyy 2026-04-01 8/400 2026-04-03 19:41 by lijunpoly
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +9 哇呼哼呼哼 2026-04-02 9/450 2026-04-03 12:05 by 1753564080
[考研] 260求调剂 +3 朱芷琳 2026-04-02 3/150 2026-04-03 08:44 by yulian1987
[考研] 372分材料与化工(085600)一志愿湖南大学求调剂 +5 蓝笺片 2026-04-02 6/300 2026-04-02 21:37 by dongzh2009
[考研] 318求调剂,计算材料方向 +10 吸喵有害笙命 2026-04-01 11/550 2026-04-02 16:29 by oooqiao
[考研] 348求调剂 +6 吴彦祖24k 2026-04-02 6/300 2026-04-02 14:07 by 给你你注意休息
[考研] 348环境工程调剂 +3 吴彦祖24k 2026-04-01 3/150 2026-04-02 09:14 by nanaliuyun
信息提示
请填处理意见