计算化学公社

标题: 求助，跑锕系元素IRC时，报错l502 [打印本页]

作者
Author: Xx_DerEnTh 时间: 2024-12-23 07:15
标题: 求助，跑锕系元素IRC时，报错l502
本帖最后由 Xx_DerEnTh 于 2024-12-23 07:43 编辑

各位老师好，又来麻烦你们了。
我模拟一篇文献，计算含Pu体系的过渡态和IRC，唯一的改动是将文献中的6-311++G(d,p)基组改为了6-311G(d,p)。
其过渡态能够成功跑出，并且通过虚频以及振动方向初步判断该过渡态是正确的。但是跑IRC时始终报错l502。
我已认真阅读了Sob老师发布的http://sobereva.com/61,http://sobereva.com/581,http://sobereva.com/400帖子，考虑到IRC的级别要和过渡态一致，已经尝试了多种符合条件的解决SCF不收敛的方法仍然没有解决问题。
在跑IRC的过程中，生成物的方向能够正常跑出，但对于反应物方向，在能量-610.27左右(图1)，就会出现报错，本认为是步长设置得比较大造成的，于是尝试减小了步长，但问题仍然存在，并且尝试的每一种解决SCF不收敛的方法大多都是在这里出现报错。
过渡态和IRC的输入输出文件已上传，其中IRC的输入和输出文件是最后一次减小步长进行尝试的文件，IRC的输出文件因过大进行了压缩。想请教各位老师，这个问题有没有好的解决办法？

作者
Author: wal 时间: 2024-12-23 09:42
XQC试了吗

作者
Author: sobereva 时间: 2024-12-23 16:16
guess=read读取TS任务的chk里的波函数当初猜

作者
Author: Xx_DerEnTh 时间: 2024-12-23 18:50

wal 发表于 2024-12-23 09:42
XQC试了吗

试过了

作者
Author: Xx_DerEnTh 时间: 2024-12-23 18:51

sobereva 发表于 2024-12-23 16:16
guess=read读取TS任务的chk里的波函数当初猜

感谢Sob老师回复，请问这样操作和把TS.out当成新的输入文件算IRC有什么区别吗？

作者
Author: wal 时间: 2024-12-23 19:22

Xx_DerEnTh 发表于 2024-12-23 18:51
感谢Sob老师回复，请问这样操作和把TS.out当成新的输入文件算IRC有什么区别吗？

保存另算的话有可能没有跟TS任务收敛到同一个波函数上

作者
Author: Xx_DerEnTh 时间: 2024-12-23 19:28

wal 发表于 2024-12-23 19:22
保存另算的话有可能没有跟TS任务收敛到同一个波函数上

好的，谢谢老师，我试试

作者
Author: Xx_DerEnTh 时间: 2024-12-23 21:07

wal 发表于 2024-12-23 19:22
保存另算的话有可能没有跟TS任务收敛到同一个波函数上

老师您好，请问一下用guess=read具体怎么写呢？是在最开始算过渡态的输入文件的下面加Link1读取产生的chk文件写IRC参数重新算过渡态和IRC，还是说新创一个gjf文件，直接读取过渡态产生的chk文件呢？如果是第二种的话，创建gjf文件时，还需要用GaussView构造分子结构吗？还需要写原子坐标吗？没用过guess=read，有点不太懂，麻烦老师指教。

作者
Author: wal 时间: 2024-12-23 21:44
本帖最后由 wal 于 2024-12-23 23:39 编辑

Xx_DerEnTh 发表于 2024-12-23 21:07
老师您好，请问一下用guess=read具体怎么写呢？是在最开始算过渡态的输入文件的下面加Link1读取产生的chk ...

我一般用第二种，要用第一种你不得重算一遍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文件。

作者
Author: Xx_DerEnTh 时间: 2024-12-23 22:41

wal 发表于 2024-12-23 21:44
我一般用第二种，要用第一种你不得重算一遍TS么，那多麻烦。几何写不写取决于加不加geom=check，guess=re ...

好的，谢谢老师！万分感谢！！

作者
Author: Xx_DerEnTh 时间: 2024-12-24 18:08

sobereva 发表于 2024-12-23 16:16
guess=read读取TS任务的chk里的波函数当初猜

Sob老师您好，我通过你的方法尝试了，最开始我只设置的原来会报错一边，即reverse，能够成功跑出。但我设置两边一起跑准备得到完整的IRC曲线时，原来能够正常跑出的那一边，即forward，报错l502了。曲线和最初的报错曲线很像(图1)，仍然像是突然的陡降，请教老师这种情况是为什么？以及可能解决的办法，谢谢老师。

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