24小时热门版块排行榜

返回列表

brothers

银虫 (小有名气)

应助: 4 (幼儿园)
金币: 220.1
散金: 100
帖子: 71
在线: 192.9小时
虫号: 900870
注册: 2009-11-12
性别: GG
专业: 金属有机化学

[求助] （已找到原因竟然是windows版gamess的问题......) $LMOEDA--DFT分解能量的？

小白求问，$LMOEDA理论上可以在DFT理论下分解能量，这是比MONOKUMA方法更为优越的，但是，请问Gamess输入文件里用什么命令来保证$LMOEDA在DFT下分解能量呢？
在Gamess 提供的exam28（LMOEDA运行成功，原程序如下）基础上，发现：
1、在exam28的$control中添加dfttyp=b3lyp（或b3lyp1，或SVWN等），算出来的本征值及能量会有一片一片的NaN。
2、在exam28的$control中添加mplevl=2，计算成功，结果中会出现 MP2 DISPERSION ENERGY 项。
感觉自己在理论及程序应用的哪个地方没搞明白，但周围实在找不到人咨询了，手册也看了不少，但始终未能解决。奉上90金币（几乎全部了），烦请各位大神不吝赐教。
注意，希望能从理论及程序两个方面解释清楚

Gamess原输入文件如下：
$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

回复此楼