24小时热门版块排行榜    

查看: 1769  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +30 医学老男孩 2026-08-13 65/3250 2026-08-15 22:31 by zxy7909
[基金申请] 时间戳又变了8-15 +12 archvillain 2026-08-15 23/1150 2026-08-15 20:47 by foyo404
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见