计算化学公社

标题: Tn和S1之间的SOC耦合常数计算求助 [打印本页]

作者
Author: 幸运兔 时间: 2021-1-5 19:07
标题: 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之间耦合常数，改如何更改输入文件

作者
Author: pwzhou 时间: 2021-1-6 10:08
**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

作者
Author: 幸运兔 时间: 2021-1-6 16:34

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的输出，不知道该读取哪些数据，非常感谢

作者
Author: pwzhou 时间: 2021-1-6 20:01

幸运兔发表于 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依次类推即可。

作者
Author: 幸运兔 时间: 2021-1-7 08:45

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

感激不尽！终于搞明白了

作者
Author: 幸运兔 时间: 2021-1-7 09:39

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是在根号下还是根号外，非常感谢你解决了我这么多天的烦恼

作者
Author: pwzhou 时间: 2021-1-7 10:14

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

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

作者
Author: 幸运兔 时间: 2021-1-7 11:24
嗯嗯，确实是哈，继续啃

作者
Author: wzkchem5 时间: 2021-1-7 11:49

幸运兔发表于 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]

作者
Author: 幸运兔 时间: 2021-1-7 12:32

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

谢谢

作者
Author: 幸运兔 时间: 2021-1-7 12:33

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

赞！

作者
Author: rabbit0602 时间: 2021-5-27 15:25

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

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

作者
Author: wzkchem5 时间: 2021-5-27 20:11

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

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

作者
Author: rabbit0602 时间: 2021-5-28 13:12

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

谢谢您！

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