计算化学公社

 找回密码 Forget password
 注册 Register
Views: 414|回复 Reply: 8
打印 Print 上一主题 Last thread 下一主题 Next thread

[辅助/分析程序] 使用fcclasses3计算IC时出现WARNING: IC rate could not be computed

[复制链接 Copy URL]

6

帖子

0

威望

161

eV
积分
167

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
刚开始学习fcclasses3,计算IC时出现WARNING: IC rate could not be computed. Change TD settings ,将input里面TD改为TI后还是不行,请问是哪里出了问题呢,应该如何修改呢
附件是计算的文件
                                            



IC.zip

1.23 MB, 下载次数 Times of downloads: 15

449

帖子

0

威望

1561

eV
积分
2010

Level 5 (御坂)

9#
发表于 Post on 2024-12-30 22:58:33 | 只看该作者 Only view this author
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的发光。

6

帖子

0

威望

161

eV
积分
167

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2024-12-30 12:35:18 | 只看该作者 Only view this author
cokie 发表于 2024-12-30 11:54
别用FC301了,改用最新版304

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

449

帖子

0

威望

1561

eV
积分
2010

Level 5 (御坂)

7#
发表于 Post on 2024-12-30 11:54:13 | 只看该作者 Only view this author
20x_deng02 发表于 2024-12-29 20:15
我该用INTERNAL后还是有两个分子算出来是负的IC速率,可以请教一下怎么用momap算IC吗

别用FC301了,改用最新版304

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
wal + 1 正解

查看全部评分 View all ratings

420

帖子

5

威望

1676

eV
积分
2196

Level 5 (御坂)

鸩羽

6#
发表于 Post on 2024-12-30 08:44:43 | 只看该作者 Only view this author
本帖最后由 wal 于 2024-12-30 12:25 编辑
20x_deng02 发表于 2024-12-29 20:15
我该用INTERNAL后还是有两个分子算出来是负的IC速率,可以请教一下怎么用momap算IC吗

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

我没看到你的版号,楼下说的对,301版本负IC速率的问题在304已经修了
某不知名实验组从苞米地里长出来的计算选手

6

帖子

0

威望

161

eV
积分
167

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2024-12-29 20:15:34 | 只看该作者 Only view this author
wal 发表于 2024-12-28 23:30
CARTESIAN算内转换有时候是会这样,我的经验来讲一般是duschinsky混合显著的容易出负IC速率。改用INTERNA ...

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

420

帖子

5

威望

1676

eV
积分
2196

Level 5 (御坂)

鸩羽

4#
发表于 Post on 2024-12-28 23:30:16 | 只看该作者 Only view this author
本帖最后由 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算了一下能算出来正常的数据,就没管
某不知名实验组从苞米地里长出来的计算选手

6

帖子

0

威望

161

eV
积分
167

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2024-12-28 23:26:17 | 只看该作者 Only view this author
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吗

420

帖子

5

威望

1676

eV
积分
2196

Level 5 (御坂)

鸩羽

2#
发表于 Post on 2024-12-28 20:28:47 | 只看该作者 Only view this author
把state1的fcc文件和state2的交换一下位置试试
某不知名实验组从苞米地里长出来的计算选手

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-17 12:38 , Processed in 0.211177 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list