sobereva 发表于 2024-12-23 16:16 Sob老师您好,我通过你的方法尝试了,最开始我只设置的原来会报错一边,即reverse,能够成功跑出。但我设置两边一起跑准备得到完整的IRC曲线时,原来能够正常跑出的那一边,即forward,报错l502了。曲线和最初的报错曲线很像(图1),仍然像是突然的陡降,请教老师这种情况是为什么?以及可能解决的办法,谢谢老师。 |
wal 发表于 2024-12-23 21:44 好的,谢谢老师!万分感谢!! |
本帖最后由 wal 于 2024-12-23 23:39 编辑 Xx_DerEnTh 发表于 2024-12-23 21:07 我一般用第二种,要用第一种你不得重算一遍TS么,那多麻烦。几何写不写取决于加不加geom=check,guess=read只读取初猜不管几何,如果不写geom=check也可以把TS结构保存出来然后改IRC的关键词。我复制一个我以前做的给你参考一下。 %oldchk=TS.chk %chk=IRC.chk(这个好像可以没有?) # irc=(LQA,rcfc,maxpoint=40) um062x/6-31g(d) nosymm guess=read geom=check title 1 1 注:oldchk的路径要指向你TS的chk文件。 |
wal 发表于 2024-12-23 19:22 老师您好,请问一下用guess=read具体怎么写呢?是在最开始算过渡态的输入文件的下面加Link1读取产生的chk文件写IRC参数重新算过渡态和IRC,还是说新创一个gjf文件,直接读取过渡态产生的chk文件呢?如果是第二种的话,创建gjf文件时,还需要用GaussView构造分子结构吗?还需要写原子坐标吗?没用过guess=read,有点不太懂,麻烦老师指教。 |
wal 发表于 2024-12-23 19:22 好的,谢谢老师,我试试 |
Xx_DerEnTh 发表于 2024-12-23 18:51 保存另算的话有可能没有跟TS任务收敛到同一个波函数上 |
sobereva 发表于 2024-12-23 16:16 感谢Sob老师回复,请问这样操作和把TS.out当成新的输入文件算IRC有什么区别吗? |
wal 发表于 2024-12-23 09:42 试过了 |
guess=read读取TS任务的chk里的波函数当初猜 |
XQC试了吗 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-14 01:01 , Processed in 0.258702 second(s), 25 queries , Gzip On.