24小时热门版块排行榜    

查看: 557  |  回复: 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

新虫 (初入文坛)

这么久了没人回答我的问题,看来用这个软件的人不多啊,但这个软件真的不错,自己赞一下。
我遇到的这个问题现在已经解决了,还是我自己问的开发者才知道怎么回事的。
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的回帖
查看全部 3 个回答

hwp199

新虫 (初入文坛)

没人回答啊,自己顶起。。。
2楼2014-08-23 12:47:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +24 医学老男孩 2026-08-13 50/2500 2026-08-15 06:38 by 学员8dgXkO
[基金申请] filecode +7 cratir 2026-08-14 11/550 2026-08-15 06:19 by 学员8dgXkO
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 04:45 by 4wMiSEwB6436
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 7lpolszZVXgi 2026-08-14 5/250 2026-08-15 04:41 by 4wMiSEwB6436
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 7lpolszZVXgi 2026-08-14 6/300 2026-08-15 04:21 by 4wMiSEwB6436
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 4/200 2026-08-15 04:17 by 4wMiSEwB6436
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 02:21 by 4wMiSEwB6436
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-15 00:16 by 4wMiSEwB6436
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-14 22:32 by 4wMiSEwB6436
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见