高斯如何输出非绝热耦合矩阵元？

WhatAYoumu · 发表于 Post on 2024-9-6 18:37:39

本帖最后由 WhatAYoumu 于 2024-9-6 18:41 编辑

各位老师好，

最近想复现一下这篇文章（J. Phys. Chem. A 2021, 125, 9000−9010）中的速率计算公式。在Internal Conversion中需要计算Vm这一项，里面需要非绝热耦合矩阵元(NACME)<k|∇(1/r)|l>（图1， k和l为电子波函数）

我自己用mokit读取高斯的波函数到pyscf，用水分子的S0和S1态算了一下Vm，但是Vm的数值很小（1E-16），不确定是不是bug，想要对照一下中间变量，想请问各位老师高斯该如何输出这些矩阵元？

我自己写的代码在压缩包里

WhatAYoumu · 发表于 Post on 2024-9-7 01:58:23

风起~ 发表于 2024-9-6 20:13
NonAdiabaticCoupling
Requests that the ground-to-excited-state non-adiabatic coupling be computed. ...

感谢！

风起~ · 发表于 Post on 2024-9-6 20:13:36

NonAdiabaticCoupling
Requests that the ground-to-excited-state non-adiabatic coupling be computed. NAC is a synonym for
this option. NoNonAdiabaticCoupling and NoNAC suppress this behavior. The default is NoNAC when
computing energies or energies+gradients because the extra cost is non-trivial. The default is NAC during
frequency calculations where the extra cost is negligible.

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

[Gaussian/gview] 高斯如何输出非绝热耦合矩阵元？