24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1715  |  回复: 2

slay3

铜虫 (初入文坛)

[求助] Gamess 能量分解的输入文件

各位大虾,本人想用Gamess做一个二聚体的能量分解,由于没有接触过此软件,不知如何编写输入文件,还有能量分解的方法也不了解,请大家赐教,在此急求能量分解的输入文件和输出文件的例子,及能量分解方法的种类和选取办法
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

liuhuisfp

金虫 (小有名气)

slay3: 回帖置顶 2012-09-14 13:53:05
引用回帖:
2楼: Originally posted by liuhuisfp at 2012-07-15 19:29:00
不清楚你体系是什么……我在这里给一个用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 ...

上面的例子是比较复杂的,考虑的PCM效应的,现在GAMESS官方程序中还没有,我们老板刚开发完(相应文章已被JCP接收)。所以你计算的时候不需要$PCM和$PCMCAV参数。

最为重要的是$CONTRL中的        runtyp=eda,用于说明你是进行EDA(能量分解)计算

然后是    $LMOEDA, 这一项用于设定EDA相应的参数:
matom(1)=9 1
   说明将下面的体系分为两个block,前一个block包含9个原子(前9个);后一个block包含一个原子(即Cu),计算这两个block间的相互作用能,并进行分解。需要注意的是一定不能让block间的原子交错。不错结果一定是错的。

mcharg(1)=0 1
用于定义两个block所带的电荷。此例子中,第一个不带电,第二个CU为+1价;

mmult(1)=1 1
定义两个block的多重度。此例子中两个BLOCK均为闭壳层,所以设定为1。

supbas=.t.
进行BSSE消除计算,这个默认选项,可以不写。

其它的一此参数与$BASIS, $SCF, SSYSTEM都是根据实际需要调整的
3楼2012-07-15 19:41:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

liuhuisfp

金虫 (小有名气)

【答案】应助回帖

不清楚你体系是什么……我在这里给一个用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

先回,下面解释参数
2楼2012-07-15 19:29:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 slay3 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 3/150 2026-05-15 12:03 by onwj4wpxp2
[基金申请] 精华III评审感受-评审感受-评审感受 +14 ferrarichen 2026-05-11 18/900 2026-05-15 11:12 by cmhchen
[基金申请] 这年头没有找到涵评专家,还有中面上的可能吗 +9 dd921ww 2026-05-12 10/500 2026-05-15 10:41 by muyiliuhui
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 3/150 2026-05-15 09:23 by onwj4wpxp2
[基金申请] 青C资助名额大幅增加! +10 西葫芦炒鸡蛋 2026-05-13 14/700 2026-05-15 09:07 by gy116024
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 cjf4bx70cj 2026-05-14 4/200 2026-05-15 09:03 by gagyerk94e
[文学芳草园] 风把牡丹吹跑了 +4 myrtle 2026-05-12 7/350 2026-05-14 23:58 by myrtle
[教师之家] 教学课件你会给同学吗 +8 硕士研究生吗 2026-05-13 8/400 2026-05-14 22:23 by 常规沥青
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
[基金申请] 重磅!青年科学基金项目(C类)资助增幅预计超过50% +5 水和泥不是水泥 2026-05-13 7/350 2026-05-14 20:57 by 水和泥不是水泥
[高分子] 本人最近太闲了,谁有问题可以提,每天会统一回复 +8 一切都是空工 2026-05-12 19/950 2026-05-14 20:03 by 一切都是空工
[考博] 申博自荐 +4 食品的橙子 2026-05-09 6/300 2026-05-14 16:05 by great1919
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
[基金申请] 请问大佬b0816评完了吗 +3 市民华南虎 2026-05-12 7/350 2026-05-14 07:41 by 市民华南虎
[论文投稿] 有带发论文的吗 +3 山楂之术 2026-05-09 3/150 2026-05-13 17:56 by Cyhcl2629
[硕博家园] 导师各种操作恶心咋办 +11 苍白的小青天 2026-05-09 13/650 2026-05-13 17:11 by 六两废铜
[论文投稿] 求助大佬sci投稿哪个好中 +3 江沅188 2026-05-12 4/200 2026-05-13 14:35 by 江沅188
[考博] 西南大学考核制博士 +3 lijunjie84 2026-05-11 6/300 2026-05-12 18:09 by lijunjie84
[文学芳草园] 窗边初夏的小雨 +7 阿美_Lml888 2026-05-09 10/500 2026-05-12 15:27 by 阿美_Lml888
[考博] 现在不知道怎么办,感觉很痛苦 +4 qweww 2026-05-11 5/250 2026-05-11 20:23 by Oversize
信息提示
请填处理意见