当前位置: 首页 > 第一原理 >【求助】castep与vasp的能量计算值有可比性吗?

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

作者 chenyang17
来源: 小木虫 650 13 举报帖子
+关注

分比用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
有什么区别? 返回小木虫查看更多

今日热帖
  • 精华评论
  • xirainbow

    我觉得,任何两个软件的计算结果都是不可比的:)
    “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.”

  • chenyang17

    castep和VASP用的方法基本是一样的 都是基于赝势 无论如何 总能不会相差这么多吧

  • topten-gg

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

    话虽这么说 但是我没深入研究过 不过常规说来 两种计算软件 计算方法 计算水平得出的结果是不比较的  只是说 这种方法和那种已经报道的方法对于这个体系得到的结果近似 (如果没有实验值 对照的话)

  • 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,参数基本都相同,为什么算出来的能量差别那么大,大家看 ...

    个人感觉不应该差别很大!如果很大的话,相信谁呢?

  • 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
    有什么区别?

    同求这个问题,帮你顶帖子了,

  • xirainbow

    引用回帖:
    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 ]

猜你喜欢
应助之星
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓