24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1600  |  回复: 14
当前主题已经存档。

xgquanjun

木虫 (正式写手)

[交流] 【求助】建立好晶面后怎样计算表面能

我用的TiO2 的体系,我切割了晶面(101),想用Castep来计算(101)晶面的 表面能,请问如何操作呀?
我怎么 在哪里都 找不到这样的 按钮或选项呀?

或者说在哪里能计算表面能,怎样计算呢 ?

[ Last edited by SHY31 on 2008-12-14 at 15:54 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shrek826

金虫 (小有名气)

★ ★ ★
zdhlover(金币+3,VIP+0):谢谢,欢迎常来
CASTEP上并没有直接计算表面能的选项
首先你需要计算TiO2晶体的能量
其次你需要切一个面,加上真空层,一般12埃,计算其能量
第三步就是拿块体的能量减去slab的能量,除以2,就是你所切面的表面能
2楼2008-12-13 19:56:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yhjiangshiyu

木虫 (正式写手)

★ ★
zdhlover(金币+2,VIP+0):如果能给出公式,那就更好了,嘿嘿
还要除以面积吧,另外与原子个数也有关,可以参考文献中的计算公式
流星,虽然短暂,却迸发出最强的光芒,耀眼夺目。
3楼2008-12-13 22:13:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xgquanjun

木虫 (正式写手)

TiO2晶体的能量是指哪个呀?
我计算的TiO2结果如下:
Pseudo atomic calculation performed for O 2s2 2p4

Converged in 31 iterations to a total energy of -430.640066 eV


Pseudo atomic calculation performed for Ti 3s2 3p6 3d2 4s2

Converged in 35 iterations to a total energy of -1598.444997 eV

请问哪个是TiO2的 能量呢?
初学者请见谅!
4楼2008-12-14 15:01:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cos60

木虫 (著名写手)

★ ★ ★
zdhlover(金币+3,VIP+0):非常感谢,欢迎以后多来
MS自带的help和tutorial是最好的教科书,除了告诉你MS怎么用以外,还告诉你物理。MS的帮助不仅是软件教师还是物理教师。专研吧,那面有公式。

open MS---->
Help ----->
Help topic ------>
index------>
enter "surface, free energy"------->
u can see Eq. Th-3, it is what u want.

[ Last edited by cos60 on 2008-12-14 at 15:04 ]
5楼2008-12-14 15:01:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cos60

木虫 (著名写手)

★ ★
zdhlover(金币+2,VIP+0):多谢
引用回帖:
Originally posted by xgquanjun at 2008-12-14 15:01:
TiO2晶体的能量是指哪个呀?
我计算的TiO2结果如下:
Pseudo atomic calculation performed for O 2s2 2p4

Converged in 31 iterations to a total energy of -430.640066 eV


Pseudo atomic calcu ...

看来楼主是要恶补一下了,这两个是O和Ti原子的能量,在这个文件的后面,还有一个总能。TiO2晶体的结合能这样算:
结合能=总能-O的原子数xO的原子能量+Ti的原子数xTi的原子能量

还是VASP好,给减完了... ...
6楼2008-12-14 16:04:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xgquanjun

木虫 (正式写手)

谢谢各位楼主回答,真是感激不尽
还有一个疑问,2楼说:”其次你需要切一个面,加上真空层,一般12埃,“
我看文献都是说有多少个真空层,怎么知道的呢? 加上的真空层与后面的埃数有什么关系呢??
7楼2008-12-15 09:50:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

amynliou

金虫 (正式写手)

★ ★ ★
fegg7502(金币+3,VIP+0):thanks
引用回帖:
Originally posted by xgquanjun at 2008-12-15 09:50:
谢谢各位楼主回答,真是感激不尽
还有一个疑问,2楼说:”其次你需要切一个面,加上真空层,一般12埃,“
我看文献都是说有多少个真空层,怎么知道的呢? 加上的真空层与后面的埃数有什么关系呢??

用来模拟表面时候,需要加上一定厚度真空层,厚度单位用埃。至少要大于10埃。
8楼2008-12-15 11:00:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xgquanjun

木虫 (正式写手)

请问6楼
我按你的方式计算的TiO2晶体的表面能是-106 eV
可是我看文献一般计算的TiO2晶体的表面能是 -21eV
怎么 相差这么大呢?
请问为什么,是什么原因呀?
Pseudo atomic calculation performed for O 2s2 2p4

Converged in 31 iterations to a total energy of -430.640066 eV


Pseudo atomic calculation performed for Ti 3s2 3p6 3d2 4s2

Converged in 35 iterations to a total energy of -1598.444997 eV
...
External pressure/stress (GPa)
                          0.00000   0.00000   0.00000
                                    0.00000   0.00000
                                              0.00000
------------------------------------------------------------------------ <-- SCF
SCF loop      Energy           Fermi           Energy gain       Timer   <-- SCF
                               energy          per atom          (sec)   <-- SCF
------------------------------------------------------------------------ <-- SCF
Initial  -9.05372665E+003  5.89364312E+001                         9.59  <-- SCF
      1  -9.84772076E+003  8.85584745E+000   6.61661763E+001      34.94  <-- SCF
      2  -9.96394821E+003  7.51189967E+000   9.68562113E+000      59.42  <-- SCF
      3  -9.95181270E+003  4.90373159E+000  -1.01129273E+000      84.05  <-- SCF
      4  -9.94594846E+003  6.19582689E+000  -4.88687071E-001     109.38  <-- SCF
      5  -9.94533818E+003  5.77295445E+000  -5.08565712E-002     133.95  <-- SCF
      6  -9.94541429E+003  5.63144052E+000   6.34300201E-003     157.64  <-- SCF
      7  -9.94537784E+003  5.77534499E+000  -3.03788343E-003     180.20  <-- SCF
      8  -9.94537817E+003  5.78948518E+000   2.79255077E-005     197.52  <-- SCF
      9  -9.94537854E+003  5.80979585E+000   3.06409531E-005     212.91  <-- SCF
     10  -9.94537858E+003  5.81275461E+000   2.98111007E-006     224.97  <-- SCF
     11  -9.94537856E+003  5.81017014E+000  -1.43547204E-006     236.77  <-- SCF
     12  -9.94537856E+003  5.80970546E+000  -1.68120395E-008     248.62  <-- SCF
------------------------------------------------------------------------ <-- SCF
Final energy =  -9945.378560459     eV
9楼2008-12-15 16:52:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xgquanjun

木虫 (正式写手)

求高手的出现?
10楼2008-12-16 08:26:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xgquanjun 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见