| 24小時(shí)熱門版塊排行榜 |
| 6 | 1/1 | 返回列表 |
| 查看: 2722 | 回復(fù): 5 | |||
tanghao86金蟲 (初入文壇)
|
[求助]
CASSCF計(jì)算如何計(jì)算CCSD(T)? 已有2人參與
|
|
對(duì)于一個(gè)金屬化合物,不知道該化合物的基態(tài)是閉殼層單重態(tài),還是開竅成單重態(tài),或三重態(tài)。對(duì)這個(gè)化合物做了CASSCF(10,6)計(jì)算,請(qǐng)教大家,我因?yàn)椴恢烙性摶衔镌趩沃貞B(tài)上會(huì)有多少個(gè)激發(fā)態(tài),那么在計(jì)算時(shí) NROOT 該等于幾呢? 另外,好像CASSCF算完一個(gè)組態(tài)會(huì)出現(xiàn)一個(gè)HF能量,但是我想要CCSD(T)的能量,我該怎么做呢?是先進(jìn)行CASSCF計(jì)算,再把計(jì)算出來(lái)的每個(gè)組態(tài)進(jìn)行CCSD(T)單點(diǎn)計(jì)算 還是 可以在計(jì)算CASSCF時(shí)就加上ccsd(t)關(guān)鍵詞? 請(qǐng)大家指教!萬(wàn)分感謝! |
經(jīng)驗(yàn) |
榮譽(yù)版主 (著名寫手)
木蟲 (正式寫手)
|
1.對(duì)于單參考方法,如DFT,CCSD等,開竅層單重態(tài)不能描述(需要兩個(gè)斯萊特行列式),所以就要使用多參考CASSCF等方法來(lái)計(jì)算。(注意活性空間選取) 2.這個(gè)可以多取幾個(gè)態(tài),因?yàn)槟阒皇且獙で笞畹湍芰康膽B(tài)。對(duì)于阿貝爾群體系而言, CAS(10,6)總共有:2s+1 : 1,3,5,7 states Number :4950,6930,2310,210 計(jì)算公式可以在molcas公司提供的材料里找到?梢钥闯鲇蟹浅6嗟膽B(tài)。 3.CCSD(t)是以HF為參考的,以CASSCF為參考那就是MRCC,這個(gè)目前都沒有什么成熟的軟件(有個(gè)程序叫MRCC,里面有MRCC模塊,目前只能計(jì)算非常小的體系)。所以建議你使用CASPT2 ,具有可以參考軟件手冊(cè)。 |
木蟲 (正式寫手)
金蟲 (初入文壇)
|
非常感謝你們的提示。我現(xiàn)在用高斯做CASSCF, 用高斯做ccsd(t)單點(diǎn),但是單重態(tài)正常結(jié)束,而在進(jìn)行三重態(tài)ccsd(t)單點(diǎn)計(jì)算時(shí),出現(xiàn)了l913錯(cuò)誤,錯(cuò)誤信息如下,請(qǐng)問(wèn)有什么辦法解決嗎? 我本想用molpro做ccsd(T) 但是我不知道該如何寫基組的輸入文件,還想請(qǐng)問(wèn) 如果我想對(duì)金屬Ru 定義SDD基組和SDD贗勢(shì); 對(duì)P S原子定義SDD基組; 對(duì)C H定義6-31g基組,我該怎么寫molpro輸入文件中基組部分呢? ********************** DD1Dir will call FoFMem 1 times, MxPair= 4294 NAB= 1443 NAA= 741 NBB= 666. Norm of the A-vectors is 7.8132901D-05 conv= 1.00D-05. RLE energy= -1.4618693533 DE(Corr)= -1.4618698 E(CORR)= -2543.9889259 Delta=-1.57D-07 NORM(A)= 0.15023207D+01 ************* *MAX. CYCLES* ************* Dominant configurations: *********************** Spin Case I J A B Value BB 70 78 -0.112989D+00 BB 71 78 0.257942D+00 BB 72 78 -0.160335D+00 BB 73 78 -0.162565D+00 BB 76 79 0.288029D+00 BB 77 78 0.108302D+00 Largest amplitude= 2.88D-01 Error termination via Lnk1e in /usr/local/g09_D01/g09/l913.exe at Fri Oct 31 07:56:46 2014. Job cpu time: 5 days 5 hours 3 minutes 33.8 seconds. File lengths (MBytes): RWF= 18105 Int= 0 D2E= 0 Chk= 10 Scr= 1 |
金蟲 (初入文壇)
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑 +3 | 兔兔兔111223314 2026-03-29 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 總分293求調(diào)劑 +8 | 加一一九 2026-03-25 | 11/550 |
|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +22 | RichLi_ 2026-03-25 | 22/1100 |
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +14 | Mr. Z 2026-03-25 | 14/700 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5 | hanamiko 2026-03-29 | 5/250 |
|
|
[考研] 330分求調(diào)劑 +5 | qzenlc 2026-03-29 | 5/250 |
|
|
[考研] 生物學(xué)學(xué)碩,一志愿湖南大學(xué),初試成績(jī)338 +6 | YYYYYNNNNN 2026-03-26 | 7/350 |
|
|
[考研] 347求調(diào)劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
[考研] 311求調(diào)劑 +9 | lin0039 2026-03-26 | 9/450 |
|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
[考研] 291求調(diào)劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
[有機(jī)交流]
高溫高壓反應(yīng)求助
10+4
|
chibby 2026-03-25 | 4/200 |
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
[考研] 一志愿上海理工能源動(dòng)力(085800)310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
[考研] 化工專碩求調(diào)劑 +3 | question挽風(fēng) 2026-03-24 | 3/150 |
|
|
[考研] 材料/農(nóng)業(yè)專業(yè),07/08開頭均可,過(guò)線就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|