计算化学公社

标题: IRC的一个方向只进行一步就正常结束了 [打印本页]
作者
Author: liyuanhe211    时间: 2015-8-12 21:45
标题: IRC的一个方向只进行一步就正常结束了
本帖最后由 liyuanhe211 于 2015-8-12 21:46 编辑

用M06-2X-D3/6-311++G(d), opt=ts，得到一简单SN2反应的过渡态，正常结束，虚频振动模式符合需要，见 (, 下载次数 Times of downloads: 15) 。取该过渡态做IRC，reverse方向正常结束且符合预期，但forward方向只进行了两个点就Normal termination了。尝试过设定maxpoint、calcall、取消empiricaldispersion，均无效。这是什么原因造成的？如何解决？
(, 下载次数 Times of downloads: 52)
IRC输入文件 (, 下载次数 Times of downloads: 5) ：

  2. #p irc=(forward,maxpoints=30,calcfc) rm062x/6-311++g(d)  empiricaldispersion=gd3

  4. Title

  6. -1 1
  7. C                  0.00000000    0.96572700    0.00000000
  8. H                 -1.05028200    1.20221300    0.00000000
  9. H                  0.51218000    1.22744000    0.90923900
  10. H                  0.51218000    1.22744000   -0.90923900
  11. O                 -0.29809000    3.31940900    0.00000000
  12. H                  0.41758400    3.96800800    0.00000000
  13. Br                 0.05694500   -1.14213500    0.00000000
复制代码
IRC输出： (, 下载次数 Times of downloads: 5)

  1. PES minimum detected on this side of the pathway.
  2.     Magnitude of the gradient =      0.0000960
  3. Calculation of FORWARD path complete.
  4. Reaction path calculation complete.

  6. Energies reported relative to the TS energy of       -2689.870684
  7. --------------------------------------------------------------------------
  8.     Summary of reaction path following
  9. --------------------------------------------------------------------------
  10.                         Energy    RxCoord
  11.    1                    0.00000   0.00000
  12.    2                   -0.00001   0.18521
  13. --------------------------------------------------------------------------

  15. Total number of points:                    1
  16. Total number of gradient calculations:     3
  17. Total number of Hessian calculations:      1
  18. IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC
复制代码


作者
Author: yjcmwgk    时间: 2015-8-13 11:10
两个点之间的能量差大概是10^(-5) a.u.
高斯认为能量曲线已经“平”了
我觉得能不能增大stepsize？或者使用tight或vtight来缩小convergence criteria
作者
Author: sobereva    时间: 2015-8-15 02:14
(, 下载次数 Times of downloads: 45)



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3