计算化学公社

标题: IRC计算过渡态优化问题求助各位老师 [打印本页]

作者
Author: fresco604 时间: 2016-10-12 22:10
标题: IRC计算过渡态优化问题求助各位老师
我想算一个分子的Isomerization

从conformer A出发，越过一个TS到达conformer B, TS我用了两种办法计算QST3和 TS(Berny), 这两个方法计算出来的结果相差很小，只有几个dihedral angle略微的差别，分别叫做TS1 和 TS2

奇怪的事情是我用TS1跑IRC，得到的曲线是从conformer A开始，到·conformer A结束，对称图形
我用TS2跑IRC, 得到的曲线是从conformer B开始，到·conformer B结束，对称图形

之前经过老师指教，对称图形可能是因为TS找的不准确，于是我也加上了一些关键字
# opt=(ts,modredundant,noeigen,calcall) freq b3lyp/6-31+g(d) geom=conn
ectivity int=grid=ultrafine scf=tight empiricaldispersion=gd3bj

# b3lyp/6-31+g(d) empiricaldispersion=gd3bj opt=(qst3,tight,calcall) f
req scf=(maxconventionalcycles=700,xqc) int=grid=ultrafine

但结果还是一样，求助各位老师，这种情况下该怎么办，我能不能在两个TS之间（TS1,TS2）再找一个TS出来？

[attach]6096[/attach]

作者
Author: fresco604 时间: 2016-10-12 22:10

作者
Author: sobereva 时间: 2016-10-13 04:36
走IRC时候把步长设大点，同时用LQA

TS1和TS2之间不可能有其它TS

你加的一些关键词完全没有丝毫意义，诸如modredundant、scf=(maxconventionalcycles=700,xqc) 。用了calcall也不需要再写freq

作者
Author: fresco604 时间: 2016-10-13 09:03

sobereva 发表于 2016-10-13 04:36
走IRC时候把步长设大点，同时用LQA

TS1和TS2之间不可能有其它TS

谢谢老师指点，之前算过渡态总失败，不知道问题出在哪，就加了一些关键词。其实好像主要是noeigen有用

作者
Author: xujian 时间: 2016-10-13 09:21
本帖最后由 xujian 于 2016-10-13 09:27 编辑

sobereva 发表于 2016-10-13 04:36
走IRC时候把步长设大点，同时用LQA

TS1和TS2之间不可能有其它TS

请问sob老师，为什么用了calcall就不用freq了。

不好意思，刚刚看了高斯说明书，看到解释了，Opt=CalcAll 关键字在几何优化的每一点进行解析的二阶导数计算。一旦完成优化，可以得到频率分析的所有必需信息。于是完成频率分析，并把计算结果作为频率任务的存档。

作者
Author: sobereva 时间: 2016-10-13 18:22

fresco604 发表于 2016-10-13 09:03
谢谢老师指点，之前算过渡态总失败，不知道问题出在哪，就加了一些关键词。其实好像主要是noeigen有用

opt=ts的时候才有必要noeigen，否则自己很难摆不出来那么好的只有一个虚频的初猜。QST用不着noeigen，因为根本不做本征值检测。

作者
Author: fresco604 时间: 2016-10-13 20:56

sobereva 发表于 2016-10-13 18:22
opt=ts的时候才有必要noeigen，否则自己很难摆不出来那么好的只有一个虚频的初猜。QST用不着noeigen，因 ...

谢谢老师指教！我把步长分别增加到了10/20/30/40/50/60/70 希望这次能算出来连接A和B的IRC, 要不然明明知道这个TS和TS1/TS2相差无几，但就是算不出来，真是令人着急啊

作者
Author: wyy015020 时间: 2019-1-15 09:23

fresco604 发表于 2016-10-13 20:56
谢谢老师指教！我把步长分别增加到了10/20/30/40/50/60/70 希望这次能算出来连接A和B的IRC, 要不然明明知 ...

请问您最后怎么解决的这个问题啊？谢谢

作者
Author: 宋程程 时间: 2023-4-18 22:58

wyy015020 发表于 2019-1-15 09:23
请问您最后怎么解决的这个问题啊？谢谢

求助您最后怎么解决的啊感谢

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