Kalinite 发表于 2022-5-9 13:24
这篇文献我以前看过。不管怎么定义，总是绕不开x1和x2（也就是文献里的g和h），但是不知道怎么让程序输出 ...

请问老师会画h和g矢量了吗？

Kalinite 发表于 2022-5-9 13:24
这篇文献我以前看过。不管怎么定义，总是绕不开x1和x2（也就是文献里的g和h），但是不知道怎么让程序输出 ...

OpenMolcas的CI结构优化收敛后会输出x1和x2. 手动计算的话就计算两个态的梯度g1和g2，得到梯度差值g。再计算两个态的非绝热耦合h。让后通过正交化把g和h投影到x1和x2。推导在J. Chem. Theory Comput. 2016, 12, 3636−3653的公式38–41。

要注意的是，CI的拓扑性质只反应附近的势能面情况，不一定能严谨地预测产物的方向和结构。因为很多激发态反应不一定要经过CI。比如FSSH模拟非绝热跃迁时，两个态的能量差可能有0.5eV左右。

wxhwbh 发表于 2022-5-8 18:48
看看一些做非绝热体系研究的文章如何定义The branching plane, 比如: J. Phys. Chem. B 2016, 120, 1940& ...

这篇文献我以前看过。不管怎么定义，总是绕不开x1和x2（也就是文献里的g和h），但是不知道怎么让程序输出出来，原理上应该是计算了这个内容的。

Kalinite 发表于 2022-5-8 16:39
请问sob老师，如何获得x1和x2呢？x1就是梯度之差，这个好办，但是x2似乎没有现成的输出。
具体到投影， ...

看看一些做非绝热体系研究的文章如何定义The branching plane, 比如: J. Phys. Chem. B 2016, 120, 1940−1949, DOI: 10.1021/acs.jpcb.5b09838

sobereva 发表于 2022-5-8 09:22
对CI点计算两个态的原子的受力，投影到x1和x2空间再对比

请问sob老师，如何获得x1和x2呢？x1就是梯度之差，这个好办，但是x2似乎没有现成的输出。
具体到投影，取x1和x2的单位法向量n，投影结果就是grad-(grad·n)n吗？

我目前是想用SF-TDDFT优化S0-S1之间的交叉，试了Q-Chem和GAMESS。Q-Chem算的时候自旋污染严重，在最低5个激发态里无法指认出单重态；GAMESS支持ROKS (BHandHLYP)参考态做SF-TDDFT，能成功优化，自旋污染不是特别严重，结构和能量上看起来也比较合理。

对CI点计算两个态的原子的受力，投影到x1和x2空间再对比