计算化学公社

标题: 使用fcclasses3计算IC时出现WARNING: IC rate could not be computed [打印本页]

作者
Author: 20x_deng02 时间: 2024-12-28 17:04
标题: 使用fcclasses3计算IC时出现WARNING: IC rate could not be computed
刚开始学习fcclasses3，计算IC时出现WARNING: IC rate could not be computed. Change TD settings ，将input里面TD改为TI后还是不行，请问是哪里出了问题呢，应该如何修改呢
附件是计算的文件

作者
Author: wal 时间: 2024-12-28 20:28
把state1的fcc文件和state2的交换一下位置试试

作者
Author: 20x_deng02 时间: 2024-12-28 23:26

wal 发表于 2024-12-28 20:28
把state1的fcc文件和state2的交换一下位置试试

换了之后就可以了，非常感谢
但是有一个分子计算后的IC trate很奇怪
IC rate constant (s-1) -3.503-191
(for Ead = 1.746 eV)
kic_vs_Ead_TD.dat文件里的数据也很奇怪
# Ead(eV) k_ic(s-1)
-0.43800 -0.45984025-191
-0.43700 -0.46161939-191
-0.43600 -0.46339843-191
-0.43500 -0.46517737-191
-0.43400 -0.46695621-191
-0.43300 -0.46873495-191
或许是-191表示E09吗

作者
Author: wal 时间: 2024-12-28 23:30
本帖最后由 wal 于 2024-12-28 23:33 编辑

20x_deng02 发表于 2024-12-28 23:26
换了之后就可以了，非常感谢
但是有一个分子计算后的IC trate很奇怪
IC rate constant (s-1) -3.503-1 ...

CARTESIAN算内转换有时候是会这样，我的经验来讲一般是duschinsky混合显著的容易出负IC速率。改用INTERNAL模式基本都能解决问题。

没留意看你的输入，你如果用了INTERNAL还是这样的话我也不知道咋搞了，上次遇到的时候借了别人的MOMAP算了一下能算出来正常的数据，就没管

作者
Author: 20x_deng02 时间: 2024-12-29 20:15

wal 发表于 2024-12-28 23:30
CARTESIAN算内转换有时候是会这样，我的经验来讲一般是duschinsky混合显著的容易出负IC速率。改用INTERNA ...

我该用INTERNAL后还是有两个分子算出来是负的IC速率，可以请教一下怎么用momap算IC吗

作者
Author: wal 时间: 2024-12-30 08:44
本帖最后由 wal 于 2024-12-30 12:25 编辑

20x_deng02 发表于 2024-12-29 20:15
我该用INTERNAL后还是有两个分子算出来是负的IC速率，可以请教一下怎么用momap算IC吗

这个如果你有版权的话直接问客服吧没版权算了发不得

我没看到你的版号，楼下说的对，301版本负IC速率的问题在304已经修了

作者
Author: cokie 时间: 2024-12-30 11:54

20x_deng02 发表于 2024-12-29 20:15
我该用INTERNAL后还是有两个分子算出来是负的IC速率，可以请教一下怎么用momap算IC吗

别用FC301了，改用最新版304

作者
Author: 20x_deng02 时间: 2024-12-30 12:35

cokie 发表于 2024-12-30 11:54
别用FC301了，改用最新版304

你好，我用FC304算了一下，但是这个分子还是负IC

作者
Author: cokie 时间: 2024-12-30 22:58

20x_deng02 发表于 2024-12-30 12:35
你好，我用FC304算了一下，但是这个分子还是负IC

还是负的话，FCclasses上应该没有什么解决方案了。
我采用你的结构在3.0.3上也试了，Internal和Cartesian的结果都是NaN，FC主页上的更新日志上也一直在提；改了用TD方法算IC出bug的问题，303改了304也还在改。就目前你的结构来看，Cartesian和Internal算出的重组能都正常，说明KK值不大，谐振模型能用。大概率还是bug吧。

临时想到：你试试把你材料的链子去掉，重新优化一个S0和S1，再试一遍用304（目前来看304应该是算IC最稳健的）。我曾经遇到过一次，删了链子IC就正常了，反正链子也不影响你计算S1的发光。

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