应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 【求助】castep与vasp的能量计算值有可比性吗?
141/2返回列表12下一页
查看: 2233  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

chenyang17

木虫 (小有名气)

[交流] 【求助】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
有什么区别？
回复此楼

» 猜你喜欢
1楼 2009-07-21 09:52:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

hgye

金虫 (小有名气)
★ ★ ★ ★
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的回帖
普通回帖

xirainbow

木虫 (正式写手)
★ ★ ★
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的回帖

chenyang17

木虫 (小有名气)
castep和VASP用的方法基本是一样的 都是基于赝势 无论如何 总能不会相差这么多吧
赞一下 回复此楼
3楼2009-07-21 10:30:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

topten-gg

银虫 (正式写手)
★ ★
chenyang17(金币+1,VIP+0): 7-21 13:45
zxzj05(金币+1,VIP+0):3Q! ^_^ 鼓励讨论！ 7-21 16:20
引用回帖:
Originally posted by chenyang17 at 2009-7-21 10:30:
castep和VASP用的方法基本是一样的 都是基于赝势 无论如何 总能不会相差这么多吧

话虽这么说 但是我没深入研究过 不过常规说来 两种计算软件 计算方法 计算水平得出的结果是不比较的  只是说 这种方法和那种已经报道的方法对于这个体系得到的结果近似 （如果没有实验值 对照的话）
赞一下(12人) 回复此楼
4楼2009-07-21 10:41:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

woshilaoer8381

至尊木虫 (知名作家)
引用回帖:
Originally posted by chenyang17 at 2009-7-21 09:52:
分比用CASTEP和VASP同样计算一个Ti单胞的能量，采用相同的结构，a=2.9506， c=4.6788，其中VASP计算时用的poscar是从.cell文件里面导出来的，K点都取9*9*2，参数基本都相同，为什么算出来的能量差别那么大，大家看 ...

个人感觉不应该差别很大！如果很大的话，相信谁呢？
赞一下 回复此楼
5楼2009-07-21 14:13:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dali2006

木虫 (著名写手)

小木虫酱油店店主
没有
回复此楼
6楼2009-07-21 14:38:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Wanghui6383

银虫 (正式写手)
还有一个问题就是一直搞不清楚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
有什么区别？

同求这个问题，帮你顶帖子了
赞一下 回复此楼
愿意结识从事第一原理研究的学者！
7楼2009-07-21 14:51:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)
★ ★ ★
zxzj05(金币+1,VIP+0):3Q! ^_^ 鼓励讨论！ 7-21 16:20
chenyang17(金币+2,VIP+0):我多研究研究这个问替 看看两种算法结合能的差别 7-21 20:30
引用回帖:
Originally posted by chenyang17 at 2009-7-21 09:52:
分比用CASTEP和VASP同样计算一个Ti单胞的能量，采用相同的结构，a=2.9506， c=4.6788，其中VASP计算时用的poscar是从.cell文件里面导出来的，K点都取9*9*2，参数基本都相同，为什么算出来的能量差别那么大，大家看 ...

任何计算给出的绝对能量都是没有可比性的。
通常孤立原子的赝势中给出的能量都要几十甚至是上百个eV。
关于Ti单胞的能量，castep中给出的结合能是多少？
VASP给出的能量的参考点是多少？

[ Last edited by xirainbow on 2009-7-21 at 16:18 ]
赞一下(12人) 回复此楼
8楼2009-07-21 16:15:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shenxianshifu

至尊木虫 (知名作家)

宗师
★
zxzj05(金币+1,VIP+0):3Q! ^_^ 鼓励讨论！ 7-21 16:20
没有可比性，因为用的是不同的近似机制，但如果有实验数据，两者具有可比性
赞一下(10人) 回复此楼
9楼2009-07-21 16:15:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxzj05

荣誉版主 (著名写手)
★
wuli8(金币+1,VIP+0):3ks 7-22 18:19
没有可比性，深入了解了程序运行计算近似原理就明白了。
只有做了差的形成能之类才有可比性。
赞一下(20人) 回复此楼
储氢家族欢迎储氢研究者！
10楼2009-07-21 16:22:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chenyang17 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭