| 3 | 1/1 | 返回列表 |
| 查看: 1687 | 回復(fù): 2 | |||
[求助]
Gamess 能量分解的輸入文件
|
| 各位大蝦,本人想用Gamess做一個二聚體的能量分解,由于沒有接觸過此軟件,不知如何編寫輸入文件,還有能量分解的方法也不了解,請大家賜教,在此急求能量分解的輸入文件和輸出文件的例子,及能量分解方法的種類和選取辦法 |
金蟲 (小有名氣)
|
不清楚你體系是什么……我在這里給一個用LMOEDA做能量分解的例子: $contrl scftyp=rhf runtyp=eda dfttyp=b3lyp1 $end $contrl icharg=1 mult=1 ispher=1 $end $lmoeda matom(1)=9 1 mcharg(1)=0 1 mmult(1)=1 1 supbas=.t. $end $basis gbasis=cct $end $pcm solvnt=water $end $pcmcav alpha=1.1 $end $pcmcav rin(1)=1.926 rin(2)=2.018 rin(3)=1.926 rin(4)=1.443 $end $pcmcav rin(5)=1.443 rin(6)=1.443 rin(7)=1.443 rin(8)=1.443 $end $pcmcav rin(9)=1.443 rin(10)=1.748 $end $tescav ntsall=240 mthall=4 $end $scf diis=.t. soscf=.f. dirscf=.t. fdiff=.f. npunch=1 ethrsh=100 $end $system mwords=100 memddi=600 $end $data itle: 1 C 6.0 10.1942801565 48.6863306161 32.3481462225 S 16.0 10.3578593780 50.3645238138 33.0032994531 C 6.0 8.7849632260 51.0813218166 32.4700178815 H 1.0 9.3667191827 48.2198062710 32.8775002151 H 1.0 11.1144700756 48.1540385817 32.5734553514 H 1.0 8.6596670994 50.9919698765 31.3954506977 H 1.0 7.9988865414 50.5450972368 32.9965472141 H 1.0 8.7786391502 52.1244074795 32.7745939847 H 1.0 10.0084600746 48.7010866606 31.2787325421 CU 29.0 11.9147019946 51.3366703824 31.8673297512 $end 先回,下面解釋參數(shù) |
金蟲 (小有名氣)
|
上面的例子是比較復(fù)雜的,考慮的PCM效應(yīng)的,現(xiàn)在GAMESS官方程序中還沒有,我們老板剛開發(fā)完(相應(yīng)文章已被JCP接收)。所以你計算的時候不需要$PCM和$PCMCAV參數(shù)。 最為重要的是$CONTRL中的 runtyp=eda,用于說明你是進(jìn)行EDA(能量分解)計算 然后是 $LMOEDA, 這一項(xiàng)用于設(shè)定EDA相應(yīng)的參數(shù): matom(1)=9 1 說明將下面的體系分為兩個block,前一個block包含9個原子(前9個);后一個block包含一個原子(即Cu),計算這兩個block間的相互作用能,并進(jìn)行分解。需要注意的是一定不能讓block間的原子交錯。不錯結(jié)果一定是錯的。 mcharg(1)=0 1 用于定義兩個block所帶的電荷。此例子中,第一個不帶電,第二個CU為+1價; mmult(1)=1 1 定義兩個block的多重度。此例子中兩個BLOCK均為閉殼層,所以設(shè)定為1。 supbas=.t. 進(jìn)行BSSE消除計算,這個默認(rèn)選項(xiàng),可以不寫。 其它的一此參數(shù)與$BASIS, $SCF, SSYSTEM都是根據(jù)實(shí)際需要調(diào)整的 |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +10 | 申子申申 2026-03-19 | 16/800 |
|
|---|---|---|---|---|
|
[考研] 319求調(diào)劑 +4 | 小力氣珂珂 2026-03-20 | 4/200 |
|
|
[考研] 305分求調(diào)劑(食品工程) +4 | Sxy112 2026-03-21 | 6/300 |
|
|
[考研] 291 求調(diào)劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 3/150 |
|
|
[考研] 303求調(diào)劑 +5 | 安憶靈 2026-03-22 | 6/300 |
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
[考研] 326求調(diào)劑 +5 | 諾貝爾化學(xué)獎覬?/a> 2026-03-15 | 8/400 |
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨; 2026-03-16 | 3/150 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|