| 查看: 1697 | 回复: 2 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
Gamess 能量分解的输入文件
|
|||
| 各位大虾，本人想用Gamess做一个二聚体的能量分解，由于没有接触过此软件，不知如何编写输入文件，还有能量分解的方法也不了解，请大家赐教，在此急求能量分解的输入文件和输出文件的例子，及能量分解方法的种类和选取办法 |
回复此楼» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
Gaussian混合基组求助
已经有13人回复
- 如何在GAMESS输入文件中指定某个原子带电？
已经有7人回复
- 求助磷酸铵，磷酸二氢铵，磷酸一氢铵分解能量顺序
已经有17人回复
- gamess 能量分解，读取gaussian初始猜测
已经有6人回复
- 回答一些关于Multiwfn的疑问以及未来Multiwfn的发展打算
已经有37人回复
- Gamess-NEDA结果求助...
已经有7人回复
- 【求助】gamess中如何做限定电荷
已经有12人回复
- 【求助】winGAMESS运行错误，请教大家
已经有5人回复
- 【求助】ORCA, Turbomole, Gamess这些软件比GAUSSIAN有什么优势吗
已经有4人回复
- 【求助】gaussian怎么做PED分析？或有其他免费软件可以做吗？
已经有6人回复
slay3: 回帖置顶 2012-09-14 13:53:05
|
上面的例子是比较复杂的，考虑的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为+１价； mmult(1)=1 1 定义两个block的多重度。此例子中两个BLOCK均为闭壳层，所以设定为1。 supbas=.t. 进行BSSE消除计算，这个默认选项，可以不写。 其它的一此参数与$BASIS, $SCF, SSYSTEM都是根据实际需要调整的 |
3楼2012-07-15 19:41:17
【答案】应助回帖
|
不清楚你体系是什么……我在这里给一个用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
