24小时热门版块排行榜    

查看: 1847  |  回复: 15

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

[求助] 求H2,O2和H2O的bonding energies计算方法和参数(新虫参与符合条件奖励20金币)

我是这样计算的:

GGA+PW91, ENCUT = 480, Gamma point, EDIFFG = 1E-4, 把小分子放在10*10*10 cubic box中结构优化然后计算总能。再计算对应单原子总能。然后用减法算binding energy.如H2: E_bind=2*E_Hatom -E_H2.  但是我的H2, O2,H2O总能绝对值偏大,分别为-6.8,-9.81,-14.27eV。这样算出的binding energy一看就大于下面的文献值(ZPE不会超过1eV)。

现在求结果,看看大家怎么算的,结果怎么样?

请新虫拿这个做练习,算的结果与下面文献接近(个位数相同)的,跟帖后我自掏腰包奖励20金币。


[ Last edited by cenwanglai on 2011-10-28 at 13:45 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

谢谢光临,如果能回复一二就更感谢了哈!
2楼2011-10-31 14:47:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)

【答案】应助回帖


cenwanglai(金币+2): 谢谢。我算完后把能量差贴出来。 2011-11-01 07:21:43
贺仪(金币+1): 多谢指教! 2011-11-02 21:07:02
能量绝对值没有意义,只有能量之差才有意义。你继续算单个原子能量,不会有问题的。
我以前算过H2分子的解离能,就是2个H原子的能量减去H2分子的能量,和实验很接近的。
3楼2011-11-01 04:56:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gougou531

禁虫 (正式写手)

cenwanglai(金币+1): 谢谢,呵呵 2011-11-02 18:11:53
本帖内容被屏蔽

4楼2011-11-02 17:48:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

引用回帖:
4楼: Originally posted by gougou531 at 2011-11-02 17:48:34:
计算单个原子要考虑自旋的话,计算出来的结合能要更接近实验值。区长试试

能否讲一下理由,为什么要考虑呢?谢谢!
嘟啦啦嘟啦啦嘟
5楼2011-11-02 21:08:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gougou531

禁虫 (正式写手)

本帖内容被屏蔽

6楼2011-11-03 08:37:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主

不知道这篇文章中的binding energy是怎么定义的,其单位是eV/mol,难道是计算出来的能量值再乘以阿佛加德罗常数?
7楼2011-11-03 09:10:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

引用回帖:
7楼: Originally posted by uuv2010 at 2011-11-03 09:10:40:
不知道这篇文章中的binding energy是怎么定义的,其单位是eV/mol,难道是计算出来的能量值再乘以阿佛加德罗常数?

是eV/molecule。

呵呵!
8楼2011-11-03 09:29:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Illusionist

银虫 (正式写手)

果然,加了自选后O2的小了很多,我算了一下-9.861268-()=-6.204434。
但是还是相差比较远。
H2算的        -6.729082-2*(-1.109629)=-4.509824
9楼2011-11-03 16:56:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

引用回帖:
9楼: Originally posted by Illusionist at 2011-11-03 16:56:11:
果然,加了自选后O2的小了很多,我算了一下-9.861268-()=-6.204434。
但是还是相差比较远。
H2算的        -6.729082-2*(-1.109629)=-4.509824

请选择应助 然后随便发点什么,跟你金币。
10楼2011-11-03 17:10:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cenwanglai 的主题更新
信息提示
请填处理意见