| 查看: 3145 | 回復: 18 | ||
[求助]
multi mrci計算問題,態(tài)平均?
|
||
|
計算一雙原子分子體系的1Delta態(tài)或3Pi態(tài),在C2v中分別是A1+A2和B1+B2. 若計算Delta態(tài)是否需要分別計算各個不同的對稱性然后取平均? 先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 } 計算分子核間距在r=1.5埃時沒有問題,但是當r增加到1.55時,出現(xiàn)下列錯誤 SYYMETRY= 4 MS2= 0 ? Error ? Open-shell singlet not implemented ? The problem occurs in scfpro ERROR EXIT CURRENT STACK: MAIN 為什么計算r=1.5埃時沒有出錯?是輸入哪里有問題?rhf中需要加入state,2么? 上述rhf,multi,mrci計算的是不是就是態(tài)平均的結果,是不是兩個對稱性分別為A1和A2的態(tài)各自計算后取其能量平均值? 請高手指點! |
|
|
|
鐵桿木蟲 (著名寫手)
|
樓主是單單要算1Delta態(tài)吧? 有兩個問題... 首先,為什么multi選這個space? 16,6,6,1比起rhf的12,6,6,1來,只是在A1上增加了4個軌道,而A2的軌道只是1...那你算symmetry 4 時不出錯才怪...1.5埃估計是碰巧算過去的... 建議:把每個軌道具體表示的是什么軌道搞清楚,在從中挑出需要的軌道組成你multi用的space。例如,如果你的計算需要用到3d軌道,那么必須把A1上的3d0,3d2+,B1上的3d1+,B2上的3d1-,以及A2上的3d2-都包含在你的active space才行,這樣你的A2至少也要是2. 其次,你把兩個態(tài)一起寫,這就是在算態(tài)平均...得到的兩個能量是已經態(tài)平均(默認50:50,比例可以用參數調)過的了,不需要你手動再平均一次。 rhf不能算態(tài)平均。 需要注意的是:active space必須要選好,要不然本來簡并的兩個態(tài)有可能算到不同能量;本來不簡并的則沒什么可能算到簡并。 |
|
rhf和multi還都是按照一篇文章里的設置的,看來有些發(fā)表了的東西也不可靠!本人目前還在摸索階段,所以先拿本人文章的例子練習。 我是要算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 } 這樣計算的A1和A2的態(tài)平均結果是不是就是我要計算的1Delta態(tài)? 這個結果與我分別計算 {multi;occ,16,6,6,2;wf,50,1,0;lquant,2} {multi;occ,16,6,6,2;wf,50,4,0;lquant,2} 得到的單點能能量的平均值應該是一樣的吧?這就是所謂的態(tài)平均吧 |
鐵桿木蟲 (著名寫手)
|
首先請把active space搞清楚,為什么你要選16,6,6,2而不是15,6,6,2或14,7,7,2或其他的呢?也許16,6,6,2是合理的,只要你搞清楚每個軌道所代表的是哪個具體軌道,能夠說服自己就行... 然后,“態(tài)平均”和“分開算再平均”是兩回事...舉例說明:假設你用“態(tài)平均”算了兩個態(tài),能量分別是-100.100000和-100.200000,這兩個能量就是每個態(tài)的“態(tài)平均”能量,不需要再加起來除以2... |
鐵桿木蟲 (著名寫手)
新蟲 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調劑院校信息 +4 | CX 330 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 材料學碩301分求調劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研] 333求調劑 +5 | 87639 2026-03-21 | 7/350 |
|
|
[考研]
|
Grand777 2026-03-21 | 3/150 |
|
|
[考研] 326求調劑 +4 | mlpqaz03 2026-03-15 | 4/200 |
|
|
[考研] 材料學碩333求調劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
[考研] 南昌大學材料專碩311分求調劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 301求調劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
[考研] 初始318分求調劑(有工作經驗) +3 | 1911236844 2026-03-17 | 3/150 |
|
|
[考研] 296求調劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 南京大學化學376求調劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
[考研] 304求調劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
[考研] 317求調劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 287求調劑 +7 | 晨昏線與星海 2026-03-19 | 8/400 |
|
|
[考研] A區(qū)線材料學調劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
[考研] 0817 化學工程 299分求調劑 有科研經歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 一志愿 南京航空航天大學大學 ,080500材料科學與工程學碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
[考研] 一志愿中國海洋大學,生物學,301分,求調劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 070303 總分349求調劑 +3 | LJY9966 2026-03-15 | 5/250 |
|