求助：Dushin计算HR因子出错

ChunLinX

各位老师好，我在用dushin计算电荷分离过程的内重组能和Hung-Rhys因子时，用D*A与D+A-两个态的平衡点几何结构计算频率，再将Gaussian09频率计算的log文件和fchk文件用来作为dushin的输入文件。参照之前的博文，我将泛函基组改成了B3LYP/6-31G*。奇怪的是，计算一个几十个原子的分子不会报错，但是当计算两个分子放在一起的体系的时候就会出错。之前有人认为是因为有虚频，但是我计算了有虚频的单体发现也没有报错，且用Windows和Linux版本计算均是同样的错误。报错内容请看下图。想向老师们请教出错的原因和解决方法，万分感谢。

ChunLinX

顶一下，不要沉了

破晓时的曙光

您好，请问您解决了这个问题了吗，我也遇到相同的问题

ChunLinX

破晓时的曙光 发表于 2019-12-17 13:50
您好，请问您解决了这个问题了吗，我也遇到相同的问题

不好意思，现在还是没有解决

破晓时的曙光

ChunLinX 发表于 2020-2-12 22:28
不好意思，现在还是没有解决

那请问有没有什么替代的方法呢

ChunLinX

破晓时的曙光 发表于 2020-2-13 12:26
那请问有没有什么替代的方法呢

如果是计算相同分子间的电荷转移的话，可以直接用单体来计算。如果是计算电荷分离（两个不同分子）也有文献用单体来算，比如 J. Phys. Chem. C 2011, 115, 2406–2415，不过对于这种情况，黄因子大部分是由低频振动（两分子相向振动）贡献的，如  J. Phys. Chem. C 2014, 118, 9780−9789，因此可以直接用高温极限而不考虑隧道效应。即使要考虑也可以像以前的文献那样分别计算给体和受体的黄因子，不必计算双体体系。简单来说就是直接计算单体的黄因子，这样应该不会报错。