24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

chenyang17

木虫 (小有名气)

应助: 2 (幼儿园)
金币: 3523.7
散金: 12
红花: 2
帖子: 205
在线: 249.6小时
虫号: 440378
注册: 2007-10-27
性别: GG
专业: 其他无机非金属材料

[交流] 【求助】castep与vasp的能量计算值有可比性吗?

分比用CASTEP和VASP同样计算一个Ti单胞的能量，采用相同的结构，a=2.9506， c=4.6788，其中VASP计算时用的poscar是从.cell文件里面导出来的，K点都取9*9*2，参数基本都相同，为什么算出来的能量差别那么大，大家看看：
这是我的castep里面算出来的结果：
部分参数：
basis set accuracy                            : FINE
plane wave basis set cut-off                : 310.0000 eV
size of standard grid                         :    2.0000
largest prime factor in FFT                   :       5
结构：
Lattice parameters(A)    Cell Angles
                  a = 2.950600       alpha = 90.000000
                  b = 2.950600       beta  = 90.000000
                  c = 4.678800       gamma =  120.000000

         Current cell volume = 35.276524    A**3
能量计算结果：
Final energy, E          =  -3206.122656156    eV
Final free energy (E-TS) =  -3206.132108357    eV
(energies not corrected for finite basis set)

NB est. 0K energy (E-0.5TS)    =  -3206.127382256    eV

这是VASP里面采取相同结构的计算结果：
   direct lattice vectors
   2.555294556 -1.475300000  0.000000000
   0.000000000  2.950600000  0.000000000
   0.000000000  0.000000000  4.678800000
  energy-cutoff  :    178.33
  volume of cell :    35.28

  FREE ENERGIE OF THE ION-ELECTRON SYSTEM (eV)
  ---------------------------------------------------
  free  energy TOTEN  =    -15.739127 eV

  energy  without entropy=    -15.380196  energy(sigma->0) =    -15.559661
为什么用CASTEP和vasp的计算总能结果会差那么多，原来以为是单位错了，可是单位都是eV啊？请教高手？
2.还有一个问题就是一直搞不清楚castep输出的能量：
Final energy, E          =  -3206.122656156    eV
Final free energy (E-TS) =  -3206.132108357    eV
(energies not corrected for finite basis set)
NB est. 0K energy (E-0.5TS)    =  -3206.127382256    eV
有什么区别总能的话应该是那个？

还有VASP的三个输出能量
  ---------------------------------------------------
  free  energy TOTEN  =    -15.739127 eV

  energy  without entropy=    -15.380196  energy(sigma->0) =    -15.559661
有什么区别？

回复此楼

» 猜你喜欢

面上项目申报已经有3人回复
有时候真觉得大城市人没有县城人甚至个体户幸福已经有5人回复
酰胺脱乙酰基已经有9人回复
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优已经有5人回复
有70后还继续奋斗在职场上的吗？已经有6人回复
博士延得我，科研能力直往上蹿已经有7人回复
退学或坚持读已经有27人回复
面上基金申报没有其他的参与者成吗已经有5人回复
遇见不省心的家人很难过已经有22人回复

1楼 2009-07-21 09:52:18

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

hgye

金虫 (小有名气)

应助: 1 (幼儿园)
金币: 1406
帖子: 180
在线: 46.1小时
虫号: 461024
注册: 2007-11-18
性别: GG
专业: 凝聚态物性 II ：电子结构

★ ★ ★ ★
chenyang17(金币+2,VIP+0):这个回答很好，我回去试试对兄弟加我QQ 229954517 我想多请教一些 7-21 20:29
qasd(金币+2,VIP+0):xiexie 7-22 17:15

引用回帖:

Originally posted by chenyang17 at 2009-7-21 10:30:
castep和VASP用的方法基本是一样的都是基于赝势无论如何总能不会相差这么多吧

两者给出的能量差很多主要是因为VASP的能量已经减去了每个原子的能量，CASTEP没有减，如果给vasp的能量加上每个原子的能量，就会发现和castep的很接近。
即使这样，两者的结果还是不可以比的。

[ Last edited by hgye on 2009-7-21 at 17:35 ]

赞一下(13人)

回复此楼

12楼2009-07-21 17:34:03

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 14 个回答

xirainbow

木虫 (正式写手)

1ST强帖: 2
应助: 0 (幼儿园)
贵宾: 0.03
金币: 2132.6
红花: 2
帖子: 775
在线: 6.1小时
虫号: 431211
注册: 2007-08-11
专业: 光谱学

★ ★ ★
chenyang17(金币+1,VIP+0):castep和VASP用的方法基本是一样的都是基于赝势无论如何总能不会相差这么多吧 7-21 10:29
spur(金币+2,VIP+0): 3Q！！谢谢帮助！ 7-21 15:45

我觉得，任何两个软件的计算结果都是不可比的：）
“Technology of DF calculation is not as mature as that of HF calculations, and because of the variety of possible procedures used, DF calculations done with two different programs using the same Exc and the same basic set may not give quite same result.”

赞一下(12人)

回复此楼

2楼2009-07-21 10:21:53

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖