| 3 | 1/1 | 返回列表 |
| 查看: 1686 | 回復(fù): 2 | ||
[求助]
Gamess 能量分解的輸入文件
|
||
| 各位大蝦,本人想用Gamess做一個(gè)二聚體的能量分解,由于沒(méi)有接觸過(guò)此軟件,不知如何編寫(xiě)輸入文件,還有能量分解的方法也不了解,請(qǐng)大家賜教,在此急求能量分解的輸入文件和輸出文件的例子,及能量分解方法的種類和選取辦法 |
金蟲(chóng) (小有名氣)
|
上面的例子是比較復(fù)雜的,考慮的PCM效應(yīng)的,現(xiàn)在GAMESS官方程序中還沒(méi)有,我們老板剛開(kāi)發(fā)完(相應(yīng)文章已被JCP接收)。所以你計(jì)算的時(shí)候不需要$PCM和$PCMCAV參數(shù)。 最為重要的是$CONTRL中的 runtyp=eda,用于說(shuō)明你是進(jìn)行EDA(能量分解)計(jì)算 然后是 $LMOEDA, 這一項(xiàng)用于設(shè)定EDA相應(yīng)的參數(shù): matom(1)=9 1 說(shuō)明將下面的體系分為兩個(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),可以不寫(xiě)。 其它的一此參數(shù)與$BASIS, $SCF, SSYSTEM都是根據(jù)實(shí)際需要調(diào)整的 |
金蟲(chóng) (小有名氣)
|
不清楚你體系是什么……我在這里給一個(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ù) |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 305分求調(diào)劑(食品工程) +4 | Sxy112 2026-03-21 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考! +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(qǐng)(333分) +5 | holy J 2026-03-21 | 5/250 |
|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
[考研] 一志愿深大,0703化學(xué),總分302,求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
[考研] 279分求調(diào)劑 一志愿211 +14 | chaojifeixia 2026-03-19 | 15/750 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 306求調(diào)劑 +4 | chuanzhu川燭 2026-03-18 | 4/200 |
|
|
[考研] 316求調(diào)劑 +6 | 梁茜雯 2026-03-19 | 6/300 |
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +5 | tcxiaoxx 2026-03-20 | 5/250 |
|
|
[考研] 08工學(xué)調(diào)劑 +5 | 用戶573181 2026-03-20 | 5/250 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 0854可跨調(diào)劑,一作一項(xiàng)核心論文五項(xiàng)專利,省、國(guó)級(jí)證書(shū)40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|