| 24小時熱門版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 2241 | 回復(fù): 6 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
swfctheochem銅蟲 (初入文壇)
|
[求助]
一個奇怪的計(jì)算結(jié)果--大家?guī)涂纯矗?
|
|
|
我對甲二醇和甲醛的相互轉(zhuǎn)化反應(yīng)進(jìn)行了計(jì)算: CH2(OH)2 →CHO+H2O 有機(jī)化學(xué)常識告訴我:稀甲醛的水溶液中水合物甲二醇占到99%以上,而甲醛只占不到1%。我用理論方法算該反應(yīng)的平衡常數(shù):B3LYP/6-31+G**, PCM溶劑化模型。 優(yōu)化結(jié)構(gòu)后,用gaussian09輸出文件中Sum of electronic and thermal Free Energies計(jì)算反應(yīng)的吉布斯自由能變,發(fā)現(xiàn)反應(yīng)的ΔG竟然為負(fù)值(-23kJ/mol),用ΔG=-RTlnK計(jì)算得到的K是一個很大的正值,這豈不是說明甲二醇轉(zhuǎn)化為甲醛在熱力學(xué)上很有利,這明顯與實(shí)驗(yàn)事實(shí)相反。后來我發(fā)現(xiàn)甲二醇的Thermal correction to Gibbs Free Energy很大,為0.031399Hartree,而CHO和H2O的校正值很小,只有0.008633。如果不進(jìn)行Thermal correction校正,該反應(yīng)的ΔG為正的36kJ/mol,這個結(jié)果就相對合理。請問一下對反應(yīng)物和產(chǎn)物的Thermal correction為什么會差別這么大?是我的計(jì)算方法有問題嗎?gaussian輸出的熱力學(xué)量的溫度應(yīng)該是298.15K吧?對這樣一個簡單反應(yīng)gaussian都稿不定,以后還能相信gaussian嗎?我不要求定量,但如果是這樣的話,連定性都不行。 我是新蟲,只有6個金幣,全部奉上了!望蟲友們幫幫忙! |
木蟲 (正式寫手)
木蟲 (正式寫手)
銅蟲 (初入文壇)
|
謝謝樓上的 提醒! 甲醛水溶液中除甲二醇和甲醛外確實(shí)還有其他物質(zhì),比如說半縮醛多聚甲醛。但稀溶液中甲二醇所占比例確實(shí)比甲醛高很多。那么從熱力學(xué)平衡來理解,甲二醇轉(zhuǎn)化為甲醛的反應(yīng)在熱力學(xué)上是不利的,平衡常數(shù)應(yīng)小于1,ΔG應(yīng)大于零。逆反應(yīng),即甲醛水合為甲二醇是有利的,在水溶液中水合甲醛的存在形式比甲醛更為穩(wěn)定,這一點(diǎn)不知樓上是否同意? 另外,我發(fā)現(xiàn):如果將產(chǎn)物看做甲醛和水分子的氫鍵絡(luò)合物,而不要將其視為單獨(dú)的甲醛和水分子,那么算出來的ΔG就大于零。在真實(shí)溶液中即便生成了甲醛,它也是被水分子包圍的,也就是說甲二醇脫水分解為單獨(dú)的甲醛和水分子的反應(yīng)在真實(shí)的溶液中是不存在的。如果我的理解是正確的,那么就可以理解ΔG小于零的原因:甲醛和水分子的氫鍵絡(luò)合物理解為甲醛和水是一個吉布斯自由能降低的過程,而我在計(jì)算中直接將產(chǎn)物視為單獨(dú)的甲醛和水,實(shí)際上是包含了這個解離過程,導(dǎo)致整個反應(yīng)的吉布斯自由能是降低的。就計(jì)算結(jié)果來講,單獨(dú)的甲醛和水,其Thermal correction to Gibbs Free Energy很。0.031399),如果按絡(luò)合物處理這一項(xiàng)為0.022302. 除了質(zhì)子直接轉(zhuǎn)移脫水的機(jī)理,我還計(jì)算了水助質(zhì)子轉(zhuǎn)移脫水的機(jī)理,在水助機(jī)理中,如果將產(chǎn)物看做甲醛與3個水分子的絡(luò)合物,發(fā)現(xiàn)吉布斯自由能進(jìn)一步增大,這一點(diǎn)是不是印證了我的觀點(diǎn)?不知樓上是否同意以上觀點(diǎn)?歡迎指正!交流! |
|
高斯本身并無問題 有問題的往往是使用高斯的人 我看到你的計(jì)算的幾個問題 1如果我是你 對于這么小的分子 我不會去用dft的 怎么也要采用ccsdt了 2而且 水是一種強(qiáng)極性物種 而且水分子與溶質(zhì)分子之間有強(qiáng)烈相互作用 那么你所用的隱式溶劑模型就十分嚴(yán)重的不合適了 如果是我 我會采用顯隱混合溶劑模型 3你這里面必然涉及到弱相互作用(不要小看它們,弱相互作用盡管只有幾個到十幾個kj/mol,但是集腋成裘,所起的作用也很大,尤其是在化學(xué)反應(yīng)過程中),但是B3LYP早就被證明是計(jì)算弱相互作用嚴(yán)重不合適的泛函,而且你的基組也偏小。關(guān)于方法和基組的選擇,量化版各位高手早已總結(jié)過了,見http://www.gaoyang168.com/bbs/viewthread.php?tid=1566146這個帖子的第2.7節(jié) 從http://www.gaoyang168.com/bbs/viewthread.php?tid=2405321中,在關(guān)于什么是好問題,什么是壞問題的章節(jié)中,我節(jié)選幾句話: “蠢問題:我的vasp無法編譯。它怎么這么爛? 他覺得都是別人的錯,這個傲慢自大的家伙 聰明問題:我的vasp無法編譯。我讀過了FAQ,但是我仍舊沒能解決。這是我編譯過程的記錄,我有什么做得不對的地方嗎? 他講明了環(huán)境,也讀過了FAQ,還指明了錯誤,并且他沒有把問題的責(zé)任推到別人頭上,這個家伙值得留意! 基本上,恕我直言,從計(jì)算方法、計(jì)算思路,樓主無一可取,請莫怪我直言為好。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 本人考085602 化學(xué)工程 專碩 +16 | 不知道叫什么! 2026-03-15 | 18/900 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項(xiàng)國獎兩個實(shí)用型專利一項(xiàng)省級立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 0854可跨調(diào)劑,一作一項(xiàng)核心論文五項(xiàng)專利,省、國級證書40+數(shù)一英一287 +3 | 小李0854 2026-03-16 | 3/150 |
|
|
[考研] 211本,11408一志愿中科院277分,曾在中科院自動化所實(shí)習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 機(jī)械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
[考研]
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
[考研] 331求調(diào)劑(0703有機(jī)化學(xué) +5 | ZY-05 2026-03-13 | 6/300 |
|
|
[考研] 四川大學(xué)085601材料工程專碩 初試294求調(diào)劑 +4 | 祝我們好在冬天 2026-03-11 | 4/200 |
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
[考研] 314求調(diào)劑 +7 | 無懈可擊的巨人 2026-03-12 | 7/350 |
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|