 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 3146 | 回復(fù): 18 | |||
[求助]
multi mrci計(jì)算問題，態(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í)沒有問題，但是當(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í)沒有出錯(cuò)？是輸入哪里有問題？rhf中需要加入state,2么？ 上述rhf，multi，mrci計(jì)算的是不是就是態(tài)平均的結(jié)果，是不是兩個(gè)對(duì)稱性分別為A1和A2的態(tài)各自計(jì)算后取其能量平均值？ 請(qǐng)高手指點(diǎn)！ |
回復(fù)此樓» 本帖已獲得的紅花（最新10朵）
18766315007» 猜你喜歡
深圳理工大學(xué)招聘AI4Science與分子模擬方向科研人員
已經(jīng)有0人回復(fù)
-
全獎(jiǎng)博士 英國(guó)利物浦大學(xué) × 臺(tái)灣清華大學(xué) 聯(lián)合培養(yǎng)
已經(jīng)有0人回復(fù)
-
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有243人回復(fù)
-
尋博導(dǎo)
已經(jīng)有13人回復(fù)
-
基金正文30頁(yè)指的是報(bào)告正文還是整個(gè)申請(qǐng)書
已經(jīng)有9人回復(fù)
-
哈爾濱工業(yè)大學(xué)（深圳）綠色化學(xué)團(tuán)隊(duì)招收2026年秋季入學(xué)博士生
已經(jīng)有11人回復(fù)
-
澳大利亞昆士蘭科技大學(xué)PhD博士獎(jiǎng)學(xué)金
已經(jīng)有0人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- molpro說明書的例子運(yùn)行不了？？
已經(jīng)有4人回復(fù)
- molpro試做旋軌耦合,怎么在輸出中找omega及其對(duì)應(yīng)的能量
已經(jīng)有5人回復(fù)
- 【求助】Molpro如何執(zhí)行旋軌耦合計(jì)算 如何執(zhí)行態(tài)平均計(jì)算?
已經(jīng)有10人回復(fù)
14樓2017-12-11 21:59:53
 |
2樓2013-07-05 21:42:06
|
|
3樓2013-07-11 22:33:37
|
|
4樓2013-07-12 11:08:34
fatpig8832
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 16 (小學(xué)生)
- 金幣: 7420.2
- 散金: 31
- 紅花: 28
- 帖子: 1438
- 在線: 610.9小時(shí)
- 蟲號(hào): 770403
- 注冊(cè): 2009-05-14
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
kent1022: 金幣+30 2013-07-13 22:30:34
gmy1990: 金幣+3 2013-07-14 07:21:27
kent1022: 金幣+30 2013-07-13 22:30:34
gmy1990: 金幣+3 2013-07-14 07:21:27
|
樓主是單單要算1Delta態(tài)吧？ 有兩個(gè)問題... 首先，為什么multi選這個(gè)space? 16,6,6,1比起rhf的12,6,6,1來，只是在A1上增加了4個(gè)軌道，而A2的軌道只是1...那你算symmetry 4 時(shí)不出錯(cuò)才怪...1.5埃估計(jì)是碰巧算過去的... 建議：把每個(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))過的了，不需要你手動(dòng)再平均一次。 rhf不能算態(tài)平均。 需要注意的是：active space必須要選好，要不然本來簡(jiǎn)并的兩個(gè)態(tài)有可能算到不同能量；本來不簡(jiǎn)并的則沒什么可能算到簡(jiǎn)并。 |
5樓2013-07-12 16:42:37
|
rhf和multi還都是按照一篇文章里的設(shè)置的，看來有些發(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)平均吧 |
6樓2013-07-13 22:19:43
7樓2013-07-13 22:25:43
fatpig8832
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 16 (小學(xué)生)
- 金幣: 7420.2
- 散金: 31
- 紅花: 28
- 帖子: 1438
- 在線: 610.9小時(shí)
- 蟲號(hào): 770403
- 注冊(cè): 2009-05-14
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
kent1022: 金幣+20 2013-07-16 10:36:13
kent1022: 金幣+20 2013-07-16 10:36:13
|
首先請(qǐng)把a(bǔ)ctive space搞清楚，為什么你要選16,6,6,2而不是15,6,6,2或14,7,7,2或其他的呢？也許16,6,6,2是合理的，只要你搞清楚每個(gè)軌道所代表的是哪個(gè)具體軌道，能夠說服自己就行... 然后，“態(tài)平均”和“分開算再平均”是兩回事...舉例說明：假設(shè)你用“態(tài)平均”算了兩個(gè)態(tài)，能量分別是-100.100000和-100.200000，這兩個(gè)能量就是每個(gè)態(tài)的“態(tài)平均”能量，不需要再加起來除以2... |
8樓2013-07-14 13:45:22
fatpig8832
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 16 (小學(xué)生)
- 金幣: 7420.2
- 散金: 31
- 紅花: 28
- 帖子: 1438
- 在線: 610.9小時(shí)
- 蟲號(hào): 770403
- 注冊(cè): 2009-05-14
9樓2013-07-14 13:50:27
18766315007
新蟲 (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 6318.5
- 散金: 100
- 紅花: 6
- 帖子: 1261
- 在線: 263.5小時(shí)
- 蟲號(hào): 2979660
- 注冊(cè): 2014-02-20
- 專業(yè): 原子和分子物理
10樓2017-12-11 16:45:20
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 kent1022 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過 +5 | 幸運(yùn)哩哩 2026-03-22 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +3 | niko- 2026-03-22 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +10 | 楓橋ZL 2026-03-18 | 12/600 |
|
|
|
[考研] 考研調(diào)劑 +3 | 呼呼？~+123456 2026-03-21 | 3/150 |
|
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|
|
|
[考研] 303求調(diào)劑 +5 | 睿08 2026-03-17 | 7/350 |
|
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|