24小时热门版块排行榜    

查看: 1784  |  回复: 2

slay3

铜虫 (初入文坛)

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

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

» 猜你喜欢

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

已阅   回复此楼   关注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的回帖

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的回帖
相关版块跳转 我要订阅楼主 slay3 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 jCd0dEvKHShX 2026-08-29 7/350 2026-08-30 19:45 by hZiFeudZyoGR
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 19:37 by hZiFeudZyoGR
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 19:35 by hZiFeudZyoGR
[博后之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 18:28 by hZiFeudZyoGR
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +7 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 18:07 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 17:14 by hZiFeudZyoGR
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +8 ASdOkHsho7FD 2026-08-28 12/600 2026-08-30 17:11 by hZiFeudZyoGR
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:48 by ZPa0EcMwuECS
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见