| 查看: 3297 | 回复: 8 | |||
[求助]
VASP计算TiO2表面能问题 求助 已有1人参与
|
|
构建了一个4*2的TiO2 SLAB模型 看到表面能计算公式里有E(bulk) 和E(slab) 想请教这两个能量具体都怎么算啊 如何从输出文件里得到呢？ 一次计算能同时得到这两个能量么，还是要计算两次。 我目前是这样理解的，是不是设置了让表面原子弛豫厚计算得到的总能就是E(slab)，同样的模型固定所有原子计算得到的总能就是E(bulk)？？ 还请各位前辈指教，菜鸟一只，没有师兄带，自己乱撞，还望有表面计算经验的前辈们多多指点。 |
回复此楼» 猜你喜欢
22408 调剂材料
已经有6人回复
-
285求调剂
已经有12人回复
-
求调剂
已经有6人回复
-
085600材料与化工301分求调剂院校
已经有19人回复
-
277工科求调剂
已经有11人回复
-
277求调剂 数一104分
已经有13人回复
-
304求调剂
已经有10人回复
-
336求调剂，一志愿中科大
已经有6人回复
-
071000生物学，一志愿深圳大学296分，求调剂
已经有3人回复
-
一志愿北京化工085600 310分求调剂
已经有19人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 极性面应该如何计算表面能？
已经有23人回复
- 表面能计算的困惑
已经有27人回复
- 用VASP计算单点能的时候总是出错，怎么办啊？大侠来给知道啊？
已经有8人回复
- vasp计算中的费米能问题
已经有20人回复
- 如何计算表面覆盖率？
已经有19人回复
- 关于vasp计算表面能中的Ebulk
已经有7人回复
- vasp计算单个氧原子和分子的能量
已经有20人回复
- 用castep计算TiO2表面能和吸附能
已经有4人回复
- ms中吸附能的计算
已经有40人回复
- TiO2(110)的面密度？
已经有3人回复
- 计算出极化强度矢量后，怎么判断它是沿哪个晶向？（如[111]，[001]）
已经有5人回复
- VASP一次性优化锐钛矿TiO2晶格常数时出错！！
已经有6人回复
- 【求助】单胞与原包能量不一样？
已经有18人回复
- 【求助】新手求助vasp计算氧分子结合能的问题
已经有12人回复
- 【求助】算VASP弹性常数时的能量-应变出错
已经有18人回复
- 【求助】如何用VASP 计算体系的表面能
已经有11人回复
- 【求助】Castep 表面能的计算
已经有13人回复
米色空白
金虫 (小有名气)
- 应助: 15 (小学生)
- 金币: 1925.4
- 红花: 3
- 帖子: 73
- 在线: 30小时
- 虫号: 2132467
- 注册: 2012-11-17
- 性别: MM
- 专业: 高分子材料的加工与成型
★ ★
identation: 金币+2, 谢谢交流 2013-11-12 13:55:01
identation: 金币+2, 谢谢交流 2013-11-12 13:55:01
|
我没做过，但是看到过一点儿资料，你不妨去拿来读一读：德国亚琛工业大学无机化学研究所在JPCC上发表的两篇文章，一篇是GeTe表面的，另一篇是GeSb2Te4、Ge2Sb2Te5、GeSb4Te7的。另外John Wiley有一本书：Density Functional Theory：A Practical Introduction里面也有相应的一章内容介绍，论坛上有，可以拿来读。 至于具体到你的问题，E(bulk)应当是不加真空层的模型充分弛豫以后的静态计算的结果，可以grep TOTEN OUTCAR来查看；同理E(slab)是加上真空层以后的值。我觉着要分别优化，分别静态计算得来。 |
» 本帖已获得的红花（最新10朵）
zjucy2楼2013-11-05 00:21:44
3楼2013-11-05 08:18:39
|
刚刚看了一下Wiley的那本书，有这样一段话 “In the case of the surface, one would typically be using a comparatively large supercell, including a vacuum space, and using comparatively few k points. In the case of the bulk, the opposite would be true. ” 意思是算E(SLAB)的时候要用大的超胞，加真空层，较少的K点? 算E(BULK)的时候就用单胞，不加真空层，较多的K点？ |
4楼2013-11-05 09:12:44
米色空白
金虫 (小有名气)
- 应助: 15 (小学生)
- 金币: 1925.4
- 红花: 3
- 帖子: 73
- 在线: 30小时
- 虫号: 2132467
- 注册: 2012-11-17
- 性别: MM
- 专业: 高分子材料的加工与成型
5楼2013-11-05 16:16:47
6楼2013-11-12 10:24:50
米色空白
金虫 (小有名气)
- 应助: 15 (小学生)
- 金币: 1925.4
- 红花: 3
- 帖子: 73
- 在线: 30小时
- 虫号: 2132467
- 注册: 2012-11-17
- 性别: MM
- 专业: 高分子材料的加工与成型
7楼2013-11-13 00:11:51
8楼2013-11-13 08:57:23
【答案】应助回帖
|
你好，我比较推荐Fiorentini and Methfesse方法 In the limit of large N one can rewrite Eq. 1 as Eslab = 2Sigma + N*Ebulk If the total energy of the slab depends linearly on slab thickness N, the bulk energy term Ebulk can be taken as the slope and used in Eq. 1, also avoiding a calculation on a separate bulk system. Eq.1就是表面能的定义式 |
9楼2014-01-04 04:44:10
10