| 9 | 1/1 | 返回列表 |
| 查看: 2622 | 回復(fù): 8 | |||
luqing6879木蟲 (著名寫手)
|
[求助]
請問mrci 和mrci+q的區(qū)別
|
|
我用molpro計算勢能面,當(dāng)用mrci的能量值進(jìn)行擬合時,得到的誤差符合所需精度,但是用mrciq進(jìn)行擬合時,擬合誤差就很大。 請問mrci 和mrci+q 的區(qū)別是什么?mrci+q“錯誤”地考慮了什么因素造成了這么大的誤差? 謝謝。^_^ |
我的收藏 |

專家顧問 (職業(yè)作家)
地溝油冶煉專家
![]() |
專家經(jīng)驗: +458 |

木蟲 (著名寫手)
|
謝謝回復(fù)。 1 我看到過評論說“+Q部分大致對應(yīng)于CISD與CCSD兩種方法之間的差別”。 但是不是很理解為什么要進(jìn)行這種校正。一個是單參考方法,一個是多參考方法。如果體系有multi-reference性質(zhì),為什么要關(guān)心ccsd呢?一個令我滿意的解釋是來自維基百科的: It allows one to estimate the value of the full configuration interaction energy from a limited configuration interaction expansion result, although more precisely it estimates the energy of configuration interaction up to quadruple excitations (CISDTQ) from the energy of configuration interaction up to double excitations (CISD). 請指教。 另外搭車問一下,molpro是不是不支持邏輯運(yùn)算符?例如: if (k.eq10 .or. k.eq.20) then ... endif 總是報錯。但是把上面的if分開來寫成兩個if就沒事。 |

木蟲 (正式寫手)
搞電器工業(yè)化的化學(xué)打工仔

木蟲 (著名寫手)

木蟲 (正式寫手)
搞電器工業(yè)化的化學(xué)打工仔

木蟲 (著名寫手)
|
+Q就是davidson correction,是一個非常簡單快捷的微擾理論估計四階激發(fā)對能量貢獻(xiàn)的方法。不需要額外的計算,只用參考態(tài)(scf或mcscf)能量,ci能量,和ci波函數(shù)與參考態(tài)波函數(shù)的交疊三個數(shù)字,就可以直接外推出來。 對mrci,一般有三種算法 (1-c0^2)*(E[ci]-E[ref]) (1-c0^2)*(E[ci]-E[ref])/(c0^2) (1-c0^2)*(E[ci]-E[ref])/(2*c0^2-1) 這不是一個很干凈的矯正,但是經(jīng)常有效,而且又不花時間,所以比較常用。當(dāng)然,因為不干凈,不準(zhǔn)的情況也有發(fā)生。 如果矯正結(jié)果不連續(xù),有幾個可能性。首先,這個矯正是基于微擾假定的,所以如果你的mcscf和mrci的解已經(jīng)出現(xiàn)定性差別,這個矯正就完全不成立了。也就是如果c0比較小,結(jié)果就不可信。一般c0應(yīng)該在0.8以上 再有可能是c0不連續(xù)造成的。有可能是你的mcscf不連續(xù),出現(xiàn)了某個點軌道忽然翻轉(zhuǎn)的情況;同時因為要把mcscf的波函數(shù)和mrci的波函數(shù)一一對應(yīng),如果對應(yīng)找錯了就很容易出現(xiàn)問題。你可以在輸出文件中找到每個mrci解和所有參考態(tài)的交疊,看一下是不是選錯了。 |

木蟲 (著名寫手)

| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 305求調(diào)劑 +7 | RuiFairyrui 2026-03-28 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 學(xué)碩274求調(diào)劑 +9 | Li李魚 2026-03-26 | 9/450 |
|
|
[考研] 生物學(xué)學(xué)碩,一志愿湖南大學(xué),初試成績338 +6 | YYYYYNNNNN 2026-03-26 | 7/350 |
|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4 | 小張zxy張 2026-03-26 | 8/400 |
|
|
[考研] 275求調(diào)劑 +10 | Micky11223 2026-03-25 | 14/700 |
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
[考研] 085404求調(diào)劑,總分309,本科經(jīng)歷較為豐富 +4 | 來財aa 2026-03-25 | 4/200 |
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +4 | 手機(jī)用戶 2026-03-26 | 4/200 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 材料292調(diào)劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|
[考研] 0703化學(xué)338求調(diào)劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
[考研] 生物學(xué) 296 求調(diào)劑 +4 | 朵朵- 2026-03-26 | 6/300 |
|
|
[考研] 352求調(diào)劑 +4 | 大米飯! 2026-03-22 | 4/200 |
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|