计算化学公社

标题: ORCA走IRC时可以输出每一步的具体信息吗？ [打印本页]

作者
Author: 霜晨月 时间: 2023-10-5 10:07
标题: ORCA走IRC时可以输出每一步的具体信息吗？
用ORCA走IRC，输出文件中只有每一步的能量，另外有一个trj.xyz文件记录每一步的坐标。
能否查看每一步的S**2、电子激发等信息？是否有关键词可以实现？
谢谢

作者
Author: ahxb 时间: 2023-10-5 10:44
本帖最后由 ahxb 于 2023-10-5 10:47 编辑

在RKS计算中将%irc中的PrintLevel设为2，可以输出SCF迭代过程、SCF能量、轨道能量和梯度。如果是UKS计算，应当能输出<S**2>。也可以尝试将PrintLevel设为更大的值，手册9.25.1节详细介绍了%irc块，但没有说明PrintLevel设为不同值时的作用。
我不清楚你说的电子激发指的是什么。

作者
Author: 霜晨月 时间: 2023-10-5 12:42

ahxb 发表于 2023-10-5 10:44
在RKS计算中将%irc中的PrintLevel设为2，可以输出SCF迭代过程、SCF能量、轨道能量和梯度。如果是UKS计算， ...

谢谢老师，我试试。
我说的电子激发信息指的是激发能、轨道跃迁等信息，是不是这样叫不合适

作者
Author: ahxb 时间: 2023-10-5 15:27

霜晨月发表于 2023-10-5 12:42
谢谢老师，我试试。
我说的电子激发信息指的是激发能、轨道跃迁等信息，是不是这样叫不合适

IRC计算过程中只计算能量和梯度，不做任何关于激发态的计算，也就没有激发能和轨道跃迁等。想得到的话只能每个结构单独计算激发态。

作者
Author: 霜晨月 时间: 2023-10-5 18:05

ahxb 发表于 2023-10-5 15:27
IRC计算过程中只计算能量和梯度，不做任何关于激发态的计算，也就没有激发能和轨道跃迁等。想得到的话只 ...

知道了，谢谢老师。
感觉很奇怪，TDDFT指定某个iroot做IRC也不做激发态计算么？

作者
Author: ahxb 时间: 2023-10-5 23:16
本帖最后由 ahxb 于 2023-10-5 23:21 编辑

霜晨月发表于 2023-10-5 18:05
知道了，谢谢老师。
感觉很奇怪，TDDFT指定某个iroot做IRC也不做激发态计算么？

原来是在激发态势能面上跑IRC。我没做过这样的计算，不过可以介绍我在基态势能面上跑IRC的计算经验。
ORCA计算像IRC，数值频率等有很多几乎完全相同的计算时，会对每个计算提供单独的basename。比如计算IRC时，如果以file.inp为输入文件，那么在计算IRC的forward方向第1个点时，会以file_IRC_F_001作为basename执行一个计算，出现很多以file_IRC_F_001开头的文件，其中有一些以.last*结尾的文件是单独计算时的输出内容，比如.lastscf中含有scf迭代和收敛部分，.lastgrad含有梯度计算部分。如果是在激发态势能面上计算IRC，可能也有对应的输出文件。可惜的是通常这一步计算结束后对应的.last*文件会被立即删除。你可以想想办法怎么保存这些文件。比如每个点应该至少按顺序做SCF、TD-DFT、梯度计算（可能还有别的），那么你可以以一定的时间间隔备份文件，趁TD-DFT已经算完，对应部分已经写入某个文件，但梯度还在计算的时候备份这个文件。

作者
Author: 霜晨月 时间: 2023-10-6 09:19

ahxb 发表于 2023-10-5 23:16
原来是在激发态势能面上跑IRC。我没做过这样的计算，不过可以介绍我在基态势能面上跑IRC的计算经验。
OR ...

好的，谢谢老师指导

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