24小时热门版块排行榜    

查看: 1688  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 311求调剂 +6 冬十三 2026-03-18 6/300 2026-03-22 20:18 by edmund7
[考研] 315分,诚求调剂,材料与化工085600 +3 13756423260 2026-03-22 3/150 2026-03-22 20:11 by edmund7
[考研] 293求调剂 +12 zjl的号 2026-03-16 17/850 2026-03-22 16:51 by i_cooler
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-20 3/150 2026-03-22 16:00 by ColorlessPI
[考研] 298求调剂一志愿211 +3 上岸6666@ 2026-03-20 3/150 2026-03-22 15:50 by ColorlessPI
[考研] 305分求调剂(食品工程) +4 Sxy112 2026-03-21 6/300 2026-03-22 15:26 by 无懈可击111
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 4/200 2026-03-22 10:25 by QGZDSYS
[考研] 生物学一志愿985,分数349求调剂 +4 zxts12 2026-03-21 7/350 2026-03-22 09:57 by zxts12
[基金申请] 山东省面上项目限额评审 +4 石瑞0426 2026-03-19 4/200 2026-03-22 08:50 by Wei_ren
[考研] 广西大学材料导师推荐 +3 夏夏夏小正 2026-03-17 5/250 2026-03-21 22:20 by 金昊ML
[考研] 工科0856求调剂 +3 沐析汀汀 2026-03-21 3/150 2026-03-21 18:30 by 学员8dgXkO
[考研] 316求调剂 +6 梁茜雯 2026-03-19 6/300 2026-03-21 06:32 by Ecowxq666!
[考研] 070300化学319求调剂 +7 锦鲤0909 2026-03-17 7/350 2026-03-21 03:46 by JourneyLucky
[考研] 299求调剂 +6 △小透明* 2026-03-17 6/300 2026-03-21 02:42 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 材料专硕英一数二306 +7 z1z2z3879 2026-03-18 7/350 2026-03-20 23:48 by JourneyLucky
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
信息提示
请填处理意见