24小时热门版块排行榜    

查看: 1768  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

slay3

铜虫 (初入文坛)

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

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

» 猜你喜欢

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

已阅   回复此楼   关注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的回帖
查看全部 3 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +8 majunge000 2026-08-11 10/500 2026-08-15 18:55 by Equinoxhua
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 关于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
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +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
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见