应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Multiwfn » (已解决)Multiwfn的能量分解与Gamess的LMOEDA
11/1返回列表
查看: 1776  |  回复: 0
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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了
回复此楼

» 收录本帖的淘帖专辑推荐

能量分解

» 猜你喜欢
1楼 2015-01-20 22:31:42
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

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

科研从小木虫开始，人人为我，我为人人
相关版块跳转 我要订阅楼主 brothers 的主题更新
11/1返回列表
信息提示
关闭
请填处理意见
关闭