计算化学公社

标题: 关于IRC能量不在最高点且两端均为反应物的问题 [打印本页]

作者
Author: Jiangxue 时间: 2024-11-15 15:33
标题: 关于IRC能量不在最高点且两端均为反应物的问题
最近在找过渡态后跑IRC时总是遇到过渡态虚频，震荡趋势都正确，但是跑irc时能量都不在最高点，且两端都是反应物的情况。如下我算的具体内容及我采取的措施。
步骤如下：
1、过渡态搜索，得到了一个虚频为-885，且震动模式正确的过渡态。
2、进行irc计算时发现过渡态不在最高点，且两端为反应物。
对此进行了如下尝试：
1、加大了步长，由原来的10，变为了20，加了关键词guess=read读取过渡态chk文件中读取最后的波函数。得到的irc两端仍然是反应物。
2、加了关键词tight，int=ultrafine重新计算过渡态，按1重新跑irc，结果仍然是两边是反应物（irc计算时也加了tight，int=ultrafine）
3、优化过渡态时用calcall，结果仍然是两边是反应物。
过渡态和irc图如下：
(, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0)

过渡态和irc输入文件如下：

%chk=3-23-2s-ts-7-2-calcall.chk
%nprocshared=48
Will use up to 48 processors via shared memory.
%mem=30GB
----------------------------------------------------------------------
#p opt=(calcall,ts,gdiis,noeigen) freq ub3lyp scrf=(smd,solvent=chloro
benzene) nosymm def2svp em=gd3bj

******************************************
%chk=3-23-2s-ts-7-2-calcall.chk
%nprocshared=46
Will use up to 46 processors via shared memory.
%mem=30GB
----------------------------------------------------------------------
#p irc=(calcall,maxpoints=30,stepsize=10,lqa) ub3lyp scrf=(smd,solvent
=chlorobenzene) nosymm def2svp em=gd3bj guess=read
----------------------------------------------------------------------

后来又在计算另外一个H提取过渡态时遇到了同样的问题，虚频震荡趋势都很符合，跑irc同样能量不在最高点，且两端均为反应物。计算体系如下：
(, 下载次数 Times of downloads: 1)
这次我尝试了加大步长，guess=read读取过渡态的最后一帧，用forward,backword两边跑时，算backword时，交进去任务就报错。输入文件和报错内容如下：
%chk=1-HAT-step2-2s-ts-5-1.chk
%nprocshared=24
%mem=16GB
# irc=(backward,calcfc,maxpoints=30,stepsize=20,lqa) ub3lyp
scrf=(smd,solvent=chlorobenzene) nosymm def2svp em=gd3bj guess=read

报错如下：
%chk=1-HAT-step2-2s-ts-5-1.chk
%nprocshared=24
Will use up to 24 processors via shared memory.
%mem=16GB
----------------------------------------------------------------------
# irc=(backward,calcfc,maxpoints=30,stepsize=20,lqa) ub3lyp scrf=(smd,
solvent=chlorobenzene) nosymm def2svp em=gd3bj guess=read
----------------------------------------------------------------------
QPErr --- A syntax error was detected in the input line.
# irc=(backward,calcfc,maxpoints=30,step
      '
Last state= "IRC2"
TCursr= 158349 LCursr=       7
Error termination via Lnk1e in /home/test/soft/g09/g09/l1.exe at Thu Nov 14 21:02:31 2024.
Job cpu time:    0 days  0 hours  0 minutes  0.9 seconds.
File lengths (MBytes):  RWF=    5 Int=    0 D2E=    0 Chk=    85 Scr=    1

取forward计算出的最后一步进行opt优化，结果像产物但又没连上去。
直接过渡态最后一帧加tight优化，H连回N上。
请教老师，这种情况要怎么才能解决呢？我已经困在这问题里很久了，恳请老师指点。

作者
Author: sai77 时间: 2024-11-15 19:42
报错是因为语法错误，语法错误是因为往逆方向的关键词是reverse而不是backward/backword，forward可以不要想当然认为另一边就是backward

作者
Author: sobereva 时间: 7 day ago
去掉溶剂模型再试。本来氯苯那么低极性溶剂对于opt freq IRC的影响就微乎其微
stepsize改成30再试
IRC的tight和opt的tight完全是两回事，找过渡态时候用tight并不需要IRC时候用tight

内存别给得那么抠，尤其是涉及到算大体系Hessian的情况，仔细看下文
Gaussian的安装方法及运行时的相关问题
http://sobereva.com/439（http://bbs.keinsci.com/thread-10814-1-1.html）

作者
Author: Jiangxue 时间: 7 day ago

sai77 发表于 2024-11-15 19:42
报错是因为语法错误，语法错误是因为往逆方向的关键词是reverse而不是backward/backword，forward可以不要 ...

好的，谢谢

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