|
发布时间: 2025-9-24 18:23
正文摘要:想请问大家关于Multiwfn3.8(dev)是否有振动分析功能，从而可以计算Huang-Rhys。 |
本帖最后由 清明雨上 于 2025-10-28 11:47 编辑 St_Maxwell 发表于 2025-10-24 14:27 您好，我用四点法算出我的体系总的重组能λ(I)＋λ(II)为0.41eV。我在分解程序中的计算：第一次计算initial_state设置为s0, final_state设置为s1；第二次计算initial_state设置为s1, final_state设置为s0，将两次计算总的结果相加得到λ(I)＋λ(II)（不知道这种做法是否正确）。Lyra中以末态为参考算出来总的重组能为1.77eV，Gaussian算出来为0.52eV。 如果计算方式没问题的话，对我而言，还是以初态为参考更合适了。 |
本帖最后由 清明雨上 于 2025-10-28 11:12 编辑 St_Maxwell 发表于 2025-10-24 14:27 感谢感谢，行家一出手，就知有没有  以末态作为参考更有物理意义一些，对于我的课题也明显解释的更合理。 |
清明雨上 发表于 2025-10-23 15:27 两者的差别来自于计算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朝向下（即惯量主轴）。虽然不太影响结果，但是以后会加上。 |
| 参与人数Participants 2 | eV +10 | 收起 理由Reason |
|---|---|---|
|
| + 5 | 谢谢 |
|
| + 5 |
St_Maxwell 发表于 2025-10-23 00:11 您好，hessian与梯度计算的fchk文件文件有点大了（超过了9.8MB无法在公社上传），请问您是否可以私发一下邮箱或者其他联系方式，我私发您。 |
清明雨上 发表于 2025-10-20 11:22 你发一下相应的hessian与梯度计算的fchk文件，还有对应的gaussian FCHT计算结果，我仔细比对比对。 如果需要引用，可以引用程序的github仓库网页（不引用也行，只要审稿人或读者没意见）。 |
本帖最后由 清明雨上 于 2025-10-20 11:34 编辑 St_Maxwell 发表于 2025-9-24 23:39 您好，您的程序非常好用！但是对于我的结构，得到的重组能为7859 (cm-1)，而基于Gaussian FCHT（http://bbs.keinsci.com/forum.php ... ighlight=huang-rhys）结果为2367 (cm-1)，有不小的差距。 对于我的体系而言，对比下来您的程序的计算结果更加合理，请问这是为什么呢？ 另外，打扰您一下，如果我在文章中使用了您的程序，我应当如何正确引用呢？ |
qeheikj 发表于 2025-10-8 11:07 北京科音高级量子化学培训班（http://www.keinsci.com/KAQC）讲振动电子光谱部分专门说了 
|
wal 发表于 2025-10-8 16:19 好的，谢谢 |
qeheikj 发表于 2025-10-8 11:07 这种基础问题自行搜索吧。指路ORCA手册有详细解释。 |
wal 发表于 2025-9-30 16:27 VG和Ah是什么代表什么？ |
qeheikj 发表于 2025-9-30 15:56 如果你要用VG，那不是这么用的，你要在同一个结构做一个gs的freq和一个td的freq。如果你要用AH，那你这个freq没有任何意义，结构优化的时候对激发态哪个root做的，freq也必须对激发态哪个root做 |
wal 发表于 2025-9-29 21:42 你好，我打算计算Huang-Rhys ，我是分开优化结构和计算频率的，两个压缩包分别是基态和三线态计算频率的输出文件。您说的TD是指优化没写root吗？我检查优化输入文件是写了的。 |
本帖最后由 wal 于 2025-9-29 22:06 编辑 qeheikj 发表于 2025-9-29 21:32 我问个题外话 为啥你这两个文件计算都是S0的freq 这样分解重组能是要算啥性质呢
结合这个抽象受力我猜你忘写TD了 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-1-24 01:30 , Processed in 0.222377 second(s), 27 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites