24小时热门版块排行榜    

查看: 2381  |  回复: 36
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

skywyy2011

新虫 (正式写手)

[求助] 团簇吸附氢分子的吸附能居然是正的?

我在用WB97XD/6-311++G(d,p)算硼团簇物理吸附单个氢分子的时候,用物理吸附能公式:ΔEAE(Bn-H2) = E(Bn-H2) – E(Bn) – E(H2) 算得的吸附能居然都是正的(此处的E都是经过零点能校正的),自发的物理吸附能不应该都是负的,也即是放热么?同时,我在一篇计算团簇吸氢的文献中也发现其中的物理吸附能是正的,但是作者使用了障眼法,将上面的公式顺序掉了一下,结果就变成负的了,作者居然还说这个负号代表放热,明显的错误最后审稿人居然没发现,还发表了。 此外,我还发现一个现象,如果不经零点能校正,则用单点能计算得到的物理吸附能都是负的(这也就解释了为什么优化可以进行),但是一经过零点能校正,则物理吸附能都变成正的了。我的问题是:
1,为什么物理吸附能是正的?
2,同一个分子它的单点能势能面和经零点能校正的单点能势能面是一样的么?(我感觉不是,因为我发现有时同一分子两个不同构型,一个单点能比另一个大,但是经过零点能校正后,他的能量反而比另一个校正的能量小了。)
ps:为了防止有人说我的方法精度不够,我又用了G3的方法算一些小团簇的吸附,结果发现最后还是正的。
求各位大虾不吝赐教,在此先行谢过了
回复此楼
天道酬勤
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

skywyy2011

新虫 (正式写手)

引用回帖:
20楼: Originally posted by smthnb at 2013-09-22 16:45:59
你把吸附构型中的硼团簇再拿出来优化看看,看看有没有变化啊...

这个肯定会变化啊,变化到未吸附氢的最稳定构型啊,只不过这种变化很小
天道酬勤
22楼2013-09-22 16:56:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 37 个回答

smthnb

禁虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
skywyy2011: 金币+6, 有帮助 2013-09-22 14:43:41
gmy1990: 金币+2 2013-09-22 18:02:32
本帖内容被屏蔽

2楼2013-09-22 11:34:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yypcnbe

木虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
你考虑过自旋多重态没?比如铁吸附氢气分子就要考虑,不然得出的结果就是吸附能是正值。
3楼2013-09-22 13:36:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

smthnb

禁虫 (正式写手)

本帖内容被屏蔽

4楼2013-09-22 13:41:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见