| 5 | 1/1 | 返回列表 |
| 查看: 3157 | 回復(fù): 18 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
[求助]
multi mrci計(jì)算問(wèn)題,態(tài)平均?
|
||
|
計(jì)算一雙原子分子體系的1Delta態(tài)或3Pi態(tài),在C2v中分別是A1+A2和B1+B2. 若計(jì)算Delta態(tài)是否需要分別計(jì)算各個(gè)不同的對(duì)稱性然后取平均? 先rhf,再multi,然后icMRCI,輸入按以下 {rHF occ,12,6,6,1 wf,50, ,0 } {multi occ,16,6,6,1 wf,50,1,0 !A1 lquant,2 wf,50,4,0 !A2 lquant,2 } {ci occ,16,6,6,1 wf,50,1,0 !A1 wf,50,4,0 !A2 ! 1Delta } 計(jì)算分子核間距在r=1.5埃時(shí)沒(méi)有問(wèn)題,但是當(dāng)r增加到1.55時(shí),出現(xiàn)下列錯(cuò)誤 SYYMETRY= 4 MS2= 0 ? Error ? Open-shell singlet not implemented ? The problem occurs in scfpro ERROR EXIT CURRENT STACK: MAIN 為什么計(jì)算r=1.5埃時(shí)沒(méi)有出錯(cuò)?是輸入哪里有問(wèn)題?rhf中需要加入state,2么? 上述rhf,multi,mrci計(jì)算的是不是就是態(tài)平均的結(jié)果,是不是兩個(gè)對(duì)稱性分別為A1和A2的態(tài)各自計(jì)算后取其能量平均值? 請(qǐng)高手指點(diǎn)! |
鐵桿木蟲 (著名寫手)
鐵桿木蟲 (著名寫手)
|
樓主是單單要算1Delta態(tài)吧? 有兩個(gè)問(wèn)題... 首先,為什么multi選這個(gè)space? 16,6,6,1比起rhf的12,6,6,1來(lái),只是在A1上增加了4個(gè)軌道,而A2的軌道只是1...那你算symmetry 4 時(shí)不出錯(cuò)才怪...1.5埃估計(jì)是碰巧算過(guò)去的... 建議:把每個(gè)軌道具體表示的是什么軌道搞清楚,在從中挑出需要的軌道組成你multi用的space。例如,如果你的計(jì)算需要用到3d軌道,那么必須把A1上的3d0,3d2+,B1上的3d1+,B2上的3d1-,以及A2上的3d2-都包含在你的active space才行,這樣你的A2至少也要是2. 其次,你把兩個(gè)態(tài)一起寫,這就是在算態(tài)平均...得到的兩個(gè)能量是已經(jīng)態(tài)平均(默認(rèn)50:50,比例可以用參數(shù)調(diào))過(guò)的了,不需要你手動(dòng)再平均一次。 rhf不能算態(tài)平均。 需要注意的是:active space必須要選好,要不然本來(lái)簡(jiǎn)并的兩個(gè)態(tài)有可能算到不同能量;本來(lái)不簡(jiǎn)并的則沒(méi)什么可能算到簡(jiǎn)并。 |
|
rhf和multi還都是按照一篇文章里的設(shè)置的,看來(lái)有些發(fā)表了的東西也不可靠!本人目前還在摸索階段,所以先拿本人文章的例子練習(xí)。 我是要算1Delta態(tài),如果把A2加到2,輸入改為 {multi occ,16,6,6,2 wf,50,1,0 !A1 lquant,2 wf,50,4,0 !A2 lquant,2 } 這樣計(jì)算的A1和A2的態(tài)平均結(jié)果是不是就是我要計(jì)算的1Delta態(tài)? 這個(gè)結(jié)果與我分別計(jì)算 {multi;occ,16,6,6,2;wf,50,1,0;lquant,2} {multi;occ,16,6,6,2;wf,50,4,0;lquant,2} 得到的單點(diǎn)能能量的平均值應(yīng)該是一樣的吧?這就是所謂的態(tài)平均吧 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 抱歉 +4 | 田洪有 2026-03-30 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 材料工程專碩求調(diào)劑 +8 | hyl3153942 2026-03-29 | 8/400 |
|
|
[考研] 311求調(diào)劑 +7 | 藍(lán)月亮亮 2026-03-30 | 7/350 |
|
|
[考研] 289求調(diào)劑 +16 | 新時(shí)代材料 2026-03-27 | 16/800 |
|
|
[考研] 26考研-291分-廈門大學(xué)(085601)-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +5 | min3 2026-03-24 | 6/300 |
|
|
[考研] 332求調(diào)劑 +6 | @MZB382400 2026-03-28 | 6/300 |
|
|
[考研] 298求調(diào)劑 +3 | 什么是胖頭魚 2026-03-30 | 5/250 |
|
|
[考研] 070300一志愿211,312分求調(diào)劑院校 +4 | 小黃鴨寶 2026-03-30 | 4/200 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +6 | @taotao 2026-03-26 | 7/350 |
|
|
[考研] 085404求調(diào)劑,總分309,本科經(jīng)歷較為豐富 +6 | 來(lái)財(cái)aa 2026-03-25 | 6/300 |
|
|
[考研] 考研調(diào)劑 +7 | 小蠟新筆 2026-03-29 | 7/350 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5 | hanamiko 2026-03-29 | 5/250 |
|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +12 | c297914 2026-03-28 | 12/600 |
|
|
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
[考研] 考研調(diào)劑 +10 | 呼呼?~+123456 2026-03-24 | 10/500 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[考研] 一志愿鄭大085600,310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|