应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gamess/Q-Chem » 求助gamess的hessian错误纠正
31/1返回列表
查看: 468  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-08-06 09:53:59
已阅   回复此楼   关注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 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭