求助ORCA如何计算vibronic SOC？

严珍珍 · 发表于 Post on 2026-1-13 01:19:03

大家好，最近在看一篇怎么计算量子产率的文章(https://pubs.acs.org/doi/10.1021/acs.jctc.3c00931)，下面有文章的附件。其中SOC的计算我不知道我理解的正不正确，所以想问问大家。

modeling-the-fluorescence-quantum-yields-of-aromatic-compounds-benchmarking-the-.pdf (3.31 MB, 下载次数 Times of downloads: 4)
文章中vibronic SOC的公式如下：

其一阶项用有限差分近似处理：

文章中d取的是0.01。

按照我的理解，大概是这样：先用优化好的S0结构算一次SOC(0阶项)，然后看频率计算的输出，下图是高斯的频率输出（文章中的频率计算是用的高斯），我截的振动模式1:

取下面的坐标乘以d (也就是0.01)，用优化的S0坐标加上它得到0+d的坐标，再用优化的S0坐标减去它得到0-d的坐标，在这两个坐标下分别计算它们的SOC得到差值后除以2d，就得到了振动模式1下的一阶项的具体数值。以此类推，得到所有的振动模式下的一阶项，将它们相加，最后再和0阶项相加，得到vibronic SOC。

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

sobereva · 发表于 Post on 2026-1-13 01:42:58

可以那么算
顺带一提，https://github.com/masaya0222/PyGraSO可以结合Gaussian算的TDDFT的Hessian得到解析的SOCME对核坐标的导数，因此不是必须靠有限差分得到

严珍珍 · 发表于 Post on 2026-1-13 19:29:55

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

谢谢sob老师！

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

[ORCA] 求助ORCA如何计算vibronic SOC？

评分 Rate

浏览过的版块