IRC的一个方向只进行一步就正常结束了

liyuanhe211 · 发表于 Post on 2015-8-12 21:45:51

本帖最后由 liyuanhe211 于 2015-8-12 21:46 编辑

用M06-2X-D3/6-311++G(d), opt=ts，得到一简单SN2反应的过渡态，正常结束，虚频振动模式符合需要，见

HCH3BR_TS.LOG (142.45 KB, 下载次数 Times of downloads: 15) 。取该过渡态做IRC，reverse方向正常结束且符合预期，但forward方向只进行了两个点就Normal termination了。尝试过设定maxpoint、calcall、取消empiricaldispersion，均无效。这是什么原因造成的？如何解决？

IRC输入文件

HCH3BRIRC_GD3.gjf (699 Bytes, 下载次数 Times of downloads: 5) ：

#p irc=(forward,maxpoints=30,calcfc) rm062x/6-311++g(d) empiricaldispersion=gd3
Title
-1 1
C 0.00000000 0.96572700 0.00000000
H -1.05028200 1.20221300 0.00000000
H 0.51218000 1.22744000 0.90923900
H 0.51218000 1.22744000 -0.90923900
O -0.29809000 3.31940900 0.00000000
H 0.41758400 3.96800800 0.00000000
Br 0.05694500 -1.14213500 0.00000000

复制代码

IRC输出：

HCH3BRIRC_Forward.out (84.65 KB, 下载次数 Times of downloads: 5)

PES minimum detected on this side of the pathway.
Magnitude of the gradient = 0.0000960
Calculation of FORWARD path complete.
Reaction path calculation complete.
Energies reported relative to the TS energy of -2689.870684
--------------------------------------------------------------------------
Summary of reaction path following
--------------------------------------------------------------------------
Energy RxCoord
1 0.00000 0.00000
2 -0.00001 0.18521
--------------------------------------------------------------------------
Total number of points: 1
Total number of gradient calculations: 3
Total number of Hessian calculations: 1
IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC

复制代码

yjcmwgk · 发表于 Post on 2015-8-13 11:10:34

两个点之间的能量差大概是10^(-5) a.u.
高斯认为能量曲线已经“平”了
我觉得能不能增大stepsize？或者使用tight或vtight来缩小convergence criteria

sobereva · 发表于 Post on 2015-8-15 02:14:07

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] IRC的一个方向只进行一步就正常结束了

评分 Rate