计算化学公社

标题: 通过relaxed scan再TS后，IRC失败了，求分析 [打印本页]

作者
Author: alan 时间: 2015-12-28 14:26
标题: 通过relaxed scan再TS后，IRC失败了，求分析
大家好，我做这个过渡态的过程用了很多方法，最后还是遇到问题，希望高人指点迷津。我的过程如下：1、自己寻找过渡态，得到虚频-47的结果，然后做IRC（stepsize=3,maxpoints=100,reverse/forward,calcfc）,得到的结果做opt(maxstep=5,notrust)都回到反应物(虽然还没优化成功)。
2、先做断键的relaxed scan(opt=modredundant freq=noraman nosymm)，然后过渡态优化Opt(ts,noeigentest,maxcyc=100,calcall)得到虚频-38的结果，再做IRC（stepsize=3,maxpoints=100,reverse/forward,calcfc），结果如下：a、只跑了2~3个点；b、换LQA算法，跑了三个点正常结束，提示PES minimum detected.；c、换GS2算法，非正常结束，提示GS2 failed.; d、增加步长stepsize=30,并加上calcall，结果是forward跑3个点提示PES minimum detected，reverse跑2个点，l123错误。
综上，会不会是过渡态不对？如果要重新relaxed scan的话，请问需不需要考虑扫描键角（考虑到键长扫描时已经是柔性扫描，包含了优化的过程）？
在线等，以上尝试过的一些方法可能没有对症下药，求指导~谢谢~

作者
Author: liyuanhe211 时间: 2015-12-28 15:53
我有点怀疑这么小的虚频不是正确的过渡态，应检查过渡态震动模式是否正确。

可以的话上传各输出（大的话压缩），看体系长什么样子

不要添加如下关键词：freq=noraman, irc=(reverse/forward)；无特殊情况不建议在calcall的情况下opt=(maxcyc=100)

作者
Author: alan 时间: 2015-12-28 16:28

liyuanhe211 发表于 2015-12-28 15:53
我有点怀疑这么小的虚频不是正确的过渡态，应检查过渡态震动模式是否正确。

可以的话上传各输出（大的话 ...

您好，谢谢你的意见，我上传了输出文件，“scan”文件夹为通过scan后再求过渡态，继而IRC等等，“no scan”为自己寻找的过渡态，烦请查看。
另外有疑问：
1、你说不加freq=noraman，意思是在过渡态计算时不加吗？如果不加的话，怎么看是否虚频？
2、我写的irc(reverse/forward)里意思是两个方向分开做的，不是一起加上去的，reverse或forward加到irc里指定方向应该是要的吧？
3、“无特殊情况不建议在calcall的情况下opt=(maxcyc=100)”，这样耗时会更大？
谢谢~

作者
Author: liyuanhe211 时间: 2015-12-28 16:39
本帖最后由 liyuanhe211 于 2015-12-28 16:47 编辑

alan 发表于 2015-12-28 16:28
您好，谢谢你的意见，我上传了输出文件，“scan”文件夹为通过scan后再求过渡态，继而IRC等等，“no scan ...

去看手册里freq=noraman什么意思，不要乱写不知道什么意思的关键词，和你说的作用毫无关系。另外这关键词加没加完全一样，是默认的。

既然你都要做，为什么还要两个方向分开。

程序默认给出的maxcyc基本是合理的，没遇到与其相关的问题没必要用。
以及这一大串iop用的。。。iop(2/16=1,5/13=1,2/16=2,1/8=2)
吐槽都不知道该从何处下口。。。知道每一个什么意思再用

（其他的回去看，这电脑上没GaussView）

作者
Author: alan 时间: 2015-12-28 16:50

liyuanhe211 发表于 2015-12-28 16:39
去看手册里freq=noraman什么意思，不要乱写不知道什么意思的关键词，和你说的作用毫无关系。另外这关键词 ...

谢谢指正~我好好看看关键词存在的问题

作者
Author: alan 时间: 2015-12-28 21:00

liyuanhe211 发表于 2015-12-28 16:39
去看手册里freq=noraman什么意思，不要乱写不知道什么意思的关键词，和你说的作用毫无关系。另外这关键词 ...

你好，由于之前沿用的关键词是参考别人和网上的，没有细究每个的意思，现在我去高斯手册看了我所用关键词的意思，有一些大致的理解，希望批评指正：
1、IRC要求用内反应坐标，那应该把internal（跟踪内坐标Z-矩阵）添加进irc里吧？
2、一般情况下不添加IOp内部选项，因为它们对计算的影响是不确定的？

作者
Author: liyuanhe211 时间: 2015-12-28 22:47

alan 发表于 2015-12-28 21:00
你好，由于之前沿用的关键词是参考别人和网上的，没有细究每个的意思，现在我去高斯手册看了我所用关键词 ...

不用
不是，iop和关键词一样都是明确的去查iop手册，但显然要先了解它是干啥的，这些iop产生的作用都不是你想要的。具体的去看iop手册，http://sobereva.com/61这里提到了5/13=1，这也不用考虑对称性的问题，还给俩矛盾的。

作者
Author: alan 时间: 2015-12-28 23:15

liyuanhe211 发表于 2015-12-28 22:47
不用
不是，iop和关键词一样都是明确的去查iop手册，但显然要先了解它是干啥的，这些iop产生的作用都不 ...

看了关于5/13=1的批判，哆嗦了一下。。。
请问SCF收敛的自查，是不是查找那四个item是否yes，或者看后两项位移与收敛标准的差距？

作者
Author: alan 时间: 2015-12-28 23:21

liyuanhe211 发表于 2015-12-28 22:47
不用
不是，iop和关键词一样都是明确的去查iop手册，但显然要先了解它是干啥的，这些iop产生的作用都不 ...

请问势能面曲率的检测1/11=2，或eigentest对TS和IRC有必要吗？

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