| 查看: 2247 | 回復(fù): 10 | |||
gmy1990榮譽(yù)版主 (著名寫手)
|
[求助]
H2鍵解離和Ar二聚體兩個(gè)原子間的弱作用計(jì)算
|
|
采用Gaussian 03 H2的鍵解離過(guò)程應(yīng)該用非限制性方法計(jì)算,但是為什么得到的和限制性方法的結(jié)果是一樣的呢?MP2計(jì)算出來(lái)的曲線怎么又會(huì)如此奇葩? 對(duì)于Ar dimer的計(jì)算的話,勢(shì)能曲線應(yīng)該是在4左右的位置有一個(gè)最低點(diǎn)(Figure 3),而不是現(xiàn)在的曲線形狀(Figure 2)。 @sobereva @zhou2009 @yjcmwgk @beefly Figure1.jpg Figure2.jpg Figure3.jpg |
專家顧問(wèn) (職業(yè)作家)
地溝油冶煉專家
![]() |
專家經(jīng)驗(yàn): +458 |

榮譽(yù)版主 (著名寫手)
|
也不行 按道理是hf的結(jié)果就是隨著兩原子距離增加而能量一直遞減。還是無(wú)法得到限制性和非限制性的差別。 關(guān)于Ar 二聚體的情況是這樣的,我換了He二聚體可以得到滿意的結(jié)果,就是換成Ar原子后就出現(xiàn)這樣的問(wèn)題,基組已經(jīng)用到aug-cc-pvtz. 這樣一來(lái),所以,我就減小了步長(zhǎng)到0.005,方才可以解決了這個(gè)問(wèn)題。是一個(gè)極小的energy minimum point [ Last edited by gmy1990 on 2013-9-11 at 19:45 ] |
木蟲 (著名寫手)
|
樓上的輸入應(yīng)該沒(méi)問(wèn)題的;你再仔細(xì)查看下是否少抄了一行。不行的話你還可以再試試guess=mix。 在解離部分你的rmp2和ump2實(shí)際上都是基于rhf參考態(tài),參考態(tài)本身有定性問(wèn)題,大概微擾理論發(fā)散了吧。 如果你算Ar2我覺(jué)得限制和非限制就是應(yīng)該沒(méi)有差別的。你用hf和mp2的話本身就不能準(zhǔn)確描述色散相互作用,即使找到極小也沒(méi)有什么參考價(jià)值。盡量要用ccsd或者casscf吧 |
榮譽(yù)版主 (著名寫手)
木蟲 (著名寫手)
|
H2用各種方法計(jì)算出來(lái)的勢(shì)能面你可以在szabo&ostlund的書上看到,uhf和rhf是在227頁(yè),ump2和ump3是在375頁(yè)。 H2在解離時(shí)相關(guān)能應(yīng)該是0,uhf,ump2,casscf和ccsd都應(yīng)該給出精確解-1a.u.;你可以以此來(lái)查看你uhf的解是否正確,如果離-1很遠(yuǎn)說(shuō)明可能有問(wèn)題。 對(duì)于Ar2,即使找到一個(gè)極小點(diǎn),這個(gè)極小點(diǎn)并沒(méi)有意義;不論mp2和hf都不描述色散相互作用,只能描述排斥力。你得到的極小很可能只是基組重疊誤差造成的,不反映實(shí)際物理相互作用。要得到有物理意義的勢(shì)能面還是需要換方法 |
榮譽(yù)版主 (著名寫手)
木蟲 (著名寫手)
|
我把這兩幅圖掃下來(lái)了,只有紙板書,手機(jī)照相不太清楚 uhf,rhf和精確解 http://www.gaoyang168.com/bbs/wpic.php?id=12650 uhr,ump2,ump3和精確解 http://www.gaoyang168.com/bbs/wpic.php?id=12651 UHF大概從bohr2.3左右開始就應(yīng)該開始低于RHF,然后逐漸收斂到和ccsd一樣的能量。 |

榮譽(yù)版主 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +5 | tcxiaoxx 2026-03-20 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
[考研] 350求調(diào)劑 +4 | weudhdk 2026-03-19 | 4/200 |
|
|
[考研] 279分求調(diào)劑 一志愿211 +10 | chaojifeixia 2026-03-19 | 11/550 |
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|
[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考! +3 | lemonzzn 2026-03-17 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽(yáng)花天天開心 2026-03-16 | 4/200 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂(lè) 2026-03-14 | 5/250 |
|