计算化学公社

标题: gamess做lmoeda能量分析的两个问题 [打印本页]

作者
Author: daogege 时间: 2017-5-3 09:20
标题: gamess做lmoeda能量分析的两个问题
  本来之前想用molpro分析作用能的，但是由于体系太大只能转用gamess。
  第一个问题就是pol项过大，体系的总能量也就是3340a.u.左右的总能量
(, 下载次数 Times of downloads: 46)
  第二个问题是看文献中用molpro分析结果中有ee（静电）,es（交换）,ind（感应或诱导） disp（色散）这四项，应该如何同gamess中ES（静电），EX（交换），REP（排斥），POL（极化）和DISP（色散）对应？

作者
Author: sobereva 时间: 2017-5-3 10:24
仔细看输出文件，看是否有诸如收敛失败等等错误提示。

molpro算SAPT0远远慢于PSI4，可以试试PSI4

SAPT的物理意义更清楚，计算方式更严格。LMOEDA的色散怎么去解释也看具体用什么理论方法，有可能完全无物理意义，其它项的解释一看原文便知。

作者
Author: daogege 时间: 2017-5-3 14:13

sobereva 发表于 2017-5-3 10:24
仔细看输出文件，看是否有诸如收敛失败等等错误提示。

molpro算SAPT0远远慢于PSI4，可以试试PSI4

谢谢老师！我先学习一下PSI4吧。
用gamess重复了一下老师那篇氢气氮气二聚体的文章，数据大多能对上，感觉Gamess的排斥+交换约等于molpro的交换，Gamess的极化约等于molpro的诱导。

作者
Author: sobereva 时间: 2017-5-4 11:17

daogege 发表于 2017-5-3 14:13
谢谢老师！我先学习一下PSI4吧。
用gamess重复了一下老师那篇氢气氮气二聚体的文章，数据大多能对上，感 ...

这么去对基本没问题

作者
Author: wei 时间: 2017-5-4 15:50
没有收敛

作者
Author: shengyuhu 时间: 2017-5-11 16:31

wei 发表于 2017-5-4 15:50
没有收敛

您好，我也在用gamess做能量分解的时候也出发现最后的能量很大，看收敛是看DENSITY CONVERGED这一项吗？

欢迎光临计算化学公社 (http://bbs.keinsci.com/)