计算化学公社

标题: 两个氢键之间过渡态问题 [打印本页]

作者
Author: 茶新味 时间: 2017-10-9 09:26
标题: 两个氢键之间过渡态问题
Sob老师，您好。我最初研究的是不同弱相互作用之前的过渡态。构型一和中间体之间的过渡态找到，IRC也能走通。
但是中间体IM与构型二（即两个不同位置接的氢键）之间能找到过渡态TS2，趋势也正确，就是走不通IRC。步长由10到30每隔5步，加LQA。
关键字：# mp2/genecp irc=(calcall,maxpoints=300,stepsize=30,LQA)
但是都只走几步就normal，请问sob老师，TS2趋势可靠，怎么才能走出合理的IRC？

作者
Author: sobereva 时间: 2017-10-9 09:34
IM和TS2基本都看不出区别，TS2应当不是真正、准确的过渡态结构
做个8-1-6角度的柔性扫描或许有帮助确认过渡态合适的初猜

作者
Author: 茶新味 时间: 2017-10-9 11:39
sob老师，您好。做了柔性扫描，输入文件如图，结果出现图二出现断点，要取最高点为TS初猜吗？

作者
Author: sobereva 时间: 2017-10-9 16:12

茶新味发表于 2017-10-9 11:39
sob老师，您好。做了柔性扫描，输入文件如图，结果出现图二出现断点，要取最高点为TS初猜吗？

如果那个断点之后都优化成产物结构了，就用这个当初猜
找过渡态步长上限设小点，最好用calcall

作者
Author: 茶新味 时间: 2017-10-9 19:11
sob老师，您说的“找过渡态步长上限设最小值”不太明白？指的是找到过渡态以后走IRC步长为stepsize=5或者maxstep=5,notrustupdate吗？
您说的加calcall指的是在找过渡态的时候opt里面加calcall吗？

作者
Author: sobereva 时间: 2017-10-9 21:50

茶新味发表于 2017-10-9 19:11
sob老师，您说的“找过渡态步长上限设最小值”不太明白？指的是找到过渡态以后走IRC步长为stepsize=5或者ma ...

opt里用maxstep=N和notrust，N设诸如3~10。这些在这里都提了
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

opt里写calcall

作者
Author: 茶新味 时间: 2017-10-10 14:54
sob老师，您好。根据初猜找到ts，如图一振动方向跟两端构型相似，但是走IRC的时候，步长由5-30每隔5设置一步，irc直走4步就结束，是不是应为弱相互作用过渡态比较敏感，用高斯默认的走irc走不通呢？

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