24小时热门版块排行榜    

查看: 2717  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 06:47 by 9fOTqj074tIg
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 06:36 by 9fOTqj074tIg
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 czAiolzR0IJC 2026-11-10 5/250 2026-11-12 06:36 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 06:25 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-12 06:06 by 9fOTqj074tIg
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 C6CRn66RYYXS 2026-11-10 7/350 2026-11-12 06:03 by 9fOTqj074tIg
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 C6CRn66RYYXS 2026-11-10 5/250 2026-11-12 06:01 by 9fOTqj074tIg
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 YHiw1gou1Eda 2026-11-10 5/250 2026-11-12 05:44 by 9fOTqj074tIg
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 05:30 by 9fOTqj074tIg
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 FChKX8meu1ds 2026-11-09 8/400 2026-11-12 05:28 by 9fOTqj074tIg
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 FChKX8meu1ds 2026-11-09 9/450 2026-11-12 05:19 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 05:19 by 9fOTqj074tIg
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +7 cqv1zKnzd7hr 2026-11-09 11/550 2026-11-12 05:08 by 9fOTqj074tIg
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 zSbvVX9JFl9J 2026-11-09 7/350 2026-11-12 04:49 by 9fOTqj074tIg
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-12 02:23 by 9fOTqj074tIg
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 C6CRn66RYYXS 2026-11-10 4/200 2026-11-12 01:28 by 9fOTqj074tIg
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-11 17:58 by TgvkbZZl9tqB
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 zSbvVX9JFl9J 2026-11-09 7/350 2026-11-11 16:41 by TgvkbZZl9tqB
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
信息提示
请填处理意见