24小时热门版块排行榜    

查看: 853  |  回复: 0

joannewanghf

铁虫 (小有名气)

[求助] 用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。希望熟悉该程序的人帮忙解决一下。多谢。
回复此楼
没有哪个港口是永远的停留。。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 joannewanghf 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 12:13 by l0VvVHGBGRLv
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:03 by l0VvVHGBGRLv
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:02 by l0VvVHGBGRLv
[考博] 找导师 +5 yuanjiabao 2026-08-29 6/300 2026-08-30 11:52 by zhouyanli11
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 11:51 by l0VvVHGBGRLv
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:41 by l0VvVHGBGRLv
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 11:32 by l0VvVHGBGRLv
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:19 by l0VvVHGBGRLv
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 11:18 by l0VvVHGBGRLv
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:17 by ZPa0EcMwuECS
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 06:10 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 02:09 by ZPa0EcMwuECS
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见