| 5 | 1/1 | 返回列表 |
| 查看: 1692 | 回復(fù): 2 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
Gamess 能量分解的輸入文件
|
|||
| 各位大蝦,本人想用Gamess做一個(gè)二聚體的能量分解,由于沒有接觸過此軟件,不知如何編寫輸入文件,還有能量分解的方法也不了解,請(qǐng)大家賜教,在此急求能量分解的輸入文件和輸出文件的例子,及能量分解方法的種類和選取辦法 |
金蟲 (小有名氣)
|
上面的例子是比較復(fù)雜的,考慮的PCM效應(yīng)的,現(xiàn)在GAMESS官方程序中還沒有,我們老板剛開發(fā)完(相應(yīng)文章已被JCP接收)。所以你計(jì)算的時(shí)候不需要$PCM和$PCMCAV參數(shù)。 最為重要的是$CONTRL中的 runtyp=eda,用于說明你是進(jìn)行EDA(能量分解)計(jì)算 然后是 $LMOEDA, 這一項(xiàng)用于設(shè)定EDA相應(yīng)的參數(shù): matom(1)=9 1 說明將下面的體系分為兩個(gè)block,前一個(gè)block包含9個(gè)原子(前9個(gè));后一個(gè)block包含一個(gè)原子(即Cu),計(jì)算這兩個(gè)block間的相互作用能,并進(jìn)行分解。需要注意的是一定不能讓block間的原子交錯(cuò)。不錯(cuò)結(jié)果一定是錯(cuò)的。 mcharg(1)=0 1 用于定義兩個(gè)block所帶的電荷。此例子中,第一個(gè)不帶電,第二個(gè)CU為+1價(jià); mmult(1)=1 1 定義兩個(gè)block的多重度。此例子中兩個(gè)BLOCK均為閉殼層,所以設(shè)定為1。 supbas=.t. 進(jìn)行BSSE消除計(jì)算,這個(gè)默認(rèn)選項(xiàng),可以不寫。 其它的一此參數(shù)與$BASIS, $SCF, SSYSTEM都是根據(jù)實(shí)際需要調(diào)整的 |
金蟲 (小有名氣)
|
不清楚你體系是什么……我在這里給一個(gè)用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ā)表 | |
|---|---|---|---|---|
|
[考研] 考研調(diào)劑 +4 | 呼呼?~+123456 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 學(xué)碩274求調(diào)劑 +6 | Li李魚 2026-03-26 | 6/300 |
|
|
[考研] 311求調(diào)劑 +8 | lin0039 2026-03-26 | 8/400 |
|
|
[考研]
|
y7czhao 2026-03-26 | 9/450 |
|
|
[考研] 085404求調(diào)劑,總分309,本科經(jīng)歷較為豐富 +4 | 來財(cái)aa 2026-03-25 | 4/200 |
|
|
[考研] 291求調(diào)劑 +7 | 孅華 2026-03-22 | 7/350 |
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +3 | 小張zxy張 2026-03-26 | 5/250 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
[考研] 0703化學(xué)338求調(diào)劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 材料調(diào)劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
[考研] 324求調(diào)劑 +8 | hanamiko 2026-03-26 | 10/500 |
|
|
[考研]
|
nnnnnnn5 2026-03-25 | 5/250 |
|
|
[考研] 一志愿華理,數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 10/500 |
|
|
[考研] 352求調(diào)劑 +4 | 大米飯! 2026-03-22 | 4/200 |
|
|
[考研] 生物技術(shù)與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑,求導(dǎo)師收 +7 | 天天好運(yùn)來上岸?/a> 2026-03-24 | 7/350 |
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|