两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版后一修订版两侧同时换到之后的修订记录 |
adf:socmatrix [2019/06/16 00:08] – [进行自旋-轨道耦合矩阵元的计算] liu.jun | adf:socmatrix [2020/08/13 21:33] – liu.jun |
---|
====== 如何计算自旋-轨道耦合矩阵元SOCMEs====== | ====== TDDFT计算S→T与T→S跃迁、自旋-轨道耦合矩阵元SOCMEs====== |
| |
**前言:** | **前言:** |
有时候我们需要关心某个特定几何结构下(例如研究$S_i$态到$T_j$的窜跃,则采用$S_i$态结构优化之后的结构),$S_i$态与$T_j$态之间自旋轨道耦合。用算符来表示即:<$S_i$|//SOC//|$T_j$>,也就是自旋-轨道耦合算符,左边乘以$S_i$态、右边乘以$T_j$态,然后在全空间积分得到的一个实数(包括实部和虚部)。而我们关心的,则是这个实数的模。而$T_j$实际上是三个态,因此实际上对应着三个实数。文献中一般关心的是这三个实数的均方根。 | 有时候我们需要关心某个特定几何结构下(例如研究$S_i$态到$T_j$的窜跃,则采用$S_i$态结构优化之后的结构),$S_i$态与$T_j$态之间自旋轨道耦合。用算符来表示即:<$S_i$|//SOC//|$T_j$>,也就是自旋-轨道耦合算符,左边乘以$S_i$态、右边乘以$T_j$态,然后在全空间积分得到的一个实数(包括实部和虚部)。而我们关心的,则是这个实数的模。而$T_j$实际上是三个态,因此实际上对应着三个实数。文献中一般关心的是这三个实数的均方根。 |
| |
| <color blue>另外,计算的是S→T跃迁还是T→S跃迁,主要取决于分子结构。例如需要计算$S_1$→$T_n$的跃迁,分子结构应该采用$S_1$的,计算结果中就是$S_1$→各个$T_n$跃迁的数据。如果是$T_1$→$S_n$跃迁,分子结构采用T1的分子结构,计算结果中,就是$T_1$→各个$S_n$跃迁的数据。</color> |
| |
**步骤:** | **步骤:** |
| |
**此处以<chem>CH4</chem>举例(C1群分子输出结果更简单)** | **此处以<chem>CH4</chem>举例(C1群分子输出结果更简单),本教程仅仅是为了演示软件功能,与体系无关,对所有体系均适用。** |
=====优化分子结构===== | =====优化分子结构===== |
根据用户关心的分子状态,可能是基态的结构,参考:[[adf:geoopt|优化分子结构(详情请点击)]]; | 根据用户关心的分子状态,可能是基态的结构,参考:[[adf:geoopt|优化分子结构(详情请点击)]]; |
</code> | </code> |
//<S|Hso|T>//实际上给出的是单重态与三重态的三个分态之间的自旋轨道耦合常数的实部的平方与虚部的平方之和,总共6项加和,之后开方,再将单位从Hatree转换为$cm^{-1}$(这也是一般文献中所使用的数据),例如<$S_1$|//$H_{so}$//|$T_3$>=4239.35$cm^{-1}$ | //<S|Hso|T>//实际上给出的是单重态与三重态的三个分态之间的自旋轨道耦合常数的实部的平方与虚部的平方之和,总共6项加和,之后开方,再将单位从Hatree转换为$cm^{-1}$(这也是一般文献中所使用的数据),例如<$S_1$|//$H_{so}$//|$T_3$>=4239.35$cm^{-1}$ |
| |
| 文献中,计算系间窜跃速率的公式中,使用的$|V_{soc}|^2$=|//<S|Hso|T>//|$^2/3$ |
| ====$S_O$与其他态之间的SOCME呢?==== |
| |
| 这种情况,并没有将$S_0$考虑进去。如果要计算$S_0$与其他态之间的SOCME,可以在Input - Details - User Input输入: |
| <code> |
| GSCORR |
| </code> |
| 输出文件的SOC列表中就会包含基态。如果体系有对称性,最好将对称性取消(Details - Symmetry - Symbol : NOSYM),否则结果中不列出该数据。 |
| =====辐射速率常数===== |
| |
| 某个辐射的速率常数k,就是该激发态的寿命τ(即上文中tau)的倒数。 |
| |
| //k//<sub>i</sub>=1/τ<sub>i</sub> |
| |
| 文献重现,参考:https://www.scm.com/adf-modeling-suite/wizard/organic-electronics/modeling-phosphorescent-lifetimes-of-oled-emitters/ |
| |
| <color blue>三重态在考虑旋轨耦合之后,劈裂成三个态。三个态的寿命和辐射跃迁速率是不一样的。如果没有其他的跃迁机理,三个三重态都直接跃迁到基态S0,那么观察到的辐射寿命就是从三个态的平均跃迁速率得到的:</color> |
| |
| τ<sub>av</sub>=3/(//k//<sub>1</sub>+//k//<sub>2</sub>+//k//<sub>3</sub>) (参考Phys.Chem.Chem.Phys., 2014, 16, 14523) |
| |
| |
| <color blue>而辐射跃迁速率,则速率最大的那个态的值更重要,例如其中一个值特别大,另外两个小几十倍,那么这个最大值就代表了这个三重态的辐射跃迁速率。ADF给出的是“辐射跃迁”的速率常数,不含非辐射跃迁。 |
| </color> |