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

qeheikj · 发表于 Post on 2025-9-29 21:32:15

本帖最后由 qeheikj 于 2025-9-29 21:37 编辑

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

你好，可以麻烦您有空在帮我看一下吗？我把固定的原子解冻，重新优化了一遍，并计算频率，查看频率fchk文件，c搜索artesian force contants没有为0的情况，还是出现了以下错误。文件有点大，分为两个压缩包上传的。

wal · 发表于 Post on 2025-9-29 21:42:20

本帖最后由 wal 于 2025-9-29 22:06 编辑

qeheikj 发表于 2025-9-29 21:32
你好，可以麻烦您有空在帮我看一下吗？我把固定的原子解冻，重新优化了一遍，并计算频率，查看频率fchk文 ...

我问个题外话为啥你这两个文件计算都是S0的freq 这样分解重组能是要算啥性质呢

结合这个抽象受力我猜你忘写TD了

qeheikj · 发表于 Post on 2025-9-30 15:56:12

wal 发表于 2025-9-29 21:42
我问个题外话为啥你这两个文件计算都是S0的freq 这样分解重组能是要算啥性质呢

你好，我打算计算Huang-Rhys ，我是分开优化结构和计算频率的，两个压缩包分别是基态和三线态计算频率的输出文件。您说的TD是指优化没写root吗？我检查优化输入文件是写了的。

wal · 发表于 Post on 2025-9-30 16:27:42

qeheikj 发表于 2025-9-30 15:56
你好，我打算计算Huang-Rhys ，我是分开优化结构和计算频率的，两个压缩包分别是基态和三线态计算频率的 ...

如果你要用VG，那不是这么用的，你要在同一个结构做一个gs的freq和一个td的freq。如果你要用AH，那你这个freq没有任何意义，结构优化的时候对激发态哪个root做的，freq也必须对激发态哪个root做

qeheikj · 发表于 Post on 2025-10-8 11:07:25

wal 发表于 2025-9-30 16:27
如果你要用VG，那不是这么用的，你要在同一个结构做一个gs的freq和一个td的freq。如果你要用AH，那你这个 ...

VG和Ah是什么代表什么？

wal · 发表于 Post on 2025-10-8 16:19:28

qeheikj 发表于 2025-10-8 11:07
VG和Ah是什么代表什么？

这种基础问题自行搜索吧。指路ORCA手册有详细解释。

qeheikj · 发表于 Post on 2025-10-8 17:56:00

wal 发表于 2025-10-8 16:19
这种基础问题自行搜索吧。指路ORCA手册有详细解释。

好的，谢谢

sobereva · 发表于 Post on 2025-10-9 14:36:09

qeheikj 发表于 2025-10-8 11:07
VG和Ah是什么代表什么？

北京科音高级量子化学培训班（http://www.keinsci.com/KAQC）讲振动电子光谱部分专门说了

清明雨上 · 发表于 Post on 2025-10-20 11:22:34

本帖最后由清明雨上于 2025-10-20 11:34 编辑

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

您好，您的程序非常好用！但是对于我的结构，得到的重组能为7859 (cm-1)，而基于Gaussian FCHT（http://bbs.keinsci.com/forum.php ... ighlight=huang-rhys）结果为2367 (cm-1)，有不小的差距。
对于我的体系而言，对比下来您的程序的计算结果更加合理，请问这是为什么呢？
另外，打扰您一下，如果我在文章中使用了您的程序，我应当如何正确引用呢？

St_Maxwell · 发表于 Post on 2025-10-23 00:11:58

清明雨上发表于 2025-10-20 11:22
您好，您的程序非常好用！但是对于我的结构，得到的重组能为7859 (cm-1)，而基于Gaussian FCHT（http://b ...

你发一下相应的hessian与梯度计算的fchk文件，还有对应的gaussian FCHT计算结果，我仔细比对比对。
如果需要引用，可以引用程序的github仓库网页（不引用也行，只要审稿人或读者没意见）。

清明雨上 · 发表于 Post on 2025-10-23 15:27:20

St_Maxwell 发表于 2025-10-23 00:11
你发一下相应的hessian与梯度计算的fchk文件，还有对应的gaussian FCHT计算结果，我仔细比对比对。
如果 ...

您好，hessian与梯度计算的fchk文件文件有点大了（超过了9.8MB无法在公社上传），请问您是否可以私发一下邮箱或者其他联系方式，我私发您。

St_Maxwell · 发表于 Post on 2025-10-24 14:27:57

清明雨上发表于 2025-10-23 15:27
您好，hessian与梯度计算的fchk文件文件有点大了（超过了9.8MB无法在公社上传），请问您是否可以私发一下 ...

两者的差别来自于计算Duschinsky转动关系时以哪一个态作为参照，以及后续计算Huang-Rhys因子和重组能时是在哪一个态的势能面上。
我看了一下fcclasses的手册（Gaussian的FCHT计算模块来自于这个程序），它计算位移是以初态势能面为参考的。（Shuai等人的文章也是这样定义的，见J. Phys. Chem. A 2010, 114, 7817–7831的(7)式）
但是我倾向于基于末态势能面计算Huang-Rhy因子和重组能，这样和跃迁过程是对应的。
例如初态是S0、末态是S1的吸收过程，垂直激发之后的结构弛豫发生在末态S1的势能面上；而垂直发射过程之后的弛豫发生在末态S0势能面。见Dushin程序的文章J. Chem. Phys. 115, 9103–9109 (2001)的(18)式。
目前我改了一下Lyra，AH分析时会分别输出以初态作为参考和以末态作为参考的结果。对于你的例子，执行python lyra.py evc S0.fchk --AH S1.fchk之后，会先以S0（初态）为参考态输出计算结果，然后以S1（末态）为参考态输出结果。前者和Gaussian结果基本一致，后者则是之前用Lyra得到的结果。你可以自己决定选择哪一个。

P.S. Lyra在计算初态频率时只会将其笛卡尔坐标平移至原点，没有投影到Eckart朝向下（即惯量主轴）。虽然不太影响结果，但是以后会加上。

清明雨上 · 发表于 Post on 2025-10-24 18:13:04

本帖最后由清明雨上于 2025-10-28 11:12 编辑

St_Maxwell 发表于 2025-10-24 14:27
两者的差别来自于计算Duschinsky转动关系时以哪一个态作为参照，以及后续计算Huang-Rhys因子和重组能时是 ...

感谢感谢，行家一出手，就知有没有

以末态作为参考更有物理意义一些，对于我的课题也明显解释的更合理。

清明雨上 · 发表于 Post on 2025-10-28 11:09:31

本帖最后由清明雨上于 2025-10-28 11:47 编辑

St_Maxwell 发表于 2025-10-24 14:27
两者的差别来自于计算Duschinsky转动关系时以哪一个态作为参照，以及后续计算Huang-Rhys因子和重组能时是 ...

您好，我用四点法算出我的体系总的重组能λ(I)＋λ(II)为0.41eV。我在分解程序中的计算：第一次计算initial_state设置为s0, final_state设置为s1；第二次计算initial_state设置为s1, final_state设置为s0，将两次计算总的结果相加得到λ(I)＋λ(II)（不知道这种做法是否正确）。Lyra中以末态为参考算出来总的重组能为1.77eV，Gaussian算出来为0.52eV。

如果计算方式没问题的话，对我而言，还是以初态为参考更合适了。

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

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

评分 Rate

浏览过的版块