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

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

木蟲 (著名寫手)
|
謝謝回復(fù)。 1 我看到過評(píng)論說“+Q部分大致對(duì)應(yīng)于CISD與CCSD兩種方法之間的差別”。 但是不是很理解為什么要進(jìn)行這種校正。一個(gè)是單參考方法,一個(gè)是多參考方法。如果體系有multi-reference性質(zhì),為什么要關(guān)心ccsd呢?一個(gè)令我滿意的解釋是來自維基百科的: 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). 請(qǐng)指教。 另外搭車問一下,molpro是不是不支持邏輯運(yùn)算符?例如: if (k.eq10 .or. k.eq.20) then ... endif 總是報(bào)錯(cuò)。但是把上面的if分開來寫成兩個(gè)if就沒事。 |

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

木蟲 (著名寫手)

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

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

木蟲 (著名寫手)

| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 299求調(diào)劑 +7 | 15188958825 2026-03-25 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
[考研] 一志愿鄭州大學(xué),080500學(xué)碩,總分317分求調(diào)劑 +5 | 舉個(gè)栗子oi 2026-03-24 | 6/300 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
[考研] 291求調(diào)劑 +7 | 孅華 2026-03-22 | 7/350 |
|
|
[考研] 085600,材料與化工321分,求調(diào)劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
[考研] 333求調(diào)劑 +6 | wfh030413@ 2026-03-23 | 6/300 |
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
[考研] 332求調(diào)劑 +6 | 032500 2026-03-25 | 6/300 |
|
|
[考研] 26考研-291分-廈門大學(xué)(085601)-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3 | min3 2026-03-24 | 4/200 |
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
[考研] 289材料與化工(085600)B區(qū)求調(diào)劑 +4 | 這么名字咋樣 2026-03-22 | 5/250 |
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
[考研] 292求調(diào)劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|