计算化学公社

标题: 求助ORCA如何计算vibronic SOC？ [打印本页]

作者
Author: 严珍珍 时间: 2026-1-13 01:19
标题: 求助ORCA如何计算vibronic SOC？
大家好，最近在看一篇怎么计算量子产率的文章(https://pubs.acs.org/doi/10.1021/acs.jctc.3c00931)，下面有文章的附件。其中SOC的计算我不知道我理解的正不正确，所以想问问大家。
(, 下载次数 Times of downloads: 4)
文章中vibronic SOC的公式如下：
(, 下载次数 Times of downloads: 0)

其一阶项用有限差分近似处理：
(, 下载次数 Times of downloads: 0)
文章中d取的是0.01。

按照我的理解，大概是这样：先用优化好的S0结构算一次SOC(0阶项)，然后看频率计算的输出，下图是高斯的频率输出（文章中的频率计算是用的高斯），我截的振动模式1:
(, 下载次数 Times of downloads: 0)
取下面的坐标乘以d (也就是0.01)，用优化的S0坐标加上它得到0+d的坐标，再用优化的S0坐标减去它得到0-d的坐标，在这两个坐标下分别计算它们的SOC得到差值后除以2d，就得到了振动模式1下的一阶项的具体数值。以此类推，得到所有的振动模式下的一阶项，将它们相加，最后再和0阶项相加，得到vibronic SOC。

不知道我的理解对不对，对ORCA和量子产率的计算都不是很熟悉，还请大家多多赐教，谢谢！

作者
Author: sobereva 时间: 2026-1-13 01:42
可以那么算
顺带一提，https://github.com/masaya0222/PyGraSO可以结合Gaussian算的TDDFT的Hessian得到解析的SOCME对核坐标的导数，因此不是必须靠有限差分得到

作者
Author: 严珍珍 时间: 2026-1-13 19:29

sobereva 发表于 2026-1-13 01:42
可以那么算
顺带一提，https://github.com/masaya0222/PyGraSO可以结合Gaussian算的TDDFT的Hessian得到解 ...

谢谢sob老师！

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