请教下使用Dushin分解重组能及HR因子的输出结果为0

zbc201011

各位老师们好，我刚开始学习使用dushin分解重组能和计算H-R因子，想请教下：


我根据sob老师的帖子（使用Dushin分解重组能和计算Huang-Rhys因子，http://bbs.keinsci.com/forum.php ... 07&highlight=Dushin）想得到带一个正电的分子的S0和S1态之间的重组能，但是输出的结果为0.并且我的supplem.*等计算文件都是空的。

我的结构有105个原子，计算方法是B3LYP/def2svp。我在dushin.dat中的首行设置为2 0，使用的计算软件是Gaussian16，我把输出文件中所有的g16均改成了g09。

想请教下各位老师：
① 为什么输出的结果为0，该做什么调整；
② 如果得到了正确的结构，我想将键长、键角和二面角所贡献的重组能投射到内坐标上，以得到各个分量对重组能的贡献的占比，我应该在哪个文件中分析这部分数据？


由于附件太大，我放在了夸克链接中：https://pan.quark.cn/s/7021920806ac
敬请各位老师批评指正！

Donkhk

同问，我也遇到了一样的问题。输出文件最后写的是 Displacement: in terms of nc of 1 THEN of nc of 2