Tn和S1之间的SOC耦合常数计算求助

幸运兔

想计算Tn和S1之间的SOC， 输入文件如下
**DALTON INPUT
.RUN RESPONSE
.NEWTRA
**INTEGRALS
.SPIN-ORBIT
**WAVE FUNCTIONS
.DFT
CAM-B3LYP
**RESPONSE
.TRPFLG
*LINEAR
.SINGLE RESIDUE
.PROPRT
X1SPNORB
.PROPRT
Y1SPNORB
.PROPRT
Z1SPNORB
.PROPRT
X2SPNORB
.PROPRT
Y2SPNORB
.PROPRT
Z2SPNORB
.ROOTS
3
**END OF DALTON INPUT
请问，现在计算得到的SOC是 Tn与S0之间的还是和S1之间的耦合常数，如果想计算Tn和S1之间耦合常数，改如何更改输入文件

pwzhou

**DALTON INPUT
.RUN RESPONSE
**INTEGRALS
.MNF-SO
**WAVE FUNCTIONS
.DFT
B3LYP
**RESPONSE
*QUADRATIC
.DOUBLE RESIDUE
.ISPABC
1 0 1
.PROPRT
X1MNF-SO
.PROPRT
Y1MNF-SO
.PROPRT
Z1MNF-SO
.ROOTS
3
**END OF INPUT

幸运兔

pwzhou 发表于 2021-1-6 10:08
**DALTON INPUT
.RUN RESPONSE
**INTEGRALS

非常感谢，我用这个输入文件进了计算，请问是不是应该读取下面最后一列数据？（输出文件已经上传）
@ B and C excitation energies, moment:   0.19511812   0.12866393   0.00000005
@ B and C excitation energies, moment:   0.19511812   0.12866393  -0.00000005
@ B and C excitation energies, moment:   0.19511812   0.12866393   0.00000041
@ B and C excitation energies, moment:   0.22592739   0.12866393   0.00000013
@ B and C excitation energies, moment:   0.22592739   0.12866393   0.00000013

但是三个态，每个态对应着三个X1MNF-SO,三个Y1MNF-SO，三个Z1MNF-SO的输出，不知道该读取哪些数据，非常感谢

pwzhou

幸运兔 发表于 2021-1-6 16:34
非常感谢，我用这个输入文件进了计算，请问是不是应该读取下面最后一列数据？（输出文件已经上传）
@ B  ...

@ Transition moment <B | A - <A> | C> in a.u. for
@ A operator label,    symmetry, spin:      X1MNF-SO    1    1
@ B excited state no., symmetry, spin:             1    1    0 (spin=0，单重态，B为S1）
@ C excited state no., symmetry, spin:             1    1    1 (spin=1，三重态，C为T1）

@ B and C excitation energies, moment:   0.19511812   0.12866393   0.00000005 （三个数字依次为S1的激发能，T1的激发能，S1和T1的SOC的X分量）

取得S1和T1的SOC的X，Y，Z三个方向分量以后，利用如下公式求出SOC，sqrt[(X^2+Y^2+Z^2)/3]
即三者的平方加和以后除以3，然后在开根号。
最后，注意单位，Dalton给出的是数值是原子单位

其他几个态之间的SOC依次类推即可。

幸运兔

pwzhou 发表于 2021-1-6 20:01
@ Transition moment  in a.u. for
@ A operator label,    symmetry, spin:      X1MNF-SO    1    1
...

感激不尽！终于搞明白了

幸运兔

pwzhou 发表于 2021-1-6 20:01
@ Transition moment  in a.u. for
@ A operator label,    symmetry, spin:      X1MNF-SO    1    1
...

我注意到这篇文献SOC的计算是sqrt[(X^2+Y^2+Z^2)]/3,现在有点迷惑，这个3是在根号下还是根号外，非常感谢你解决了我这么多天的烦恼

pwzhou

幸运兔 发表于 2021-1-7 09:39
我注意到这篇文献SOC的计算是sqrt[(X^2+Y^2+Z^2)]/3,现在有点迷惑，这个3是在根号下还是根号外，非常感谢 ...

我也有点迷糊了，我之前看的是这篇文献：
https://pubs.rsc.org/en/content/ ... 858e/unauth#tab2fna

幸运兔

嗯嗯，确实是哈，继续啃

wzkchem5

幸运兔 发表于 2021-1-7 09:39
我注意到这篇文献SOC的计算是sqrt[(X^2+Y^2+Z^2)]/3,现在有点迷惑，这个3是在根号下还是根号外，非常感谢 ...

是在里面，这个结论其实很好记，因为费米黄金规则里面，SOC矩阵元只以模平方的形式出现，所以只有对SOC矩阵元的模平方取平均才是有意义的，也就是应该取方均根。
有的程序（比如BDF）给的不是xyz分量，而是m=-1,0,1的分量，得到的矩阵元是复数，这时盲目套sqrt[(X^2+Y^2+Z^2)/3]就错了，应该用模平方，sqrt[(|Hsoc(m=-1)|^2+|Hsoc(m=0)|^2+|Hsoc(m=1)|^2)/3]

幸运兔

wzkchem5 发表于 2021-1-7 11:49
是在里面，这个结论其实很好记，因为费米黄金规则里面，SOC矩阵元只以模平方的形式出现，所以只有对SOC矩 ...

谢谢

幸运兔

wzkchem5 发表于 2021-1-7 11:49
是在里面，这个结论其实很好记，因为费米黄金规则里面，SOC矩阵元只以模平方的形式出现，所以只有对SOC矩 ...

赞！

rabbit0602

幸运兔 发表于 2021-1-6 16:34
非常感谢，我用这个输入文件进了计算，请问是不是应该读取下面最后一列数据？（输出文件已经上传）
@ B  ...

请问如果想要计算S1与T1之间的SOC耦合常数，输入文件mol是要用S1的结构吗？

wzkchem5

rabbit0602 发表于 2021-5-27 08:25
请问如果想要计算S1与T1之间的SOC耦合常数，输入文件mol是要用S1的结构吗？

对，一般用初态结构。例外情况是如果你知道这个体系的ISC主要在MECP点附近发生，此时应该用MECP结构。

rabbit0602

wzkchem5 发表于 2021-5-27 20:11
对，一般用初态结构。例外情况是如果你知道这个体系的ISC主要在MECP点附近发生，此时应该用MECP结构。

谢谢您！