想请问大家Multiwfn3.8(dev)版本，关于计算Huang-Rhys

qeheikj · 发表于 Post on 2025-9-24 18:23:54

想请问大家关于Multiwfn3.8(dev)是否有振动分析功能，从而可以计算Huang-Rhys。

wal · 发表于 Post on 2025-9-24 18:29:42

完全不属于波函数分析范畴

qeheikj · 发表于 Post on 2025-9-24 18:33:55

wal 发表于 2025-9-24 18:29
完全不属于波函数分析范畴

那就是没有吧，我老师说可以算，但是我在手册也没找到。以防漏看，我求证一下哈哈

wal · 发表于 Post on 2025-9-24 18:41:25

qeheikj 发表于 2025-9-24 18:33
那就是没有吧，我老师说可以算，但是我在手册也没找到。以防漏看，我求证一下哈哈

肯定不能的找专业程序做吧

qeheikj · 发表于 Post on 2025-9-24 18:50:15

wal 发表于 2025-9-24 18:41
肯定不能的找专业程序做吧

好的好的谢谢

SharkYYX2025 · 发表于 Post on 2025-9-24 23:08:53

在这个帖子有计算方法，可以参考
使用Dushin分解重组能和计算Huang-Rhys因子 - 思想家公社的门口：量子化学·分子模拟·二次元 http://sobereva.com/330

St_Maxwell · 发表于 Post on 2025-9-24 23:39:44

我有一个用python写的小程序，支持绝热Hessian（AH）以及垂直梯度近似（VG） https://github.com/St-Maxwell/Lyra-Py

sobereva · 发表于 Post on 2025-9-25 02:54:46

Gaussian和ORCA都可以算Huang-Rhys因子，北京科音高级量子化学培训班（http://www.keinsci.com/KAQC）讲振动电子光谱的部分专门给了详细例子。第三方程序也可以用dushin算。和Multiwfn没任何直接联系

qeheikj · 发表于 Post on 2025-9-25 10:34:20

SharkYYX2025 发表于 2025-9-24 23:08
在这个帖子有计算方法，可以参考
使用Dushin分解重组能和计算Huang-Rhys因子 - 思想家公社的门口：量子化 ...

好的，谢谢，我有用Dushin试着算，但总是报错。

qeheikj · 发表于 Post on 2025-9-25 10:35:18

St_Maxwell 发表于 2025-9-24 23:39
我有一个用python写的小程序，支持绝热Hessian（AH）以及垂直梯度近似（VG） https://github.com/St-Maxwel ...

好的，谢谢。我来学习一下

qeheikj · 发表于 Post on 2025-9-25 10:47:38

sobereva 发表于 2025-9-25 02:54
Gaussian和ORCA都可以算Huang-Rhys因子，北京科音高级量子化学培训班（http://www.keinsci.com ...

好的，谢谢，我有用dushin进行计算，总是报错，如何解决呢？

St_Maxwell · 发表于 Post on 2025-9-25 14:51:50

本帖最后由 St_Maxwell 于 2025-9-25 15:15 编辑

感觉 RPy2+BP2-S0S2T4-opt-cam3b-free-freq.fchk 这个文件中的Hessian矩阵信息有问题，不知为何有好多为零的值。

---
原来S2或者T4的频率计算有冻结原子坐标，前18个模式是要舍弃的。这导致了初末两个态的简正模式的个数不一致。

qeheikj · 发表于 Post on 2025-9-25 17:11:13

本帖最后由 qeheikj 于 2025-9-25 17:13 编辑

St_Maxwell 发表于 2025-9-25 14:51
感觉 RPy2+BP2-S0S2T4-opt-cam3b-free-freq.fchk 这个文件中的Hessian矩阵信息有问题，不知为何有好多为零 ...

是的，我在优化时固定了部分原子，您说的舍弃前18个模式，是指这些吗，？

St_Maxwell · 发表于 Post on 2025-9-25 21:54:08

本帖最后由 St_Maxwell 于 2025-9-25 22:21 编辑

qeheikj 发表于 2025-9-25 17:11
是的，我在优化时固定了部分原子，您说的舍弃前18个模式，是指这些吗，？

我拿Lyra根据fchk里的Hessian算了下频率，结果得到了好多接近零的值。
一开始我以为自己的代码出毛病了，不过看了Gaussian输出文件发现这个体系有58个原子，但是只打印了156个简正模式（3*58-18）。
然后看到Gaussian还打印了 "6 atoms frozen in the vibrational analysis" 这样一句话，所以我估计舍弃的18个模式就来自这6×3个冻结的自由度。
总之我建议放开冻结原子优化后再算频率，初末态简正模式数不一致的话，计算结果应该是没意义的。
要是末态的结构优化比较困难，也可以考虑采用垂直梯度近似，只需要在S0极小点结构上算末态的能量梯度。

我觉得初末态的结构差异不算很大，垂直梯度近似还算可靠。

qeheikj · 发表于 Post on 2025-9-26 11:16:42

St_Maxwell 发表于 2025-9-25 21:54
我拿Lyra根据fchk里的Hessian算了下频率，结果得到了好多接近零的值。
一开始我以为自己的代码出毛病了 ...

好的谢谢您，优化是芘-嘧啶之间有点弯曲，所以我在优化三线态时候固定了部分原子。我不固定试试计算吧

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

[Multiwfn使用咨询] 想请问大家Multiwfn3.8(dev)版本，关于计算Huang-Rhys

浏览过的版块