【求助】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
有什么区别?
返回小木虫查看更多
京公网安备 11010802022153号
我觉得,任何两个软件的计算结果都是不可比的:)
“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.”
castep和VASP用的方法基本是一样的 都是基于赝势 无论如何 总能不会相差这么多吧
话虽这么说 但是我没深入研究过 不过常规说来 两种计算软件 计算方法 计算水平得出的结果是不比较的 只是说 这种方法和那种已经报道的方法对于这个体系得到的结果近似 (如果没有实验值 对照的话)
个人感觉不应该差别很大!如果很大的话,相信谁呢?
还有一个问题就是一直搞不清楚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
有什么区别?
同求这个问题,帮你顶帖子了,
任何计算给出的绝对能量都是没有可比性的。
通常孤立原子的赝势中给出的能量都要几十甚至是上百个eV。
关于Ti单胞的能量,castep中给出的结合能是多少?
VASP给出的能量的参考点是多少?
[ Last edited by xirainbow on 2009-7-21 at 16:18 ]