计算化学公社

标题: IRC只能向一边跑，如何解决？ [打印本页]

作者
Author: 一直在路上 时间: 2017-12-12 10:42
标题: IRC只能向一边跑，如何解决？
对过渡态进行IRC验证，结果IRC只能跑一半，无论加reverse还是forward都不行？请问该如何解决？

作者
Author: xujian 时间: 2017-12-12 14:18
1.先尝试用calcfc或calcall、opt=tight. 2.增大步长

作者
Author: 一直在路上 时间: 2017-12-12 14:30

xujian 发表于 2017-12-12 14:18
1.先尝试用calcfc或calcall、opt=tight. 2.增大步长

我用ccsd方法，用上述关键词的话没有二阶解析函数。增大步长一般不起作用，不过我还是再加到200试一下吧。谢谢你

作者
Author: xujian 时间: 2017-12-12 15:20

一直在路上发表于 2017-12-12 14:30
我用ccsd方法，用上述关键词的话没有二阶解析函数。增大步长一般不起作用，不过我还是再加到200试一下吧 ...

服务器好牛逼，用CCSD走IRC

作者
Author: 一直在路上 时间: 2017-12-12 15:37

xujian 发表于 2017-12-12 15:20
服务器好牛逼，用CCSD走IRC

算的慢

作者
Author: xylz6188 时间: 2017-12-12 19:00
ccsd优化输入方法是怎么样的啊，我之前一直想用它来做优化，就是一直报错

作者
Author: sobereva 时间: 2017-12-12 20:01

xylz6188 发表于 2017-12-12 19:00
ccsd优化输入方法是怎么样的啊，我之前一直想用它来做优化，就是一直报错

# opt CCSD/cc-pVTZ

作者
Author: sobereva 时间: 2017-12-12 20:02
优化过渡态时用tight收敛限，之后走IRC再试
步长不可能增加到200，绝对死掉。适当增大步长有可能解决当前问题。另外也可以尝试改用GS2

作者
Author: 一直在路上 时间: 2017-12-12 21:23

sobereva 发表于 2017-12-12 20:02
优化过渡态时用tight收敛限，之后走IRC再试
步长不可能增加到200，绝对死掉。适当增大步长有可能解决当前 ...

谢谢sob老师！
过渡态我用tight重新计算一下，再跑IRC看看结果。
步长增加到200果断死掉

作者
Author: xylz6188 时间: 2017-12-12 21:26

sobereva 发表于 2017-12-12 20:01
# opt CCSD/cc-pVTZ

谢谢sobereva老师！我试一试。另外前次在群里请教你计算偶极矩的方法，你推荐用b3lyp/Def2-QZVPD。我用这个方法和基组来做反应机理，优化反应复合物，最后变成产物了，过渡态有很多虚频，这是怎么回事sobereva老师，再次请教一下！

作者
Author: sobereva 时间: 2017-12-12 21:50

xylz6188 发表于 2017-12-12 21:26
谢谢sobereva老师！我试一试。另外前次在群里请教你计算偶极矩的方法，你推荐用b3lyp/Def2-QZVPD。我用这 ...

我推荐的是def2-TZVPD。没必要QZ
算反应机理用def2-TZVP就已经足够了，不要加D，也不要用QZ，否则白费时间。
你当前的问题不是基组的事

作者
Author: lip 时间: 2017-12-13 15:54
这一个问题Sob老师之前也解释过，可以在站内搜索一下。

看你的问题应该是过渡态位置不正确，找到的不是s=0的位置，而是在“坡”上。
先检查过渡态振动模式是否正确，在看虚频大小是否合适

作者
Author: 一直在路上 时间: 2017-12-13 20:38

lip 发表于 2017-12-13 15:54
这一个问题Sob老师之前也解释过，可以在站内搜索一下。

看你的问题应该是过渡态位置不正确，找到的不是s ...

过渡态振动模式正确，通过其他方法基组已经优化出TS并且进行了IRC的验证，现在换成CCSD方法重新进行计算，TS也优化出来，且振动模式正确，但就是IRC前半部分跑不出来

作者
Author: 一直在路上 时间: 2017-12-14 08:13

sobereva 发表于 2017-12-12 20:02
优化过渡态时用tight收敛限，之后走IRC再试
步长不可能增加到200，绝对死掉。适当增大步长有可能解决当前 ...

sob老师，用tight重新计算了过渡态，再跑irc仍然只有一半。

作者
Author: sobereva 时间: 2017-12-14 16:35

一直在路上发表于 2017-12-14 08:13
sob老师，用tight重新计算了过渡态，再跑irc仍然只有一半。

步长比默认增大一倍看看
其实用双杂化怎么也够了，此时用calcfc、recalc都方便

作者
Author: 一直在路上 时间: 2017-12-14 19:49

sobereva 发表于 2017-12-14 16:35
步长比默认增大一倍看看
其实用双杂化怎么也够了，此时用calcfc、recalc都方便

我用CCSD跑的irc，加calcfc的话没有解析二阶倒数。目前我用的是rcfc

作者
Author: lijiayisjtu 时间: 2018-4-8 14:23

一直在路上发表于 2017-12-14 08:13
sob老师，用tight重新计算了过渡态，再跑irc仍然只有一半。

应该不是计算方法和基组的问题，想问下您，这问题后来是怎么解决的，IRC 只跑了一半

作者
Author: sobereva 时间: 2018-4-9 02:34

lijiayisjtu 发表于 2018-4-8 14:23
应该不是计算方法和基组的问题，想问下您，这问题后来是怎么解决的，IRC 只跑了一半

在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400

作者
Author: 自制奶酪 时间: 2019-4-27 14:39

一直在路上发表于 2017-12-14 08:13
sob老师，用tight重新计算了过渡态，再跑irc仍然只有一半。

你好。我对过渡态做IRC验证的时候也是只跑出一半，请问你怎么解决的呢，能告诉我一下吗

作者
Author: sobereva 时间: 2019-4-28 04:25

自制奶酪发表于 2019-4-27 14:39
你好。我对过渡态做IRC验证的时候也是只跑出一半，请问你怎么解决的呢，能告诉我一下吗

这里说得极为清楚
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com/thread-7736-1-1.html）

作者
Author: sun20201028 时间: 2020-11-30 09:27

一直在路上发表于 2017-12-14 08:13
sob老师，用tight重新计算了过渡态，再跑irc仍然只有一半。

楼主，请问

这个问题解决了吗？

作者
Author: 小西瓜123 时间: 2022-1-22 11:02

xujian 发表于 2017-12-12 15:20
服务器好牛逼，用CCSD走IRC

你好，请问您用ccsd做IRC的输入文件怎么写的呢？我先做了半数值的freq, 然后读取chk,readfc,还是没算就报错了，您可以分享一下您的输入文件吗？我的输入文件如下
%chk=r3-ts-irc-test
%nprocshared=8
%mem=16gb
# ccsd(maxcycle=500)/aug-cc-pvtz freq

freq

0 1
O                0.05120100 1.24334900 0.00000000
H                0.06339300 0.07509000 0.00000000
O                0.05120100 -1.08448300 0.00000000
H                -0.88260300 -1.34602500 0.00000000

--Link1--
%chk=r3-ts-irc-test
# CCSD/aug-cc-pVtZ Geom=Check Guess=Read IRC=(Maxpoints=40,readFC,lqa)

irc

0 1

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