24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1756  |  回复: 0

brothers

银虫 (小有名气)

[求助] (已解决)Multiwfn的能量分解与Gamess的LMOEDA

一、用Multiwfn的EDA功能作了简单的能量分解,发现只能按照zigler的方式分出E(steric),不能分出其中的静电势能。因此没办法与轨道作用能相比较;此外,不能分出fragments之间的sigma和pai的能量,不能深入讨论DCD模型里的那四种对称形式的键的贡献。真希望有朝一日能加入更详细的EDA功能。
二、此外,请问,用Gamess的LMOEDA进行能量分解,如何确保是在DFT下来分解能量的?
自带的示例中有LMOEDA分解的例子,但是不清楚怎么样写才是用DFT的而不是HF的啊?
如果在原示例文件的$control中添加dfttyp=b3lyp,算出来的结果就是一片一片的NaN。感觉似乎是自己在理论的哪个地方没搞明白,但周围实在找不到人咨询了,烦请各位大神不吝赐教。
原输入文件如下:
$contrl scftyp=rhf runtyp=eda coord=zmt $end
$system timlim=1 $end
$basis  gbasis=n31 ngauss=4 $end
$guess  guess=huckel $end

x$morokm iatm(1)=3 $end

$lmoeda matom(1)=3,4 mcharg(1)=0,0 mmult(1)=1,1 $end
$contrl nosym=1 $end
$system mwords=8 $end

$data
water-ammonia dimer...4-31G basis set
Cs

H
O 1 rOH
H 2 rOH 1 aHOH
N 2  R  1 aHOH    3    0.0
H 4 rNH 3 aHNaxis 1  180.0
H 4 rNH 3 aHNaxis 5 +120.0
H 4 rNH 3 aHNaxis 5 -120.0

rOH=0.956
aHOH=105.2
rNH=1.0124
aHNaxis=112.1451  ! makes HNH=106.67
R=2.93
$end

原来是软件问题,换了64位的就OK了
回复此楼
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 brothers 的主题更新
信息提示
请填处理意见