FCclasses3.0计算IC速率求助

xezhang · 发表于 Post on 2023-5-15 15:59:14

各位老师好，最近在学习用FCclasses3.0计算S1-S0的内转换速率Kic，需要非绝热耦合的结果文件。看了使用手册，还是不太明白。想请教大家这个文件用高斯怎么得到呢？谢谢大家

xuhj199508 · 发表于 Post on 2023-5-15 17:21:44

用g16算一下激发态得到S1.fchk，然后用gen_fcc_dipfile -i S1.fchk -nac就可以得到所需nac文件了。

xezhang · 发表于 Post on 2023-5-15 19:37:17

xuhj199508 发表于 2023-5-15 17:21
用g16算一下激发态得到S1.fchk，然后用gen_fcc_dipfile -i S1.fchk -nac就可以得到所需nac文件了。

哇塞！真的可以！太感谢啦！！！

xezhang · 发表于 Post on 2023-5-15 20:42:54

xuhj199508 发表于 2023-5-15 17:21
用g16算一下激发态得到S1.fchk，然后用gen_fcc_dipfile -i S1.fchk -nac就可以得到所需nac文件了。

你好~刚刚用FCclasses3算了一个分子的IC速率，都正常结束了，但是在fcc.out文件里显示IC值为0，请教一下这是什么原因呀？相关的输入输出文件已上传到附件中。麻烦你啦~谢谢你喔

xuhj199508 · 发表于 Post on 2023-5-16 08:38:38

xezhang 发表于 2023-5-15 20:42
你好~刚刚用FCclasses3算了一个分子的IC速率，都正常结束了，但是在fcc.out文件里显示IC值为0，请教一下 ...

你的nac文件中的元素都为0，我用你的激发态文件生成的nac文件见附件。

xezhang · 发表于 Post on 2023-5-16 09:46:33

xuhj199508 发表于 2023-5-16 08:38
你的nac文件中的元素都为0，我用你的激发态文件生成的nac文件见附件。

我刚刚用你给的nac文件试着算了一下IC速率，得到的数值是正常的耶！太感谢啦！可以再麻烦您给一个用g16计算nac的激发态S1的输入文件示例吗？

xuhj199508 · 发表于 Post on 2023-5-16 13:26:51

xezhang 发表于 2023-5-16 09:46
我刚刚用你给的nac文件试着算了一下IC速率，得到的数值是正常的耶！太感谢啦！可以再麻烦您给一个用g16计 ...

就按你现在这么算就可以了。

xm.liu · 发表于 Post on 2023-9-7 13:29:04

xuhj199508 发表于 2023-5-16 08:38
你的nac文件中的元素都为0，我用你的激发态文件生成的nac文件见附件。

您好，我用FCclass3.0计算好几个分子的IC，但是出现有一个分子IC是负值，请问您遇到过这种情况吗？

xuhj199508 · 发表于 Post on 2023-9-7 13:40:57

xm.liu 发表于 2023-9-7 13:29
您好，我用FCclass3.0计算好几个分子的IC，但是出现有一个分子IC是负值，请问您遇到过这种情况吗？

如果当前COORDS选择的是CARTESIAN，可以换成INTERNAL试试。

xm.liu · 发表于 Post on 2023-9-7 13:45:29

本帖最后由 xm.liu 于 2023-9-7 13:55 编辑

xuhj199508 发表于 2023-9-7 13:40
如果当前COORDS选择的是CARTESIAN，可以换成INTERNAL试试。

好的好的，感谢您的回复，就是我计算前几个都是用的CARTESIAN（我看手册的绝热过程使用的都是CARTESIAN）,算出来的数值很正常，就是这个是负数，这两种坐标形式对计算结果有啥影响呢

xuhj199508 · 发表于 Post on 2023-9-7 15:21:11

xm.liu 发表于 2023-9-7 13:45
好的好的，感谢您的回复，就是我计算前几个都是用的CARTESIAN（我看手册的绝热过程使用的都是CARTESIAN） ...

对于AH和VH模型，如果初态和末态结构变化较大，使用CARTESIAN可能没法很好的描述分子的扭转和旋转。

xm.liu · 发表于 Post on 2023-9-7 17:15:00

xuhj199508 发表于 2023-9-7 15:21
对于AH和VH模型，如果初态和末态结构变化较大，使用CARTESIAN可能没法很好的描述分子的扭转和旋转。

非常感谢您的回复，我尝试了更换坐标类型，得到了正的IC值，但是out文件里输出的重组能值两种坐标之下算出来的相差巨大，这是什么原因呢？因为我之前用别的方法计算了黄昆因子（计算结果与CARTESIAN坐标系计算结果一致），所以我的计算结果还是有点问题，您有什么建议吗？

wzkchem5 · 发表于 Post on 2023-9-7 17:32:03

xm.liu 发表于 2023-9-7 10:15
非常感谢您的回复，我尝试了更换坐标类型，得到了正的IC值，但是out文件里输出的重组能值两种坐标之下算 ...

看Phys. Chem. Chem. Phys., 2012, 14, 13549–13563

xm.liu · 发表于 Post on 2023-9-7 19:27:29

wzkchem5 发表于 2023-9-7 17:32
看Phys. Chem. Chem. Phys., 2012, 14, 13549–13563

谢谢老师的回复

xuhj199508 · 发表于 Post on 2023-9-8 08:43:48

xm.liu 发表于 2023-9-7 17:15
非常感谢您的回复，我尝试了更换坐标类型，得到了正的IC值，但是out文件里输出的重组能值两种坐标之下算 ...

一般两种COORDS得到的重组能不会相差很多，所以直接采用internal的就可以了。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[辅助/分析程序] FCclasses3.0计算IC速率求助