24小时热门版块排行榜

返回列表

joannewanghf

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 497.6
散金: 30
帖子: 127
在线: 81.5小时
虫号: 941966
注册: 2010-01-13
性别: MM
专业: 理论和计算化学

[求助] 用firefly做morakuma能量分解分析求助

我有高斯优化的结构，想读取高斯的轨道信息，用gamess做能量分解分析。我做了一个gamess的输入文件，结果出错。请教大神们帮我看一眼。
文件的前三行是我自己写的，后面的是用chemcraft从高斯的*chk文件中转换过来的。
$contrl scftyp=rhf runtyp=morokuma coord=cart DFTTYP=BLYP icharg=1 mult=1 $end
$system memory=1000000 timlim=5 $end
$morokm iatm(1)=2 ICHM=-1 $end
$GUESS GUESS=MOREAD $END
$DATA
CoC12N4H28O2
C1
OXYGEN                   8.0 -0.000003    -0.719549    2.160780
S 6
  1    5484.671660000000    0.001831074430
  2    825.234946000000    0.013950172200
  3    188.046958000000    0.068445078100
  4       52.964500000000    0.232714336000
  5       16.897570400000    0.470192898000
  6       5.799635340000    0.358520853000
L 3
  1       15.539616200000    -0.110777549523    0.070874268235
  2       3.599933590000    -0.148026262363    0.339752839168
  3       1.013761750000    1.130767015132    0.727158577359
L 1
  1       0.270005823000    1.000000000000    1.000000000000
D 1
  1       0.800000000000    1.000000000000
。
。
。
HYDROGEN                1.0 -0.884356    2.868426    0.858971
S 3
  1       18.731137000000    0.033494604317
  2       2.825394360000    0.234726953116
  3       0.640121692000    0.813757326402
S 1
  1       0.161277759000    1.000000000000

$END

---  80 occupied orbitals are printed (391 AO coefficients in each orbital)  ---
$VEC
1  1 7.04172930E-01 2.34363045E-02 5.33577372E-08 1.22420124E-03-5.53465003E-04
1  2-2.09598275E-02-3.10633177E-07-1.41858416E-03 4.15200790E-03 3.71883517E-04
1  3 7.89072408E-05-4.50790758E-04-8.46154524E-09 2.07272382E-08 3.09049438E-04
1  4 7.03301507E-01 2.34084957E-02-5.32409437E-08-1.22356892E-03-5.52703761E-04
1  5-2.09355310E-02 3.11938160E-07 1.41341985E-03 4.14804713E-03 3.71299014E-04
1  6 7.94895473E-05-4.50788561E-04-7.90720273E-09-1.93805526E-08-3.08805001E-04
1  7 1.00960775E-04 7.26615414E-04-2.03105945E-03 9.43352469E-03-1.67975214E-02
1  8 7.89341915E-04-2.99897083E-11 2.47962193E-09-8.42177483E-06-2.03922667E-10
。
。
。
80 77 1.91662066E-03 1.09936419E-06-2.43275026E-03 1.06607291E-02-5.98053645E-02
80 78 6.58375343E-06 3.08212352E-05-1.06816305E-02 5.97265589E-02-3.96476149E-03
80 79 2.56671129E-02
$END

我按照这个运行了一下，结果出错：
*** SCFTYP=RHF MUST HAVE MULT=1 ***

ADDRESS 0x0059618B HAS INITIATED PROGRAM ABORT BECAUSE OF FATAL ERROR(S)

CPU       TIME: STEP =    0.03 ,  TOTAL =       0.0 SECONDS ( 0.0 MIN)
WALL CLOCK TIME: STEP =    0.04 ,  TOTAL =       0.0 SECONDS ( 0.0 MIN)
CPU UTILIZATION: STEP =    79.11%,  TOTAL =    118.66%

   397122 WORDS OF DYNAMIC MEMORY USED
   81280 BYTES OF HEAP MEMORY USED,    24384 BYTES REMAIN IN USE
EXECUTION OF FIREFLY TERMINATED ABNORMALLY AT 12:46:15  4-DEC-2015

疑问：我已经写了icharg=1，可以结果文件中说分子电荷为0；写了muti=1，结果输出却说muti不能等于1。希望熟悉该程序的人帮忙解决一下。多谢。

回复此楼

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

能量分解

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

没有哪个港口是永远的停留。。

1楼 2015-12-07 10:31:20

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 joannewanghf 的主题更新

返回列表