24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

skywyy2011

新虫 (正式写手)

应助: 1 (幼儿园)
金币: 364.7
红花: 3
帖子: 585
在线: 196.3小时
虫号: 1847776
注册: 2012-06-05
性别: GG
专业: 宇宙学

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

我在用WB97XD/6-311++G(d,p)算硼团簇物理吸附单个氢分子的时候，用物理吸附能公式：ΔEAE(Bn-H2) = E(Bn-H2) – E(Bn) – E(H2) 算得的吸附能居然都是正的（此处的E都是经过零点能校正的），自发的物理吸附能不应该都是负的，也即是放热么？同时，我在一篇计算团簇吸氢的文献中也发现其中的物理吸附能是正的，但是作者使用了障眼法，将上面的公式顺序掉了一下，结果就变成负的了，作者居然还说这个负号代表放热

，明显的错误最后审稿人居然没发现，还发表了。此外，我还发现一个现象，如果不经零点能校正，则用单点能计算得到的物理吸附能都是负的（这也就解释了为什么优化可以进行），但是一经过零点能校正，则物理吸附能都变成正的了。我的问题是：
1，为什么物理吸附能是正的？
2，同一个分子它的单点能势能面和经零点能校正的单点能势能面是一样的么？（我感觉不是，因为我发现有时同一分子两个不同构型，一个单点能比另一个大，但是经过零点能校正后，他的能量反而比另一个校正的能量小了。）
ps:为了防止有人说我的方法精度不够，我又用了G3的方法算一些小团簇的吸附，结果发现最后还是正的。
求各位大虾不吝赐教，在此先行谢过了

回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

天道酬勤

1楼 2013-09-22 10:04:17

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

专家经验: +224
应助: 47 (小学生)
金币: 3489.8
散金: 14
红花: 223
帖子: 2477
在线: 1253.7小时
虫号: 1369423
注册: 2011-08-16
性别: GG
专业: 金属材料的磨损与磨蚀
管辖: 第一性原理

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
skywyy2011: 金币+6, ★有帮助 2013-09-22 15:56:07

引用回帖:

13楼: Originally posted by skywyy2011 at 2013-09-22 15:41:36
我之前先用那个泛函算过硼的二聚体的一些参数，比如键长，键能啥的都和实验值贴合的比较好。如果用其他泛函的话，还有哪些泛函呢？...

具体哪个Functional比较好，这要看计算的体系和要计算的性质，我不做你这个体系我不是很清楚。看了小卒的东西，也是引用的http://hi.baidu.com/qidft/item/64a303b1bb092a452aebe3d6
tpssh处理硼团簇也许还不错，但是我没测试过。
你仅仅只能参考一下。

如果你的体系不大，有必要多测试一下。