24小时热门版块排行榜    

查看: 3789  |  回复: 6

rik203

金虫 (小有名气)

[求助] 求用振动频率算零点能

MgH2零点能其他人计算值一般是0.39ev,我计算的四个文件和outcar都在压缩包里.零点能公式E=(∑hωj)/2,然后就是那个ωj我不知道取哪一个频率值,是1f对应的那个吗?还是1f,2f...等对应的频率值求和?怎么算都不对。求解。
回复此楼
岂能尽如人意,但求竭尽所能
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)

【答案】应助回帖

★ ★
cenwanglai(金币+2): 谢谢回复呵呵,早啊! 2011-11-01 07:26:10
rik203(金币+4): 18个频率算出来的零点能就den打咯。氢气是怎么算零点能的:http://emuch.net/html/201103/2880954.html 2011-11-02 09:41:36
rik203(金币+4): 18个频率算出来的零点能就den打咯。氢气是怎么算零点能的:http://emuch.net/html/201103/2880954.html 2011-11-02 09:41:46
原胞6个原子,共有18个本征频率。OUTCAR也给出了18个频率。最后三个是声学波对于的频率,在gamma点为0。之所以不严格为0,是数值误差造成的。
零点能就是这18个频率相加再乘以 h/2.
2楼2011-11-01 05:03:49
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

【答案】应助回帖

rik203(金币+3): 还是不太懂怎么算最后的零点能。。。。 2011-11-02 09:48:56
3楼2011-11-01 07:25:53
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)

干嘛这么费劲啊,直接加能量不可以么
4楼2011-11-01 15:16:47
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

【答案】应助回帖

rik203(金币+1): 很细心 2011-12-20 11:31:53
引用回帖:
4楼: Originally posted by hakuna at 2011-11-01 15:16:47:
干嘛这么费劲啊,直接加能量不可以么

对啊。不过要注意,OUTCAR中给出的能量需要除以2。把所有实频对应的能量除以2然后加和。
5楼2011-11-02 09:51:37
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jiaonahappy

金虫 (小有名气)

引用回帖:
3240672楼: Originally posted by rik203 at 2011-10-31 09:52:09:
MgH2零点能其他人计算值一般是0.39ev,我计算的四个文件和outcar都在压缩包里.零点能公式E=(∑hωj)/2,然后就是那个ωj我不知道取哪一个频率值,是1f对应的那个吗?还是1f,2f...等对应的频率值求和?怎么算都不对 ...

我想请问楼主,你的INCAR设置是进行relax吧?可是零点能计算不是应该得到了体系的波函数(即体系的基态)之后再进行计算的吗?(可能我理解有问题,如有不对之处请不吝赐教,谢谢。)
踏踏实实的走好每一步
6楼2012-04-12 09:46:17
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

刘小伟110

木虫 (小有名气)

引用回帖:
2楼: Originally posted by stractor at 2011-11-01 05:03:49
原胞6个原子,共有18个本征频率。OUTCAR也给出了18个频率。最后三个是声学波对于的频率,在gamma点为0。之所以不严格为0,是数值误差造成的。
零点能就是这18个频率相加再乘以 h/2.

stractor好,
   我想问一下在计算零点能时,不是有个零点能的计算公式ZPE还有个3N-6什么的,为什么不用那个算啊,难道vasp中计算的都是振动频率吗
  还望您不吝赐教,谢谢啦
慢慢的思考。。。
7楼2013-11-07 09:28:22
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 rik203 的主题更新
信息提示
请填处理意见