| 5 | 1/1 | 返回列表 |
| 查看: 2721 | 回復(fù): 5 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
CASSCF計(jì)算如何計(jì)算CCSD(T)? 已有2人參與
|
|||
|
對于一個(gè)金屬化合物,不知道該化合物的基態(tài)是閉殼層單重態(tài),還是開竅成單重態(tài),或三重態(tài)。對這個(gè)化合物做了CASSCF(10,6)計(jì)算,請教大家,我因?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ì)算出來的每個(gè)組態(tài)進(jìn)行CCSD(T)單點(diǎn)計(jì)算 還是 可以在計(jì)算CASSCF時(shí)就加上ccsd(t)關(guān)鍵詞? 請大家指教!萬分感謝! |
經(jīng)驗(yàn) |
木蟲 (正式寫手)
|
1.對于單參考方法,如DFT,CCSD等,開竅層單重態(tài)不能描述(需要兩個(gè)斯萊特行列式),所以就要使用多參考CASSCF等方法來計(jì)算。(注意活性空間選取) 2.這個(gè)可以多取幾個(gè)態(tài),因?yàn)槟阒皇且獙で笞畹湍芰康膽B(tài)。對于阿貝爾群體系而言, 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 ,具有可以參考軟件手冊。 |
榮譽(yù)版主 (著名寫手)
木蟲 (正式寫手)
|
非常感謝你們的提示。我現(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ò)誤信息如下,請問有什么辦法解決嗎? 我本想用molpro做ccsd(T) 但是我不知道該如何寫基組的輸入文件,還想請問 如果我想對金屬Ru 定義SDD基組和SDD贗勢; 對P S原子定義SDD基組; 對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 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +21 | RichLi_ 2026-03-25 | 21/1050 |
|
|---|---|---|---|---|
|
[考研] 調(diào)劑310 +11 | 溫柔的晚安 2026-03-25 | 12/600 |
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
[考研] 329求調(diào)劑 +7 | 星野? 2026-03-26 | 7/350 |
|
|
[考研] 321求調(diào)劑 +7 | 璞玉~~ 2026-03-25 | 8/400 |
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動(dòng)285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
[考研] 266求調(diào)劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
[考研] 07化學(xué)280分求調(diào)劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
[考研] 復(fù)試調(diào)劑,一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
[考研] 085600,材料與化工321分,求調(diào)劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
[考研] 一志愿華東理工大學(xué)081700,初試分?jǐn)?shù)271 +6 | kotoko_ik 2026-03-23 | 7/350 |
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
[考研] 各位老師您好:本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|