24小时热门版块排行榜    

查看: 572  |  回复: 2

hwp199

新虫 (初入文坛)

[交流] 求助gamess的hessian错误纠正

使用Firefly运行:
INPUT CARD>!   File created by MacMolPlt 7.3                                               
INPUT CARD> $CONTRL SCFTYP=RHF RUNTYP=HESSIAN MAXIT=30 MULT=1 COORD=CART FRIEND=GAMESSUK   
INPUT CARD>    MOLPLT=.TRUE. $END                                                         
INPUT CARD> $SYSTEM TIMLIM=525600 MEMORY=30000000 $END                                    
INPUT CARD> $BASIS GBASIS=N31 NGAUSS=6 NDFUNC=1 $END                                       
INPUT CARD> $SCF DIRSCF=.TRUE. $END                                                        
INPUT CARD> $FORCE METHOD=ANALYTIC VIBANL=.TRUE. $END                                      
INPUT CARD> $DATA                                                                          
……
.out文件:
              E         H         G         CV        CP        S
           KJ/MOL    KJ/MOL    KJ/MOL   J/MOL-K   J/MOL-K   J/MOL-K
ELEC.      0.000     0.000     0.000     0.000     0.000     0.000
TRANS.     3.718     6.197   -45.189    12.472    20.786   172.351
ROT.       3.718     3.718   -33.905    12.472    12.472   126.190
VIB.     848.605   848.605   729.408    27.613    27.613   399.790
TOTAL    856.042   858.521   650.313    52.556    60.871   698.331

              E         H         G         CV        CP        S
         KCAL/MOL  KCAL/MOL  KCAL/MOL CAL/MOL-K CAL/MOL-K CAL/MOL-K
ELEC.      0.000     0.000     0.000     0.000     0.000     0.000
TRANS.     0.889     1.481   -10.800     2.981     4.968    41.193
ROT.       0.889     0.889    -8.104     2.981     2.981    30.160
VIB.     202.821   202.821   174.333     6.600     6.600    95.552
TOTAL    204.599   205.191   155.429    12.561    14.548   166.905

     *******************************************************
     * THIS IS NOT A STATIONARY POINT ON THE MOLECULAR PES *
     *     THE VIBRATIONAL ANALYSIS IS NOT VALID !!!       *
     *******************************************************
(1)这里提示:
* THIS IS NOT A STATIONARY POINT ON THE MOLECULAR PES *
     *     THE VIBRATIONAL ANALYSIS IS NOT VALID !!!       *
怎么解决?
(2)计算结果和gaussian同样方法计算出来的差别很大,这是guassian计算出来的(单位cal/mol或kcal/mol):
E        CV        S
72.602        34.835        96.678
有什么方法可以使得gamess计算的结果与gaussian一致(至少Cv、S应该一致)?
(3)gamess算出的H是正的,但gaussian算出来的是负的,它们两者的H是同一个含义(标准焓)吗?
本人刚学习gamess-us(Firefly),诸多问题不能明白,希望各位不胜赐教!致敬!
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : Ben1.out
  • 2014-08-06 09:48:42, 107.49 K

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

hwp199

新虫 (初入文坛)

没人回答啊,自己顶起。。。
2楼2014-08-23 12:47:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hwp199

新虫 (初入文坛)

这么久了没人回答我的问题,看来用这个软件的人不多啊,但这个软件真的不错,自己赞一下。
我遇到的这个问题现在已经解决了,还是我自己问的开发者才知道怎么回事的。
INPUT CARD> $CONTRL SCFTYP=RHF RUNTYP=HESSIAN MAXIT=30 MULT=1 COORD=CART FRIEND=GAMESSUK   这里的 COORD=CART 用COORD=UNIQUE 就好了,否则软件只是检查并没有真正运行
gamess算出的H是正的,对于气态物质,可以验证这里H=E+PV=E+R
3楼2014-09-04 01:46:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hwp199 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 14:03 by GkCVyODlx4mR
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 13:54 by GkCVyODlx4mR
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:14 by l0VvVHGBGRLv
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:03 by l0VvVHGBGRLv
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:41 by l0VvVHGBGRLv
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:30 by l0VvVHGBGRLv
[论文投稿] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 11:07 by l0VvVHGBGRLv
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:39 by ZPa0EcMwuECS
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +7 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 05:40 by ZPa0EcMwuECS
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:45 by ZPa0EcMwuECS
[基金申请] 中青基了要发朋友圈吗? +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 基金未中,这种答复是模板吗? +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
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
信息提示
请填处理意见