| 24小時(shí)熱門版塊排行榜 |
| 查看: 4174 | 回復(fù): 13 | |||
20928140木蟲 (正式寫手)
|
[求助]
高斯計(jì)算過(guò)渡態(tài)時(shí)的能壘總是負(fù)值
|
||
|
自己現(xiàn)在正在計(jì)算CF2ClCFClH和HO^-1發(fā)生E2消除反應(yīng)的過(guò)渡態(tài)(此反應(yīng)過(guò)程應(yīng)該是基元反應(yīng),可以認(rèn)為是強(qiáng)堿拔氫)。可是自己計(jì)算出來(lái)的能壘一直是負(fù)值,不知道該如何處理。過(guò)渡態(tài)有唯一虛頻且較大,用IRC掃描確認(rèn)無(wú)誤。計(jì)算方法如下: 剛開始用的是HF/6-31+G(d,p),發(fā)現(xiàn)能壘是-122.8Kj/mol,采用B3LYP/6-31+G(d,p),也是負(fù)的。計(jì)算方法是采用E(TS)-E(CF2ClCFClH)-E(HO^-1)得到的。 哪位前輩指點(diǎn)一下。繛槭裁词秦(fù)值? 謝謝! |
計(jì)算化學(xué) |

鐵桿木蟲 (著名寫手)
|
我汗,樓上某位專家貌似沒(méi)看清楚樓主的帖子啊 樓主說(shuō)的是用 E(TS)-E(CF2ClCFClH)-E(HO^-1) 算出來(lái)的能壘,這樣算顯然是不對(duì)的... 應(yīng)該用 E(TS)-E[(CF2ClCFClH)+(HO^-1)] 一般來(lái)說(shuō) E(CF2ClCFClH)-E(HO^-1) 是 大于 E[(CF2ClCFClH)+(HO^-1)] 的。 簡(jiǎn)單點(diǎn)的做法:把你TS的IRC,朝向反應(yīng)物的那邊,的結(jié)果,拿出來(lái)再opt一次;最后看看是不是 (CF2ClCFClH)在(HO^-1)附近的樣子,這個(gè)能量就是 E[(CF2ClCFClH)+(HO^-1)] 。 順便吐槽下某專家,我看了你不少帖子了,一上來(lái)不管樓主什么問(wèn)題,就是一大堆專業(yè)名詞刷出...具體問(wèn)題具體分析好吧... |
至尊木蟲 (著名寫手)
我叫MT

至尊木蟲 (著名寫手)
我叫MT
|
說(shuō)得挺好的,也提醒了我忘了BSSE校正也可能引起樓主的問(wèn)題。那就用counterpoise做一下BSSE校正好了,不過(guò)一般BSSE不會(huì)太大,2-3kcal/mol就很了不起了吧。樓主雖然沒(méi)說(shuō)能壘是多少,但我估計(jì)不會(huì)這么小的。但是給我提醒一個(gè)思路,還是感謝你的。 我學(xué)gaussian才不過(guò)2年半,這期間曾得到版面上好多人的指導(dǎo),重要的不僅僅是個(gè)別問(wèn)題的答案,而是怎樣思考問(wèn)題。所以有人問(wèn)1+1=幾的時(shí)候,我愿意給出怎么解題的思路,而不是直接回答2.在高手看來(lái)是顯擺,在我自身學(xué)習(xí)的過(guò)程體會(huì),這比答案重要。何況我所學(xué)不精,我不把我的錯(cuò)誤擺出來(lái),怎么有機(jī)會(huì)得到高手的糾正呢。。他人遇到的問(wèn)題,我也會(huì)遇到,我如果一直用錯(cuò)誤的方式解決而不暴露出來(lái),怎么能得到教訓(xùn)呢。我也反感那個(gè)專家顧問(wèn)的頭銜,我只是勇于展示自己的錯(cuò)誤想法,期望得到高手指正的機(jī)會(huì)。你的吐槽不會(huì)讓我因此而不敢露怯,我會(huì)繼續(xù)宣揚(yáng)我的思路,以期盼有高手看不過(guò)去,進(jìn)而來(lái)糾正我。 |

木蟲 (正式寫手)

至尊木蟲 (著名寫手)
我叫MT
|
1、上文第一條意指你可能找到的不是全局最小點(diǎn),而且你這里還有一個(gè)手性C(手性C不影響能量),這個(gè)可能占80% 2、基組的選擇,以及方法的選擇是否恰當(dāng),這個(gè)可能占10% 3、第三條是指你是否考慮了熱力學(xué)內(nèi)能,是否做了零點(diǎn)能校正,這個(gè)可能占10% 4、極少情況下,實(shí)驗(yàn)測(cè)的活化能可以是負(fù)值(多發(fā)生在生化反應(yīng)中),估計(jì)你這個(gè)反應(yīng)不在此列中。 |

木蟲 (正式寫手)

木蟲 (正式寫手)

至尊木蟲 (著名寫手)
我叫MT
|
原子數(shù)多?8個(gè)原子也叫多,那80個(gè)原子怎么辦? 先說(shuō)你這個(gè)基組下的基函數(shù)一共多少個(gè)吧。如果只有300-500個(gè),那就太少了。不超過(guò)2000個(gè),應(yīng)該是比較合適的——計(jì)算量是按照基函數(shù)的4次方遞增的。 再說(shuō)校正,我說(shuō)的是零點(diǎn)能校正,你說(shuō)的校正因子是什么?是頻率的校正因子么? 再說(shuō)方法,為什么用HF不用DFT呢?HF不考慮電子相關(guān),誤差偏大。 再說(shuō)能量,這是個(gè)非常模糊的概念,你必須清楚地指明是什么能量,如你所說(shuō)的E(RHF),僅僅是HF方法算出的電子的勢(shì)能,你沒(méi)有考慮熱力學(xué)能——平動(dòng)、振動(dòng)、轉(zhuǎn)動(dòng)的能量,沒(méi)有做零點(diǎn)能校正(ZPE),是不正確的,一般應(yīng)該以包括零點(diǎn)能校正的焓來(lái)計(jì)算。 再說(shuō)體系,你以O(shè)H-為反應(yīng)物之一,顯然應(yīng)該考慮溶劑化效應(yīng),你的反應(yīng)物分子極性都偏大,而產(chǎn)物極性偏小,在極性溶劑中,反應(yīng)物占優(yōu)勢(shì),在非極性溶劑中,產(chǎn)物占優(yōu)勢(shì)。你的計(jì)算需要在溶劑環(huán)境下重新做,至于采用顯性模型還是隱性模型,尚需根據(jù)經(jīng)驗(yàn)探討。 再說(shuō)初始構(gòu)型以及全局與局域最小點(diǎn),這個(gè)也是化學(xué)的經(jīng)驗(yàn),初始構(gòu)象合理即可找到合理的最小點(diǎn),不合理的初始構(gòu)象是找不到全局最小點(diǎn)的。這么簡(jiǎn)單的分子,能存在的構(gòu)象并不多,你多試幾個(gè)構(gòu)象,很快會(huì)找到全局最小點(diǎn)的。復(fù)雜的分子就很難了,要憑經(jīng)驗(yàn)。。。 |

鐵桿木蟲 (著名寫手)
木蟲 (正式寫手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +10 | Faune 2026-03-21 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +3 | vv迷 2026-03-21 | 4/200 |
|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分,一志愿四川大學(xué),誠(chéng)心求調(diào)劑 +11 | 吃吃吃才有意義 2026-03-19 | 11/550 |
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
[基金申請(qǐng)]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +6 | 梨花珞晚風(fēng) 2026-03-17 | 6/300 |
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +5 | tcxiaoxx 2026-03-20 | 5/250 |
|
|
[考研] 295材料求調(diào)劑,一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語(yǔ) 2026-03-15 | 5/250 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|