计算化学公社

标题: 高斯如何输出非绝热耦合矩阵元？ [打印本页]

作者
Author: WhatAYoumu 时间: 2024-9-6 18:37
标题: 高斯如何输出非绝热耦合矩阵元？
本帖最后由 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为电子波函数） (, 下载次数 Times of downloads: 0) 我自己用mokit读取高斯的波函数到pyscf，用水分子的S0和S1态算了一下Vm，但是Vm的数值很小（1E-16），不确定是不是bug，想要对照一下中间变量，想请问各位老师高斯该如何输出这些矩阵元？

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

作者
Author: 风起~ 时间: 2024-9-6 20:13
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.

作者
Author: WhatAYoumu 时间: 2024-9-7 01:58

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

感谢！

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