计算化学公社

标题: 求助 dushin 分解重组能 [打印本页]

作者
Author: 玛卡胖不拉叽 时间: 2025-11-5 22:56
标题: 求助 dushin 分解重组能
本帖最后由玛卡胖不拉叽于 2025-12-8 16:09 编辑

您好，我根据使用Dushin分解重组能和计算Huang-Rhys因子 - 思想家公社的门口：量子化学·分子模拟·二次元这个帖子，计算了重组能但出现了下面问题。

① 首先，我们用#P，在g16下计算了小分子（和萘类似的结构），在将log文件里面的g16全部替换成g09，再用dushin可以进行分解；
   但是  dushin.out 输出文件中没有这三个信息：Dushinsky matrix, ncs 2 in terms of ncs 1  、Dushinsky matrix, ncs 1 in terms of ncs 2 、Displacement: in terms of nc of 1 THEN of nc of 2；
   重组能计算的结果也是零；
   尝试了两种泛函基组#P opt freq b3lyp/6-31g(d)和#P opt freq cam-b3lyp/6-31g(d,p) scrf=(solvent=dichloromethane,pcm)；dushin.out输出文件都是一样的。
   很奇怪这里问题出现在哪里？文件是1-dushin.out    ——感觉dushin还没进行分解，没有正则坐标的计算，也没有振动信息结构的输出；

(, 下载次数 Times of downloads: 0)
(, 下载次数 Times of downloads: 2)
(, 下载次数 Times of downloads: 2)
(, 下载次数 Times of downloads: 2)

② 后面，我们也用同样的方法算了，计算了我的分子（一共50个原子），出现报错：cut= 1.d-7, nber non-redundant internal coords=       257 nv is incorrect
   用的函数基组是：#P opt freq cam-b3lyp/6-31g(d,p) scrf=(solvent=dichloromethane,pcm)
   dushin.out文件是这样的：文件是2-dushin.out
(, 下载次数 Times of downloads: 1)
(, 下载次数 Times of downloads: 0)

log文件太大了，上传不了；

③ 同时，我们也用同样的方法算了一个分子（一共100个原子），出现报错：unconnected atom       72          1    ； unconnected atom in bmatred ——这里是苯环上的C-H键
   用的函数基组是：#P opt freq cam-b3lyp/6-31g(d,p) scrf=(solvent=dichloromethane,pcm)

是我的分子太大了，导致算不出来吗，是需要再去尝试MOMAP这个吗；但是感觉这个dushin还没开始计算；

④ ☆☆☆另一个问题是，分子的重组能应该是S0的正则坐标对重组能贡献+S1的正则坐标对重组能贡献的和嘛？还有一个问题，不同模式对重组能贡献的图，是需要把S0振动模式对λ(I)的贡献和S1振动模式对λ(I)的贡献绘制在一张图上吗，因为看文献分子的重组能只有一个，图也只有一个；——这个问题不懂
(, 下载次数 Times of downloads: 0)

求助求助，非常感谢！

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