24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1697  |  回复: 13
当前主题已经存档。

majun04

木虫 (著名写手)

majun

[交流] 【讨论】计算反应路径HF方法和DFT方法比较。

想通过计算来解释反应机理,描述催化过程。
首先计算空白反应(不加催化剂)路径的能垒,然后计算引入催化剂以后的反应路径的能垒。
使用方法和基组分别是HF/6-311++G(d,p)和B3LYP/6-311++G(d,p),对后得到结果相比较,发现用HF方法得到的结构和B3LYP方法得到的结构是差不多的,但是能量上确相差好大。在空白反应路径和加了催化以后的反应路径上,HF方法得到的能垒要明显高于B3LYP得到的能垒。请问这个什么原因?
认为,催化剂的作用在于把能垒降低了多少。那么就算了一下,催化剂的引入使得能垒下降的数值。结果显示:用HF方法计算,和用B3LYP方法计算,能垒下降的数值是差不多的。请问这样能说明本反应体系选择HF方法和DFT方法得到的结果是相同的吗?
回复此楼
乐观、执着、勤奋
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

majun04

木虫 (著名写手)

majun

★ ★ ★
lei0736(金币+3,VIP+0):谢谢 鼓励发起讨论 5-13 20:20
在一篇文献中看到这样的话:
"Energetically, the calculated activation energy of the hydration of 1 via transition structure 2 to from 3, is 69.43 kcal/mol at HF, but noly 49.00 kcal/mol at MP2. This large difference between the calculated HF and MP2 results indicateds that the correlation energy was enhanced in the calculation for the structure involving new bond fomation or destruction." J. Phys. Chem. A 2001, 105, 6543-6551
请问这个MP2方法和HF方法为什么在计算有成键断键的过渡态,结果差别这么大?现在好多计算HF方法都不用了,这个是不是就是HF被淘汰的原因?

我的体系有弱的分子相互作用,但是用MP2方法太昂贵了,根本就玩儿不起,所以用DFT方法计算的(这个方法在计算范德华力作用体系有缺陷)。补充HF方法的结果后,也是发现HF方法计算出来的能垒是明显的高于DFT的结果。
请问这个是为什么?

[ Last edited by majun04 on 2009-5-13 at 19:36 ]
乐观、执着、勤奋
2楼2009-05-13 19:34:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

313588908

木虫 (正式写手)

★ ★ ★
majun04(金币+1,VIP+0):谢谢你 5-13 22:24
lei0736(金币+2,VIP+0):谢谢 5-14 20:42
好像是因为HF考虑没有电子相关能,过渡态的能量一般不用HF算了。DFT方法还是比较可信的
3楼2009-05-13 22:22:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

majun04

木虫 (著名写手)

majun


loovfnd(金币+1):谢谢参与!
引用回帖:
Originally posted by 313588908 at 2009-5-13 22:22:
好像是因为HF考虑没有电子相关能,过渡态的能量一般不用HF算了。DFT方法还是比较可信的

投稿的时候是用DFT方法算得,但是审稿人提出质疑,说DFT在算弱相关作用的时候是有缺陷的,叫补充HF和MP2方法的结果。
HF的是算了,差别比较大,MP2的太耗用机时了,就算了个别的点,感觉这个MP2方法的结果和DFT的还比较接近。
挺郁闷的,好多文献都在用DFT方法算反应路径了,也是多个分子反应的,有的还有氢键作用,怎么我这个就不行了,咳。
乐观、执着、勤奋
4楼2009-05-13 22:27:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shalene

至尊木虫 (著名写手)

★ ★
loovfnd(金币+1):谢谢参与!
majun04(金币+1,VIP+0):谢谢你,但是对于大的体系用MP2算不起,怎么办?我现在面临的困难就是这个,MP2太昂贵了! 5-13 23:25
DFT在算长程弱相关作用的时候是有缺陷的,短程时可以接受
5楼2009-05-13 23:01:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

313588908

木虫 (正式写手)

★ ★
loovfnd(金币+1):谢谢参与!
majun04(金币+1,VIP+0):谢谢 5-13 23:25
发文章算弱相关作用还是MP2比较站得住脚,只能多花点时间了。HF都没怎么见到了。
6楼2009-05-13 23:02:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

313588908

木虫 (正式写手)


lei0736(金币+1,VIP+0):谢谢 5-14 20:42
而且HF与DFT,MP2能量相差都很大的,只做一定参考吧
7楼2009-05-13 23:04:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

majun04

木虫 (著名写手)

majun

引用回帖:
Originally posted by 313588908 at 2009-5-13 23:04:
而且HF与DFT,MP2能量相差都很大的,只做一定参考吧

能讲解一下为什么差别会比较大吗?本质上是什么原因,为什么他算的能垒数值就相对高?
乐观、执着、勤奋
8楼2009-05-13 23:24:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shalene

至尊木虫 (著名写手)

★ ★ ★
majun04(金币+1,VIP+0):谢谢你,不想换软件了,改动太大了,文章会重新送审的,已经延期了,想着就发愁。 5-13 23:32
lei0736(金币+2,VIP+0):谢谢 5-14 20:43
建议用turbomole里的ri-mp2,快很多。
9楼2009-05-13 23:28:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shalene

至尊木虫 (著名写手)

★ ★ ★ ★
majun04(金币+1,VIP+0):咳! 5-13 23:47
lei0736(金币+3,VIP+0):谢谢 5-14 20:43
已经延期了,还愁啥?

用turbomole里的ri-dft和ri-mp2把你的体系全部重算一遍,计算的花费也不会比仅用gaussian里的dft算大。

从实际计算的角度看,ri是个好的近似方法,虽然缺少物理意义,但是能在保持足够精度情况下,提升好几倍甚至一个数量级的速度!

不过turbomole是商业软件,不好找。替代品有免费的ORCA,也包含ri方法,不过我没用过。
10楼2009-05-13 23:45:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 majun04 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 318求调剂 +12 ykyhsa 2026-04-05 14/700 2026-04-06 17:46 by fuyu_
[考研] 生物学求调剂 +5 15064154688 2026-04-03 5/250 2026-04-06 11:56 by lijunpoly
[考研] 一志愿武汉理工大学080200机械工程308分,求调剂 +4 终不似从前 2026-04-05 4/200 2026-04-06 11:46 by 考研学校招点人
[考研] 求调剂,一志愿厦门大学,生物与医药,总分272,本科211 +4 Electron1cc 2026-04-01 5/250 2026-04-06 10:45 by barlinike
[考研] 0703调剂 +13 拾玖壹 2026-04-04 15/750 2026-04-06 10:35 by 蓝云思雨
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +10 大火山小火山 2026-04-05 10/500 2026-04-05 19:58 by 小蜗牛小蜗牛
[考研] 283分求调剂 +9 试试看呗 2026-04-04 9/450 2026-04-05 10:27 by 果冻大王
[考研] 286求调剂 +3 草木不言 2026-04-04 3/150 2026-04-04 22:40 by lbsjt
[考研] 321求调剂 +6 认真求上学 2026-04-03 6/300 2026-04-04 19:51 by dongzh2009
[考研] 342求调剂 +3 Liang7111 2026-04-04 5/250 2026-04-04 19:47 by dongzh2009
[考研] 一志愿华南师范361分,化学求调剂 +7 Nicole88888 2026-04-01 7/350 2026-04-04 18:28 by macy2011
[考研] 考研调剂 +3 Draa 2026-04-03 3/150 2026-04-03 17:37 by hgwz7468
[考研] 专硕085601求调剂 +7 suyifei 2026-04-03 8/400 2026-04-03 14:00 by 欣喜777
[考研] 338求调剂,一志愿能源动力,外语是日语203 +5 zzz,,r 2026-04-02 5/250 2026-04-03 09:45 by 蓝云思雨
[考研] 309求调剂 +14 呆菇不是戴夫 2026-04-02 14/700 2026-04-03 09:42 by 蓝云思雨
[考研] 346求调剂 +5 郑诚乐 2026-04-02 5/250 2026-04-02 16:38 by SZW_UJN
[考研] 279求调剂 +6 学而思兮知 2026-04-01 6/300 2026-04-02 09:16 by vgtyfty
[考研] 材料调剂 +12 一样YWY 2026-04-01 12/600 2026-04-02 00:21 by 百秒光年
[考研] 材料调剂 +11 一样YWY 2026-03-31 11/550 2026-04-01 22:25 by zhouyuwinner
[考研] 292求调剂 +17 木虫er12138 2026-04-01 17/850 2026-04-01 21:37 by 七度不信任
信息提示
请填处理意见