关于利用pysoc拆分分子中原子 SOC 贡献的方法请教

youngzq

各位老师好，我最近想定量拆分分子里各个原子（比如重原子碘和硫）对 S1​→T1​ 的自旋轨道耦合（SOC）贡献。

我目前在 PySOC 和 MolSOC 中采用的方法是：在单电子 Breit-Pauli 算符下，通过修改soc.py文件，让其生成的输入文件（molsoc.inp）把非目标原子的核电荷缩放因子设为 0.0，只保留目标原子为1.0（相当于一个数值开关），并且把init.py中的soc_key设置为['ANG', 'ONE']从而直接在积分级别把其他原子屏蔽掉。依次计算各原子后，最后和分子整体进行对比。

有几个问题想向各位老师请教：

1. 这种通过“屏蔽其他原子电荷”来计算单原子贡献的方案，在理论上和发表论文时合不合适？

2. 我发现有的状态下，可能存在复数相位的相消干涉，单独重原子的 SOC 贡献算出来甚至比分子总 SOC 还要大一点，这是否是正常的？

3. 有没有其他更推荐的、也是基于 TDDFT 级别计算原子 SOC 贡献的方法？

youngzq

谢谢Sob老师！我现在就去试试

sobereva

把原本的结果，和只把特定那个原子的有效核电荷设0，得到的SOC相关结果求差，衡量它的影响，有可能更为合理。至少可以和你当前得到的结果进行对照，不管绝对值如何，若得到的不同原子的贡献趋势与之前的做法相同，会更有说服力。