xujian 发表于 2017-12-12 15:20 你好，请问您用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 |
一直在路上 发表于 2017-12-14 08:13 楼主，请问  这个问题解决了吗？ |
自制奶酪 发表于 2019-4-27 14:39 这里说得极为清楚 在Gaussian中计算IRC的方法和常见问题 http://sobereva.com/400（http://bbs.keinsci.com/thread-7736-1-1.html） |
一直在路上 发表于 2017-12-14 08:13 你好。我对过渡态做IRC验证的时候也是只跑出一半，请问你怎么解决的呢，能告诉我一下吗 |
lijiayisjtu 发表于 2018-4-8 14:23 在Gaussian中计算IRC的方法和常见问题 http://sobereva.com/400 |
一直在路上 发表于 2017-12-14 08:13 应该不是计算方法和基组的问题，想问下您，这问题后来是怎么解决的，IRC 只跑了一半 |
sobereva 发表于 2017-12-14 16:35 我用CCSD跑的irc，加calcfc的话没有解析二阶倒数。目前我用的是rcfc |
一直在路上 发表于 2017-12-14 08:13 步长比默认增大一倍看看 其实用双杂化怎么也够了，此时用calcfc、recalc都方便 |
sobereva 发表于 2017-12-12 20:02 sob老师，用tight重新计算了过渡态，再跑irc仍然只有一半。 |
lip 发表于 2017-12-13 15:54 过渡态振动模式正确，通过其他方法基组已经优化出TS并且进行了IRC的验证，现在换成CCSD方法重新进行计算，TS也优化出来，且振动模式正确，但就是IRC前半部分跑不出来 |
|
这一个问题Sob老师之前也解释过，可以在站内搜索一下。 看你的问题应该是过渡态位置不正确，找到的不是s=0的位置，而是在“坡”上。 先检查过渡态振动模式是否正确，在看虚频大小是否合适 |
xylz6188 发表于 2017-12-12 21:26 我推荐的是def2-TZVPD。没必要QZ 算反应机理用def2-TZVP就已经足够了，不要加D，也不要用QZ，否则白费时间。 你当前的问题不是基组的事 |
sobereva 发表于 2017-12-12 20:01 谢谢sobereva老师！我试一试。另外前次在群里请教你计算偶极矩的方法，你推荐用b3lyp/Def2-QZVPD。我用这个方法和基组来做反应机理，优化反应复合物，最后变成产物了，过渡态有很多虚频，这是怎么回事sobereva老师，再次请教一下！ |
sobereva 发表于 2017-12-12 20:02 谢谢sob老师！ 过渡态我用tight重新计算一下，再跑IRC看看结果。 步长增加到200果断死掉 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-21 09:16 , Processed in 0.196840 second(s), 31 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites